static int CVDiagSolve(CVodeMem cv_mem, N_Vector b, N_Vector weight,
N_Vector ycur, N_Vector fcur)
{
booleantype invOK;
realtype r;
CVDiagMem cvdiag_mem;
cvdiag_mem = (CVDiagMem) lmem;
/* If gamma has changed, update factor in M, and save gamma value */
if (gammasv != gamma) {
r = gamma / gammasv;
N_VInv(M, M);
N_VAddConst(M, -ONE, M);
N_VScale(r, M, M);
N_VAddConst(M, ONE, M);
invOK = N_VInvTest(M, M);
if (!invOK) {
last_flag = CVDIAG_INV_FAIL;
return (1);
}
gammasv = gamma;
}
/* Apply M-inverse to b */
N_VProd(b, M, b);
last_flag = CVDIAG_SUCCESS;
return(0);
}
static real KINScSteplength(KINMem kin_mem, N_Vector ucur,
N_Vector ss, N_Vector usc)
{
N_VInv(usc, vtemp1);
N_VAbs(ucur, vtemp2);
N_VLinearSum(ONE, vtemp1, ONE, vtemp2, vtemp1);
N_VDiv(ss, vtemp1, vtemp1);
return(N_VMaxNorm(vtemp1));
}
void N_VInv_SensWrapper(N_Vector x, N_Vector z)
{
int i;
for (i=0; i < NV_NVECS_SW(x); i++)
N_VInv(NV_VEC_SW(x,i), NV_VEC_SW(z,i));
return;
}
