Changeset 16768

Timestamp:

11/14/13 16:14:35 (11 years ago)

Author:

Eric.Larour

Message:

CHG: new fetch multiple data for integers.

Location:

issm/trunk-jpl/src/c/classes

Files:

• IoModel.cpp (modified) (1 diff)
• IoModel.h (modified) (1 diff)

issm/trunk-jpl/src/c/classes/IoModel.cpp

@@ -1300 +1300 @@
 }
 /*}}}*/
+/*FUNCTION IoModel::FetchMultipleData(int*** pmatrices,int** pmdims,int** pndims, int* pM,int data_enum){{{*/
+void  IoModel::FetchMultipleData(int*** pmatrices,int** pmdims,int** pndims, int* pnumrecords,int data_enum){
+
+        int     num_instances;
+        fpos_t* file_positions=NULL;
+
+        /*output: */
+        int        **matrices = NULL;
+        int         *mdims    = NULL;
+        int         *ndims    = NULL;
+
+        /*intermediary: */
+        int          M, N;
+        IssmPDouble *pmatrix = NULL;
+        IssmDouble  *matrix  = NULL;
+        int         *integer_matrix=NULL;
+        int         *codes   = NULL;
+        int          code;
+
+        /*recover my_rank:*/
+        int my_rank=IssmComm::GetRank();
+        
+        /*Get file pointers to beginning of the data (multiple instances of it): */
+        file_positions=this->SetFilePointersToData(&codes,NULL,&num_instances,data_enum);
+
+        if(num_instances){
+
+                /*Allocate matrices :*/
+                matrices=xNew<int*>(num_instances);
+                mdims=xNew<int>(num_instances);
+                ndims=xNew<int>(num_instances);
+
+                for(int i=0;i<num_instances;i++){
+
+                        if(my_rank==0){
+                                code=codes[i];
+
+                                if((code!=5) && (code!=6) && (code!=7))_error_("expecting a IssmDouble, integer or boolean matrix for enum " << EnumToStringx(data_enum));
+
+                                /*We have to read a matrix from disk. First read the dimensions of the matrix, then the whole matrix: */
+                                /*numberofelements: */
+                                fsetpos(fid,file_positions+i);
+                                if(fread(&M,sizeof(int),1,fid)!=1) _error_("could not read number of rows for matrix ");
+                        }
+                        ISSM_MPI_Bcast(&M,1,ISSM_MPI_INT,0,IssmComm::GetComm()); 
+
+                        if(my_rank==0){  
+                                if(fread(&N,sizeof(int),1,fid)!=1) _error_("could not read number of columns for matrix ");
+                        }
+                        ISSM_MPI_Bcast(&N,1,ISSM_MPI_INT,0,IssmComm::GetComm()); 
+
+
+                        /*Now allocate matrix: */
+                        if(M*N){
+                                pmatrix=xNew<IssmPDouble>(M*N);
+                                integer_matrix=xNew<int>(M*N);
+
+                                /*Read matrix on node 0, then broadcast: */
+                                if(my_rank==0){  
+                                        if(fread(pmatrix,M*N*sizeof(IssmPDouble),1,fid)!=1) _error_("could not read matrix ");
+                                }
+                                ISSM_MPI_Bcast(pmatrix,M*N,ISSM_MPI_PDOUBLE,0,IssmComm::GetComm()); 
+
+                                _assert_(this->independents);
+                                if(this->independents[data_enum]){
+                                        /*this data has already been checked out! So cancel all that we've done here, and return 
+                                         * the data[data_enum] directly: */
+                                        matrix=this->data[data_enum];
+                                        for (int i=0;i<M*N;++i) integer_matrixmatrix[i]=reCast<int>(matrix[i]);
+                                }
+                                else{
+                                        for (int i=0;i<M*N;++i) integer_matrix[i]=pmatrix[i];
+                                }
+                                xDelete<IssmPDouble>(pmatrix);
+                        }
+                        else
+                                integer_matrix=NULL;
+                        
+                        
+                        /*Assign: */
+                        mdims[i]=M;
+                        matrices[i]=integer_matrix;
+                        ndims[i]=N;
+                }
+        }
+                        
+        /*Free ressources:*/
+        xDelete<fpos_t>(file_positions);
+        xDelete<int>(codes);
+
+        /*Assign output pointers: */
+        *pmatrices=matrices;
+        if(pmdims){
+                *pmdims=mdims;
+        }
+        else{
+                xDelete<int>(mdims);
+        }
+        if(pndims){
+                *pndims=ndims;
+        }
+        else{
+                xDelete<int>(ndims);
+        }
+        *pnumrecords=num_instances;
+}
+/*}}}*/
 /*FUNCTION IoModel::FetchDataToInput(Elements* elements,int vector_enum,IssmDouble default_value){{{*/
 void IoModel::FetchDataToInput(Elements* elements,int vector_enum,IssmDouble default_value){

issm/trunk-jpl/src/c/classes/IoModel.h

@@ -79 +79 @@
                 void        FetchMultipleData(char***   pstringarray,int* pnumstrings,int data_enum);
                 void        FetchMultipleData(IssmDouble*** pmatrixarray,int** pmdims,int** pndims, int* pnumrecords,int data_enum);
+                void        FetchMultipleData(int*** pmatrices,int** pmdims,int** pndims, int* pnumrecords,int data_enum);
                 void        FetchData(Option **poption,int data_enum);
                 void        FetchData(int num,...);