/[escript]/trunk/ripley/src/Rectangle.cpp
ViewVC logotype

Diff of /trunk/ripley/src/Rectangle.cpp

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

revision 4154 by jfenwick, Tue Jan 22 09:30:23 2013 UTC revision 4174 by caltinay, Wed Jan 30 03:21:27 2013 UTC
# Line 243  void Rectangle::readNcGrid(escript::Data Line 243  void Rectangle::readNcGrid(escript::Data
243          for (index_t x=0; x<num0; x++) {          for (index_t x=0; x<num0; x++) {
244              const int dataIndex = (first1+y)*myN0+first0+x;              const int dataIndex = (first1+y)*myN0+first0+x;
245              const int srcIndex=y*num0+x;              const int srcIndex=y*num0+x;
246              double* dest = out.getSampleDataRW(dataIndex);              if (!isnan(values[srcIndex])) {
247              for (index_t q=0; q<dpp; q++) {                  double* dest = out.getSampleDataRW(dataIndex);
248                  *dest++ = values[srcIndex];                  for (index_t q=0; q<dpp; q++) {
249                        *dest++ = values[srcIndex];
250                    }
251              }              }
252          }          }
253      }      }
# Line 313  void Rectangle::readBinaryGrid(escript:: Line 315  void Rectangle::readBinaryGrid(escript::
315          for (index_t x=0; x<num0; x++) {          for (index_t x=0; x<num0; x++) {
316              double* dest = out.getSampleDataRW(first0+x+(first1+y)*myN0);              double* dest = out.getSampleDataRW(first0+x+(first1+y)*myN0);
317              for (index_t c=0; c<numComp; c++) {              for (index_t c=0; c<numComp; c++) {
318                  for (index_t q=0; q<dpp; q++) {                  if (!isnan(values[x*numComp+c])) {
319                      *dest++ = static_cast<double>(values[x*numComp+c]);                      for (index_t q=0; q<dpp; q++) {
320                            *dest++ = static_cast<double>(values[x*numComp+c]);
321                        }
322                  }                  }
323              }              }
324          }          }

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

  ViewVC Help
Powered by ViewVC 1.1.26