diff options
Diffstat (limited to 'src/patch/patch.cc')
-rw-r--r-- | src/patch/patch.cc | 109 |
1 files changed, 2 insertions, 107 deletions
diff --git a/src/patch/patch.cc b/src/patch/patch.cc index 956f85d..d67569a 100644 --- a/src/patch/patch.cc +++ b/src/patch/patch.cc @@ -15,24 +15,13 @@ // patch::setup_patch_frontier // patch::border_on_edge // -// patch::write_grid3d_gridfn_at_iwr -// patch::write_grid3d_gridfn__adat -// patch::radially_extrapolate_gridfn_vector -// *** template instantiations for patch::radially_extrapolate_gridfn_vector -// #include <stdio.h> #include <assert.h> -#include <limits.h> #include <math.h> -#include <float.h> - -#include <vector.h> -#include <jt/stdc.h> -#include <jt/util.h> -#include <jt/util++.hh> -#include <jt/misc++.hh> +#include "jt/stdc.h" +#include "jt/util.hh" #include "policy.hh" #include "fp.hh" @@ -311,97 +300,3 @@ return edge_in.edge_is_min() : ( edge_in.edge_is_rho() ? max_rho_ghost_zone() : max_sigma_ghost_zone() ); } - -//***************************************************************************** -//***************************************************************************** -//***************************************************************************** - -// -// This function writes the values of a grid3d gridfn at a specified iwr -// to the already-open stdio stream output_fp (which is *not* closed -// when we're through), using the format -// # dpx dpy gridfn -// Each set of dpx = const lines is followed by a blank line. -// Each patch is followed by two blank lines. -// -void patch::write_grid3d_gridfn__adat(FILE *output_fp, - int gfn, - int iwr, bool want_borders) - const -{ -if (output_fp == NULL) - then error_exit(PANIC_EXIT, -"***** patch::write_grid3d_gridfn__adat:\n" -" got NULL output stream for gfn=%d iwr=%d!\n" -, - gfn, iwr); /*NOTREACHED*/ - -fprintf(output_fp, - "# %s patch: (dpx,dpy) = (%s,%s)\n", - name(), defn_of_dpx(), defn_of_dpy()); -fprintf(output_fp, - "# dpx\tdpy\tgridfn\n"); - - for (int irho = effective_min_irho(want_borders) ; - irho <= effective_max_irho(want_borders) ; - ++irho) - { - for (int isigma = effective_min_isigma(want_borders) ; - isigma <= effective_max_isigma(want_borders) ; - ++isigma) - { - fprintf(output_fp, - "%g\t%g\t%.20e\n", - double(dpx_of_irho_isigma(irho,isigma)), - double(dpy_of_irho_isigma(irho,isigma)), - double(gridfn__grid3d(gfn, iwr, irho, isigma))); - } - fprintf(output_fp, "\n"); - } - -fprintf(output_fp, "\n"); -} - -//***************************************************************************** - -// -// This function does radial extrapolation at a specified radial border, -// for all the gridfns in a specified gridfn vector. -// The loop nesting is cache-optimized for interleaved gridfns. -// -template <int pm> -void patch::radially_extrapolate_gridfn_vector(radial_extrapolation_type ret, - const gridfn_vector& gv) -{ - for (int irho = bordered_min_irho() ; - irho <= bordered_max_irho() ; - ++irho) - { - for (int isigma = bordered_min_isigma() ; - isigma <= bordered_max_isigma() ; - ++isigma) - { - for (int ii = 0 ; ii < gv.N_gridfns() ; ++ii) - { - const int gfn = gv.ith_gfn(ii); - - radially_extrapolate<pm>(ret, gfn, irho,isigma); - } - } - } -} - -//***************************************************************************** - -// -// *** template instantiations for patch::radially_extrapolate_gridfn_range -// - -template - void patch::radially_extrapolate_gridfn_vector<inner_boundary> - (radial_extrapolation_type ret, - const gridfn_vector& gv); -template - void patch::radially_extrapolate_gridfn_vector<outer_boundary> - (radial_extrapolation_type ret, - const gridfn_vector& gv); |