Skip to content

Commit

Permalink
fix tabs and format issue in Castro_react.cpp (#2707)
Browse files Browse the repository at this point in the history
  • Loading branch information
zhichen3 authored Jan 10, 2024
1 parent b14f494 commit 50f4266
Showing 1 changed file with 24 additions and 24 deletions.
48 changes: 24 additions & 24 deletions Source/reactions/Castro_react.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -214,10 +214,10 @@ Castro::react_state(MultiFab& s, MultiFab& r, Real time, Real dt, const int stra

burn_t burn_state;
#ifdef NSE_NET
burn_state.mu_p = U(i,j,k,UMUP);
burn_state.mu_n = U(i,j,k,UMUN);
burn_state.mu_p = U(i,j,k,UMUP);
burn_state.mu_n = U(i,j,k,UMUN);

burn_state.y_e = -1.0_rt;
burn_state.y_e = -1.0_rt;
#endif

#if AMREX_SPACEDIM == 1
Expand Down Expand Up @@ -251,9 +251,9 @@ Castro::react_state(MultiFab& s, MultiFab& r, Real time, Real dt, const int stra

burn_state.rho = U(i,j,k,URHO);

// e is used as an input for some NSE solvers
// e is used as an input for some NSE solvers

burn_state.e = U(i,j,k,UEINT) * rhoInv;
burn_state.e = U(i,j,k,UEINT) * rhoInv;

// this T is consistent with UEINT because we did an EOS call before
// calling this function
Expand Down Expand Up @@ -363,19 +363,19 @@ Castro::react_state(MultiFab& s, MultiFab& r, Real time, Real dt, const int stra
}
}
#ifdef NSE
if (store_omegadot == 1) {
reactions(i,j,k,NumSpec+NumAux+1) = burn_state.nse;
}
else {
reactions(i,j,k,1) = burn_state.nse;
}
if (store_omegadot == 1) {
reactions(i,j,k,NumSpec+NumAux+1) = burn_state.nse;
}
else {
reactions(i,j,k,1) = burn_state.nse;
}
#endif
}

// update the state
#ifdef NSE_NET
U(i,j,k,UMUP) = burn_state.mu_p;
U(i,j,k,UMUN) = burn_state.mu_n;
U(i,j,k,UMUP) = burn_state.mu_p;
U(i,j,k,UMUN) = burn_state.mu_n;
#endif
for (int n = 0; n < NumSpec; ++n) {
U(i,j,k,UFS+n) = U(i,j,k,URHO) * burn_state.xn[n];
Expand Down Expand Up @@ -554,10 +554,10 @@ Castro::react_state(Real time, Real dt)
#endif

#ifdef NSE_NET
burn_state.mu_p = U_old(i,j,k,UMUP);
burn_state.mu_n = U_old(i,j,k,UMUN);
burn_state.mu_p = U_old(i,j,k,UMUP);
burn_state.mu_n = U_old(i,j,k,UMUN);

burn_state.y_e = -1.0_rt;
burn_state.y_e = -1.0_rt;
#endif
// Initialize some data for later.

Expand Down Expand Up @@ -692,8 +692,8 @@ Castro::react_state(Real time, Real dt)

// update the state data.
#ifdef NSE_NET
U_new(i,j,k,UMUP) = burn_state.mu_p;
U_new(i,j,k,UMUN) = burn_state.mu_n;
U_new(i,j,k,UMUP) = burn_state.mu_p;
U_new(i,j,k,UMUN) = burn_state.mu_n;
#endif
U_new(i,j,k,UEDEN) = burn_state.y[SEDEN];
U_new(i,j,k,UEINT) = burn_state.y[SEINT];
Expand Down Expand Up @@ -740,12 +740,12 @@ Castro::react_state(Real time, Real dt)
}
}
#ifdef NSE
if (store_omegadot == 1) {
react_src(i,j,k,NumSpec+NumAux+1) = burn_state.nse;
}
else {
react_src(i,j,k,1) = burn_state.nse;
}
if (store_omegadot == 1) {
react_src(i,j,k,NumSpec+NumAux+1) = burn_state.nse;
}
else {
react_src(i,j,k,1) = burn_state.nse;
}
#endif
}

Expand Down

0 comments on commit 50f4266

Please sign in to comment.