From 5dc8ef1a85a6e54093cbeda55763c378d19d305d Mon Sep 17 00:00:00 2001
From: Tim Behrens <behrens@fmrib.ox.ac.uk>
Date: Thu, 29 Jul 2004 14:21:56 +0000
Subject: [PATCH] *** empty log message ***
---
diff_pvm.cc | 12 ++++++++----
1 file changed, 8 insertions(+), 4 deletions(-)
diff --git a/diff_pvm.cc b/diff_pvm.cc
index 02575f6..834b50b 100644
--- a/diff_pvm.cc
+++ b/diff_pvm.cc
@@ -307,10 +307,10 @@ void Diff_pvmModel::initialise(const ColumnVector& S){
else fsquared=0;
if(fsquared>0){f=sqrt(fsquared);}
else{f=0;}
- if(f>=0.95) f=0.95;
- if(f<=0.001) f=0.001;
- cout<<"S0 "<<S0<<endl;
- cout<<"S1 "<<S(1)<<endl;
+ if(f>=0.95) f=0.95;
+ if(f<=0.001) f=0.001;
+ //cout<<"S0 "<<S0<<endl;
+ //cout<<"S1 "<<S(1)<<endl;
getparam(0).setinitvalue(th);
getparam(1).setinitvalue(ph);
getparam(2).setinitvalue(f);
@@ -339,8 +339,12 @@ int main(int argc, char *argv[])
Tracer_Plus::settimingon();
// read data
+
VolumeSeries data;
data.read(opts.datafile.value());
+ data.writeAsFloat(LogSingleton::getInstance().appendDir("data"));
+ cout<<"done"<<endl;
+ return 0;
int ntpts = data.tsize();
Matrix bvecs = read_ascii_matrix(opts.bvecsfile.value());
Matrix bvals = read_ascii_matrix(opts.bvalsfile.value());
--
GitLab