[cig-commits] r13979 - cs/cigma/trunk/src

[luis at geodynamics.org]

Author: luis
Date: 2009-01-28 09:05:06 -0800 (Wed, 28 Jan 2009)
New Revision: 13979

Modified:
   cs/cigma/trunk/src/Quadrature.cpp
Log:
Adjust default integration rule for Cell::TRI3

Modified: cs/cigma/trunk/src/Quadrature.cpp
===================================================================
--- cs/cigma/trunk/src/Quadrature.cpp	2009-01-28 17:05:05 UTC (rev 13978)
+++ cs/cigma/trunk/src/Quadrature.cpp	2009-01-28 17:05:06 UTC (rev 13979)
@@ -147,6 +147,19 @@
         double qwts[nno] = {
             0.63608276,  0.36391724,  0.63608276,  0.36391724
         };
+
+        for (int i = 0; i < nno; i++)
+        {   
+            for (int j = 0; j < celldim; j++)
+            {   
+                int n = celldim * i + j;
+                qpts[n] += 1;
+                qpts[n] *= 0.5;
+            }
+            // don't forget to adjust the integration weights
+            // to account for the reference domain transformation!
+            qwts[i] *= 0.25;
+        }
         
         Q.reset(new Quadrature(nno, celldim));
         Q->setData(qpts, qwts);