1 /* -------------------------------------------------------------------------- */
2 /* ssmult mexFunction */
3 /* -------------------------------------------------------------------------- */
4
5 #include "ssmult.h"
6
7 /* C = A*B and variants. Both A and B must be sparse */
8
mexFunction(int nargout,mxArray * pargout[],int nargin,const mxArray * pargin[])9 void mexFunction
10 (
11 int nargout,
12 mxArray *pargout [ ],
13 int nargin,
14 const mxArray *pargin [ ]
15 )
16 {
17 const mxArray *A, *B ;
18 int at, ac, bt, bc, ct, cc ;
19
20 /* ---------------------------------------------------------------------- */
21 /* get inputs */
22 /* ---------------------------------------------------------------------- */
23
24 if (nargin < 2 || nargin > 8 || nargout > 1)
25 {
26 mexErrMsgTxt ("Usage: C = ssmult (A,B, at,ac, bt,bc, ct,cc)") ;
27 }
28
29 A = pargin [0] ;
30 B = pargin [1] ;
31
32 at = (nargin > 2) ? mxGetScalar (pargin [2]) : 0 ;
33 ac = (nargin > 3) ? mxGetScalar (pargin [3]) : 0 ;
34 bt = (nargin > 4) ? mxGetScalar (pargin [4]) : 0 ;
35 bc = (nargin > 5) ? mxGetScalar (pargin [5]) : 0 ;
36 ct = (nargin > 6) ? mxGetScalar (pargin [6]) : 0 ;
37 cc = (nargin > 7) ? mxGetScalar (pargin [7]) : 0 ;
38
39 /* ---------------------------------------------------------------------- */
40 /* C = A*B or variants */
41 /* ---------------------------------------------------------------------- */
42
43 pargout [0] = ssmult (A, B, at, ac, bt, bc, ct, cc) ;
44 }
45