[cig-commits] r11521 - cs/benchmark/cigma/trunk/src

[luis at geodynamics.org]

Author: luis
Date: 2008-03-24 09:27:09 -0700 (Mon, 24 Mar 2008)
New Revision: 11521

Modified:
   cs/benchmark/cigma/trunk/src/Residuals.cpp
   cs/benchmark/cigma/trunk/src/Residuals.h
Log:
Renamed ResidualField to Residuals


Modified: cs/benchmark/cigma/trunk/src/Residuals.cpp
===================================================================
--- cs/benchmark/cigma/trunk/src/Residuals.cpp	2008-03-24 16:27:08 UTC (rev 11520)
+++ cs/benchmark/cigma/trunk/src/Residuals.cpp	2008-03-24 16:27:09 UTC (rev 11521)
@@ -1,7 +1,8 @@
 #include <iostream>
+#include <cassert>
 #include <string>
 #include <cmath>
-#include "ResidualField.h"
+#include "Residuals.h"
 #include "VtkWriter.h"
 
 using namespace std;
@@ -9,21 +10,21 @@
 
 // ---------------------------------------------------------------------------
 
-ResidualField::ResidualField()
+Residuals::Residuals()
 {
     nel = 0;
     epsilon = 0;
     meshPart = 0;
 }
 
-ResidualField::~ResidualField()
+Residuals::~Residuals()
 {
     if (epsilon != 0) delete [] epsilon;
 }
 
 // ---------------------------------------------------------------------------
 
-void ResidualField::set_mesh(MeshPart *meshPart)
+void Residuals::set_mesh(MeshPart *meshPart)
 {
     assert(meshPart != 0);
     this->meshPart = meshPart;
@@ -33,7 +34,7 @@
 
 // ---------------------------------------------------------------------------
 
-void ResidualField::zero_out()
+void Residuals::zero_out()
 {
     assert(nel > 0);
     assert(epsilon != 0);
@@ -46,7 +47,7 @@
     }
 }
 
-void ResidualField::update(int e, double cell_residual)
+void Residuals::update(int e, double cell_residual)
 {
     // remember residual on this cell
     epsilon[e] = cell_residual;
@@ -56,14 +57,14 @@
 
 }
 
-double ResidualField::L2()
+double Residuals::L2()
 {
     return sqrt(global_error); // XXX: generalize to norm_pow(global_error)
 }
 
 // ---------------------------------------------------------------------------
 
-void ResidualField::write_vtk(const char *filename)
+void Residuals::write_vtk(const char *filename)
 {
     assert(meshPart != 0); // XXX: add support for dumping residuals w/o mesh
     assert(meshPart->nel == nel);
@@ -72,7 +73,7 @@
     cout << "Creating file " << output_filename << endl;
 
     VtkWriter writer;
-    writer.open(output_filename); // XXX: change signature to 'const char *'
+    writer.open(output_filename.c_str());
     //writer.write_header();
     writer.write_points(meshPart->coords, meshPart->nno, meshPart->nsd);
     writer.write_cells(meshPart->connect, meshPart->nel, meshPart->ndofs);

Modified: cs/benchmark/cigma/trunk/src/Residuals.h
===================================================================
--- cs/benchmark/cigma/trunk/src/Residuals.h	2008-03-24 16:27:08 UTC (rev 11520)
+++ cs/benchmark/cigma/trunk/src/Residuals.h	2008-03-24 16:27:09 UTC (rev 11521)
@@ -1,19 +1,18 @@
-#ifndef __RESIDUAL_FIELD_H__
-#define __RESIDUAL_FIELD_H__
+#ifndef __RESIDUALS_H__
+#define __RESIDUALS_H__
 
-//#include "Field.h"
 #include "MeshPart.h"
 
 namespace cigma
 {
-    class ResidualField;
+    class Residuals;
 }
 
-class cigma::ResidualField
+class cigma::Residuals
 {
 public:
-    ResidualField();
-    ~ResidualField();
+    Residuals();
+    ~Residuals();
 
 public:
     void set_mesh(MeshPart *meshPart);