From 72b1b4f2a251ad157d3f1cd798c87c26ec7bab02 Mon Sep 17 00:00:00 2001
From: Jesper Andersson <jesper@fmrib.ox.ac.uk>
Date: Wed, 6 Mar 2013 15:05:16 +0000
Subject: [PATCH] removed debug printouts
---
splinterpolator.h | 15 ---------------
1 file changed, 15 deletions(-)
diff --git a/splinterpolator.h b/splinterpolator.h
index 8781cb2..49b258a 100644
--- a/splinterpolator.h
+++ b/splinterpolator.h
@@ -608,17 +608,6 @@ double Splinterpolator<T>::value_at(const double *coord) const
ni = get_start_indicies(coord,inds);
get_wgts(coord,inds,wgts);
- 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++) {
- for (int i=0; i<4; i++) {
- printf("wgts[%d] = %f\n",j*4+i,wgts[j][i]);
- }
- }
- }
-
-
double val=0.0;
for (unsigned int m=0, me=(_ndim>4)?ni:1; m<me; m++) {
for (unsigned int l=0, le=(_ndim>3)?ni:1; l<le; l++) {
@@ -630,10 +619,6 @@ 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] < -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, coef[%d] = %f\n",i,j,k,ii,wgt2,i,*iiwgt,ii,cptr[ii]);
- }
val += cptr[linear2+indx2indx(inds[0]+i,0)]*(*iiwgt)*wgt2;
}
}
--
GitLab