------------------------------------------------------------
revno: 3965
committer: Jan Stránský <jan.stran...@fsv.cvut.cz>
timestamp: Tue 2016-11-08 17:18:55 +0100
message:
Matrix_computeUnitaryPositive, changed computeThinX to computeFullX flags
should solve question #403514
modified:
lib/base/Math.hpp
--
lp:yade
https://code.launchpad.net/~yade-pkg/yade/git-trunk
Your team Yade developers is subscribed to branch lp:yade.
To unsubscribe from this branch go to
https://code.launchpad.net/~yade-pkg/yade/git-trunk/+edit-subscription
=== modified file 'lib/base/Math.hpp'
--- lib/base/Math.hpp 2015-11-09 20:26:06 +0000
+++ lib/base/Math.hpp 2016-11-08 16:18:55 +0000
@@ -177,7 +177,8 @@
template<typename MatrixT>
void Matrix_computeUnitaryPositive(const MatrixT& in, MatrixT* unitary, MatrixT* positive){
assert(unitary); assert(positive);
- Eigen::JacobiSVD<MatrixT> svd(in, Eigen::ComputeThinU | Eigen::ComputeThinV);
+ //Eigen::JacobiSVD<MatrixT> svd(in, Eigen::ComputeThinU | Eigen::ComputeThinV);
+ Eigen::JacobiSVD<MatrixT> svd(in, Eigen::ComputeFullU | Eigen::ComputeFullV);
MatrixT mU, mV, mS;
mU = svd.matrixU();
mV = svd.matrixV();
_______________________________________________
Mailing list: https://launchpad.net/~yade-dev
Post to : yade-dev@lists.launchpad.net
Unsubscribe : https://launchpad.net/~yade-dev
More help : https://help.launchpad.net/ListHelp
