rifread - load intermediate computation as saved using the 'store' options of rifsimp

rifread( filename )

Parameters

filename - (optional) name of the file containing the partial rifsimp computation.

• The rifread command loads a partial rifsimp computation that was run using the store or storeall options (see rifsimp[adv_options] ).

• The storage options are most useful for large or complex computations, where the resources required to complete the computation may exceed the capability of the machine. The state of the system at the last iteration (using store ), or at all previous iterations and/or cases (using storeall ), can be retrieved using this command.
• When called with no arguments, the file "RifStorage.m" is used. This is the default file name for use with the store option of rifsimp . If a file name was specified for the rifsimp run, or if the storeall option was used (storing all iterations and/or cases in separate files), the file name must be included in the rifread command.

> with(Rif):

Consider the following system, run with rifsimp .

> sys1 := [diff(xi(x,y),y,y),
diff(eta(x,y),y,y)-2*diff(xi(x,y),x,y),
-3*y^2*diff(xi(x,y),y)+2*diff(eta(x,y),x,y)-diff(xi(x,y),x,x),
-2*eta(x,y)*y-2*y^2*diff(xi(x,y),x)
+y^2*diff(eta(x,y),y)+diff(eta(x,y),x,x)];

> ans1:=rifsimp(sys1,[xi,eta],store);

Partial results can be obtained using rifread , even if the full computation did not succeed (notice the presence of the redundant equation eta[yy]=0 ).

> temp1 := rifread();

Note that the results obtained above were stored in the file "RifStorage.m".

The same example can have results stored under a different name. Here the temp results are stored in "tmpstore.m":

> ans1:=rifsimp(sys1,[xi,eta],store="tmpstore");

> temp1 := rifread("tmpstore");

It is also possible to store the results at the end of each iteration and for each separate case using storeall . Increasing the infolevel will display the file name as each partial calculation is saved.

> sys2 := [-diff(xi(x,y),y)*y-diff(xi(x,y),y,y)*y^2,
eta(x,y)-2*diff(xi(x,y),x,y)*y^2+diff(eta(x,y),y,y)*y^2
-diff(eta(x,y),y)*y-2*diff(xi(x,y),y)*y
+2*diff(xi(x,y),y)*a*y^2,
-2*y*diff(eta(x,y),x)-diff(xi(x,y),x,x)*y^2
+2*diff(eta(x,y),x,y)*y^2+eta(x,y)-diff(xi(x,y),x)*y
+diff(xi(x,y),x)*a*y^2+3*y^3*diff(xi(x,y),y)*a^2
-3*y^4*diff(xi(x,y),y)*b^2-3*y^5*diff(xi(x,y),y)*b^2
+3*y^2*diff(xi(x,y),y)*a^2,
diff(eta(x,y),x,x)*y^2+y^2*diff(eta(x,y),x)*a
-y*diff(eta(x,y),x)+2*y^3*diff(xi(x,y),x)*a^2
-2*y^4*diff(xi(x,y),x)*b^2-y^3*diff(eta(x,y),y)*a^2
+y^4*diff(eta(x,y),y)*b^2+eta(x,y)*y^2*a^2
-2*eta(x,y)*y^3*b^2-y^2*diff(eta(x,y),y)*a^2
-3*eta(x,y)*y^4*b^2-2*y^5*diff(xi(x,y),x)*b^2
+2*y^2*diff(xi(x,y),x)*a^2+y^5*diff(eta(x,y),y)*b^2];

> infolevel[rifsimp]:=1;

> rr:=rifsimp(sys2,[xi,eta],storeall,casesplit);

      Storing current system in RifStorage_1_1.m

      Storing current system in RifStorage_1_2.m

      Storing current system in RifStorage_1_3.m

      Storing current system in RifStorage_1_4.m

      Storing current system in RifStorage_1_5.m

      Storing current system in RifStorage_1_6.m

      Storing current system in RifStorage_1_7.m

      Storing current system in RifStorage_1_8.m

      Storing current system in RifStorage_1_9.m

      Storing current system in RifStorage_1_10.m

      Storing current system in RifStorage_1_11.m

      Storing current system in RifStorage_1_12.m

      Storing current system in RifStorage_1_13.m

      Storing current system in RifStorage_1_14.m

      Storing current system in RifStorage_2_13.m

      Storing current system in RifStorage_2_14.m

      Storing current system in RifStorage_2_15.m

      Storing current system in RifStorage_3_14.m

      Storing current system in RifStorage_3_15.m

      Storing current system in RifStorage_3_16.m

      Storing current system in RifStorage_4_15.m

> rifread("RifStorage_2_13");

A description of the meaning of each entry appearing in the output tables can be found on the rifsimp[output] page.

rifsimp , rifsimp[adv_options] , rifsimp[options] , rifsimp[output] , rifsimp[overview]