Block-Structured AMR Software Framework
AMReX_EB_Redistribution.H
Go to the documentation of this file.
1 #ifndef AMREX_EB_REDIST_H_
2 #define AMREX_EB_REDIST_H_
3 #include <AMReX_Config.H>
4 
5 #include <AMReX.H>
6 #include <AMReX_BCRec.H>
7 #include <AMReX_MultiFab.H>
8 #include <AMReX_MultiFabUtil.H>
9 #include <AMReX_MultiCutFab.H>
10 #include <AMReX_EB2.H>
11 
12 namespace amrex {
13 
14 #ifdef AMREX_USE_FLOAT
15  static constexpr amrex::Real eb_covered_val = amrex::Real(1.e20);
16 #else
17  static constexpr amrex::Real eb_covered_val = amrex::Real(1.e40);
18 #endif
19 
20  void single_level_redistribute (amrex::MultiFab& div_tmp_in, amrex::MultiFab& div_out,
21  int div_comp, int ncomp, const amrex::Geometry& geom);
22 
23  void single_level_weighted_redistribute (amrex::MultiFab& div_tmp_in, amrex::MultiFab& div_out,
24  const amrex::MultiFab& weights,
25  int div_comp, int ncomp, const amrex::Geometry& geom,
26  bool use_wts_in_divnc);
27 
28  // Flux redistribution scheme for case where C/F interface does *not* cross EB
29  void apply_flux_redistribution ( const amrex::Box& bx,
30  amrex::Array4<amrex::Real > const& div,
31  amrex::Array4<amrex::Real const> const& divc,
32  amrex::Array4<amrex::Real const> const& wt,
33  int icomp, int ncomp,
34  amrex::Array4<amrex::EBCellFlag const> const& flag_arr,
35  amrex::Array4<amrex::Real const> const& vfrac,
36  const amrex::Geometry & geom,
37  bool use_wts_in_divnc);
38 
39  // Flux redistribution scheme for case where C/F interface crosses EB
40  void amrex_flux_redistribute (const amrex::Box& bx,
41  amrex::Array4<amrex::Real > const& dqdt,
42  amrex::Array4<amrex::Real const> const& divc,
43  amrex::Array4<amrex::Real const> const& wt,
44  amrex::Array4<amrex::Real const> const& vfrac,
45  amrex::Array4<amrex::EBCellFlag const> const& flag,
46  int as_crse,
47  amrex::Array4<amrex::Real > const& rr_drho_crse,
48  amrex::Array4<int const> const& rr_flag_crse,
49  int as_fine,
50  amrex::Array4<amrex::Real > const& dm_as_fine,
51  amrex::Array4<int const> const& levmsk,
52  const amrex::Geometry & geom,
53  bool use_wts_in_divnc,
54  int level_mask_not_covered,
55  int icomp, int ncomp, amrex::Real dt);
56 
57  // Interface to redistribution schemes that only calls single-level routines
58  void ApplyRedistribution ( amrex::Box const& bx, int ncomp,
59  amrex::Array4<amrex::Real> const& dUdt_out,
60  amrex::Array4<amrex::Real> const& dUdt_in,
61  amrex::Array4<amrex::Real const> const& U_in,
62  amrex::Array4<amrex::Real> const& scratch,
63  amrex::Array4<amrex::EBCellFlag const> const& flag,
64  AMREX_D_DECL(amrex::Array4<amrex::Real const> const& apx,
65  amrex::Array4<amrex::Real const> const& apy,
66  amrex::Array4<amrex::Real const> const& apz),
67  amrex::Array4<amrex::Real const> const& vfrac,
68  AMREX_D_DECL(amrex::Array4<amrex::Real const> const& fcx,
69  amrex::Array4<amrex::Real const> const& fcy,
70  amrex::Array4<amrex::Real const> const& fcz),
71  amrex::Array4<amrex::Real const> const& ccc,
72  amrex::BCRec const* d_bcrec_ptr,
73  amrex::Geometry const& lev_geom,
74  amrex::Real dt, std::string const& redistribution_type,
75  bool use_wts_in_divnc = false,
76  int srd_max_order = 2,
77  amrex::Real target_volfrac = 0.5_rt,
78  amrex::Array4<amrex::Real const> const& update_scale={});
79 
80  // Interface to redistribution schemes that calls multi-level routines
81  void ApplyMLRedistribution (
82  amrex::Box const& bx, int ncomp,
83  amrex::Array4<amrex::Real > const& dUdt_out,
84  amrex::Array4<amrex::Real > const& dUdt_in,
85  amrex::Array4<amrex::Real const> const& U_in,
86  amrex::Array4<amrex::Real > const& scratch,
87  amrex::Array4<amrex::EBCellFlag const> const& flag,
88  AMREX_D_DECL(amrex::Array4<amrex::Real const> const& apx,
89  amrex::Array4<amrex::Real const> const& apy,
90  amrex::Array4<amrex::Real const> const& apz),
91  amrex::Array4<amrex::Real const> const& vfrac,
92  AMREX_D_DECL(amrex::Array4<amrex::Real const> const& fcx,
93  amrex::Array4<amrex::Real const> const& fcy,
94  amrex::Array4<amrex::Real const> const& fcz),
95  amrex::Array4<amrex::Real const> const& ccc,
96  amrex::BCRec const* d_bcrec_ptr,
97  amrex::Geometry const& lev_geom, amrex::Real dt,
98  std::string const& redistribution_type,
99  int as_crse,
100  amrex::Array4<amrex::Real > const& rr_drho_crse,
101  amrex::Array4<int const> const& rr_flag_crse,
102  int as_fine,
103  amrex::Array4<amrex::Real > const& dm_as_fine,
104  amrex::Array4<int const> const& levmsk,
105  int level_mask_not_covered,
106  amrex::Real fac_for_deltaR = 1.0_rt,
107  bool use_wts_in_divnc = false,
108  int icomp = 0,
109  int srd_max_order = 2,
110  amrex::Real target_volfrac = 0.5_rt,
111  amrex::Array4<amrex::Real const> const& update_scale={});
112 
113  void ApplyInitialRedistribution (
114  amrex::Box const& bx, int ncomp,
115  amrex::Array4<amrex::Real > const& U_out,
116  amrex::Array4<amrex::Real > const& U_in,
117  amrex::Array4<amrex::EBCellFlag const> const& flag,
118  AMREX_D_DECL(amrex::Array4<amrex::Real const> const& apx,
119  amrex::Array4<amrex::Real const> const& apy,
120  amrex::Array4<amrex::Real const> const& apz),
121  amrex::Array4<amrex::Real const> const& vfrac,
122  AMREX_D_DECL(amrex::Array4<amrex::Real const> const& fcx,
123  amrex::Array4<amrex::Real const> const& fcy,
124  amrex::Array4<amrex::Real const> const& fcz),
125  amrex::Array4<amrex::Real const> const& ccc,
126  amrex::BCRec const* d_bcrec_ptr,
127  amrex::Geometry const& geom, std::string const& redistribution_type,
128  int srd_max_order = 2,
129  amrex::Real target_volfrac = 0.5_rt);
130 
131  void StateRedistribute ( amrex::Box const& bx, int ncomp,
132  amrex::Array4<amrex::Real> const& U_out,
133  amrex::Array4<amrex::Real> const& U_in,
134  amrex::Array4<amrex::EBCellFlag const> const& flag,
135  amrex::Array4<amrex::Real const> const& vfrac,
136  AMREX_D_DECL(amrex::Array4<amrex::Real const> const& fcx,
137  amrex::Array4<amrex::Real const> const& fcy,
138  amrex::Array4<amrex::Real const> const& fcz),
139  amrex::Array4<amrex::Real const> const& ccent,
140  amrex::BCRec const* d_bcrec_ptr,
141  amrex::Array4<int const> const& itracker,
142  amrex::Array4<amrex::Real const> const& nrs,
143  amrex::Array4<amrex::Real const> const& alpha,
144  amrex::Array4<amrex::Real const> const& nbhd_vol,
145  amrex::Array4<amrex::Real const> const& cent_hat,
146  amrex::Geometry const& geom,
147  int max_order = 2);
148 
149  void MLStateRedistribute ( amrex::Box const& bx, int ncomp,
150  amrex::Array4<amrex::Real> const& U_out,
151  amrex::Array4<amrex::Real> const& U_in,
152  amrex::Array4<amrex::EBCellFlag const> const& flag,
153  amrex::Array4<amrex::Real const> const& vfrac,
154  AMREX_D_DECL(amrex::Array4<amrex::Real const> const& fcx,
155  amrex::Array4<amrex::Real const> const& fcy,
156  amrex::Array4<amrex::Real const> const& fcz),
157  amrex::Array4<amrex::Real const> const& ccent,
158  amrex::BCRec const* d_bcrec_ptr,
159  amrex::Array4<int const> const& itracker,
160  amrex::Array4<amrex::Real const> const& nrs,
161  amrex::Array4<amrex::Real const> const& alpha,
162  amrex::Array4<amrex::Real const> const& nbhd_vol,
163  amrex::Array4<amrex::Real const> const& cent_hat,
164  amrex::Geometry const& geom,
165  int as_crse,
166  Array4<Real > const& drho_as_crse,
167  Array4<int const> const& flag_as_crse,
168  int as_fine,
169  Array4<Real > const& dm_as_fine,
170  Array4<int const> const& levmsk,
171  int is_ghost_cell,
172  amrex::Real fac_for_deltaR,
173  int max_order = 2);
174 
175  void MakeITracker ( amrex::Box const& bx,
176  AMREX_D_DECL(amrex::Array4<amrex::Real const> const& apx,
177  amrex::Array4<amrex::Real const> const& apy,
178  amrex::Array4<amrex::Real const> const& apz),
179  amrex::Array4<amrex::Real const> const& vfrac,
180  amrex::Array4<int> const& itracker,
181  amrex::Geometry const& geom,
182  amrex::Real target_volfrac);
183 
184  void MakeStateRedistUtils ( amrex::Box const& bx,
185  amrex::Array4<amrex::EBCellFlag const> const& flag,
186  amrex::Array4<amrex::Real const> const& vfrac,
187  amrex::Array4<amrex::Real const> const& ccent,
188  amrex::Array4< int const> const& itracker,
189  amrex::Array4<amrex::Real> const& nrs,
190  amrex::Array4<amrex::Real> const& alpha,
191  amrex::Array4<amrex::Real> const& nbhd_vol,
192  amrex::Array4<amrex::Real> const& cent_hat,
193  amrex::Geometry const& geom,
194  amrex::Real target_volfrac);
195 
196 }
197 
198 #endif
AMReX_BCRec.H
AMReX_EB2.H
AMReX_MultiCutFab.H
AMReX_MultiFabUtil.H
AMReX_MultiFab.H
AMREX_D_DECL
#define AMREX_D_DECL(a, b, c)
Definition: AMReX_SPACE.H:104
amrex::BCRec
Boundary Condition Records. Necessary information and functions for computing boundary conditions.
Definition: AMReX_BCRec.H:17
amrex::Geometry
Rectangular problem domain geometry.
Definition: AMReX_Geometry.H:73
amrex::MultiFab
A collection (stored as an array) of FArrayBox objects.
Definition: AMReX_MultiFab.H:38
amrex
Definition: AMReX_Amr.cpp:49
amrex::single_level_redistribute
void single_level_redistribute(amrex::MultiFab &div_tmp_in, amrex::MultiFab &div_out, int div_comp, int ncomp, const amrex::Geometry &geom)
amrex::MakeStateRedistUtils
void MakeStateRedistUtils(amrex::Box const &bx, amrex::Array4< amrex::EBCellFlag const > const &flag, amrex::Array4< amrex::Real const > const &vfrac, amrex::Array4< amrex::Real const > const &ccent, amrex::Array4< int const > const &itracker, amrex::Array4< amrex::Real > const &nrs, amrex::Array4< amrex::Real > const &alpha, amrex::Array4< amrex::Real > const &nbhd_vol, amrex::Array4< amrex::Real > const &cent_hat, amrex::Geometry const &geom, amrex::Real target_volfrac)
amrex::amrex_flux_redistribute
void amrex_flux_redistribute(const amrex::Box &bx, amrex::Array4< amrex::Real > const &dqdt, amrex::Array4< amrex::Real const > const &divc, amrex::Array4< amrex::Real const > const &wt, amrex::Array4< amrex::Real const > const &vfrac, amrex::Array4< amrex::EBCellFlag const > const &flag, int as_crse, amrex::Array4< amrex::Real > const &rr_drho_crse, amrex::Array4< int const > const &rr_flag_crse, int as_fine, amrex::Array4< amrex::Real > const &dm_as_fine, amrex::Array4< int const > const &levmsk, const amrex::Geometry &geom, bool use_wts_in_divnc, int level_mask_not_covered, int icomp, int ncomp, amrex::Real dt)
amrex::MLStateRedistribute
void MLStateRedistribute(amrex::Box const &bx, int ncomp, amrex::Array4< amrex::Real > const &U_out, amrex::Array4< amrex::Real > const &U_in, amrex::Array4< amrex::EBCellFlag const > const &flag, amrex::Array4< amrex::Real const > const &vfrac, AMREX_D_DECL(amrex::Array4< amrex::Real const > const &fcx, amrex::Array4< amrex::Real const > const &fcy, amrex::Array4< amrex::Real const > const &fcz), amrex::Array4< amrex::Real const > const &ccent, amrex::BCRec const *d_bcrec_ptr, amrex::Array4< int const > const &itracker, amrex::Array4< amrex::Real const > const &nrs, amrex::Array4< amrex::Real const > const &alpha, amrex::Array4< amrex::Real const > const &nbhd_vol, amrex::Array4< amrex::Real const > const &cent_hat, amrex::Geometry const &geom, int as_crse, Array4< Real > const &drho_as_crse, Array4< int const > const &flag_as_crse, int as_fine, Array4< Real > const &dm_as_fine, Array4< int const > const &levmsk, int is_ghost_cell, amrex::Real fac_for_deltaR, int max_order=2)
amrex::ApplyMLRedistribution
void ApplyMLRedistribution(amrex::Box const &bx, int ncomp, amrex::Array4< amrex::Real > const &dUdt_out, amrex::Array4< amrex::Real > const &dUdt_in, amrex::Array4< amrex::Real const > const &U_in, amrex::Array4< amrex::Real > const &scratch, amrex::Array4< amrex::EBCellFlag const > const &flag, AMREX_D_DECL(amrex::Array4< amrex::Real const > const &apx, amrex::Array4< amrex::Real const > const &apy, amrex::Array4< amrex::Real const > const &apz), amrex::Array4< amrex::Real const > const &vfrac, AMREX_D_DECL(amrex::Array4< amrex::Real const > const &fcx, amrex::Array4< amrex::Real const > const &fcy, amrex::Array4< amrex::Real const > const &fcz), amrex::Array4< amrex::Real const > const &ccc, amrex::BCRec const *d_bcrec_ptr, amrex::Geometry const &lev_geom, amrex::Real dt, std::string const &redistribution_type, int as_crse, amrex::Array4< amrex::Real > const &rr_drho_crse, amrex::Array4< int const > const &rr_flag_crse, int as_fine, amrex::Array4< amrex::Real > const &dm_as_fine, amrex::Array4< int const > const &levmsk, int level_mask_not_covered, amrex::Real fac_for_deltaR=1.0_rt, bool use_wts_in_divnc=false, int icomp=0, int srd_max_order=2, amrex::Real target_volfrac=0.5_rt, amrex::Array4< amrex::Real const > const &update_scale={})
amrex::MakeITracker
void MakeITracker(amrex::Box const &bx, AMREX_D_DECL(amrex::Array4< amrex::Real const > const &apx, amrex::Array4< amrex::Real const > const &apy, amrex::Array4< amrex::Real const > const &apz), amrex::Array4< amrex::Real const > const &vfrac, amrex::Array4< int > const &itracker, amrex::Geometry const &geom, amrex::Real target_volfrac)
amrex::StateRedistribute
void StateRedistribute(amrex::Box const &bx, int ncomp, amrex::Array4< amrex::Real > const &U_out, amrex::Array4< amrex::Real > const &U_in, amrex::Array4< amrex::EBCellFlag const > const &flag, amrex::Array4< amrex::Real const > const &vfrac, AMREX_D_DECL(amrex::Array4< amrex::Real const > const &fcx, amrex::Array4< amrex::Real const > const &fcy, amrex::Array4< amrex::Real const > const &fcz), amrex::Array4< amrex::Real const > const &ccent, amrex::BCRec const *d_bcrec_ptr, amrex::Array4< int const > const &itracker, amrex::Array4< amrex::Real const > const &nrs, amrex::Array4< amrex::Real const > const &alpha, amrex::Array4< amrex::Real const > const &nbhd_vol, amrex::Array4< amrex::Real const > const &cent_hat, amrex::Geometry const &geom, int max_order=2)
amrex::eb_covered_val
static constexpr amrex::Real eb_covered_val
Definition: AMReX_EB_Redistribution.H:17
amrex::single_level_weighted_redistribute
void single_level_weighted_redistribute(amrex::MultiFab &div_tmp_in, amrex::MultiFab &div_out, const amrex::MultiFab &weights, int div_comp, int ncomp, const amrex::Geometry &geom, bool use_wts_in_divnc)
amrex::ApplyRedistribution
void ApplyRedistribution(amrex::Box const &bx, int ncomp, amrex::Array4< amrex::Real > const &dUdt_out, amrex::Array4< amrex::Real > const &dUdt_in, amrex::Array4< amrex::Real const > const &U_in, amrex::Array4< amrex::Real > const &scratch, amrex::Array4< amrex::EBCellFlag const > const &flag, AMREX_D_DECL(amrex::Array4< amrex::Real const > const &apx, amrex::Array4< amrex::Real const > const &apy, amrex::Array4< amrex::Real const > const &apz), amrex::Array4< amrex::Real const > const &vfrac, AMREX_D_DECL(amrex::Array4< amrex::Real const > const &fcx, amrex::Array4< amrex::Real const > const &fcy, amrex::Array4< amrex::Real const > const &fcz), amrex::Array4< amrex::Real const > const &ccc, amrex::BCRec const *d_bcrec_ptr, amrex::Geometry const &lev_geom, amrex::Real dt, std::string const &redistribution_type, bool use_wts_in_divnc=false, int srd_max_order=2, amrex::Real target_volfrac=0.5_rt, amrex::Array4< amrex::Real const > const &update_scale={})
amrex::ApplyInitialRedistribution
void ApplyInitialRedistribution(amrex::Box const &bx, int ncomp, amrex::Array4< amrex::Real > const &U_out, amrex::Array4< amrex::Real > const &U_in, amrex::Array4< amrex::EBCellFlag const > const &flag, AMREX_D_DECL(amrex::Array4< amrex::Real const > const &apx, amrex::Array4< amrex::Real const > const &apy, amrex::Array4< amrex::Real const > const &apz), amrex::Array4< amrex::Real const > const &vfrac, AMREX_D_DECL(amrex::Array4< amrex::Real const > const &fcx, amrex::Array4< amrex::Real const > const &fcy, amrex::Array4< amrex::Real const > const &fcz), amrex::Array4< amrex::Real const > const &ccc, amrex::BCRec const *d_bcrec_ptr, amrex::Geometry const &geom, std::string const &redistribution_type, int srd_max_order=2, amrex::Real target_volfrac=0.5_rt)
amrex::apply_flux_redistribution
void apply_flux_redistribution(const amrex::Box &bx, amrex::Array4< amrex::Real > const &div, amrex::Array4< amrex::Real const > const &divc, amrex::Array4< amrex::Real const > const &wt, int icomp, int ncomp, amrex::Array4< amrex::EBCellFlag const > const &flag_arr, amrex::Array4< amrex::Real const > const &vfrac, const amrex::Geometry &geom, bool use_wts_in_divnc)
amrex::Array4
Definition: AMReX_Array4.H:61