/[escript]/branches/domexper/dudley/src/Mesh_saveVTK.c
ViewVC logotype

Diff of /branches/domexper/dudley/src/Mesh_saveVTK.c

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

revision 2256 by jfenwick, Mon Feb 9 06:20:32 2009 UTC revision 2271 by jfenwick, Mon Feb 16 05:08:29 2009 UTC
# Line 762  void Finley_Mesh_saveVTK(const char *fil Line 762  void Finley_Mesh_saveVTK(const char *fil
762                  txtBufferInUse = 0;                  txtBufferInUse = 0;
763                  for (i=0; i<numCells; i++) {                  for (i=0; i<numCells; i++) {
764                      if (elements->Owner[i] == my_mpi_rank) {                      if (elements->Owner[i] == my_mpi_rank) {
765                          double *values = getSampleData(data_pp[dataIdx], i);              void* sampleBuffer=allocSampleBuffer(data_pp[dataIdx]);
766                            __const double *values = getSampleDataRO(data_pp[dataIdx], i,sampleBuffer);
767                          for (l = 0; l < numCellFactor; l++) {                          for (l = 0; l < numCellFactor; l++) {
768                              double sampleAvg[NCOMP_MAX];                              double sampleAvg[NCOMP_MAX];
769                              dim_t nCompUsed = MIN(nComp, NCOMP_MAX);                              dim_t nCompUsed = MIN(nComp, NCOMP_MAX);
# Line 828  void Finley_Mesh_saveVTK(const char *fil Line 829  void Finley_Mesh_saveVTK(const char *fil
829                                  fputs(tmpBuffer, fileHandle_p);                                  fputs(tmpBuffer, fileHandle_p);
830                              }                              }
831                          } /* for l (numCellFactor) */                          } /* for l (numCellFactor) */
832                freeSampleBuffer(sampleBuffer);
833                      } /* if I am the owner */                      } /* if I am the owner */
834                  } /* for i (numCells) */                  } /* for i (numCells) */
835    
# Line 943  void Finley_Mesh_saveVTK(const char *fil Line 945  void Finley_Mesh_saveVTK(const char *fil
945                  for (i=0; i<mesh_p->Nodes->numNodes; i++) {                  for (i=0; i<mesh_p->Nodes->numNodes; i++) {
946                      k = globalNodeIndex[i];                      k = globalNodeIndex[i];
947                      if ( (myFirstNode <= k) && (k < myLastNode) ) {                      if ( (myFirstNode <= k) && (k < myLastNode) ) {
948                          double *values = getSampleData(data_pp[dataIdx], nodeMapping->target[i]);              void* sampleBuffer=allocSampleBuffer(data_pp[dataIdx]);
949                            __const double *values = getSampleDataRO(data_pp[dataIdx], nodeMapping->target[i], sampleBuffer);
950                          /* if the number of mpi_required components is more than                          /* if the number of mpi_required components is more than
951                           * the number of actual components, pad with zeros.                           * the number of actual components, pad with zeros.
952                           * Probably only need to get shape of first element */                           * Probably only need to get shape of first element */
# Line 983  void Finley_Mesh_saveVTK(const char *fil Line 986  void Finley_Mesh_saveVTK(const char *fil
986                          } else {                          } else {
987                              fputs(tmpBuffer, fileHandle_p);                              fputs(tmpBuffer, fileHandle_p);
988                          }                          }
989                freeSampleBuffer(sampleBuffer);         /* no-one needs values anymore */
990                      } /* if this is my node */                      } /* if this is my node */
991                  } /* for i (numNodes) */                  } /* for i (numNodes) */
992    

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

  ViewVC Help
Powered by ViewVC 1.1.26