void mexFunction( int nlhs, mxArray *plhs[],
int nrhs, const mxArray *prhs[] )
{
int mrows, mcols;
double *V, *F, *RCT, *JVS;
/* Check for the right number and size of input arguments */
if ( nrhs != 3 ) {
mexErrMsgTxt("saprc99_Jac_SP requires 3 input vectors: V(74), F(5), RCT(211)");
}
mrows = mxGetM(prhs[0]); mcols = mxGetN(prhs[0]);
if ( ( mrows != 74 )||( mcols != 1 ) ) {
mexPrintf("First saprc99_Jac_SP input argument is of size V(%d,%d).",
mrows, mcols);
mexErrMsgTxt("First saprc99_Jac_SP input argument should be a column vector V(74,1)");
}
mrows = mxGetM(prhs[1]); mcols = mxGetN(prhs[1]);
if ( ( mrows != 5 )||( mcols != 1 ) ) {
mexPrintf("Second saprc99_Jac_SP input argument is of size F(%d,%d).",
mrows, mcols);
mexErrMsgTxt("Second saprc99_Jac_SP input argument should be a column vector F(5,1)");
}
mrows = mxGetM(prhs[2]); mcols = mxGetN(prhs[2]);
if ( ( mrows != 211 )||( mcols != 1 ) ) {
mexPrintf("Third saprc99_Jac_SP input argument is of size RCT(%d,%d).",
mrows, mcols);
mexErrMsgTxt("Third saprc99_Jac_SP input argument should be a column vector RCT(211,1)");
}
/* Check for the right number of output arguments */
if ( nlhs != 1 ) {
mexErrMsgTxt("saprc99_Jac_SP requires 1 output column vector: JVS(920)");
}
V = mxGetPr(prhs[0]);
F = mxGetPr(prhs[1]);
RCT = mxGetPr(prhs[2]);
plhs[0] = mxCreateDoubleMatrix(920,1,mxREAL);
JVS = mxGetPr(plhs[0]);
Jac_SP( V, F, RCT, JVS );
}
/*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/
void JacTemplate( KPP_REAL T, KPP_REAL Y[], KPP_REAL Jcb[] )
/*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Template for the ODE Jacobian call.
Updates the rate coefficients (and possibly the fixed species) at each call
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/
{
/*~~~> Local variables */
KPP_REAL Told;
Told = TIME;
TIME = T ;
Update_SUN();
Update_RCONST();
Jac_SP( Y, FIX, RCONST, Jcb );
TIME = Told;
Njac++;
} /* JacTemplate */
