naginterfaces.library.ode.dae_​dassl_​linalg

naginterfaces.library.ode.dae_dassl_linalg(neq, ml, mu, comm)

dae_dassl_linalg is a setup function which you must call prior to dae_dassl_gen() and after a call to dae_dassl_setup(), if the Jacobian is to be considered as having a banded structure.

For full information please refer to the NAG Library document for d02np

https://www.nag.com/numeric/nl/nagdoc_30/flhtml/d02/d02npf.html

Parameters

neqint

The number of differential-algebraic equations to be solved.

mlint

$m_L$, the number of subdiagonals in the band.

muint

$m_U$, the number of superdiagonals in the band.

commdict, communication object, modified in place

Communication structure.

This argument must have been initialized by a prior call to dae_dassl_setup().

Raises

NagValueError

(errno $1$)

On entry, $neq=⟨value⟩$.

Constraint: $neq\ge 1$.

(errno $2$)

On entry, $ml=⟨value⟩$.

Constraint: $ml\ge 0$.

(errno $2$)

On entry, $ml=⟨value⟩$ and $neq=⟨value⟩$.

Constraint: $ml\le neq-1$.

(errno $3$)

On entry, $mu=⟨value⟩$.

Constraint: $mu\ge 0$.

(errno $3$)

On entry, $mu=⟨value⟩$ and $neq=⟨value⟩$.

Constraint: $mu\le neq-1$.

(errno $4$)

Either the initialization function has not been called prior to the first call of this function or the communication array has become corrupted.

(errno $5$)

On entry, $\text{licom}$ is too small: $\text{licom}=⟨value⟩$.

Notes

A call to dae_dassl_linalg specifies that the Jacobian to be used is banded in structure. If dae_dassl_linalg is not called before a call to dae_dassl_gen() then the Jacobian is assumed to be full.