|
4 | 4 |
|
5 | 5 | #include "radiation/radiation_system.hpp" // IWYU pragma: keep |
6 | 6 |
|
7 | | - |
8 | 7 | // Compute kappaE and kappaP based on the opacity model. The result is stored in the last five arguments: alpha_P, alpha_E, kappaP, kappaE, and kappaPoverE. |
9 | 8 | template <typename problem_t> |
10 | 9 | AMREX_GPU_DEVICE auto RadSystem<problem_t>::ComputeModelDependentKappaEAndKappaP( |
@@ -786,8 +785,8 @@ void RadSystem<problem_t>::AddSourceTermsMultiGroup(array_t &consVar, arrayconst |
786 | 785 | work = updated_energy.work; |
787 | 786 |
|
788 | 787 | // Check for convergence of the work term |
789 | | - auto const Egastot1 = quokka::EOS<problem_t>::ComputeEgasFromEint( |
790 | | - rho, updated_flux.gasMomentum[0], updated_flux.gasMomentum[1], updated_flux.gasMomentum[2], Egas_guess, B_cc); |
| 788 | + auto const Egastot1 = quokka::EOS<problem_t>::ComputeEgasFromEint(rho, updated_flux.gasMomentum[0], updated_flux.gasMomentum[1], |
| 789 | + updated_flux.gasMomentum[2], Egas_guess, B_cc); |
791 | 790 | const double rel_lag_tol = 1.0e-8; |
792 | 791 | const double lag_tol = 1.0e-13; |
793 | 792 | double ref_work = rel_lag_tol * sum(abs(work)); |
@@ -831,8 +830,7 @@ void RadSystem<problem_t>::AddSourceTermsMultiGroup(array_t &consVar, arrayconst |
831 | 830 | if constexpr (gamma_ != 1.0) { |
832 | 831 | Egas_guess = Egas0 + (Egas_guess - Egas0) * gas_update_factor; |
833 | 832 | consNew(i, j, k, gasInternalEnergy_index) = Egas_guess; |
834 | | - consNew(i, j, k, gasEnergy_index) = |
835 | | - quokka::EOS<problem_t>::ComputeEgasFromEint(rho, x1GasMom1, x2GasMom1, x3GasMom1, Egas_guess, B_cc); |
| 833 | + consNew(i, j, k, gasEnergy_index) = quokka::EOS<problem_t>::ComputeEgasFromEint(rho, x1GasMom1, x2GasMom1, x3GasMom1, Egas_guess, B_cc); |
836 | 834 | } else { |
837 | 835 | amrex::ignore_unused(Egas_guess); |
838 | 836 | amrex::ignore_unused(Egas0); |
|
0 commit comments