Skip to content
Snippets Groups Projects
Commit 833488ca authored by Tim Behrens's avatar Tim Behrens
Browse files

*** empty log message ***

parent c9ca994c
No related branches found
No related tags found
No related merge requests found
...@@ -135,7 +135,6 @@ void tensorfit(DiagonalMatrix& Dd,ColumnVector& evec1,ColumnVector& evec2,Column ...@@ -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); logS(i)=(S(i)/s0)>0.01 ? log(S(i)):log(0.01*s0);
} }
Dvec = -pinv(Amat)*logS; Dvec = -pinv(Amat)*logS;
cout << Dvec;
s0=exp(-Dvec(7)); s0=exp(-Dvec(7));
if(s0<S.Sum()/S.Nrows()){ s0=S.Sum()/S.Nrows(); } if(s0<S.Sum()/S.Nrows()){ s0=S.Sum()/S.Nrows(); }
tens = vec2tens(Dvec); tens = vec2tens(Dvec);
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment