[cig-commits] r11854 - cs/benchmark/cigma/trunk/src

[luis at geodynamics.org]

Author: luis
Date: 2008-04-23 09:56:40 -0700 (Wed, 23 Apr 2008)
New Revision: 11854

Modified:
   cs/benchmark/cigma/trunk/src/Tet.cpp
   cs/benchmark/cigma/trunk/src/rules.py
Log:
Need to adjust integration weights for Tet cell!
This will be unnecessary once we switch to the same domain as FIAT


Modified: cs/benchmark/cigma/trunk/src/Tet.cpp
===================================================================
--- cs/benchmark/cigma/trunk/src/Tet.cpp	2008-04-23 16:56:37 UTC (rev 11853)
+++ cs/benchmark/cigma/trunk/src/Tet.cpp	2008-04-23 16:56:40 UTC (rev 11854)
@@ -88,6 +88,9 @@
             tet_qpts[n] += 1;
             tet_qpts[n] *= 0.5;
         }
+        // don't forget to adjust the integration weights
+        // to account for the reference domain transformation!
+        tet_qwts[i] *= 0.125;
     }
 
     *npts = tet_nno;

Modified: cs/benchmark/cigma/trunk/src/rules.py
===================================================================
--- cs/benchmark/cigma/trunk/src/rules.py	2008-04-23 16:56:37 UTC (rev 11853)
+++ cs/benchmark/cigma/trunk/src/rules.py	2008-04-23 16:56:40 UTC (rev 11854)
@@ -12,6 +12,7 @@
         x,w = qr(order)
         x += 1
         x *= 0.5
+        w *= 0.5 ** x.shape[1]
         return (x,w)
     return new_qr