/[escript]/trunk/escript/src/DataExpanded.cpp
ViewVC logotype

Diff of /trunk/escript/src/DataExpanded.cpp

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

revision 1872 by jfenwick, Mon Oct 13 00:18:55 2008 UTC revision 1962 by gross, Tue Nov 4 07:17:31 2008 UTC
# Line 439  DataExpanded::copyToDataPoint(const int Line 439  DataExpanded::copyToDataPoint(const int
439       if ((dataPointNo >= numDataPointsPerSample) || (dataPointNo < 0)) {       if ((dataPointNo >= numDataPointsPerSample) || (dataPointNo < 0)) {
440             throw DataException("Error - DataExpanded::copyDataPoint invalid dataPointNoInSample.");             throw DataException("Error - DataExpanded::copyDataPoint invalid dataPointNoInSample.");
441       }       }
442         ValueType::size_type offset = getPointOffset(sampleNo, dataPointNo);
443       ValueType& vec=getVector();       ValueType& vec=getVector();
444       if (dataPointRank==0) {       if (dataPointRank==0) {
445           vec[0]=extract<double>(value[0]);           vec[offset]=extract<double>(value[0]);
446       } else if (dataPointRank==1) {       } else if (dataPointRank==1) {
447          for (int i=0; i<shape[0]; i++) {          for (int i=0; i<shape[0]; i++) {
448              vec[i]=extract<double>(value[i]);              vec[offset+i]=extract<double>(value[i]);
449          }          }
450       } else if (dataPointRank==2) {       } else if (dataPointRank==2) {
451          for (int i=0; i<shape[0]; i++) {          for (int i=0; i<shape[0]; i++) {
452             for (int j=0; j<shape[1]; j++) {             for (int j=0; j<shape[1]; j++) {
453                vec[getRelIndex(shape,i,j)]=extract<double>(value[i][j]);                vec[offset+getRelIndex(shape,i,j)]=extract<double>(value[i][j]);
454             }             }
455          }          }
456       } else if (dataPointRank==3) {       } else if (dataPointRank==3) {
457          for (int i=0; i<shape[0]; i++) {          for (int i=0; i<shape[0]; i++) {
458             for (int j=0; j<shape[1]; j++) {             for (int j=0; j<shape[1]; j++) {
459                for (int k=0; k<shape[2]; k++) {                for (int k=0; k<shape[2]; k++) {
460                   vec[getRelIndex(shape,i,j,k)]=extract<double>(value[i][j][k]);                   vec[offset+getRelIndex(shape,i,j,k)]=extract<double>(value[i][j][k]);
461                }                }
462             }             }
463          }          }
# Line 465  DataExpanded::copyToDataPoint(const int Line 466  DataExpanded::copyToDataPoint(const int
466             for (int j=0; j<shape[1]; j++) {             for (int j=0; j<shape[1]; j++) {
467               for (int k=0; k<shape[2]; k++) {               for (int k=0; k<shape[2]; k++) {
468                 for (int l=0; l<shape[3]; l++) {                 for (int l=0; l<shape[3]; l++) {
469                    vec[getRelIndex(shape,i,j,k,l)]=extract<double>(value[i][j][k][l]);                    vec[offset+getRelIndex(shape,i,j,k,l)]=extract<double>(value[i][j][k][l]);
470                 }                 }
471               }               }
472             }             }

Legend:
Removed from v.1872  
changed lines
  Added in v.1962

  ViewVC Help
Powered by ViewVC 1.1.26