Skip to content
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

MAC projector with Robin BCs needs initialization every time. The #125

Merged
merged 2 commits into from
Mar 4, 2024

Conversation

cgilet
Copy link
Collaborator

@cgilet cgilet commented Mar 4, 2024

underlying solver, (EB)MLABecLaplacian, is not safe for reuse with Robin BCs.

underlying solver, (EB)MLABecLaplacian, is not save for reuse with
Robin BCs.
@cgilet cgilet requested a review from WeiqunZhang March 4, 2024 20:24
@@ -212,6 +212,7 @@ MacProjector::setLevelBC (int amrlev, const MultiFab* levelbcdata, const MultiFa
"setDomainBC must be called before setLevelBC");
m_linop->setLevelBC(amrlev, levelbcdata, robin_a, robin_b, robin_f);
m_needs_level_bcs[amrlev] = false;
m_needs_init = true; // Solver is not safe for reuse with Robin BC
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Maybe we should only do this when the robin_* pointers are non null.

Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Simply checking on robin_a should work right? MLABecLap will fail if only robin_a is non-null, right?

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Yes, I think so.

@WeiqunZhang WeiqunZhang merged commit 111a42c into AMReX-Fluids:development Mar 4, 2024
13 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants