diff options
author | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2002-07-18 17:40:09 +0000 |
---|---|---|
committer | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2002-07-18 17:40:09 +0000 |
commit | e01a3c4dfd2cea40f95c7c49084ccc26d68f1a81 (patch) | |
tree | 047647e67474062bf7f5e99776a42980ba6cbeb4 /src/elliptic/Jacobian.cc | |
parent | 8b734e46072a9614f8ffe12d3a10dcca558644a6 (diff) |
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
Diffstat (limited to 'src/elliptic/Jacobian.cc')
-rw-r--r-- | src/elliptic/Jacobian.cc | 81 |
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 +} |