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

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

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

revision 2670 by gross, Thu Sep 17 04:04:09 2009 UTC revision 2671 by gross, Thu Sep 17 06:37:22 2009 UTC
# Line 3107  Data::interpolateFromTable1D(const Wrapp Line 3107  Data::interpolateFromTable1D(const Wrapp
3107                          break;                          break;
3108                      }                      }
3109              // map x*Astep <= a << (x+1)*Astep to [-1,1]              // map x*Astep <= a << (x+1)*Astep to [-1,1]
3110                      double la = 2.0*(a-(x*Astep))/Astep-1;                      double la = 2.0*(a-Amin-(x*Astep))/Astep-1;
3111                      rdat[l]=((1-la)*e + (1+la)*w)/2;                      rdat[l]=((1-la)*e + (1+la)*w)/2;
3112                  }                  }
3113              }              }
# Line 3215  Data::interpolateFromTable2D(const Wrapp Line 3215  Data::interpolateFromTable2D(const Wrapp
3215                   error=2;                   error=2;
3216                   break;                   break;
3217                   }                   }
3218                   double lb = 2.0*(b-(y*Bstep))/Bstep-1;                   double lb = 2.0*(b-Bmin-(y*Bstep))/Bstep-1;
3219                   rdat[l]=((1-lb)*sw + (1+lb)*nw )/2.;                   rdat[l]=((1-lb)*sw + (1+lb)*nw )/2.;
3220    
3221                       }                       }
# Line 3228  Data::interpolateFromTable2D(const Wrapp Line 3228  Data::interpolateFromTable2D(const Wrapp
3228                   error=2;                   error=2;
3229                   break;                   break;
3230                   }                   }
3231                   double la = 2.0*(a-(x*Astep))/Astep-1;                   double la = 2.0*(a-Amin-(x*Astep))/Astep-1;
3232                   rdat[l]=((1-la)*sw + (1+la)*se )/2;                   rdat[l]=((1-la)*sw + (1+la)*se )/2;
3233    
3234                       } else {                       } else {
# Line 3243  Data::interpolateFromTable2D(const Wrapp Line 3243  Data::interpolateFromTable2D(const Wrapp
3243                   }                   }
3244                   // map x*Astep <= a << (x+1)*Astep to [-1,1]                   // map x*Astep <= a << (x+1)*Astep to [-1,1]
3245                   // same with b                   // same with b
3246                   double la = 2.0*(a-(x*Astep))/Astep-1;                   double la = 2.0*(a-Amin-(x*Astep))/Astep-1;
3247                   double lb = 2.0*(b-(y*Bstep))/Bstep-1;                   double lb = 2.0*(b-Bmin-(y*Bstep))/Bstep-1;
3248                   rdat[l]=((1-la)*(1-lb)*sw + (1-la)*(1+lb)*nw +                   rdat[l]=((1-la)*(1-lb)*sw + (1-la)*(1+lb)*nw +
3249                        (1+la)*(1-lb)*se + (1+la)*(1+lb)*ne)/4;                        (1+la)*(1-lb)*se + (1+la)*(1+lb)*ne)/4;
   
3250                       }                       }
3251    
3252                  }                  }

Legend:
Removed from v.2670  
changed lines
  Added in v.2671

  ViewVC Help
Powered by ViewVC 1.1.26