Skip to content
Snippets Groups Projects
Commit ea40f9d7 authored by Oliver Sander's avatar Oliver Sander Committed by sander@PCPOOL.MI.FU-BERLIN.DE
Browse files

remove code which calls the lapack++ svd 

[[Imported from SVN: r2420]]
parent 37b71a28
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 0 additions and 6 deletions
src/averageinterface.hh
+ 0
− 6
View file @ ea40f9d7
...@@ -791,7 +791,6 @@ void computeAverageInterface(const BoundaryPatch<GridType>& interface, ...@@ -791,7 +791,6 @@ void computeAverageInterface(const BoundaryPatch<GridType>& interface,
FieldMatrix<double,dim,dim> deformationGradientBackup = deformationGradient; FieldMatrix<double,dim,dim> deformationGradientBackup = deformationGradient;
#if 1
// returns a decomposition U W VT, where U is returned in the first argument // returns a decomposition U W VT, where U is returned in the first argument
svdcmp<double,dim,dim>(deformationGradient, W, V); svdcmp<double,dim,dim>(deformationGradient, W, V);
...@@ -800,11 +799,6 @@ void computeAverageInterface(const BoundaryPatch<GridType>& interface, ...@@ -800,11 +799,6 @@ void computeAverageInterface(const BoundaryPatch<GridType>& interface,
VT[i][j] = V[j][i]; VT[i][j] = V[j][i];
deformationGradient.rightmultiply(VT); deformationGradient.rightmultiply(VT);
#else
lapackSVD(deformationGradientBackup, U, W, VT);
deformationGradient = U;
deformationGradient.rightmultiply(VT);
#endif
//std::cout << deformationGradient << std::endl; //std::cout << deformationGradient << std::endl;
// deformationGradient now contains the orthogonal part of the polar decomposition // deformationGradient now contains the orthogonal part of the polar decomposition
......
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