int KINSpilsAtimes(void *kinsol_mem, N_Vector v, N_Vector z)
{
KINMem kin_mem;
KINSpilsMem kinspils_mem;
int ret;
kin_mem = (KINMem) kinsol_mem;
kinspils_mem = (KINSpilsMem) lmem;
ret = jtimes(v, z, uu, &new_uu, J_data);
njtimes++;
return(ret);
}
int IDASpilsAtimes(void *ida_mem, N_Vector v, N_Vector z)
{
IDAMem IDA_mem;
IDASpilsMem idaspils_mem;
int jtflag;
IDA_mem = (IDAMem) ida_mem;
idaspils_mem = (IDASpilsMem) lmem;
jtflag = jtimes(tn, ycur, ypcur, rcur, v, z, cj, jdata, ytemp, yptemp);
njtimes++;
return(jtflag);
}
int CVSpilsAtimes(void *cvode_mem, N_Vector v, N_Vector z)
{
CVodeMem cv_mem;
CVSpilsMem cvspils_mem;
int retval;
cv_mem = (CVodeMem) cvode_mem;
cvspils_mem = (CVSpilsMem) lmem;
retval = jtimes(v, z, tn, ycur, fcur, j_data, ytemp);
njtimes++;
if (retval != 0) return(retval);
N_VLinearSum(ONE, v, -gamma, z, z);
return(0);
}
static int CVSpgmrAtimes(void *cvode_mem, N_Vector v, N_Vector z)
{
CVodeMem cv_mem;
CVSpgmrMem cvspgmr_mem;
int jtflag;
cv_mem = (CVodeMem) cvode_mem;
cvspgmr_mem = (CVSpgmrMem) lmem;
jtflag = jtimes(v, z, tn, ycur, fcur, j_data, ytemp);
njtimes++;
if (jtflag != 0) return(jtflag);
N_VLinearSum(ONE, v, -gamma, z, z);
return(0);
}
