From 169bad15644dc458e8e2a73f225451ee4de05f4d Mon Sep 17 00:00:00 2001
From: Jesper Andersson <jesper@fmrib.ox.ac.uk>
Date: Wed, 6 Mar 2013 13:11:08 +0000
Subject: [PATCH] debug printouts for working on eddy
---
splinterpolator.h | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/splinterpolator.h b/splinterpolator.h
index 7f8da34..26174a9 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);
}
--
GitLab