Changed how element value is calculated from material point values in…

… FEMeshAdaptorCriterion.
febiosoftware · Oct 11, 2024 · d4174dc · d4174dc
1 parent cf0b86b
commit d4174dc
Showing 1 changed file with 8 additions and 6 deletions.
diff --git a/FECore/FEMeshAdaptorCriterion.cpp b/FECore/FEMeshAdaptorCriterion.cpp
@@ -88,18 +88,20 @@ FEMeshAdaptorSelection FEMeshAdaptorCriterion::GetElementSelection(FEElementSet*
 
 bool FEMeshAdaptorCriterion::GetElementValue(FEElement& el, double& value)
 {
-	bool bvalid = true;
 	value = 0.0;
-	int ni = el.GaussPoints();
+	int ni = el.GaussPoints(), nv = 0;
 	for (int i = 0; i < ni; ++i)
 	{
 		double vali = 0.0;
 		bool b = GetMaterialPointValue(*el.GetMaterialPoint(i), vali);
-		if (b) value += vali;
-		bvalid = (bvalid && b);
+		if (b) { value += vali; nv++; }
 	}
-	value /= (double)ni;
-	return bvalid;
+	if (nv > 0)
+	{
+		value /= (double)nv;
+		return true;
+	}
+	else return false;
 }
 
 bool FEMeshAdaptorCriterion::GetMaterialPointValue(FEMaterialPoint& mp, double& elemVal)