Procedures

this routine copies one rectangular array, a, to another, b, where a and b may have different row dimensions, nrowa and nrowb. the data copied consists of nrow rows and ncol columns.

constant times a vector plus a vector. uses unrolled loops for increments equal to one.

copies a vector, x, to a vector, y. uses unrolled loops for increments equal to one.

forms the dot product of two vectors. uses unrolled loops for increments equal to one.

Set error weight vector.

factor a band matrix using gaussian elimination.

solve the real band system a*x=b or trans(a)*x=b using the factors computed by dgbco or dgbfa.

factor a matrix using gaussian elimination.

solve the real system a*x=b or trans(a)*x=b using the factors computed by dgeco or dgefa.

Function that returns the Euclidean norm $\sqrt{ \mathbf{x}^T \mathbf{x} }$ of a vector $\mathbf{x}$.

scales a vector by a constant. uses unrolled loops for increment equal to one.

this routine computes the step size, h0, to be attempted on the first step, when the user has not supplied a value for this.

dvindy computes interpolated values of the k-th derivative of the dependent variable vector y, and stores it in dky. this routine is called within the package with k = 0 and t = tout, but may also be called by the user for any k up to the current order. (see detailed instructions in the usage documentation.)

dvjac is called by dvnlsd to compute and process the matrix p = i - h*rl1*j , where j is an approximation to the jacobian.

this subroutine adjusts the yh array on reduction of order, and also when the order is increased for the stiff option (meth = 2).

subroutine dvnlsd is a nonlinear system solver, which uses functional iteration or a chord (modified newton) method. for the chord method direct linear algebraic system solvers are used. subroutine dvnlsd then handles the corrector phase of this integration package.

weighted root-mean-square vector norm.

dvode: variable-coefficient ordinary differential equation solver, with fixed-leading-coefficient implementation.

dvset is called by dvstep and sets coefficients for use there.

This routine manages the solution of the linear system arising from a chord iteration. it is called if miter /= 0:

this routine saves or restores (depending on job) the contents of the dvode internal variables.

dvstep performs one step of the integration of an initial value problem for a system of ordinary differential equations.

Find the smallest index of that component of a vector having the maximum magnitude.

Set the function pointers. This must be called before dvode is called.

save and recall error message control parameters.

write error message with values.

reset the error print control flag.

reset the logical unit number for error messages.