From ea40f9d70deb7a69d6eec286aaca4e80c8e4f320 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Wed, 11 Jun 2008 13:59:26 +0000
Subject: [PATCH] remove code which calls the lapack++ svd
[[Imported from SVN: r2420]]
---
src/averageinterface.hh | 6 ------
1 file changed, 6 deletions(-)
diff --git a/src/averageinterface.hh b/src/averageinterface.hh
index f71fce44..36e615d4 100644
--- a/src/averageinterface.hh
+++ b/src/averageinterface.hh
@@ -791,7 +791,6 @@ void computeAverageInterface(const BoundaryPatch<GridType>& interface,
FieldMatrix<double,dim,dim> deformationGradientBackup = deformationGradient;
-#if 1
// returns a decomposition U W VT, where U is returned in the first argument
svdcmp<double,dim,dim>(deformationGradient, W, V);
@@ -800,11 +799,6 @@ void computeAverageInterface(const BoundaryPatch<GridType>& interface,
VT[i][j] = V[j][i];
deformationGradient.rightmultiply(VT);
-#else
- lapackSVD(deformationGradientBackup, U, W, VT);
- deformationGradient = U;
- deformationGradient.rightmultiply(VT);
-#endif
//std::cout << deformationGradient << std::endl;
// deformationGradient now contains the orthogonal part of the polar decomposition
--
GitLab