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

Update rhoX for nse_solver for SDC and discard strang #1607

Merged
merged 15 commits into from
Jul 11, 2024

Conversation

zhichen3
Copy link
Collaborator

@zhichen3 zhichen3 commented Jul 9, 2024

This is to address #1606 and partially addresses #1542

Also some clean-ups to nse-solver. It mainly stores state.y[SFS+n] at the end.

@zhichen3
Copy link
Collaborator Author

zhichen3 commented Jul 9, 2024

overflow in hybrj

@zhichen3 zhichen3 marked this pull request as ready for review July 9, 2024 20:43
@zhichen3 zhichen3 requested a review from zingale July 10, 2024 19:35
@zhichen3
Copy link
Collaborator Author

Mostly just cleaning. I ended up storing state.y[SFS+n] after we solved NSE.

@@ -73,11 +73,12 @@ T get_nonexponent_nse_state(const T& state) {

// find nse mass frac without the exponent term.

amrex::Real power = 2.0 * M_PI * network::mion(n+1) *
C::k_B * T_in / (C::hplanck * C::hplanck);
constexpr amrex::Real ihplanck2 = 1.0_rt / (C::hplanck * C::hplanck);
Copy link
Member

Choose a reason for hiding this comment

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

this can be moved outside of the loop

@zingale zingale merged commit 4142d4f into AMReX-Astro:development Jul 11, 2024
29 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