[cig-commits] r8954 - cs/benchmark/cigma/trunk/src/tests

[luis at geodynamics.org]

Author: luis
Date: 2007-12-19 12:07:37 -0800 (Wed, 19 Dec 2007)
New Revision: 8954

Added:
   cs/benchmark/cigma/trunk/src/tests/TestQuadrature.cpp
Log:
Quick test for Quadrature class

Added: cs/benchmark/cigma/trunk/src/tests/TestQuadrature.cpp
===================================================================
--- cs/benchmark/cigma/trunk/src/tests/TestQuadrature.cpp	2007-12-19 20:07:30 UTC (rev 8953)
+++ cs/benchmark/cigma/trunk/src/tests/TestQuadrature.cpp	2007-12-19 20:07:37 UTC (rev 8954)
@@ -0,0 +1,92 @@
+#include <iostream>
+#include <iomanip>
+#include "../Tet.h"
+#include "../Quadrature.h"
+
+int main()
+{
+    using namespace std;
+    using namespace cigma;
+
+    int i,j;
+
+    const int npts = 8;
+
+    double qwts[npts] = {
+        0.29583885,  0.12821632,  0.16925605,
+        0.07335544,  0.29583885,  0.12821632,
+        0.16925605,  0.07335544
+    };
+
+    double qpts[npts*3] = {
+        -0.68663473, -0.72789005, -0.75497035,
+        -0.83720867, -0.85864055,  0.08830369,
+        -0.86832263,  0.13186633, -0.75497035,
+        -0.93159441, -0.4120024 ,  0.08830369,
+         0.16949513, -0.72789005, -0.75497035,
+        -0.39245447, -0.85864055,  0.08830369,
+        -0.50857335,  0.13186633, -0.75497035,
+        -0.74470688, -0.4120024 ,  0.08830369
+    };
+
+    Tet tet;
+    Quadrature Q;
+    Q.set_data(&tet, qpts, qwts, npts);
+
+    /* change default numeric format */
+
+    cout << std::showpos
+         << std::fixed
+         << std::setprecision(8);
+
+
+    /* print reference points on tet */
+
+    cout << "Reference points" << endl;
+    for (i = 0; i < Q.n_points(); i++)
+    {
+        for (j = 0; j < Q.n_refdim(); j++)
+        {
+            cout << " " << Q(i,j);
+        }
+        cout << endl;
+    }
+    cout << endl;
+
+    
+    /* apply reference map to quadrature points
+     * and print out the resulting global coordinates
+     */ 
+    
+    cout << "Global coordinates" << endl;
+    
+    /*
+    double verts1[4*3] = {
+        0.0, 0.0, 0.0,
+        1.0, 0.0, 0.0,
+        0.0, 1.0, 0.0,
+        0.0, 0.0, 1.0
+    }; // */
+    double verts2[4*3] = {
+        -1, -1, -1,
+        +1, -1, -1,
+        -1, +1, -1,
+        -1, -1, +1
+    };
+
+    //tet.update_vertices(verts1, 4, 3);
+    tet.update_vertices(verts2, 4, 3);
+
+    Q.apply_refmap();
+
+    for (i = 0; i < Q.n_points(); i++)
+    {
+        for (j = 0; j < Q.n_globaldim(); j++)
+        {
+            cout << " " << Q.global_value(i,j);
+        }
+        cout << endl;
+    }
+
+    return 0;
+}