-
Notifications
You must be signed in to change notification settings - Fork 98
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
move the fourth order center <-> average stuff to C++ #884
Conversation
all SDC tests pass: http://groot.astro.sunysb.edu/Castro/test-suite/gfortran/2020-04-27-004/index.html blocked by #885 |
@@ -289,23 +288,29 @@ Castro::construct_mol_hydro_source(Real time, Real dt, MultiFab& A_update) | |||
|
|||
#if AMREX_SPACEDIM >= 2 | |||
// construct the face-center interface states q_fc | |||
AMREX_PARALLEL_FOR_4D(nbx, NQ, i, j, k, n, { | |||
bool test = (n == QGC) || (n == QTEMP); | |||
const int* lo_bc = phys_bc.lo(); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
phys_bc
is sized AMREX_SPACEDIM
so we need to avoid dereferencing it outside those bounds. There are a few other instances of that below.
Source/hydro/Castro_hydro.H
Outdated
void make_cell_center(const Box& bx, | ||
Array4<Real const> const& U, | ||
Array4<Real> const& U_cc, | ||
GpuArray<int, 3>& domlo, GpuArray<int, 3>& domhi); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
domlo and domhi should be const in all of these routines
I think I addressed all the comments |
PR summary
PR checklist
CHANGES
file has been updated