yade-dev team mailing list archive

[noreply@xxxxxxxxxxxxx]

------------------------------------------------------------
revno: 3965
committer: Jan Stránský <jan.stransky@xxxxxxxxxxx>
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();