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

wrap SDC+NSE neutrino calls with NEUTRINOS #1400

Merged
merged 2 commits into from
Nov 28, 2023

Conversation

zingale
Copy link
Member

@zingale zingale commented Nov 28, 2023

this will allow them to be disabled at compile-time just like we do for a network

this will allow them to be disabled at compile-time just like
we do for a network
@@ -75,11 +75,13 @@ void sdc_nse_burn(BurnT& state, const Real dt) {
Real abar = state.y[SFX+iabar] / rho_old;
Real zbar = abar * ye_in;

#ifdef NEUTRINOS
constexpr int do_derivatives = 0;
sneut5<do_derivatives>(T_in, rho_old, abar, zbar,
snu, dsnudt, dsnudd, dsnuda, dsnudz);

Real snu_old = snu;
Copy link
Collaborator

Choose a reason for hiding this comment

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

I think you meant to put snu_old outside just like the NSE_NET version below.

Copy link
Member Author

Choose a reason for hiding this comment

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

fixed

@zingale zingale merged commit 0c5ef6a into AMReX-Astro:development Nov 28, 2023
24 checks passed
@zingale zingale deleted the disable_nse_neutrinos branch November 28, 2023 18:13
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