-
Notifications
You must be signed in to change notification settings - Fork 2
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
Add eos information to floors, misc UserWorkBeforeOutput fix #210
Conversation
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.
I think i see why you're getting segfaults. See comments below. After making those changese everywhere relevant (I only highlighted once), give it a try on device again.
src/fluid/fluid.cpp
Outdated
auto MaxDensity = [](MeshData<Real> *md) { | ||
return ReduceOneVar<Kokkos::Max<Real>>(md, fluid_prim::density::name(), 0); | ||
}; | ||
|
||
hst_vars.emplace_back(HistoryOutputVar(HstMax, MaxDensity, "maximum density")); |
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.
Looks like this branch might be out of date with main... double check you and @mari2895 are in sync so your PRs don't undo each other.
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.
Looks like this is the same as current main, unless I missed something.
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.
Why is there a diff showing here then? Odd.
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.
Oh I see it now... not sure why it's different from main, don't think I did this myself..
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.
maybe just merge main in again... see if that fixes it?
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.
@AstroBarker Looks the same on my branch too. Don't see anything different.
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.
but why is it deleted? should not be deleted
Ready for re-review, @AstroBarker ? |
Good to go |
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.
Two tiny nitpicks. Otherwise this is good to go.
src/fluid/fluid.cpp
Outdated
auto MaxDensity = [](MeshData<Real> *md) { | ||
return ReduceOneVar<Kokkos::Max<Real>>(md, fluid_prim::density::name(), 0); | ||
}; | ||
|
||
hst_vars.emplace_back(HistoryOutputVar(HstMax, MaxDensity, "maximum density")); |
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.
Why is there a diff showing here then? Odd.
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.
LGTM, thanks for explaining!
This PR adds knowledge of eos bounds to the floors. Specifically, for tabulated EOS, demand that the floors do not extrapolate off of the table. This required making the
bounds
param mutable.Currently, for ideal gas, the option to enforce user set bounds is supported by supplying, e.g.,
<eos>/rho_min
in the input deck. If unspecified, reasonable defaults are set.Also fixed a misc bug in the
UserWorkBeforeOutput
code where a package and param were pulled out inside a device kernel that caused a crash on device.