Skip to content

Commit

Permalink
some more true-SDC simplication
Browse files Browse the repository at this point in the history
we can basically do what simplified-SDC does now
  • Loading branch information
zingale committed Sep 23, 2023
1 parent ffae4da commit d098295
Show file tree
Hide file tree
Showing 3 changed files with 2 additions and 24 deletions.
3 changes: 0 additions & 3 deletions integration/RKC/rkc_type.H
Original file line number Diff line number Diff line change
Expand Up @@ -26,9 +26,6 @@ struct rkc_t {
// Tolerances
amrex::Real rtol_spec, atol_spec;
amrex::Real rtol_enuc, atol_enuc;
#ifdef TRUE_SDC
amrex::Real rtol_dens, atol_dens;
#endif

// RKC had an info array with 4 elements

Expand Down
4 changes: 1 addition & 3 deletions integration/VODE/vode_type.H
Original file line number Diff line number Diff line change
Expand Up @@ -139,9 +139,7 @@ struct dvode_t
// Tolerances
Real rtol_spec, atol_spec;
Real rtol_enuc, atol_enuc;
#ifdef TRUE_SDC
Real rtol_dens, atol_dens;
#endif

// Local time and integration end time
amrex::Real t, tout;

Expand Down
19 changes: 1 addition & 18 deletions integration/integrator_data.H
Original file line number Diff line number Diff line change
Expand Up @@ -5,18 +5,7 @@

// Define the size of the ODE system that VODE will integrate

#ifdef STRANG
const int INT_NEQS = NumSpec + 1;
#endif

#ifdef SIMPLIFIED_SDC
const int INT_NEQS = SVAR_EVOLVE;
#endif

#ifdef TRUE_SDC
const int INT_NEQS = NumSpec + 1;
#endif


// We will use this parameter to determine if a given species
// abundance is unreasonably small or large (each X must satisfy
Expand Down Expand Up @@ -49,13 +38,7 @@ constexpr int integrator_neqs ()
if constexpr (has_energy<BurnT>::value) {
int_neqs += 1;
}
#endif

#ifdef SIMPLIFIED_SDC
int_neqs = SVAR_EVOLVE;
#endif

#ifdef TRUE_SDC
#else
int_neqs = NumSpec + 1;
#endif

Expand Down

0 comments on commit d098295

Please sign in to comment.