Fixed a bug for the extrapolation in line, quad or hex element

ufz · Jun 24, 2016 · b1b0c8c · b1b0c8c
1 parent 7f61caa
commit b1b0c8c
Showing 1 changed file with 9 additions and 7 deletions.
diff --git a/FEM/fem_ele.cpp b/FEM/fem_ele.cpp
@@ -1396,17 +1396,19 @@ void CElement::SetExtropoGaussPoints(const int i)
 double CElement::CalcXi_p()
 {
 	MshElemType::type ElementType = MeshElement->GetElementType();
-	if (ElementType == MshElemType::LINE || ElementType == MshElemType::QUAD || ElementType == MshElemType::HEXAHEDRON)
+	Xi_p = 0.0;
+	if (   ElementType == MshElemType::LINE
+		|| ElementType == MshElemType::QUAD
+		|| ElementType == MshElemType::QUAD8
+		|| ElementType == MshElemType::HEXAHEDRON)
 	{
-		double r = .0;
 		for (gp = 0; gp < nGauss; gp++)
 		{
-			r = MXPGaussPkt(nGauss, gp);
-			if (fabs(r) > Xi_p)
-				Xi_p = fabs(r);
+			const double r = fabs(MXPGaussPkt(nGauss, gp));
+			if (r > Xi_p)
+				Xi_p = r;
 		}
-		r = 1.0 / Xi_p;
-		Xi_p = r;
+		Xi_p = 1.0 / Xi_p;
 		return Xi_p;
 	}
 	else