diff options
| author | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2002-09-13 14:12:18 +0000 |
|---|---|---|
| committer | jthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5> | 2002-09-13 14:12:18 +0000 |
| commit | e77d77af39eabf9aea5ca024ed1aca9df48460dc (patch) | |
| tree | 5bad2b1762ed2231efa74ce9e764465851726427 /src/driver/Newton.cc | |
| parent | 3bccbe234a26d06513131bf185440b2cd01d9222 (diff) | |
change $Id:$ --> $Header:$
rename src/util/ --> src/patch/
add Jacobian code for computing surface integrals
add support for choosing ASCII or HDF5 data files (HDF5 not implemented yet)
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@727 f88db872-0e4f-0410-b76b-b9085cfa78c5
Diffstat (limited to 'src/driver/Newton.cc')
| -rw-r--r-- | src/driver/Newton.cc | 22 |
1 files changed, 11 insertions, 11 deletions
diff --git a/src/driver/Newton.cc b/src/driver/Newton.cc index e7ce466..7267088 100644 --- a/src/driver/Newton.cc +++ b/src/driver/Newton.cc @@ -1,5 +1,5 @@ // Newton.cc -- solve H(h) = 0 via Newton's method -// $Id$ +// $Header$ // // Newton_solve - driver to solve H(h) = 0 via Newton's method // @@ -21,14 +21,14 @@ #include "../jtutil/linear_map.hh" using jtutil::error_exit; -#include "../util/coords.hh" -#include "../util/grid.hh" -#include "../util/fd_grid.hh" -#include "../util/patch.hh" -#include "../util/patch_edge.hh" -#include "../util/patch_interp.hh" -#include "../util/ghost_zone.hh" -#include "../util/patch_system.hh" +#include "../patch/coords.hh" +#include "../patch/grid.hh" +#include "../patch/fd_grid.hh" +#include "../patch/patch.hh" +#include "../patch/patch_edge.hh" +#include "../patch/patch_interp.hh" +#include "../patch/ghost_zone.hh" +#include "../patch/patch_system.hh" #include "../elliptic/Jacobian.hh" @@ -67,7 +67,7 @@ bool Newton_solve(patch_system& ps, then CCTK_VInfo(CCTK_THORNSTRING, "Newton iteration %d", iteration); - if (solver_info.output_h_and_H_at_each_Newton_iteration) + if (IO_info.output_h_and_H_at_each_Newton_iteration) then output_gridfn(ps, gfns::gfn__h, IO_info, IO_info.h_base_file_name, hn, verbose_info.print_algorithm_details, @@ -86,7 +86,7 @@ bool Newton_solve(patch_system& ps, iteration, H_norms.rms_norm(), H_norms.infinity_norm()); - if (solver_info.output_h_and_H_at_each_Newton_iteration) + if (IO_info.output_h_and_H_at_each_Newton_iteration) then output_gridfn(ps, gfns::gfn__H, IO_info, IO_info.H_base_file_name, hn, verbose_info.print_algorithm_details, |