Changeset 13378


Ignore:
Timestamp:
09/15/12 11:29:55 (13 years ago)
Author:
Mathieu Morlighem
Message:

CHG:cosmetics

Location:
issm/trunk-jpl/src/c/matlab/io
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/c/matlab/io/FetchMatlabData.cpp

    r13375 r13378  
    6161void FetchData(double** pmatrix,int* pnumel,int* pndims,int** psize,const mxArray* dataref){
    6262
    63         double*  outmatrix=NULL;
    64         int      outmatrix_numel,outmatrix_ndims;
    65         int*     outmatrix_size=NULL;
     63        int     outmatrix_numel,outmatrix_ndims;
     64        double *outmatrix       = NULL;
     65        int    *outmatrix_size  = NULL;
    6666
    6767        if(mxIsEmpty(dataref) ){
    6868                /*Nothing to pick up. Just initialize matrix pointer to NULL: */
    69                 outmatrix_numel=0;
    70                 outmatrix_ndims=0;
    71                 outmatrix_size =NULL;
    72                 outmatrix=NULL;
     69                outmatrix_numel = 0;
     70                outmatrix_ndims = 0;
     71                outmatrix_size  = NULL;
     72                outmatrix       = NULL;
    7373        }
    7474        else if( mxIsClass(dataref,"double") ||
     
    7979
    8080                /*Check dataref is not pointing to NaN: */
    81                 if ( mxIsNaN(*(mxGetPr(dataref))) && (mxGetNumberOfElements(dataref)==1) ){
    82                         outmatrix_numel=0;
    83                         outmatrix_ndims=0;
    84                         outmatrix_size =NULL;
    85                         outmatrix=NULL;
     81                if (mxIsNaN(*(mxGetPr(dataref))) && (mxGetNumberOfElements(dataref)==1)){
     82                        outmatrix_numel = 0;
     83                        outmatrix_ndims = 0;
     84                        outmatrix_size  = NULL;
     85                        outmatrix       = NULL;
    8686                }
    8787                else{
  • issm/trunk-jpl/src/c/matlab/io/MatlabNArrayToNArray.cpp

    r12595 r13378  
    2828
    2929        /*matlab indices: */
    30         mwIndex*    ir=NULL;
    31         mwIndex*    jc=NULL;
    32         double* pr=NULL;
    33         int     count;
    34         int     nnz;
    35         int     nz;
     30        mwIndex *ir    = NULL;
     31        mwIndex *jc    = NULL;
     32        double  *pr    = NULL;
     33        int      count;
    3634
    3735        /*get Matlab matrix information: */
     
    4644
    4745                /*Dealing with sparse matrix: recover size first: */
    48                 rows=mxGetM(mxmatrix);
    49                 cols=mxGetN(mxmatrix);
    50                 nnz=mxGetNzmax(mxmatrix);
    51                 nz=(int)((double)nnz/(double)rows);
     46                rows = mxGetM(mxmatrix);
     47                cols = mxGetN(mxmatrix);
    5248
    5349                matrix=xNewZeroInit<double>(rows*cols);
    5450
    5551                /*Now, get ir,jc and pr: */
    56                 ir=mxGetIr(mxmatrix);
    57                 jc=mxGetJc(mxmatrix);
    58                 pr=mxGetPr(mxmatrix);
     52                ir = mxGetIr(mxmatrix);
     53                jc = mxGetJc(mxmatrix);
     54                pr = mxGetPr(mxmatrix);
    5955
    6056                /*Now, start inserting data into double* matrix: */
     
    7975                for(i=0;i<numel;i++){
    8076                        ColumnWiseDimsFromIndex(dims,i,size,ndims);
    81                         j=IndexFromRowWiseDims(dims,size,ndims);
    82                         *(matrix+j)=*(mxmatrix_ptr+i);
     77                        j = IndexFromRowWiseDims(dims,size,ndims);
     78                        matrix[j]=(double)mxmatrix_ptr[i];
    8379                }
    8480                xDelete<int>(dims);
    85                
    8681        }
    8782
    8883        /*Assign output pointer: */
    89         *pmatrix=matrix;
    90         *pmatrix_numel=numel;
    91         *pmatrix_ndims=ndims;
    92         *pmatrix_size=size;
     84        *pmatrix       = matrix;
     85        *pmatrix_numel = numel;
     86        *pmatrix_ndims = ndims;
     87        *pmatrix_size  = size;
    9388
    9489        return 1;
     
    108103
    109104        /*matlab indices: */
    110         mwIndex*    ir=NULL;
    111         mwIndex*    jc=NULL;
    112         bool*   pm=NULL;
    113         int     count;
    114         int     nnz;
    115         int     nz;
     105        mwIndex *ir    = NULL;
     106        mwIndex *jc    = NULL;
     107        bool    *pm    = NULL;
     108        int      count;
    116109
    117110        /*get Matlab matrix information: */
    118         numel=mxGetNumberOfElements(mxmatrix);
    119         ndims=mxGetNumberOfDimensions(mxmatrix);
    120         ipt  =mxGetDimensions(mxmatrix);
    121         size =xNew<int>(ndims);
     111        numel = mxGetNumberOfElements(mxmatrix);
     112        ndims = mxGetNumberOfDimensions(mxmatrix);
     113        ipt   = mxGetDimensions(mxmatrix);
     114        size  = xNew<int>(ndims);
    122115        for (i=0;i<ndims;i++) size[i]=(int)ipt[i];
    123116
     
    128121                rows=mxGetM(mxmatrix);
    129122                cols=mxGetN(mxmatrix);
    130                 nnz=mxGetNzmax(mxmatrix);
    131                 nz=(int)((double)nnz/(double)rows);
    132 
    133123                matrix=xNewZeroInit<bool>(rows*cols);
    134124
     
    142132                for(i=0;i<cols;i++){
    143133                        for(j=0;j<(jc[i+1]-jc[i]);j++){
    144                                 *(matrix+rows*ir[count]+i)=pm[count];
     134                                matrix[rows*ir[count]+i]=pm[count];
    145135                                count++;
    146136                        }
    147137                }
    148 
    149138        }
    150139        else{
     
    159148                        ColumnWiseDimsFromIndex(dims,i,size,ndims);
    160149                        j=IndexFromRowWiseDims(dims,size,ndims);
    161                         *(matrix+j)=(bool)*(mxmatrix_ptr+i);
     150                        matrix[j]=(bool)mxmatrix_ptr[i];
    162151                }
    163152                xDelete<int>(dims);
     
    165154
    166155        /*Assign output pointer: */
    167         *pmatrix=matrix;
    168         *pmatrix_numel=numel;
    169         *pmatrix_ndims=ndims;
    170         *pmatrix_size=size;
     156        *pmatrix       = matrix;
     157        *pmatrix_numel = numel;
     158        *pmatrix_ndims = ndims;
     159        *pmatrix_size  = size;
    171160
    172161        return 1;
     
    190179        char    *pm    = NULL;
    191180        int      count;
    192         int      nnz;
    193         int      nz;
    194181
    195182        /*get Matlab matrix information: */
    196         numel=mxGetNumberOfElements(mxmatrix);
    197         ndims=mxGetNumberOfDimensions(mxmatrix);
    198         ipt  =mxGetDimensions(mxmatrix);
    199         size =xNew<int>(ndims);
     183        numel = mxGetNumberOfElements(mxmatrix);
     184        ndims = mxGetNumberOfDimensions(mxmatrix);
     185        ipt   = mxGetDimensions(mxmatrix);
     186        size  = xNew<int>(ndims);
    200187        for (i=0;i<ndims;i++) size[i]=(int)ipt[i];
    201188
     
    204191
    205192                /*Dealing with sparse matrix: recover size first: */
    206                 rows=mxGetM(mxmatrix);
    207                 cols=mxGetN(mxmatrix);
    208                 nnz=mxGetNzmax(mxmatrix);
    209                 nz=(int)((double)nnz/(double)rows);
    210 
     193                rows = mxGetM(mxmatrix);
     194                cols = mxGetN(mxmatrix);
    211195                matrix=xNew<char>(rows*cols);
    212196
    213197                /*Now, get ir,jc and pm: */
    214                 ir=mxGetIr(mxmatrix);
    215                 jc=mxGetJc(mxmatrix);
    216                 pm=(char*)mxGetData(mxmatrix);
     198                ir = mxGetIr(mxmatrix);
     199                jc = mxGetJc(mxmatrix);
     200                pm = (char*)mxGetData(mxmatrix);
    217201
    218202                /*Now, start inserting data into char* matrix: */
     
    220204                for(i=0;i<cols;i++){
    221205                        for(j=0;j<(jc[i+1]-jc[i]);j++){
    222                                 *(matrix+rows*ir[count]+i)=(char)pm[count];
     206                                matrix[rows*ir[count]+i]=(char)pm[count];
    223207                                count++;
    224208                        }
    225209                }
    226 
    227210        }
    228211        else{
    229 
    230212                /*Dealing with dense matrix: recover pointer and size: */
    231213                mxmatrix_ptr=mxGetChars(mxmatrix);
     
    261243
    262244        /*Assign output pointer: */
    263         *pmatrix=matrix;
    264         *pmatrix_numel=numel;
    265         *pmatrix_ndims=ndims;
    266         *pmatrix_size=size;
     245        *pmatrix       = matrix;
     246        *pmatrix_numel = numel;
     247        *pmatrix_ndims = ndims;
     248        *pmatrix_size  = size;
    267249
    268250        return 1;
Note: See TracChangeset for help on using the changeset viewer.