diff options
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 +} |