Lower constant tet_volume multiplication below

accumulation loop.
2012-06-26 20:59:19 +02:00 · 2012-06-26 20:59:19 +02:00 · 5d5404aebf
commit 5d5404aebf
parent 0e308c3be1
1 changed files with 2 additions and 1 deletions
--- a/geometry.c
+++ b/geometry.c
@ -200,8 +200,9 @@ compute_cell_geometry(int ndims, double *coords,

            tet_volume = 0;
            for(i=0; i<ndims; ++i){
-                tet_volume += 0.5/3 * w[i]*(x[i]-xcell[i]);
+                tet_volume += w[i] * (x[i] - xcell[i]);
            }
+	    tet_volume *= 0.5 / 3;

            subnormal_sign=0.0;
            for(i=0; i<ndims; ++i){