DEtools
rifread
load intermediate computation as saved using the 'store' options of rifsimp
Calling Sequence
Parameters
Description
Examples
rifread(filename)
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.
Please note that the system obtained using rifread is not in final form and may have redundant equations or unresolved integrability conditions.
with⁡DEtools:
Consider the following system, run with rifsimp.
sys1≔diff⁡ξ⁡x,y,y,y,diff⁡η⁡x,y,y,y−2⁢diff⁡ξ⁡x,y,x,y,−3⁢y2⁢diff⁡ξ⁡x,y,y+2⁢diff⁡η⁡x,y,x,y−diff⁡ξ⁡x,y,x,x,−2⁢η⁡x,y⁢y−2⁢y2⁢diff⁡ξ⁡x,y,x+y2⁢diff⁡η⁡x,y,y+diff⁡η⁡x,y,x,x
sys1≔∂2∂y2ξ⁡x,y,∂2∂y2η⁡x,y−2⁢∂2∂x∂yξ⁡x,y,−3⁢y2⁢∂∂yξ⁡x,y+2⁢∂2∂x∂yη⁡x,y−∂2∂x2ξ⁡x,y,−2⁢η⁡x,y⁢y−2⁢y2⁢∂∂xξ⁡x,y+y2⁢∂∂yη⁡x,y+∂2∂x2η⁡x,y
ans1≔rifsimp⁡sys1,ξ,η,store
ans1≔table⁡Solved=∂∂xξ⁡x,y=−η⁡x,y2⁢y,∂∂xη⁡x,y=0,∂∂yξ⁡x,y=0,∂∂yη⁡x,y=η⁡x,yy
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≔table⁡Solved=∂∂yξ⁡x,y=0,∂2∂y2η⁡x,y=0,∂∂xη⁡x,y=0,∂∂xξ⁡x,y=−η⁡x,y2⁢y,∂∂yη⁡x,y=η⁡x,yy
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,ξ,η,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⁡ξ⁡x,y,y⁢y−diff⁡ξ⁡x,y,y,y⁢y2,η⁡x,y−2⁢diff⁡ξ⁡x,y,x,y⁢y2+diff⁡η⁡x,y,y,y⁢y2−diff⁡η⁡x,y,y⁢y−2⁢diff⁡ξ⁡x,y,y⁢y+2⁢diff⁡ξ⁡x,y,y⁢a⁢y2,−2⁢y⁢diff⁡η⁡x,y,x−diff⁡ξ⁡x,y,x,x⁢y2+2⁢diff⁡η⁡x,y,x,y⁢y2+η⁡x,y−diff⁡ξ⁡x,y,x⁢y+diff⁡ξ⁡x,y,x⁢a⁢y2+3⁢y3⁢diff⁡ξ⁡x,y,y⁢a2−3⁢y4⁢diff⁡ξ⁡x,y,y⁢b2−3⁢y5⁢diff⁡ξ⁡x,y,y⁢b2+3⁢y2⁢diff⁡ξ⁡x,y,y⁢a2,diff⁡η⁡x,y,x,x⁢y2+y2⁢diff⁡η⁡x,y,x⁢a−y⁢diff⁡η⁡x,y,x+2⁢y3⁢diff⁡ξ⁡x,y,x⁢a2−2⁢y4⁢diff⁡ξ⁡x,y,x⁢b2−y3⁢diff⁡η⁡x,y,y⁢a2+y4⁢diff⁡η⁡x,y,y⁢b2+η⁡x,y⁢y2⁢a2−2⁢η⁡x,y⁢y3⁢b2−y2⁢diff⁡η⁡x,y,y⁢a2−3⁢η⁡x,y⁢y4⁢b2−2⁢y5⁢diff⁡ξ⁡x,y,x⁢b2+2⁢y2⁢diff⁡ξ⁡x,y,x⁢a2+y5⁢diff⁡η⁡x,y,y⁢b2
sys2≔−∂∂yξ⁡x,y⁢y−∂2∂y2ξ⁡x,y⁢y2,η⁡x,y−2⁢∂2∂x∂yξ⁡x,y⁢y2+∂2∂y2η⁡x,y⁢y2−∂∂yη⁡x,y⁢y−2⁢∂∂yξ⁡x,y⁢y+2⁢∂∂yξ⁡x,y⁢a⁢y2,−2⁢y⁢∂∂xη⁡x,y−∂2∂x2ξ⁡x,y⁢y2+2⁢∂2∂x∂yη⁡x,y⁢y2+η⁡x,y−∂∂xξ⁡x,y⁢y+∂∂xξ⁡x,y⁢a⁢y2+3⁢y3⁢∂∂yξ⁡x,y⁢a2−3⁢y4⁢∂∂yξ⁡x,y⁢b2−3⁢y5⁢∂∂yξ⁡x,y⁢b2+3⁢y2⁢∂∂yξ⁡x,y⁢a2,∂2∂x2η⁡x,y⁢y2+y2⁢∂∂xη⁡x,y⁢a−y⁢∂∂xη⁡x,y+2⁢y3⁢∂∂xξ⁡x,y⁢a2−2⁢y4⁢∂∂xξ⁡x,y⁢b2−y3⁢∂∂yη⁡x,y⁢a2+y4⁢∂∂yη⁡x,y⁢b2+η⁡x,y⁢y2⁢a2−2⁢η⁡x,y⁢y3⁢b2−y2⁢∂∂yη⁡x,y⁢a2−3⁢η⁡x,y⁢y4⁢b2−2⁢y5⁢∂∂xξ⁡x,y⁢b2+2⁢y2⁢∂∂xξ⁡x,y⁢a2+y5⁢∂∂yη⁡x,y⁢b2
If infolevel is set to a greater integer (possible settings are 1 through 5), more detailed information about the computation method is displayed.
infolevelrifsimp≔1
rr≔rifsimp⁡sys2,ξ,η,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≔table⁡1=table⁡Solved=∂∂xξ⁡x,y=0,∂∂yξ⁡x,y=0,η⁡x,y=0,Case=b≠0,η⁡x,y,Pivots=b≠0,2=table⁡Solved=∂∂xξ⁡x,y=0,∂∂yξ⁡x,y=0,η⁡x,y=0,b=0,Case=b=0,η⁡x,y,a⁢a−1≠0,η⁡x,y,Pivots=a≠0,a−1≠0,3=table⁡Solved=∂∂xξ⁡x,y=0,∂∂yξ⁡x,y=0,η⁡x,y=0,a=1,b=0,Case=b=0,η⁡x,y,a⁢a−1=0,η⁡x,y,a≠0,η⁡x,y,4=table⁡Solved=∂∂xξ⁡x,y=η⁡x,yy,∂∂xη⁡x,y=0,∂∂yξ⁡x,y=0,∂∂yη⁡x,y=η⁡x,yy,a=0,b=0,Case=b=0,η⁡x,y,a⁢a−1=0,η⁡x,y,a=0,η⁡x,y,casecount=4
rifread⁡RifStorage_2_13
table⁡Solved=∂∂yξ⁡x,y=0,b=0,∂∂xξ⁡x,y=η⁡x,yy,∂∂xη⁡x,y=a⁢3⁢a−2⁢η⁡x,y,∂∂yη⁡x,y=η⁡x,yy,DiffConstraint=0=η⁡x,y⁢a⁢a−1,0=η⁡x,y⁢a⁢a−1,0=η⁡x,y⁢a⁢a−1,0=η⁡x,y⁢a⁢9⁢a3⁢y−9⁢a2⁢y+4⁢a⁢y−2⁢a+2
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