From 3184a38ffab57e06782606997727d2472f3f3940 Mon Sep 17 00:00:00 2001
From: Innocent Niyonzima <in2199@columbia.edu>
Date: Thu, 16 Dec 2010 10:45:52 +0000
Subject: [PATCH] fix bug in MElement::getVolume (uninitialized variable!)
---
Geo/MElement.cpp | 12 ++++++------
1 file changed, 6 insertions(+), 6 deletions(-)
diff --git a/Geo/MElement.cpp b/Geo/MElement.cpp
index 3150b52e67..bad7035122 100644
--- a/Geo/MElement.cpp
+++ b/Geo/MElement.cpp
@@ -146,19 +146,19 @@ SPoint3 MElement::barycenter()
return p;
}
-double MElement::getVolume(){
+double MElement::getVolume()
+{
int npts;
IntPt *pts;
- getIntegrationPoints(getDim()*(getPolynomialOrder()-1), &npts, &pts);
- double vol;
- for (int i=0;i<npts;i++){
- vol += getJacobianDeterminant(pts[i].pt[0],pts[i].pt[1],pts[i].pt[2])
+ getIntegrationPoints(getDim() * (getPolynomialOrder() - 1), &npts, &pts);
+ double vol = 0.;
+ for (int i = 0; i < npts; i++){
+ vol += getJacobianDeterminant(pts[i].pt[0], pts[i].pt[1], pts[i].pt[2])
* pts[i].weight;
}
return vol;
}
-
int MElement::getVolumeSign()
{
double v = getVolume();
--
GitLab