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

Online Help

All Products    Maple    MapleSim


Home : Support : Online Help : Mathematics : Differential Equations : Rif : DEtools/rifread

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

Calling Sequence

rifread(filename)

Parameters

filename

-

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

Description

• 

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.

• 

Please note that the system obtained using rifread is not in final form and may have redundant equations or unresolved integrability conditions.

Examples

withDEtools:

Consider the following system, run with rifsimp.

sys1:=2y2ξx,y,2y2ηx,y22yxξx,y,3y2yξx,y+22yxηx,y2x2ξx,y,2ηx,yy2y2xξx,y+y2yηx,y+2x2ηx,y

sys1:=2y2ξx,y,2y2ηx,y22yxξx,y,3y2yξx,y+22yxηx,y2x2ξx,y,2ηx,yy2y2xξx,y+y2yηx,y+2x2ηx,y

(1)

ans1:=rifsimpsys1,ξ,η,store

ans1:=tableSolved=xξx,y=12ηx,yy,xηx,y=0,yξx,y=0,yηx,y=ηx,yy

(2)

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

temp1:=rifread

temp1:=tableSolved=yξx,y=0,2y2ηx,y=0,xηx,y=0,xξx,y=12ηx,yy,yηx,y=ηx,yy

(3)

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:=rifsimpsys1,ξ,η,store=tmpstore

ans1:=tableSolved=xξx,y=12ηx,yy,xηx,y=0,yξx,y=0,yηx,y=ηx,yy

(4)

temp1:=rifreadtmpstore

temp1:=tableSolved=yξx,y=0,2y2ηx,y=0,xηx,y=0,xξx,y=12ηx,yy,yηx,y=ηx,yy

(5)

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:=yξx,yy2y2ξx,yy2,ηx,y22yxξx,yy2+2y2ηx,yy2yηx,yy2yξx,yy+2yξx,yay2,2yxηx,y2x2ξx,yy2+22yxηx,yy2+ηx,yxξx,yy+xξx,yay2+3y3yξx,ya23y4yξx,yb23y5yξx,yb2+3y2yξx,ya2,2x2ηx,yy2+y2xηx,yayxηx,y+2y3xξx,ya22y4xξx,yb2y3yηx,ya2+y4yηx,yb2+ηx,yy2a22ηx,yy3b2y2yηx,ya23ηx,yy4b22y5xξx,yb2+2y2xξx,ya2+y5yηx,yb2

sys2:=yξx,yy2y2ξx,yy2,ηx,y22yxξx,yy2+2y2ηx,yy2yηx,yy2yξx,yy+2yξx,yay2,2yxηx,y2x2ξx,yy2+22yxηx,yy2+ηx,yxξx,yy+xξx,yay2+3y3yξx,ya23y4yξx,yb23y5yξx,yb2+3y2yξx,ya2,2x2ηx,yy2+y2xηx,yayxηx,y+2y3xξx,ya22y4xξx,yb2y3yηx,ya2+y4yηx,yb2+ηx,yy2a22ηx,yy3b2y2yηx,ya23ηx,yy4b22y5xξx,yb2+2y2xξx,ya2+y5yηx,yb2

(6)

If infolevel is set to a greater integer (possible settings are 1 through 5), more detailed information about the computation method is displayed.

infolevel[rifsimp]:=1

infolevelrifsimp:=1

(7)

rr:=rifsimpsys2,ξ,η,storeall,casesplit

      Storing current system in RifStorage_1_1.mStoring 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

The system has been identified as follows:

The system has the following dependent variables:

ξ,η

The system has the following independent variables:

x,y

The system has the following constants:

a,b

The following are to be treated as solve variables:

ξ,η,a,b

rr:=table1=tableCase=b0,ηx,y,Pivots=b0,Solved=xξx,y=0,yξx,y=0,ηx,y=0,2=tableCase=b=0,ηx,y,aa10,ηx,y,Pivots=a0,a10,Solved=xξx,y=0,yξx,y=0,ηx,y=0,b=0,3=tableCase=b=0,ηx,y,aa1=0,ηx,y,a0,ηx,y,Solved=xξx,y=0,yξx,y=0,ηx,y=0,a=1,b=0,4=tableCase=b=0,ηx,y,aa1=0,ηx,y,a=0,ηx,y,Solved=xξx,y=ηx,yy,xηx,y=0,yξx,y=0,yηx,y=ηx,yy,a=0,b=0,casecount=4

(8)

rifreadRifStorage_2_13

tableDiffConstraint=0=aηx,ya1,0=aηx,ya1,0=aηx,ya1,0=aηx,y9a3y9a2y+4ay2a+2,Solved=yξx,y=0,b=0,xξx,y=ηx,yy,xηx,y=a3a2ηx,y,yηx,y=ηx,yy

(9)

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

See Also

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


Download Help Document

Was this information helpful?



Please add your Comment (Optional)
E-mail Address (Optional)
What is ? This question helps us to combat spam