Click or drag to resize

ESTRADRun Method

Computes the real-valued local error estimate for the selected `Radau5` linear algebra mode.

Namespace: Altaxo.Calc.Ode.Obsolete.Radau5
Assembly: AltaxoCore (in AltaxoCore.dll) Version: 4.8.3572.0 (4.8.3572.0)
Syntax
C#
Copy
public void Run(
	int N,
	double[] FJAC,
	int offset_fjac,
	int LDJAC,
	int MLJAC,
	int MUJAC,
	double[] FMAS,
	int offset_fmas,
	int LDMAS,
	int MLMAS,
	int MUMAS,
	double H,
	double DD1,
	double DD2,
	double DD3,
	IFVPOL FCN,
	ref int NFCN,
	double[] Y0,
	int offset_y0,
	double[] Y,
	int offset_y,
	int IJOB,
	double X,
	int M1,
	int M2,
	int NM1,
	double[] E1,
	int offset_e1,
	int LDE1,
	double[] Z1,
	int offset_z1,
	double[] Z2,
	int offset_z2,
	double[] Z3,
	int offset_z3,
	ref double[] CONT,
	int offset_cont,
	ref double[] F1,
	int offset_f1,
	ref double[] F2,
	int offset_f2,
	int[] IP1,
	int offset_ip1,
	int[] IPHES,
	int offset_iphes,
	double[] SCAL,
	int offset_scal,
	ref double ERR,
	bool FIRST,
	bool REJECT,
	double FAC1,
	double[] RPAR,
	int offset_rpar,
	int[] IPAR,
	int offset_ipar
)

Parameters

N  Int32
The system dimension.
FJAC  Double
The Jacobian storage array.
offset_fjac  Int32
The starting offset in FJAC.
LDJAC  Int32
The leading dimension of FJAC.
MLJAC  Int32
The lower bandwidth of the Jacobian.
MUJAC  Int32
The upper bandwidth of the Jacobian.
FMAS  Double
The mass matrix storage array.
offset_fmas  Int32
The starting offset in FMAS.
LDMAS  Int32
The leading dimension of FMAS.
MLMAS  Int32
The lower bandwidth of the mass matrix.
MUMAS  Int32
The upper bandwidth of the mass matrix.
H  Double
The step size.
DD1  Double
The first error coefficient.
DD2  Double
The second error coefficient.
DD3  Double
The third error coefficient.
FCN  IFVPOL
The differential equation evaluator.
NFCN  Int32
The function evaluation counter.
Y0  Double
The base solution vector.
offset_y0  Int32
The starting offset in Y0.
Y  Double
The current solution vector.
offset_y  Int32
The starting offset in Y.
IJOB  Int32
The linear algebra job selector.
X  Double
The current integration point.
M1  Int32
The size of the first block in the second-order formulation.
M2  Int32
The block stride used by the second-order formulation.
NM1  Int32
The reduced dimension used by the second-order formulation.
E1  Double
The real factorized system matrix.
offset_e1  Int32
The starting offset in E1.
LDE1  Int32
The leading dimension of E1.
Z1  Double
The first stage vector.
offset_z1  Int32
The starting offset in Z1.
Z2  Double
The second stage vector.
offset_z2  Int32
The starting offset in Z2.
Z3  Double
The third stage vector.
offset_z3  Int32
The starting offset in Z3.
CONT  Double
The continuation and error workspace.
offset_cont  Int32
The starting offset in CONT.
F1  Double
The first real work vector.
offset_f1  Int32
The starting offset in F1.
F2  Double
The second real work vector.
offset_f2  Int32
The starting offset in F2.
IP1  Int32
The pivot array for the real factorization.
offset_ip1  Int32
The starting offset in IP1.
IPHES  Int32
The pivot information for the Hessenberg transformation.
offset_iphes  Int32
The starting offset in IPHES.
SCAL  Double
The scaling vector used for the error norm.
offset_scal  Int32
The starting offset in SCAL.
ERR  Double
The resulting error estimate.
FIRST  Boolean
Indicates whether the current step is the first step.
REJECT  Boolean
Indicates whether the previous step was rejected.
FAC1  Double
The scalar factor used by the transformed system.
RPAR  Double
The user-provided real parameter array.
offset_rpar  Int32
The starting offset in RPAR.
IPAR  Int32
The user-provided integer parameter array.
offset_ipar  Int32
The starting offset in IPAR.
See Also