/[escript]/branches/diaplayground/ripley/src/Brick.cpp
ViewVC logotype

Diff of /branches/diaplayground/ripley/src/Brick.cpp

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

revision 4173 by jfenwick, Tue Jan 22 09:30:23 2013 UTC revision 4174 by caltinay, Wed Jan 30 03:21:27 2013 UTC
# Line 314  void Brick::readBinaryGrid(escript::Data Line 314  void Brick::readBinaryGrid(escript::Data
314              for (index_t x=0; x<num0; x++) {              for (index_t x=0; x<num0; x++) {
315                  double* dest = out.getSampleDataRW(dataIndex+x);                  double* dest = out.getSampleDataRW(dataIndex+x);
316                  for (index_t c=0; c<numComp; c++) {                  for (index_t c=0; c<numComp; c++) {
317                      for (index_t q=0; q<dpp; q++) {                      if (!isnan(values[x*numComp+c])) {
318                          *dest++ = static_cast<double>(values[x*numComp+c]);                          for (index_t q=0; q<dpp; q++) {
319                                *dest++ = static_cast<double>(values[x*numComp+c]);
320                            }
321                      }                      }
322                  }                  }
323              }              }
# Line 418  void Brick::readNcGrid(escript::Data& ou Line 420  void Brick::readNcGrid(escript::Data& ou
420              for (index_t x=0; x<num0; x++) {              for (index_t x=0; x<num0; x++) {
421                  const int dataIndex = first0+(first1+y)*myN0+(first2+z)*myN0*myN1+x;                  const int dataIndex = first0+(first1+y)*myN0+(first2+z)*myN0*myN1+x;
422                  const int srcIndex=z*num1*num0+y*num0+x;                  const int srcIndex=z*num1*num0+y*num0+x;
423                  double* dest = out.getSampleDataRW(dataIndex);                  if (!isnan(values[srcIndex])) {
424                  for (index_t q=0; q<dpp; q++) {                      double* dest = out.getSampleDataRW(dataIndex);
425                      *dest++ = values[srcIndex];                      for (index_t q=0; q<dpp; q++) {
426                            *dest++ = values[srcIndex];
427                        }
428                  }                  }
429              }              }
430          }          }

Legend:
Removed from v.4173  
changed lines
  Added in v.4174

  ViewVC Help
Powered by ViewVC 1.1.26