|
3998
|
1 CLASS = "ODESSA" |
|
|
2 |
|
4044
|
3 INCLUDE = "ODES.h" |
|
|
4 |
|
3998
|
5 OPTION |
|
|
6 NAME = "absolute tolerance" |
|
4051
|
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. |
|
|
10 END_DOC_ITEM |
|
3998
|
11 TYPE = "Array<double>" |
|
|
12 SET_ARG_TYPE = "const $TYPE&" |
|
|
13 INIT_BODY |
|
|
14 $OPTVAR.resize (1); |
|
|
15 $OPTVAR(0) = ::sqrt (DBL_EPSILON); |
|
|
16 END_INIT_BODY |
|
|
17 SET_CODE |
|
|
18 void set_$OPT (double val) |
|
|
19 { |
|
|
20 $OPTVAR.resize (1); |
|
|
21 $OPTVAR(0) = (val > 0.0) ? val : ::sqrt (DBL_EPSILON); |
|
4049
|
22 reset = true; |
|
3998
|
23 } |
|
|
24 |
|
|
25 void set_$OPT (const $TYPE& val) |
|
4049
|
26 { $OPTVAR = val; reset = true; } |
|
3998
|
27 END_SET_CODE |
|
|
28 END_OPTION |
|
|
29 |
|
|
30 OPTION |
|
4051
|
31 NAME = "relative tolerance" |
|
|
32 DOC_ITEM |
|
|
33 Relative tolerance parameter. Unlike the absolute tolerance, this |
|
|
34 parameter may only be a scalar. |
|
|
35 |
|
|
36 The local error test applied at each integration step is |
|
|
37 |
|
|
38 @example |
|
5016
|
39 abs (local error in x(i)) |
|
|
40 <= rtol * abs (y(i)) + atol(i) |
|
4051
|
41 @end example |
|
|
42 END_DOC_ITEM |
|
4231
|
43 TYPE = "double" |
|
|
44 INIT_VALUE = "::sqrt (DBL_EPSILON)" |
|
|
45 SET_EXPR = "(val > 0.0) ? val : ::sqrt (DBL_EPSILON)" |
|
|
46 END_OPTION |
|
|
47 |
|
|
48 OPTION |
|
|
49 NAME = "integration method" |
|
|
50 DOC_ITEM |
|
|
51 A string specifing the method of integration to use to solve the ODE |
|
|
52 system. Valid values are |
|
|
53 |
|
|
54 @table @asis |
|
|
55 @item \"adams\" |
|
|
56 @itemx \"non-stiff\" |
|
|
57 No Jacobian used (even if it is available). |
|
|
58 @item \"bdf\" |
|
|
59 @item \"stiff\" |
|
|
60 Use stiff backward differentiation formula (BDF) method. If a |
|
|
61 function to compute the Jacobian is not supplied, @code{lsode} will |
|
|
62 compute a finite difference approximation of the Jacobian matrix. |
|
|
63 @end table |
|
|
64 END_DOC_ITEM |
|
3998
|
65 TYPE = "std::string" |
|
|
66 SET_ARG_TYPE = "const $TYPE&" |
|
|
67 INIT_VALUE = {"stiff"} |
|
|
68 SET_BODY |
|
|
69 if (val == "stiff" || val == "bdf") |
|
|
70 $OPTVAR = "stiff"; |
|
|
71 else if (val == "non-stiff" || val == "adams") |
|
|
72 $OPTVAR = "non-stiff"; |
|
|
73 else |
|
|
74 (*current_liboctave_error_handler) |
|
|
75 ("lsode_options: method must be \"stiff\", \"bdf\", \"non-stiff\", or \"adams\""); |
|
|
76 END_SET_BODY |
|
|
77 END_OPTION |
|
|
78 |
|
|
79 OPTION |
|
|
80 NAME = "initial step size" |
|
4051
|
81 DOC_ITEM |
|
|
82 The step size to be attempted on the first step (default is determined |
|
|
83 automatically). |
|
|
84 END_DOC_ITEM |
|
3998
|
85 TYPE = "double" |
|
|
86 INIT_VALUE = "-1.0" |
|
|
87 SET_EXPR = "(val >= 0.0) ? val : -1.0" |
|
|
88 END_OPTION |
|
|
89 |
|
|
90 OPTION |
|
4231
|
91 NAME = "maximum order" |
|
|
92 DOC_ITEM |
|
|
93 Restrict the maximum order of the solution method. If using the Adams |
|
|
94 method, this option must be between 1 and 12. Otherwise, it must be |
|
|
95 between 1 and 5, inclusive. |
|
|
96 END_DOC_ITEM |
|
|
97 TYPE = "int" |
|
|
98 INIT_VALUE = "-1" |
|
|
99 SET_EXPR = "val" |
|
|
100 END_OPTION |
|
|
101 |
|
|
102 OPTION |
|
3998
|
103 NAME = "maximum step size" |
|
4051
|
104 DOC_ITEM |
|
|
105 Setting the maximum stepsize will avoid passing over very large |
|
|
106 regions (default is not specified). |
|
|
107 END_DOC_ITEM |
|
3998
|
108 TYPE = "double" |
|
|
109 INIT_VALUE = "-1.0" |
|
|
110 SET_EXPR = "(val >= 0.0) ? val : -1.0" |
|
|
111 END_OPTION |
|
|
112 |
|
|
113 OPTION |
|
|
114 NAME = "minimum step size" |
|
4051
|
115 DOC_ITEM |
|
|
116 The minimum absolute step size allowed (default is 0). |
|
|
117 END_DOC_ITEM |
|
3998
|
118 TYPE = "double" |
|
|
119 INIT_VALUE = "0.0" |
|
|
120 SET_EXPR = "(val >= 0.0) ? val : 0.0" |
|
|
121 END_OPTION |
|
|
122 |
|
|
123 OPTION |
|
|
124 NAME = "step limit" |
|
4051
|
125 DOC_ITEM |
|
|
126 Maximum number of steps allowed (default is 100000). |
|
|
127 END_DOC_ITEM |
|
3998
|
128 TYPE = "int" |
|
|
129 INIT_VALUE = "100000" |
|
|
130 SET_EXPR = "val" |
|
|
131 END_OPTION |
|
|
132 |
