storage for Jacobian matrix which knows ../util/ sparsity structure

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@631 f88db872-0e4f-0410-b76b-b9085cfa78c5

1 files changed, 81 insertions, 0 deletions

diff --git a/src/elliptic/Jacobian.cc b/src/elliptic/Jacobian.cc
new file mode 100644
index 0000000..1bd2b5b
--- /dev/null
+++ b/src/elliptic/Jacobian.cc
@@ -0,0 +1,81 @@
+// Jacobian.cc -- data structures for the Jacobian matrix
+// $Id$
+
+//
+// Jacobian::Jacobian
+//
+// dense_Jacobian::dense_Jacobian
+// dense_Jacobian::zero
+// dense_Jacobian::zero_row
+//
+
+#include <stdio.h>
+#include <assert.h>
+#include <math.h>
+#include <vector>
+
+#include "util_Table.h"
+#include "cctk.h"
+#include "cctk_Arguments.h"
+
+#include "stdc.h"
+#include "config.hh"
+#include "../jtutil/util.hh"
+#include "../jtutil/array.hh"
+#include "../jtutil/cpm_map.hh"
+#include "../jtutil/linear_map.hh"
+using jtutil::error_exit;
+
+#include "coords.hh"
+#include "grid.hh"
+#include "fd_grid.hh"
+#include "patch.hh"
+#include "patch_edge.hh"
+#include "patch_interp.hh"
+#include "ghost_zone.hh"
+#include "patch_system.hh"
+#include "Jacobian.hh"
+
+//******************************************************************************
+//******************************************************************************
+//******************************************************************************
+
+//
+// This function constructs a  Jacobian  object.
+//
+Jacobian::Jacobian(const patch_system& ps)
+	: ps_(ps)
+{ }
+
+//******************************************************************************
+//******************************************************************************
+//******************************************************************************
+
+//
+// This function constructs a  dense_Jacobian  object.
+//
+dense_Jacobian::dense_Jacobian(const patch_system& ps)
+	: Jacobian(ps),
+	  matrix_(0, ps.N_grid_points(),
+		  0, ps.N_grid_points())
+{ }
+
+//******************************************************************************
+
+//
+// This function zeros a  dense_Jacobian  object.
+//
+void dense_Jacobian::zero()
+{
+// FIXME FIXME
+}
+
+//******************************************************************************
+
+//
+// This function zeros a single row of a  dense_Jacobian  object.
+//
+void dense_Jacobian::zero_row(int II)
+{
+// FIXME FIXME
+}