chem/KPP/kpp/kpp-2.1/util/Mex_Jac_SP.c

   1 /*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   2             Matlab Gateway for the Jacobian Jac_SP
   3 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/
   4
   5 #include "mex.h"
   6 #define min( x, y ) (x) < (y) ? (x) : (y)
   7 #define max( x, y ) (x) > (y) ? (x) : (y)
   8
   9 void mexFunction( int nlhs, mxArray *plhs[],
  10                 int nrhs, const mxArray *prhs[] )
  11 {
  12  int mrows, mcols;
  13  KPP_REAL *V, *F, *RCT, *JVS;
  14
  15 /* Check for the right number and size of input arguments */
  16  if ( nrhs != 3 ) {
  17    mexErrMsgTxt("KPP_ROOT_Jac_SP requires 3 input vectors: V(KPP_NVAR), F(KPP_NFIX), RCT(KPP_NREACT)");
  18  }
  19  mrows =  mxGetM(prhs[0]); mcols = mxGetN(prhs[0]);
  20  if ( ( mrows != KPP_NVAR )||( mcols != 1 ) ) {
  21    mexPrintf("First KPP_ROOT_Jac_SP input argument is of size V(%d,%d).",
  22                mrows, mcols);
  23    mexErrMsgTxt("First KPP_ROOT_Jac_SP input argument should be a column vector V(KPP_NVAR,1)");
  24  }
  25  mrows =  mxGetM(prhs[1]); mcols = mxGetN(prhs[1]);
  26  if ( ( mrows != KPP_NFIX )||( mcols != 1 ) ) {
  27    mexPrintf("Second KPP_ROOT_Jac_SP input argument is of size F(%d,%d).",
  28                mrows, mcols);
  29    mexErrMsgTxt("Second KPP_ROOT_Jac_SP input argument should be a column vector F(KPP_NFIX,1)");
  30  }
  31  mrows =  mxGetM(prhs[2]); mcols = mxGetN(prhs[2]);
  32  if ( (  mrows != KPP_NREACT )||( mcols != 1 ) ) {
  33    mexPrintf("Third KPP_ROOT_Jac_SP input argument is of size RCT(%d,%d).",
  34                mrows, mcols);
  35    mexErrMsgTxt("Third KPP_ROOT_Jac_SP input argument should be a column vector RCT(KPP_NREACT,1)");
  36  }
  37
  38 /* Check for the right number of output arguments */
  39  if ( nlhs != 1 ) {
  40    mexErrMsgTxt("KPP_ROOT_Jac_SP requires 1 output column vector: JVS(KPP_LU_NONZERO)");
  41  }
  42
  43
  44  V   = mxGetPr(prhs[0]);
  45  F   = mxGetPr(prhs[1]);
  46  RCT = mxGetPr(prhs[2]);
  47
  48  plhs[0] = mxCreateDoubleMatrix(KPP_LU_NONZERO,1,mxREAL);
  49  JVS = mxGetPr(plhs[0]);
  50
  51  Jac_SP( V, F, RCT, JVS );
  52
  53 }