Changeset 16770

Timestamp:

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

Author:

Eric.Larour

Message:

CHG: added multiple integer vector fetch.

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

@@ -1195 +1195 @@
 }
 /*}}}*/
+/*FUNCTION IoModel::FetchMultipleData(int** pvector, int* pM,int data_enum){{{*/
+void  IoModel::FetchMultipleData(int** pvector, int* pnum_instances,int data_enum){
+
+        int     num_instances;
+        fpos_t* file_positions=NULL;
+
+        /*output: */
+        int* vector=NULL;
+        int num
+
+        /*intermediary: */
+        int          M, N;
+        IssmPDouble *pmatrix = NULL;
+        IssmDouble  *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 vector :*/
+                vector=xNew<int>(num_instances);
+
+                for(int i=0;i<num_instances;i++){
+
+                        if(my_rank==0){
+                                code=codes[i];
+
+                                if(code!=2)_error_("expecting an integer for enum " << EnumToStringx(data_enum));
+                                
+                                /*We have to read a integer from disk. First read the dimensions of the integer, then the integer: */
+                                fsetpos(fid,file_positions+i);
+                                if(my_rank==0){  
+                                        if(fread(&integer,sizeof(int),1,fid)!=1) _error_("could not read integer ");
+                                }
+                        }
+                        ISSM_MPI_Bcast(&integer,1,ISSM_MPI_INT,0,IssmComm::GetComm()); 
+
+                        /*Assign: */
+                        vector[i]=integer;
+                }
+        }
+                        
+        /*Free ressources:*/
+        xDelete<fpos_t>(file_positions);
+        xDelete<int>(codes);
+
+        /*Assign output pointers: */
+        *pvector=vector;
+        *pnum_instances=num_instances;
+}
+/*}}}*/
 /*FUNCTION IoModel::FetchMultipleData(IssmDouble*** pmatrices,int** pmdims,int** pndims, int* pM,int data_enum){{{*/
 void  IoModel::FetchMultipleData(IssmDouble*** pmatrices,int** pmdims,int** pndims, int* pnumrecords,int data_enum){

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

@@ -80 +80 @@
                 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        FetchMultipleData(int** pvector, int* pnum_instances,int data_enum);
                 void        FetchData(Option **poption,int data_enum);
                 void        FetchData(int num,...);