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

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

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

revision 2754 by gross, Wed Nov 18 07:44:26 2009 UTC revision 2770 by jfenwick, Wed Nov 25 01:24:51 2009 UTC
# Line 803  void Finley_Mesh_saveVTK(const char *fil Line 803  void Finley_Mesh_saveVTK(const char *fil
803                  txtBufferInUse = 0;                  txtBufferInUse = 0;
804                  for (i=0; i<numCells; i++) {                  for (i=0; i<numCells; i++) {
805                      if (elements->Owner[i] == my_mpi_rank) {                      if (elements->Owner[i] == my_mpi_rank) {
806                          void* sampleBuffer=allocSampleBuffer(data_pp[dataIdx]);                          __const double *values = getSampleDataRO(data_pp[dataIdx], i);
                         __const double *values = getSampleDataRO(data_pp[dataIdx], i,sampleBuffer);  
807                          for (l = 0; l < numCellFactor; l++) {                          for (l = 0; l < numCellFactor; l++) {
808                              double sampleAvg[NCOMP_MAX];                              double sampleAvg[NCOMP_MAX];
809                              dim_t nCompUsed = MIN(nComp, NCOMP_MAX);                              dim_t nCompUsed = MIN(nComp, NCOMP_MAX);
# Line 869  void Finley_Mesh_saveVTK(const char *fil Line 868  void Finley_Mesh_saveVTK(const char *fil
868                                  fputs(tmpBuffer, fileHandle_p);                                  fputs(tmpBuffer, fileHandle_p);
869                              }                              }
870                          } /* for l (numCellFactor) */                          } /* for l (numCellFactor) */
                         freeSampleBuffer(sampleBuffer);  
871                      } /* if I am the owner */                      } /* if I am the owner */
872                  } /* for i (numCells) */                  } /* for i (numCells) */
873    
# Line 985  void Finley_Mesh_saveVTK(const char *fil Line 983  void Finley_Mesh_saveVTK(const char *fil
983                  for (i=0; i<mesh_p->Nodes->numNodes; i++) {                  for (i=0; i<mesh_p->Nodes->numNodes; i++) {
984                      k = globalNodeIndex[i];                      k = globalNodeIndex[i];
985                      if ( (myFirstNode <= k) && (k < myLastNode) ) {                      if ( (myFirstNode <= k) && (k < myLastNode) ) {
986                          void* sampleBuffer=allocSampleBuffer(data_pp[dataIdx]);                          __const double *values = getSampleDataRO(data_pp[dataIdx], nodeMapping->target[i]);
                         __const double *values = getSampleDataRO(data_pp[dataIdx], nodeMapping->target[i], sampleBuffer);  
987                          /* if the number of mpi_required components is more than                          /* if the number of mpi_required components is more than
988                           * the number of actual components, pad with zeros.                           * the number of actual components, pad with zeros.
989                           * Probably only need to get shape of first element */                           * Probably only need to get shape of first element */
# Line 1026  void Finley_Mesh_saveVTK(const char *fil Line 1023  void Finley_Mesh_saveVTK(const char *fil
1023                          } else {                          } else {
1024                              fputs(tmpBuffer, fileHandle_p);                              fputs(tmpBuffer, fileHandle_p);
1025                          }                          }
                         freeSampleBuffer(sampleBuffer);                 /* no-one needs values anymore */  
1026                      } /* if this is my node */                      } /* if this is my node */
1027                  } /* for i (numNodes) */                  } /* for i (numNodes) */
1028    

Legend:
Removed from v.2754  
changed lines
  Added in v.2770

  ViewVC Help
Powered by ViewVC 1.1.26