diff --git a/dtifit.cc b/dtifit.cc index ae84c19015d96c221bd24df0e69fea70662c2d13..73f7532fe5d67c77775398f09095260bec0bfe45 100644 --- a/dtifit.cc +++ b/dtifit.cc @@ -135,7 +135,6 @@ void tensorfit(DiagonalMatrix& Dd,ColumnVector& evec1,ColumnVector& evec2,Column logS(i)=(S(i)/s0)>0.01 ? log(S(i)):log(0.01*s0); } Dvec = -pinv(Amat)*logS; - cout << Dvec; s0=exp(-Dvec(7)); if(s0<S.Sum()/S.Nrows()){ s0=S.Sum()/S.Nrows(); } tens = vec2tens(Dvec);