Changeset 18806

Timestamp:

11/18/14 14:03:44 (10 years ago)

Author:

Eric.Larour

Message:

CHG: adding fetch multiple data for double vectors.

File:

• issm/trunk-jpl/src/c/classes/IoModel.cpp (modified) (1 diff)

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

@@ -1229 +1229 @@
 }
 /*}}}*/
+void  IoModel::FetchMultipleData(IssmDouble** pvector, int* pnum_instances,int data_enum){/*{{{*/
+
+        int     num_instances;
+        fpos_t* file_positions=NULL;
+
+        /*output: */
+        IssmDouble* vector=NULL;
+
+        /*intermediary: */
+        IssmPDouble          scalar;
+        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<IssmDouble>(num_instances);
+
+                for(int i=0;i<num_instances;i++){
+
+                        if(my_rank==0){
+                                code=codes[i];
+
+                                if(code!=3)_error_("expecting a double for enum " << EnumToStringx(data_enum));
+                                
+                                /*We have to read a double from disk: */
+                                fsetpos(fid,file_positions+i);
+                                if(my_rank==0){  
+                                        if(fread(&scalar,sizeof(IssmPDouble),1,fid)!=1) _error_("could not read scalar ");
+                                }
+                        }
+                        ISSM_MPI_Bcast(&scalar,1,ISSM_MPI_PDOUBLE,0,IssmComm::GetComm()); 
+
+                        /*Assign: */
+                        vector[i]=scalar;
+                }
+        }
+                        
+        /*Free ressources:*/
+        xDelete<fpos_t>(file_positions);
+        xDelete<int>(codes);
+
+        /*Assign output pointers: */
+        *pvector=vector;
+        *pnum_instances=num_instances;
+}
+/*}}}*/
 void  IoModel::FetchMultipleData(IssmDouble*** pmatrices,int** pmdims,int** pndims, int* pnumrecords,int data_enum){/*{{{*/