From 6595e9dc0dfaf0e331d658e13b87d2cf3625f94a Mon Sep 17 00:00:00 2001
From: Stamatios Sotiropoulos <stam@fmrib.ox.ac.uk>
Date: Fri, 3 Aug 2012 12:00:12 +0000
Subject: [PATCH] Minor changes in BIC, f0 flags
---
xfibres.cc | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/xfibres.cc b/xfibres.cc
index c68a14c..3bdc9fc 100644
--- a/xfibres.cc
+++ b/xfibres.cc
@@ -594,7 +594,7 @@ public:
}
}
else{ //Do constrained optimization
- PVM_single_c pvm(m_data,m_bvecs,m_bvals,opts.nfibres.value(),opts.f0.value());
+ PVM_single_c pvm(m_data,m_bvecs,m_bvals,opts.nfibres.value(),false,opts.f0.value());
pvm.fit(); // this will give th,ph,f in the correct order
pvmf = pvm.get_f();
@@ -610,7 +610,7 @@ public:
//If the full model gives values that are considered implausible, or we are in a CSF voxel (f1<0.05)
//then fit a model without the f0 and drive f0_init to almost zero
if ((opts.nfibres.value()>0 && pvmf(1)<0.05) || pvmd>0.007 || pvmf0>0.4){
- PVM_single_c pvm2(m_data,m_bvecs,m_bvals,opts.nfibres.value(),false);
+ PVM_single_c pvm2(m_data,m_bvecs,m_bvals,opts.nfibres.value(),false,false);
pvm2.fit(); // this will give th,ph,f in the correct order
pvmf0=0.001;
pvmS0=pvm2.get_s0();
--
GitLab