From 60e42e68bb0ee596bf434e96ebf302319f91ad16 Mon Sep 17 00:00:00 2001
From: Saad Jbabdi <saad@fmrib.ox.ac.uk>
Date: Mon, 5 Nov 2012 09:24:25 +0000
Subject: [PATCH] bug fix
---
dtigen.cc | 6 +++---
1 file changed, 3 insertions(+), 3 deletions(-)
diff --git a/dtigen.cc b/dtigen.cc
index 82db386..9594f35 100644
--- a/dtigen.cc
+++ b/dtigen.cc
@@ -78,10 +78,10 @@ int do_dtigen(){
for(int t=1;t<=data.tsize();t++){
data(x,y,z,t-1) = S0(x,y,z);
data(x,y,z,t-1) *= exp( -b(1,t) * r(1,t) * r(1,t) * tensor(x,y,z,0) );
- data(x,y,z,t-1) *= exp( -b(1,t) * r(1,t) * r(2,t) * tensor(x,y,z,1) );
- data(x,y,z,t-1) *= exp( -b(1,t) * r(1,t) * r(3,t) * tensor(x,y,z,2) );
+ data(x,y,z,t-1) *= exp( -2*b(1,t) * r(1,t) * r(2,t) * tensor(x,y,z,1) );
+ data(x,y,z,t-1) *= exp( -2*b(1,t) * r(1,t) * r(3,t) * tensor(x,y,z,2) );
data(x,y,z,t-1) *= exp( -b(1,t) * r(2,t) * r(2,t) * tensor(x,y,z,3) );
- data(x,y,z,t-1) *= exp( -b(1,t) * r(2,t) * r(3,t) * tensor(x,y,z,4) );
+ data(x,y,z,t-1) *= exp( -2*b(1,t) * r(2,t) * r(3,t) * tensor(x,y,z,4) );
data(x,y,z,t-1) *= exp( -b(1,t) * r(3,t) * r(3,t) * tensor(x,y,z,5) );
}
--
GitLab