Changeset 17145

Timestamp:

01/20/14 20:46:03 (11 years ago)

Author:

Eric.Larour

Message:

CHG: missing instantiation.

File:

• issm/trunk-jpl/src/c/classes/ExternalResults/GenericExternalResult.h (modified) (1 diff)

issm/trunk-jpl/src/c/classes/ExternalResults/GenericExternalResult.h

@@ -368 +368 @@
 } /*}}}*/
 
+/*Specific instantiations for IssmDouble*: */
+#if defined(_HAVE_ADOLC_) && !defined(_WRAPPERS_)  //We hook off this specific specialization when not running ADOLC, otherwise we get a redeclaration with the next specialization. 
+template <> inline void GenericExternalResult<IssmDouble*>::WriteData(FILE* fid,bool io_gather){ /*{{{*/
+
+        int     i;
+        int     my_rank;
+        int     type;
+        int     rows,cols;
+        char   *name    = NULL;
+        IssmPDouble passiveDouble;
+        IssmPDouble* passiveDoubles;
+
+        /*recover my_rank:*/
+        my_rank=IssmComm::GetRank();
+
+        if(io_gather){
+                /*we are gathering the data on cpu 0, don't write on other cpus: */
+                if(my_rank) return;
+        }
+
+        /*First write enum: */
+        int length=(strlen(this->result_name)+1)*sizeof(char);
+        fwrite(&length,sizeof(int),1,fid);
+        fwrite(this->result_name,length,1,fid);
+
+        /*Now write time and step: */
+        passiveDouble=reCast<IssmPDouble>(time);
+        fwrite(&passiveDouble,sizeof(IssmPDouble),1,fid);
+        fwrite(&step,sizeof(int),1,fid);
+
+        /*writing a IssmDouble array, type is 3:*/
+        type=3;
+        fwrite(&type,sizeof(int),1,fid);
+        rows=this->M;
+        fwrite(&rows,sizeof(int),1,fid);
+        cols=this->N;
+        fwrite(&cols,sizeof(int),1,fid);
+
+        passiveDoubles=xNew<IssmPDouble>(this->M*this->N);
+        for (i=0;i<this->M*this->N;i++)passiveDoubles[i]=reCast<IssmPDouble>(value[i]);
+        fwrite(passiveDoubles,cols*rows*sizeof(IssmPDouble),1,fid);
+        xDelete<IssmPDouble>(passiveDoubles);
+
+}
+/*}}}*/
+#endif
+
 /*Specifics instantiations for Vector*/
 template <> inline GenericExternalResult<Vector<IssmPDouble>*>::GenericExternalResult(int in_id, int in_enum_type,Vector<IssmPDouble>* in_values, int in_M,int in_N,int in_step,IssmDouble in_time){/*{{{*/