diff --git a/Geo/MElement.cpp b/Geo/MElement.cpp index 3150b52e676f3ff9ff4a09f33b48d51af96b5952..bad7035122383a1a7dff22d77084c63cb9485576 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();