From 26e12abf8c67cc47d414315869d3773fdc1eceab Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Tue, 20 Nov 2007 18:21:06 +0000
Subject: [PATCH] minor cleanup
[[Imported from SVN: r1790]]
---
src/averageinterface.hh | 8 +-------
1 file changed, 1 insertion(+), 7 deletions(-)
diff --git a/src/averageinterface.hh b/src/averageinterface.hh
index 816f6279..ce94d518 100644
--- a/src/averageinterface.hh
+++ b/src/averageinterface.hh
@@ -264,9 +264,6 @@ void computeAverageInterface(const BoundaryPatch<GridType>& interface,
const typename NeighborIterator::Geometry& segmentGeometry = nIt.intersectionGlobal();
- const ReferenceElement<double,dim>& refElement = ReferenceElements<double, dim>::general(eIt->type());
- int nDofs = refElement.size(nIt.numberInSelf(),1,dim);
-
// Get quadrature rule
const QuadratureRule<double, dim-1>& quad = QuadratureRules<double, dim-1>::rule(segmentGeometry.type(), dim-1);
@@ -295,9 +292,6 @@ void computeAverageInterface(const BoundaryPatch<GridType>& interface,
const FieldMatrix<double,dim,dim>& inv = eIt->geometry().jacobianInverseTransposed(quadPos);
- /* Compute the weight of the current integration point */
- double weight = quad[ip].weight() * integrationElement;
-
/**********************************************/
/* compute gradients of the shape functions */
/**********************************************/
@@ -381,7 +375,7 @@ void computeAverageInterface(const BoundaryPatch<GridType>& interface,
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
assert( std::abs(1-deformationGradient.determinant()) < 1e-3);
--
GitLab