NAG CL Interface
d02npc (dae_​dassl_​linalg)

1 Purpose

d02npc is a setup function which you must call prior to d02nec and after a call to d02mwc, if the Jacobian is to be considered as having a banded structure.

2 Specification

#include <nag.h>
void  d02npc (Integer neq, Integer ml, Integer mu, Integer icom[], Integer licom, NagError *fail)
The function may be called by the names: d02npc, nag_ode_dae_dassl_linalg or nag_dae_ivp_dassl_linalg.

3 Description

A call to d02npc specifies that the Jacobian to be used is banded in structure. If d02npc is not called before a call to d02nec then the Jacobian is assumed to be full.

4 References


5 Arguments

1: neq Integer Input
On entry: the number of differential-algebraic equations to be solved.
Constraint: 1neq.
2: ml Integer Input
On entry: ml, the number of subdiagonals in the band.
Constraint: 0mlneq-1.
3: mu Integer Input
On entry: mu, the number of superdiagonals in the band.
Constraint: 0muneq-1.
4: icom[licom] Integer Communication Array
icom is used to communicate details of the integration from d02mwc and details of the banded structure of the Jacobian to d02nec.
5: licom Integer Input
On entry: the dimension of the array icom.
Constraint: licom50+neq.
6: fail NagError * Input/Output
The NAG error argument (see Section 7 in the Introduction to the NAG Library CL Interface).

6 Error Indicators and Warnings

Dynamic memory allocation failed.
See Section 3.1.2 in the Introduction to the NAG Library CL Interface for further information.
On entry, argument value had an illegal value.
Either the initialization function has not been called prior to the first call of this function or the communication array has become corrupted.
On entry, licom is too small: licom=value.
On entry, ml=value.
Constraint: ml0.
On entry, mu=value.
Constraint: mu0.
On entry, neq=value.
Constraint: neq1.
On entry, ml=value and neq=value.
Constraint: mlneq-1.
On entry, mu=value and neq=value.
Constraint: muneq-1.
An internal error has occurred in this function. Check the function call and any array sizes. If the call is correct then please contact NAG for assistance.
See Section 7.5 in the Introduction to the NAG Library CL Interface for further information.
Your licence key may have expired or may not have been installed correctly.
See Section 8 in the Introduction to the NAG Library CL Interface for further information.

7 Accuracy

Not applicable.

8 Parallelism and Performance

d02npc is not threaded in any implementation.

9 Further Comments


10 Example

See Section 10 in d02nec and d02mwc.