From 6e6ba9b7b6e3d9b23d3590221704be7c5d0a920f Mon Sep 17 00:00:00 2001
From: Tim Behrens <behrens@fmrib.ox.ac.uk>
Date: Wed, 22 Mar 2006 10:39:34 +0000
Subject: [PATCH] *** empty log message ***
---
dtifit.cc | 10 ++++++----
1 file changed, 6 insertions(+), 4 deletions(-)
diff --git a/dtifit.cc b/dtifit.cc
index cb10e23..491ce14 100644
--- a/dtifit.cc
+++ b/dtifit.cc
@@ -189,10 +189,12 @@ int main(int argc, char** argv)
Matrix r = read_ascii_matrix(opts.bvecsfile.value());
if(r.Nrows()>3) r=r.t();
for(int i=1;i<=r.Ncols();i++){
- float tmpsum=sqrt(r(i,1)*r(i,1)+r(i,2)*r(i,2)+r(i,3)*r(i,3));
- r(i,1)=r(i,1)/tmpsum;
- r(i,2)=r(i,2)/tmpsum;
- r(i,3)=r(i,3)/tmpsum;
+ float tmpsum=sqrt(r(1,i)*r(1,i)+r(2,i)*r(2,i)+r(3,i)*r(3,i));
+ if(tmpsum!=0){
+ r(1,i)=r(1,i)/tmpsum;
+ r(2,i)=r(2,i)/tmpsum;
+ r(3,i)=r(3,i)/tmpsum;
+ }
}
Matrix b = read_ascii_matrix(opts.bvalsfile.value());
if(b.Nrows()>1) b=b.t();
--
GitLab