/[escript]/trunk/finley/src/Mesh_saveDX.c
ViewVC logotype

Diff of /trunk/finley/src/Mesh_saveDX.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1811 by ksteube, Thu Sep 25 23:11:13 2008 UTC revision 2271 by jfenwick, Mon Feb 16 05:08:29 2009 UTC
# Line 37  void Finley_Mesh_saveDX(const char * fil Line 37  void Finley_Mesh_saveDX(const char * fil
37    FILE * fileHandle_p = NULL;    FILE * fileHandle_p = NULL;
38    int i,j,k,i_data, elementtype, numPoints = 0, nDim, *resortIndex=NULL, p,    int i,j,k,i_data, elementtype, numPoints = 0, nDim, *resortIndex=NULL, p,
39        numDXNodesPerElement=0, numCells, NN, object_count, rank, nComp, numPointsPerSample;        numDXNodesPerElement=0, numCells, NN, object_count, rank, nComp, numPointsPerSample;
40    double* values,rtmp;    __const double* values;
41      double rtmp;
42    bool_t *isCellCentered=NULL;    bool_t *isCellCentered=NULL;
43    Finley_ElementFile* elements=NULL;    Finley_ElementFile* elements=NULL;
44    ElementTypeId TypeId;    ElementTypeId TypeId;
# Line 237  void Finley_Mesh_saveDX(const char * fil Line 238  void Finley_Mesh_saveDX(const char * fil
238                  numPointsPerSample=elements->ReferenceElement->numQuadNodes;                  numPointsPerSample=elements->ReferenceElement->numQuadNodes;
239               }               }
240               if (numPointsPerSample>0) {               if (numPointsPerSample>0) {
241            void* buffer=allocSampleBuffer(data_pp[i_data]);
242                  fprintf(fileHandle_p, "items %d data follows\n", numCells);                  fprintf(fileHandle_p, "items %d data follows\n", numCells);
243                  for (i=0;i<elements->numElements;i++) {                  for (i=0;i<elements->numElements;i++) {
244                      values=getSampleData(data_pp[i_data],i);                      values=getSampleDataRO(data_pp[i_data],i,buffer);
245                      for (k=0;k<nComp;k++) {                      for (k=0;k<nComp;k++) {
246                          if ( isExpanded(data_pp[i_data]) ) {                          if ( isExpanded(data_pp[i_data]) ) {
247                              rtmp=0.;                              rtmp=0.;
# Line 251  void Finley_Mesh_saveDX(const char * fil Line 253  void Finley_Mesh_saveDX(const char * fil
253                      }                      }
254                  fprintf(fileHandle_p, "\n");                  fprintf(fileHandle_p, "\n");
255                  }                  }
256            freeSampleBuffer(buffer);
257                  fprintf(fileHandle_p, "attribute \"dep\" string \"connections\"\n");                  fprintf(fileHandle_p, "attribute \"dep\" string \"connections\"\n");
258              }              }
259           } else {           } else {
260             void* buffer=allocSampleBuffer(data_pp[i_data]);
261               fprintf(fileHandle_p, "items %d data follows\n", numPoints);               fprintf(fileHandle_p, "items %d data follows\n", numPoints);
262               for (i=0;i<numPoints;i++) {               for (i=0;i<numPoints;i++) {
263                     values=getSampleData(data_pp[i_data],i);                     values=getSampleDataRO(data_pp[i_data],i,buffer);
264                     for (k=0;k<nComp;k++) fprintf(fileHandle_p, " %g", values[k]);                     for (k=0;k<nComp;k++) fprintf(fileHandle_p, " %g", values[k]);
265                 fprintf(fileHandle_p, "\n");                 fprintf(fileHandle_p, "\n");
266               }               }
267             freeSampleBuffer(buffer);
268               fprintf(fileHandle_p, "attribute \"dep\" string \"positions\"\n");               fprintf(fileHandle_p, "attribute \"dep\" string \"positions\"\n");
269           }           }
270       }       }

Legend:
Removed from v.1811  
changed lines
  Added in v.2271

  ViewVC Help
Powered by ViewVC 1.1.26