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

[luis at geodynamics.org]

Author: luis
Date: 2008-12-09 18:13:05 -0800 (Tue, 09 Dec 2008)
New Revision: 13516

Modified:
   cs/cigma/trunk/src/FE.cpp
   cs/cigma/trunk/src/FE.h
Log:
Allocate jxw array in cigma::FE (how did this work before?)

Modified: cs/cigma/trunk/src/FE.cpp
===================================================================
--- cs/cigma/trunk/src/FE.cpp	2008-12-10 02:13:03 UTC (rev 13515)
+++ cs/cigma/trunk/src/FE.cpp	2008-12-10 02:13:05 UTC (rev 13516)
@@ -35,7 +35,9 @@
     assert(quadrature);
     const int nq = quadrature->n_points();
     const int ndofs = cell->n_nodes();
+    if (jxw != 0) { delete [] jxw; }
     if (basis_tab != 0) { delete [] basis_tab; }
+    jxw = new double[nq];
     basis_tab = new double[nq * ndofs];
     cell->shape(nq, quadrature->points, basis_tab);
 }

Modified: cs/cigma/trunk/src/FE.h
===================================================================
--- cs/cigma/trunk/src/FE.h	2008-12-10 02:13:03 UTC (rev 13515)
+++ cs/cigma/trunk/src/FE.h	2008-12-10 02:13:05 UTC (rev 13516)
@@ -19,12 +19,13 @@
     FE(const FE& other);
     ~FE();
 
-    //void init(boost::shared_ptr<Quadrature> Q, boost::shared_ptr<Cell> cell);
     void init(boost::shared_ptr<Quadrature> Q);
     void init_basis();
     void init_basis_jet();
     void update_jxw();
 
+    //void apply_refmap(Cell& cell);
+
 public:
     boost::shared_ptr<Cell> cell;
     boost::shared_ptr<Quadrature> quadrature;