nag_dae_ivp_dassl_linalg (d02npc) (PDF version)
d02 Chapter Contents
d02 Chapter Introduction
NAG C Library Manual

NAG Library Function Document

nag_dae_ivp_dassl_linalg (d02npc)

+ Contents

    1  Purpose
    7  Accuracy
    9  Example

1  Purpose

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

2  Specification

#include <nag.h>
#include <nagd02.h>
void  nag_dae_ivp_dassl_linalg (Integer neq, Integer ml, Integer mu, Integer icom[], Integer licom, NagError *fail)

3  Description

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

4  References


5  Arguments

1:     neqIntegerInput
On entry: the number of differential-algebraic equations to be solved.
Constraint: 1neq.
2:     mlIntegerInput
On entry: ml, the number of subdiagonals in the band.
Constraint: 0mlneq-1.
3:     muIntegerInput
On entry: mu, the number of superdiagonals in the band.
Constraint: 0muneq-1.
4:     icom[licom]IntegerCommunication Array
icom is used to communicate details of the integration from nag_dae_ivp_dassl_setup (d02mwc) and details of the banded structure of the Jacobian to nag_dae_ivp_dassl_gen (d02nec).
5:     licomIntegerInput
On entry: the dimension of the array icom.
Constraint: licom50+neq.
6:     failNagError *Input/Output
The NAG error argument (see Section 3.6 in the Essential Introduction).

6  Error Indicators and Warnings

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-1=value.
Constraint: mlneq-1.
On entry, mu=value and neq-1=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.

7  Accuracy

Not applicable.

8  Further Comments


9  Example

See Section 9 in nag_dae_ivp_dassl_gen (d02nec) and nag_dae_ivp_dassl_setup (d02mwc).

nag_dae_ivp_dassl_linalg (d02npc) (PDF version)
d02 Chapter Contents
d02 Chapter Introduction
NAG C Library Manual

© The Numerical Algorithms Group Ltd, Oxford, UK. 2012