Changeset 4748

Timestamp:

07/22/10 14:14:00 (15 years ago)

Author:

Eric.Larour

Message:

Preliminary version of SparseToVector mex module

Location:

issm/trunk/src/mex/SparseToVector

Files:

• SparseToVector.c (modified) (3 diffs)
• SparseToVector.h (added)

issm/trunk/src/mex/SparseToVector/SparseToVector.c

@@ +1 @@
+/*\file SparseToVector.cpp
+ *\brief: extract from a sparse matlab matrix ir, jc, and other information 
+ */
 
-#define THISFUNCTION "SparseToVector"
+#include "./SparseToVector.h"
 
-#include <stdio.h>
-#include <string.h>    /*  strcasecmp  */
-#include <time.h>      /*  clock,time,difftime  */
-#include "mex.h"
-
-
-/* Input Arguments */
-
-#define    A_IN         prhs[0]
-
-/* Output Arguments */
-
-#define    IR_OUT       plhs[0]
-#define    JC_OUT       plhs[1]
-#define    PR_OUT       plhs[2]
-#define    M_OUT        plhs[3]
-#define    N_OUT        plhs[4]
-#define    NZMAX_OUT    plhs[5]
-
-
-void SparseToVectorUsage( void );
-
-
-void mexFunction( int nlhs,
-                                  mxArray *plhs[],
-                                  int nrhs,
-                                  const mxArray *prhs[] )
-{
-        mwIndex *ir =NULL,*jc =NULL;
-        double  *ird=NULL,*jcd=NULL,*pr=NULL,*prd=NULL;
+void mexFunction( int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[] ) {
+        
         int     i;
+        
+        /*Output: */
+        mwIndex *ir =NULL;
+        mwIndex *jc =NULL;
+        double  *ird=NULL;
+        double  *jcd=NULL;
+        double   *pr=NULL;
+        double   *prd=NULL;
 
         /* Check for proper number of arguments */
@@ -48 +31 @@
 
         if (!mxIsNumeric(A_IN) || !mxIsSparse(A_IN)) {
-                mexPrintf("%s -- Input matrix must be numeric and sparse.\n",THISFUNCTION);
+                mexPrintf("%s -- Input matrix must be numeric and sparse.\n",__FUNCT__);
                 mexErrMsgTxt(" ");
         }
@@ -95 +78 @@
                 }
                 mexPrintf("%s -- Input matrix is of size %d by %d with %d non-zeroes.\n",
-                                  THISFUNCTION,mxGetM(A_IN),mxGetN(A_IN),mxGetNzmax(A_IN));
+                                  __FUNCT__,mxGetM(A_IN),mxGetN(A_IN),mxGetNzmax(A_IN));
         }
-
-        return;
 }
 
-void SparseToVectorUsage( void )
-{
+void SparseToVectorUsage( void ){
 
     mexPrintf("\n");
     mexPrintf("Usage: [ir,jc,pr,m,n,nzmax]=SparseToVector(a);\n");
     mexPrintf("\n");
-
-    return;
 }