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

[luis at geodynamics.org]

Author: luis
Date: 2008-12-09 18:12:17 -0800 (Tue, 09 Dec 2008)
New Revision: 13484

Modified:
   cs/cigma/trunk/src/Quadrature.cpp
   cs/cigma/trunk/src/Quadrature.h
Log:
Copy constructor for cigma::Quadrature

Modified: cs/cigma/trunk/src/Quadrature.cpp
===================================================================
--- cs/cigma/trunk/src/Quadrature.cpp	2008-12-10 02:12:16 UTC (rev 13483)
+++ cs/cigma/trunk/src/Quadrature.cpp	2008-12-10 02:12:17 UTC (rev 13484)
@@ -160,10 +160,23 @@
     }
 }
 
-/*
 Quadrature::Quadrature(const Quadrature& other)
 {
-} // */
+    cell_type = other.cell_type;
+    npts = other.npts;
+    ndim = other.ndim;
+    points = new double[npts * ndim];
+    weights = new double[npts];
+    for (int q = 0; q < npts; q++)
+    {
+        weights[q] = other.weights[q];
+        for (int j = 0; j < ndim; j++)
+        {
+            int n = ndim*q + j;
+            points[n] = other.points[n];
+        }
+    }
+}
 
 Quadrature::~Quadrature()
 {

Modified: cs/cigma/trunk/src/Quadrature.h
===================================================================
--- cs/cigma/trunk/src/Quadrature.h	2008-12-10 02:12:16 UTC (rev 13483)
+++ cs/cigma/trunk/src/Quadrature.h	2008-12-10 02:12:17 UTC (rev 13484)
@@ -18,7 +18,7 @@
 public:
     Quadrature();
     Quadrature(int npts, int ndim);
-    //Quadrature(const Quadrature& other);
+    Quadrature(const Quadrature& other);
     ~Quadrature();
 
     int n_points() const;