Commit 753224ca authored by David Bommes's avatar David Bommes

don't perform a unnecessary matrix copy

git-svn-id: http://www.openflipper.org/svnrepo/CoMISo/trunk@214 1355f012-dd97-4b2f-ae87-10fa9f823a57
parent bc520782
......@@ -991,13 +991,14 @@ void factored_to_quadratic( MatrixT& _F, MatrixT2& _Q, VectorT& _rhs)
gmm::resize(_Q, n-1, n-1);
gmm::resize(_rhs, n);
// set up transposed
MatrixT Ft(n,m);
gmm::copy(gmm::transposed(_F), Ft);
// // set up transposed
// MatrixT Ft(n,m);
// gmm::copy(gmm::transposed(_F), Ft);
// compute quadratic matrix
MatrixT Q(n,n);
gmm::mult(Ft,_F,Q);
// gmm::mult(Ft,_F,Q);
gmm::mult(gmm::transposed(_F),_F,Q);
// extract rhs
gmm::copy( gmm::scaled(gmm::mat_const_row( Q, n - 1),-1.0), _rhs);
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment