-
Notifications
You must be signed in to change notification settings - Fork 41
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
Generalized multi-level relaxation #1233
Conversation
…of an affect on the boundaries than expected.
…rners. The mask seems OK but the logic is hard to follow.
…will be overwritten with FIP. Need to remove G2E from fast_rhs.
@WeiqunZhang can you focus on the ERF_FillPatcher.cpp/.H files? The new code for generating the mask array for set and relax regions is here. |
com_ba.define(com_bl); | ||
com_ba.complementIn(com_bl, fba_bnd); | ||
|
||
// Fill mask based upon the com_bl BoxList |
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.
There is a more efficient way to do this. Since this pattern of setVal in a BoxArray's region can be useful in other places, I will add a function to amrex. After that is merged, we can do a new PR to use it.
Co-authored-by: Weiqun Zhang <[email protected]>
Co-authored-by: Weiqun Zhang <[email protected]>
Co-authored-by: Weiqun Zhang <[email protected]>
Building the mask array in ERFFillPatcher.cpp is a bit off a mess. The code works but it is hard to follow and not perfectly general.
The tradeoff between space and time interpolation has not been assessed. We currently do both for all the fill calls; this should be modified (currently inefficient).
Grids to evolve has been removed. With the fast integrator, we will compute an update on all interior cells; but the set region is overwritten because we call apply_bcs() at the end of a fast iteration.