Corrected oldstyle usage of eigen SVD
This commit is contained in:
parent
b95684e6e5
commit
d1b901cbb5
|
@ -187,8 +187,10 @@ public:
|
||||||
|
|
||||||
b[c] = n;
|
b[c] = n;
|
||||||
}
|
}
|
||||||
|
|
||||||
A.svd().solve(b, &sol);
|
Eigen::JacobiSVD<Eigen::MatrixXd> svd(A);
|
||||||
|
sol=svd.solve(b);
|
||||||
|
// A.svd().solve(b, &sol);
|
||||||
|
|
||||||
vector<double> r(16);
|
vector<double> r(16);
|
||||||
|
|
||||||
|
@ -435,7 +437,10 @@ public:
|
||||||
bm[c] = onedimensional[c];
|
bm[c] = onedimensional[c];
|
||||||
}
|
}
|
||||||
|
|
||||||
Am.svd().solve(bm, &sol);
|
|
||||||
|
// Am.svd().solve(bm, &sol);
|
||||||
|
Eigen::JacobiSVD<Eigen::MatrixXd> svd(Am);
|
||||||
|
sol=svd.solve(bm);
|
||||||
|
|
||||||
it->Q() = pow((double)sol[0],0.25);
|
it->Q() = pow((double)sol[0],0.25);
|
||||||
|
|
||||||
|
|
Loading…
Reference in New Issue