diff --git a/splinterpolator.h b/splinterpolator.h index 7f8da3445998756e3c4a656a41dba4bac22366cd..26174a9a958c004386f9814770673fccf3621cdc 100644 --- a/splinterpolator.h +++ b/splinterpolator.h @@ -608,7 +608,7 @@ double Splinterpolator<T>::value_at(const double *coord) const ni = get_start_indicies(coord,inds); get_wgts(coord,inds,wgts); - if (coord[0] < -4 && coord[0] > -5 && coord[1] < 143 && coord[1] > 142 && coord[2] < 92 && coord[2] > 91) { + if (coord[0] < -1 && coord[0] > -2 && coord[1] < 89 && coord[1] > 88 && coord[2] < 65 && coord[2] > 64) { for (unsigned int i=0; i<ni; i++) printf("coord[%d] = %f\n",i,coord[i]); for (unsigned int i=0; i<ni; i++) printf("inds[%d] = %d\n",i,inds[i]); for (int j=0; j<3; j++) { @@ -630,7 +630,7 @@ double Splinterpolator<T>::value_at(const double *coord) const int linear2 = add2linear(linear1,inds[1]+j); double *iiwgt=iwgt; for (unsigned int i=0; i<ni; i++, iiwgt++) { - if (coord[0] < -4 && coord[0] > -5 && coord[1] < 143 && coord[1] > 142 && coord[2] < 92 && coord[2] > 91) { + if (coord[0] < -1 && coord[0] > -2 && coord[1] < 89 && coord[1] > 88 && coord[2] < 65 && coord[2] > 64) { int ii = linear2+indx2indx(inds[0]+i,0); // printf("i = %d, j = %d, k = %d, ii = %d, wgt2 = %f, wgts[%d] = %f\n",i,j,k,ii,wgt2,i,*iiwgt); }