Mercurial > octave
annotate liboctave/numeric/DASRT-opts.in @ 31249:de6fc38c78c6
Make Jacobian types offered by dlsode.f accessible by lsode (bug #31626).
* liboctave/numeric/LSODE-opts.in: Add options "jacobian type", "lower jacobian
subdiagonals", and "upper jacobian subdiagonals".
* liboctave/numeric/LSODE.cc (file scope, lsode_j,
LSODE::do_integrate (double)): Handle new configurable Jacobian types.
* build-aux/mk-opts.pl: Don't implicitly convert to integer in condition.
| author | Olaf Till <olaf.till@uni-jena.de> |
|---|---|
| date | Fri, 12 Nov 2010 08:53:05 +0100 |
| parents | 796f54d4ddbf |
| children | 597f3ee61a48 |
| rev | line source |
|---|---|
|
27923
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
1 ######################################################################## |
|
29357
bbbe4dcc7200
use the same comment style for copyright headers in .m files and shell scripts
John W. Eaton <jwe@octave.org>
parents:
27923
diff
changeset
|
2 ## |
|
30564
796f54d4ddbf
update Octave Project Developers copyright for the new year
John W. Eaton <jwe@octave.org>
parents:
30046
diff
changeset
|
3 ## Copyright (C) 2002-2022 The Octave Project Developers |
|
29357
bbbe4dcc7200
use the same comment style for copyright headers in .m files and shell scripts
John W. Eaton <jwe@octave.org>
parents:
27923
diff
changeset
|
4 ## |
|
bbbe4dcc7200
use the same comment style for copyright headers in .m files and shell scripts
John W. Eaton <jwe@octave.org>
parents:
27923
diff
changeset
|
5 ## See the file COPYRIGHT.md in the top-level directory of this |
|
bbbe4dcc7200
use the same comment style for copyright headers in .m files and shell scripts
John W. Eaton <jwe@octave.org>
parents:
27923
diff
changeset
|
6 ## distribution or <https://octave.org/copyright/>. |
|
bbbe4dcc7200
use the same comment style for copyright headers in .m files and shell scripts
John W. Eaton <jwe@octave.org>
parents:
27923
diff
changeset
|
7 ## |
|
bbbe4dcc7200
use the same comment style for copyright headers in .m files and shell scripts
John W. Eaton <jwe@octave.org>
parents:
27923
diff
changeset
|
8 ## This file is part of Octave. |
|
bbbe4dcc7200
use the same comment style for copyright headers in .m files and shell scripts
John W. Eaton <jwe@octave.org>
parents:
27923
diff
changeset
|
9 ## |
|
bbbe4dcc7200
use the same comment style for copyright headers in .m files and shell scripts
John W. Eaton <jwe@octave.org>
parents:
27923
diff
changeset
|
10 ## Octave is free software: you can redistribute it and/or modify it |
|
bbbe4dcc7200
use the same comment style for copyright headers in .m files and shell scripts
John W. Eaton <jwe@octave.org>
parents:
27923
diff
changeset
|
11 ## under the terms of the GNU General Public License as published by |
|
bbbe4dcc7200
use the same comment style for copyright headers in .m files and shell scripts
John W. Eaton <jwe@octave.org>
parents:
27923
diff
changeset
|
12 ## the Free Software Foundation, either version 3 of the License, or |
|
bbbe4dcc7200
use the same comment style for copyright headers in .m files and shell scripts
John W. Eaton <jwe@octave.org>
parents:
27923
diff
changeset
|
13 ## (at your option) any later version. |
|
bbbe4dcc7200
use the same comment style for copyright headers in .m files and shell scripts
John W. Eaton <jwe@octave.org>
parents:
27923
diff
changeset
|
14 ## |
|
bbbe4dcc7200
use the same comment style for copyright headers in .m files and shell scripts
John W. Eaton <jwe@octave.org>
parents:
27923
diff
changeset
|
15 ## Octave is distributed in the hope that it will be useful, but |
|
bbbe4dcc7200
use the same comment style for copyright headers in .m files and shell scripts
John W. Eaton <jwe@octave.org>
parents:
27923
diff
changeset
|
16 ## WITHOUT ANY WARRANTY; without even the implied warranty of |
|
bbbe4dcc7200
use the same comment style for copyright headers in .m files and shell scripts
John W. Eaton <jwe@octave.org>
parents:
27923
diff
changeset
|
17 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
|
bbbe4dcc7200
use the same comment style for copyright headers in .m files and shell scripts
John W. Eaton <jwe@octave.org>
parents:
27923
diff
changeset
|
18 ## GNU General Public License for more details. |
|
bbbe4dcc7200
use the same comment style for copyright headers in .m files and shell scripts
John W. Eaton <jwe@octave.org>
parents:
27923
diff
changeset
|
19 ## |
|
bbbe4dcc7200
use the same comment style for copyright headers in .m files and shell scripts
John W. Eaton <jwe@octave.org>
parents:
27923
diff
changeset
|
20 ## You should have received a copy of the GNU General Public License |
|
bbbe4dcc7200
use the same comment style for copyright headers in .m files and shell scripts
John W. Eaton <jwe@octave.org>
parents:
27923
diff
changeset
|
21 ## along with Octave; see the file COPYING. If not, see |
|
bbbe4dcc7200
use the same comment style for copyright headers in .m files and shell scripts
John W. Eaton <jwe@octave.org>
parents:
27923
diff
changeset
|
22 ## <https://www.gnu.org/licenses/>. |
|
bbbe4dcc7200
use the same comment style for copyright headers in .m files and shell scripts
John W. Eaton <jwe@octave.org>
parents:
27923
diff
changeset
|
23 ## |
|
27923
bd51beb6205e
update formatting of copyright notices
John W. Eaton <jwe@octave.org>
parents:
27919
diff
changeset
|
24 ######################################################################## |
| 7017 | 25 |
| 3998 | 26 CLASS = "DASRT" |
| 27 | |
| 4044 | 28 INCLUDE = "DAERT.h" |
| 29 | |
| 3998 | 30 OPTION |
| 31 NAME = "absolute tolerance" | |
| 4050 | 32 DOC_ITEM |
| 33 Absolute tolerance. May be either vector or scalar. If a vector, it | |
| 34 must match the dimension of the state vector, and the relative | |
| 35 tolerance must also be a vector of the same length. | |
| 10840 | 36 |
| 4050 | 37 END_DOC_ITEM |
| 3998 | 38 TYPE = "Array<double>" |
| 39 SET_ARG_TYPE = "const $TYPE&" | |
| 40 INIT_BODY | |
|
11574
a83bad07f7e3
attempt better backward compatibility for Array resize functions
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
41 $OPTVAR.resize (dim_vector (1, 1)); |
|
15220
61822c866ba1
use std::numeric_limits<T>::epsilon in C++ code
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
42 $OPTVAR(0) = ::sqrt (std::numeric_limits<double>::epsilon ()); |
| 3998 | 43 END_INIT_BODY |
| 44 SET_CODE | |
| 45 void set_$OPT (double val) | |
| 46 { | |
|
11574
a83bad07f7e3
attempt better backward compatibility for Array resize functions
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
47 $OPTVAR.resize (dim_vector (1, 1)); |
|
15220
61822c866ba1
use std::numeric_limits<T>::epsilon in C++ code
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
48 $OPTVAR(0) = (val > 0.0) ? val : ::sqrt (std::numeric_limits<double>::epsilon ()); |
|
30046
b3717fd85e49
maint: use "m_" prefix for member variables ODE/DAE classes in liboctave/numeric.
Rik <rik@octave.org>
parents:
29358
diff
changeset
|
49 m_reset = true; |
| 3998 | 50 } |
| 51 | |
| 52 void set_$OPT (const $TYPE& val) | |
|
30046
b3717fd85e49
maint: use "m_" prefix for member variables ODE/DAE classes in liboctave/numeric.
Rik <rik@octave.org>
parents:
29358
diff
changeset
|
53 { $OPTVAR = val; m_reset = true; } |
| 3998 | 54 END_SET_CODE |
| 55 END_OPTION | |
| 56 | |
| 57 OPTION | |
| 58 NAME = "relative tolerance" | |
| 4050 | 59 DOC_ITEM |
| 60 Relative tolerance. May be either vector or scalar. If a vector, it | |
| 61 must match the dimension of the state vector, and the absolute | |
| 62 tolerance must also be a vector of the same length. | |
| 63 | |
| 64 The local error test applied at each integration step is | |
| 10840 | 65 |
| 4050 | 66 @example |
|
9067
8970b4b10e9f
Cleanup documentation for quad.texi and diffeq.texi
Rik <rdrider0-list@yahoo.com>
parents:
7081
diff
changeset
|
67 @group |
| 7081 | 68 abs (local error in x(i)) <= ... |
| 69 rtol(i) * abs (Y(i)) + atol(i) | |
|
9067
8970b4b10e9f
Cleanup documentation for quad.texi and diffeq.texi
Rik <rdrider0-list@yahoo.com>
parents:
7081
diff
changeset
|
70 @end group |
| 4050 | 71 @end example |
| 10840 | 72 |
| 4050 | 73 END_DOC_ITEM |
| 3998 | 74 TYPE = "Array<double>" |
| 75 SET_ARG_TYPE = "const $TYPE&" | |
| 76 INIT_BODY | |
|
11574
a83bad07f7e3
attempt better backward compatibility for Array resize functions
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
77 $OPTVAR.resize (dim_vector (1, 1)); |
|
15220
61822c866ba1
use std::numeric_limits<T>::epsilon in C++ code
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
78 $OPTVAR(0) = ::sqrt (std::numeric_limits<double>::epsilon ()); |
| 3998 | 79 END_INIT_BODY |
| 80 SET_CODE | |
| 81 void set_$OPT (double val) | |
| 82 { | |
|
11574
a83bad07f7e3
attempt better backward compatibility for Array resize functions
John W. Eaton <jwe@octave.org>
parents:
11523
diff
changeset
|
83 $OPTVAR.resize (dim_vector (1, 1)); |
|
15220
61822c866ba1
use std::numeric_limits<T>::epsilon in C++ code
John W. Eaton <jwe@octave.org>
parents:
14138
diff
changeset
|
84 $OPTVAR(0) = (val > 0.0) ? val : ::sqrt (std::numeric_limits<double>::epsilon ()); |
|
30046
b3717fd85e49
maint: use "m_" prefix for member variables ODE/DAE classes in liboctave/numeric.
Rik <rik@octave.org>
parents:
29358
diff
changeset
|
85 m_reset = true; |
| 3998 | 86 } |
| 87 | |
| 88 void set_$OPT (const $TYPE& val) | |
|
30046
b3717fd85e49
maint: use "m_" prefix for member variables ODE/DAE classes in liboctave/numeric.
Rik <rik@octave.org>
parents:
29358
diff
changeset
|
89 { $OPTVAR = val; m_reset = true; } |
| 3998 | 90 END_SET_CODE |
| 91 END_OPTION | |
| 92 | |
| 93 OPTION | |
| 94 NAME = "initial step size" | |
| 4050 | 95 DOC_ITEM |
| 7001 | 96 Differential-algebraic problems may occasionally suffer from severe |
| 4050 | 97 scaling difficulties on the first step. If you know a great deal |
| 98 about the scaling of your problem, you can help to alleviate this | |
| 99 problem by specifying an initial stepsize. | |
| 10840 | 100 |
| 4050 | 101 END_DOC_ITEM |
| 3998 | 102 TYPE = "double" |
| 103 INIT_VALUE = "-1.0" | |
| 104 SET_EXPR = "(val >= 0.0) ? val : -1.0" | |
| 105 END_OPTION | |
| 106 | |
| 107 OPTION | |
| 4049 | 108 NAME = "maximum order" |
| 4050 | 109 DOC_ITEM |
| 110 Restrict the maximum order of the solution method. This option must | |
| 111 be between 1 and 5, inclusive. | |
| 10840 | 112 |
| 4050 | 113 END_DOC_ITEM |
| 5275 | 114 TYPE = "octave_idx_type" |
| 4049 | 115 INIT_VALUE = "-1" |
| 116 SET_EXPR = "val" | |
| 117 END_OPTION | |
| 118 | |
| 119 OPTION | |
| 3998 | 120 NAME = "maximum step size" |
| 4050 | 121 DOC_ITEM |
| 122 Setting the maximum stepsize will avoid passing over very large | |
| 123 regions. | |
| 10840 | 124 |
| 4050 | 125 END_DOC_ITEM |
| 3998 | 126 TYPE = "double" |
| 127 INIT_VALUE = "-1.0" | |
| 128 SET_EXPR = "(val >= 0.0) ? val : -1.0" | |
| 129 END_OPTION | |
| 130 | |
| 131 OPTION | |
| 132 NAME = "step limit" | |
| 4050 | 133 DOC_ITEM |
| 134 Maximum number of integration steps to attempt on a single call to the | |
| 135 underlying Fortran code. | |
| 136 END_DOC_ITEM | |
| 5275 | 137 TYPE = "octave_idx_type" |
| 3998 | 138 INIT_VALUE = "-1" |
| 139 SET_EXPR = "(val >= 0) ? val : -1" | |
| 140 END_OPTION |