|
3998
|
1 CLASS = "DASRT" |
|
|
2 |
|
4044
|
3 INCLUDE = "DAERT.h" |
|
|
4 |
|
3998
|
5 OPTION |
|
|
6 NAME = "absolute tolerance" |
|
4050
|
7 DOC_ITEM |
|
|
8 Absolute tolerance. May be either vector or scalar. If a vector, it |
|
|
9 must match the dimension of the state vector, and the relative |
|
|
10 tolerance must also be a vector of the same length. |
|
|
11 END_DOC_ITEM |
|
3998
|
12 TYPE = "Array<double>" |
|
|
13 SET_ARG_TYPE = "const $TYPE&" |
|
|
14 INIT_BODY |
|
|
15 $OPTVAR.resize (1); |
|
|
16 $OPTVAR(0) = ::sqrt (DBL_EPSILON); |
|
|
17 END_INIT_BODY |
|
|
18 SET_CODE |
|
|
19 void set_$OPT (double val) |
|
|
20 { |
|
|
21 $OPTVAR.resize (1); |
|
|
22 $OPTVAR(0) = (val > 0.0) ? val : ::sqrt (DBL_EPSILON); |
|
4049
|
23 reset = true; |
|
3998
|
24 } |
|
|
25 |
|
|
26 void set_$OPT (const $TYPE& val) |
|
4049
|
27 { $OPTVAR = val; reset = true; } |
|
3998
|
28 END_SET_CODE |
|
|
29 END_OPTION |
|
|
30 |
|
|
31 OPTION |
|
|
32 NAME = "relative tolerance" |
|
4050
|
33 DOC_ITEM |
|
|
34 Relative tolerance. May be either vector or scalar. If a vector, it |
|
|
35 must match the dimension of the state vector, and the absolute |
|
|
36 tolerance must also be a vector of the same length. |
|
|
37 |
|
|
38 The local error test applied at each integration step is |
|
|
39 @example |
|
|
40 abs (local error in x(i)) <= rtol(i) * abs (Y(i)) + atol(i) |
|
|
41 @end example |
|
|
42 END_DOC_ITEM |
|
3998
|
43 TYPE = "Array<double>" |
|
|
44 SET_ARG_TYPE = "const $TYPE&" |
|
|
45 INIT_BODY |
|
|
46 $OPTVAR.resize (1); |
|
|
47 $OPTVAR(0) = ::sqrt (DBL_EPSILON); |
|
|
48 END_INIT_BODY |
|
|
49 SET_CODE |
|
|
50 void set_$OPT (double val) |
|
|
51 { |
|
|
52 $OPTVAR.resize (1); |
|
|
53 $OPTVAR(0) = (val > 0.0) ? val : ::sqrt (DBL_EPSILON); |
|
4049
|
54 reset = true; |
|
3998
|
55 } |
|
|
56 |
|
|
57 void set_$OPT (const $TYPE& val) |
|
4049
|
58 { $OPTVAR = val; reset = true; } |
|
3998
|
59 END_SET_CODE |
|
|
60 END_OPTION |
|
|
61 |
|
|
62 OPTION |
|
|
63 NAME = "initial step size" |
|
4050
|
64 DOC_ITEM |
|
|
65 Differential-algebraic problems may occaisionally suffer from severe |
|
|
66 scaling difficulties on the first step. If you know a great deal |
|
|
67 about the scaling of your problem, you can help to alleviate this |
|
|
68 problem by specifying an initial stepsize. |
|
|
69 END_DOC_ITEM |
|
3998
|
70 TYPE = "double" |
|
|
71 INIT_VALUE = "-1.0" |
|
|
72 SET_EXPR = "(val >= 0.0) ? val : -1.0" |
|
|
73 END_OPTION |
|
|
74 |
|
|
75 OPTION |
|
4049
|
76 NAME = "maximum order" |
|
4050
|
77 DOC_ITEM |
|
|
78 Restrict the maximum order of the solution method. This option must |
|
|
79 be between 1 and 5, inclusive. |
|
|
80 END_DOC_ITEM |
|
4049
|
81 TYPE = "int" |
|
|
82 INIT_VALUE = "-1" |
|
|
83 SET_EXPR = "val" |
|
|
84 END_OPTION |
|
|
85 |
|
|
86 OPTION |
|
3998
|
87 NAME = "maximum step size" |
|
4050
|
88 DOC_ITEM |
|
|
89 Setting the maximum stepsize will avoid passing over very large |
|
|
90 regions. |
|
|
91 END_DOC_ITEM |
|
3998
|
92 TYPE = "double" |
|
|
93 INIT_VALUE = "-1.0" |
|
|
94 SET_EXPR = "(val >= 0.0) ? val : -1.0" |
|
|
95 END_OPTION |
|
|
96 |
|
|
97 OPTION |
|
|
98 NAME = "step limit" |
|
4050
|
99 DOC_ITEM |
|
|
100 Maximum number of integration steps to attempt on a single call to the |
|
|
101 underlying Fortran code. |
|
|
102 END_DOC_ITEM |
|
3998
|
103 TYPE = "int" |
|
|
104 INIT_VALUE = "-1" |
|
|
105 SET_EXPR = "(val >= 0) ? val : -1" |
|
|
106 END_OPTION |
