Skip to content

Commit

Permalink
Replace tabs with spaces
Browse files Browse the repository at this point in the history
  • Loading branch information
yut23 committed Oct 5, 2023
1 parent e3e866c commit 05cbbda
Show file tree
Hide file tree
Showing 15 changed files with 429 additions and 429 deletions.
6 changes: 3 additions & 3 deletions integration/VODE/vode_dvstep.H
Original file line number Diff line number Diff line change
Expand Up @@ -225,9 +225,9 @@ int dvstep (BurnT& state, DvodeT& vstate)
bool valid_update = true;

#ifdef SDC
if (vstate.y(SEINT+1) < 0.0_rt) {
valid_update = false;
}
if (vstate.y(SEINT+1) < 0.0_rt) {
valid_update = false;
}

Real rho_current = state.rho_orig + vstate.tn * state.ydot_a[SRHO];

Expand Down
2 changes: 1 addition & 1 deletion interfaces/burner.H
Original file line number Diff line number Diff line change
Expand Up @@ -63,7 +63,7 @@ void burner (BurnT& state, Real dt)
#ifdef NSE_TABLE
if (in_nse(state, true) && state.success == false && dt_remaining > 0.0) {
#else
if (in_nse(state, nse_skip_molar) && state.success == false && dt_remaining > 0.0) {
if (in_nse(state, nse_skip_molar) && state.success == false && dt_remaining > 0.0) {
#endif

#ifndef AMREX_USE_GPU
Expand Down
38 changes: 19 additions & 19 deletions nse_solver/make_table/burn_cell.H
Original file line number Diff line number Diff line change
Expand Up @@ -31,32 +31,32 @@ void burn_cell_c()
state.rho = rho;
state.y_e = Ye;

if (state.y_e > 0.52_rt){
state.mu_p = -1.0_rt;
state.mu_n = -16.0_rt;
}
else if (state.y_e > 0.48_rt){
state.mu_p = -6.0_rt;
state.mu_n = -11.0_rt;
}
else if (state.y_e > 0.4_rt){
state.mu_p = -10.0_rt;
state.mu_n = -7.0_rt;
}
else{
state.mu_p = -18.0;
state.mu_n = -1.0;
}
if (state.y_e > 0.52_rt){
state.mu_p = -1.0_rt;
state.mu_n = -16.0_rt;
}
else if (state.y_e > 0.48_rt){
state.mu_p = -6.0_rt;
state.mu_n = -11.0_rt;
}
else if (state.y_e > 0.4_rt){
state.mu_p = -10.0_rt;
state.mu_n = -7.0_rt;
}
else{
state.mu_p = -18.0;
state.mu_n = -1.0;
}

// find the nse state

const bool assume_ye_is_valid = true;
Real eps = 1.e-10;
use_hybrid_solver = 1;
use_hybrid_solver = 1;

auto nse_state = get_actual_nse_state(state, eps, assume_ye_is_valid);

std::cout << std::scientific;
std::cout << std::scientific;
std::cout << std::setw(20) << state.rho << " "
<< std::setw(20) << state.T << " " << std::fixed
<< std::setw(20) << state.y_e << " "
Expand Down
Loading

0 comments on commit 05cbbda

Please sign in to comment.