[cig-commits] r11716 - in cs/benchmark/cigma/trunk/src: . hold

luis at geodynamics.org luis at geodynamics.org
Wed Apr 2 03:24:49 PDT 2008


Author: luis
Date: 2008-04-02 03:24:48 -0700 (Wed, 02 Apr 2008)
New Revision: 11716

Added:
   cs/benchmark/cigma/trunk/src/hold/FE.cpp
   cs/benchmark/cigma/trunk/src/hold/FE.h
Removed:
   cs/benchmark/cigma/trunk/src/FE.cpp
   cs/benchmark/cigma/trunk/src/FE.h
Log:
Put FE.{h,cpp} on hold


Deleted: cs/benchmark/cigma/trunk/src/FE.cpp
===================================================================
--- cs/benchmark/cigma/trunk/src/FE.cpp	2008-04-02 10:24:47 UTC (rev 11715)
+++ cs/benchmark/cigma/trunk/src/FE.cpp	2008-04-02 10:24:48 UTC (rev 11716)
@@ -1,40 +0,0 @@
-#include <cassert>
-#include "FE.h"
-
-// ---------------------------------------------------------------------------
-
-cigma::FE::FE()
-{
-    basis_tab = 0;
-    basis_jet = 0;
-}
-
-cigma::FE::~FE()
-{
-    if (basis_tab != 0) delete [] basis_tab;
-    if (basis_jet != 0) delete [] basis_jet;
-}
-
-// ---------------------------------------------------------------------------
-
-void cigma::FE::set_quadrature_points(QuadraturePoints *pts)
-{
-    assert(meshPart != 0);
-    assert(meshPart->cell != 0);
-
-    this->QuadratureRule::set_quadrature_points(pts);
-
-    assert(points != 0);
-    const int nq = points->n_points();
-    const int ndofs = meshPart->cell->n_nodes();
-
-    // get shape function values at known quadrature points
-    basis_tab = new double[nq * ndofs];
-    meshPart->cell->shape(nq, points->qpts, basis_tab);
-
-    // get shape function values at known quadrature points
-    //basis_jet = new double[nq * ndofs * dim];
-    //meshPart->cell->grad_shape(nq, quadrature->qpts, basis_jet);
-}
-
-// ---------------------------------------------------------------------------

Deleted: cs/benchmark/cigma/trunk/src/FE.h
===================================================================
--- cs/benchmark/cigma/trunk/src/FE.h	2008-04-02 10:24:47 UTC (rev 11715)
+++ cs/benchmark/cigma/trunk/src/FE.h	2008-04-02 10:24:48 UTC (rev 11716)
@@ -1,32 +0,0 @@
-#ifndef __FE_H__
-#define __FE_H__
-
-#include "Cell.h"
-#include "Points.h"
-#include "QuadraturePoints.h"
-#include "QuadratureRule.h"
-
-namespace cigma
-{
-    class FE;   // XXX: rename to FE_QR since now we derive from QuadratureRule
-}
-
-/**
- * @brief Tabulation for Finite Elements;
- *
- */
-class cigma::FE : public QuadratureRule
-{
-public:
-    FE();
-    ~FE();
-
-public:
-    void set_quadrature_points(QuadraturePoints *pts);
-
-public:
-    double *basis_tab;  // [nq x ndofs]
-    double *basis_jet;  // [nq x ndofs x celldim]
-};
-
-#endif

Copied: cs/benchmark/cigma/trunk/src/hold/FE.cpp (from rev 11715, cs/benchmark/cigma/trunk/src/FE.cpp)
===================================================================
--- cs/benchmark/cigma/trunk/src/hold/FE.cpp	                        (rev 0)
+++ cs/benchmark/cigma/trunk/src/hold/FE.cpp	2008-04-02 10:24:48 UTC (rev 11716)
@@ -0,0 +1,40 @@
+#include <cassert>
+#include "FE.h"
+
+// ---------------------------------------------------------------------------
+
+cigma::FE::FE()
+{
+    basis_tab = 0;
+    basis_jet = 0;
+}
+
+cigma::FE::~FE()
+{
+    if (basis_tab != 0) delete [] basis_tab;
+    if (basis_jet != 0) delete [] basis_jet;
+}
+
+// ---------------------------------------------------------------------------
+
+void cigma::FE::set_quadrature_points(QuadraturePoints *pts)
+{
+    assert(meshPart != 0);
+    assert(meshPart->cell != 0);
+
+    this->QuadratureRule::set_quadrature_points(pts);
+
+    assert(points != 0);
+    const int nq = points->n_points();
+    const int ndofs = meshPart->cell->n_nodes();
+
+    // get shape function values at known quadrature points
+    basis_tab = new double[nq * ndofs];
+    meshPart->cell->shape(nq, points->qpts, basis_tab);
+
+    // get shape function values at known quadrature points
+    //basis_jet = new double[nq * ndofs * dim];
+    //meshPart->cell->grad_shape(nq, quadrature->qpts, basis_jet);
+}
+
+// ---------------------------------------------------------------------------

Copied: cs/benchmark/cigma/trunk/src/hold/FE.h (from rev 11715, cs/benchmark/cigma/trunk/src/FE.h)
===================================================================
--- cs/benchmark/cigma/trunk/src/hold/FE.h	                        (rev 0)
+++ cs/benchmark/cigma/trunk/src/hold/FE.h	2008-04-02 10:24:48 UTC (rev 11716)
@@ -0,0 +1,32 @@
+#ifndef __FE_H__
+#define __FE_H__
+
+#include "Cell.h"
+#include "Points.h"
+#include "QuadraturePoints.h"
+#include "QuadratureRule.h"
+
+namespace cigma
+{
+    class FE;   // XXX: rename to FE_QR since now we derive from QuadratureRule
+}
+
+/**
+ * @brief Tabulation for Finite Elements;
+ *
+ */
+class cigma::FE : public QuadratureRule
+{
+public:
+    FE();
+    ~FE();
+
+public:
+    void set_quadrature_points(QuadraturePoints *pts);
+
+public:
+    double *basis_tab;  // [nq x ndofs]
+    double *basis_jet;  // [nq x ndofs x celldim]
+};
+
+#endif



More information about the cig-commits mailing list