Prepares or updates the linear system matrix used in the DVODE corrector iteration.
Namespace: Altaxo.Calc.Ode.Obsolete.DVodeAssembly: AltaxoCore (in AltaxoCore.dll) Version: 4.8.3572.0 (4.8.3572.0)
Syntaxpublic void Run(
ref double[] Y,
int offset_y,
double[] YH,
int offset_yh,
int LDYH,
double[] EWT,
int offset_ewt,
ref double[] FTEM,
int offset_ftem,
double[] SAVF,
int offset_savf,
ref double[] WM,
int offset_wm,
ref int[] IWM,
int offset_iwm,
IFEX F,
IJEX JAC,
ref int IERPJ,
double[] RPAR,
int offset_rpar,
int[] IPAR,
int offset_ipar
)
Parameters
- Y Double
-
= Vector containing predicted values on entry.
- offset_y Int32
- The starting offset in Y.
- YH Double
-
= The Nordsieck array, an LDYH by LMAX array, input.
- offset_yh Int32
- The starting offset in YH.
- LDYH Int32
-
= A constant .ge. N, the first dimension of YH, input.
- EWT Double
-
= An error weight vector of length N.
- offset_ewt Int32
- The starting offset in EWT.
- FTEM Double
- Temporary storage for function evaluations.
- offset_ftem Int32
- The starting offset in FTEM.
- SAVF Double
-
= Array containing f evaluated at predicted y, input.
- offset_savf Int32
- The starting offset in SAVF.
- WM Double
-
= Real work space for matrices. In the output, it containS
the inverse diagonal matrix if MITER = 3 and the LU
decomposition of P if MITER is 1, 2 , 4, or 5.
Storage of matrix elements starts at WM(3).
Storage of the saved Jacobian starts at WM(LOCJS).
WM also contains the following matrix-related data:
WM(1) = SQRT(UROUND), used in numerical Jacobian step.
WM(2) = H*RL1, saved for later use if MITER = 3.
- offset_wm Int32
- The starting offset in WM.
- IWM Int32
-
= Integer work space containing pivot information,
starting at IWM(31), if MITER is 1, 2, 4, or 5.
IWM also contains band parameters ML = IWM(1) and
MU = IWM(2) if MITER is 4 or 5.
- offset_iwm Int32
- The starting offset in IWM.
- F IFEX
-
= Dummy name for the user supplied subroutine for f.
- JAC IJEX
-
= Dummy name for the user supplied Jacobian subroutine.
- IERPJ Int32
-
= Output error flag, = 0 if no trouble, 1 if the P
matrix is found to be singular.
- RPAR Double
- User-defined real parameter array passed to callback routines.
- offset_rpar Int32
- The starting offset in RPAR.
- IPAR Int32
- User-defined integer parameter array passed to callback routines.
- offset_ipar Int32
- The starting offset in IPAR.
See Also