get rid of the wrapper around step_rate_gauss_seidel

Author: mabruzzo

src/clib/fortran_func_wrappers.hpp

@@ -591,27 +591,6 @@ inline void scale_fields_g(
 
 }
 
-/// Uses one linearly implicit Gauss-Seidel sweep of a backward-Euler time
-/// integrator to advance the rate equations by one (sub-)cycle (dtit).
-inline void step_rate_g(
-  double* dtit, IndexRange idx_range, gr_mask_type anydust, double* h2dust,
-  double* rhoH, double* dedot_prev, double* HIdot_prev,
-  gr_mask_type* itmask, gr_mask_type* itmask_metal, int imetal,
-  chemistry_data* my_chemistry, grackle_field_data* my_fields,
-  photo_rate_storage my_uvb_rates,
-  grackle::impl::GrainSpeciesCollection grain_growth_rates,
-  grackle::impl::SpeciesCollection species_tmpdens,
-  grackle::impl::CollisionalRxnRateCollection kcr_buf,
-  grackle::impl::PhotoRxnRateCollection kshield_buf
-) {
-  grackle::impl::step_rate_gauss_seidel(
-    dtit, anydust, h2dust, rhoH, dedot_prev, HIdot_prev, itmask, itmask_metal,
-    my_chemistry, my_fields, my_uvb_rates, idx_range, grain_growth_rates,
-    species_tmpdens, kcr_buf, kshield_buf
-  );
-}
-
-
 } // namespace grackle::impl::fortran_wrapper
 
 #endif /* FORTRAN_FUNC_WRAPPERS_HPP */

src/clib/solve_rate_cool_g-cpp.cpp

@@ -883,12 +883,12 @@ int solve_rate_cool_g(
           // Solve rate equations with one linearly implicit Gauss-Seidel
           // sweep of a backward Euler method (for all cells specified by
           // itmask_gs)
-          f_wrap::step_rate_g(
+          grackle::impl::step_rate_gauss_seidel(
             dtit.data(), idx_range, anydust, spsolvbuf.h2dust, rhoH.data(),
             spsolvbuf.dedot_prev, spsolvbuf.HIdot_prev, spsolvbuf.itmask_gs,
-            itmask_metal.data(), imetal, my_chemistry, my_fields,
-            *my_uvb_rates, spsolvbuf.grain_growth_rates,
-            spsolvbuf.species_tmpdens, spsolvbuf.kcr_buf, spsolvbuf.kshield_buf
+            itmask_metal.data(), my_chemistry, my_fields, *my_uvb_rates,
+            spsolvbuf.grain_growth_rates, spsolvbuf.species_tmpdens,
+            spsolvbuf.kcr_buf, spsolvbuf.kshield_buf
           );
 
           // Solve rate equations with one linearly implicit Gauss-Seidel

src/clib/step_rate_gauss_seidel.hpp

@@ -32,11 +32,10 @@ namespace grackle::impl {
 /// Uses one linearly implicit Gauss-Seidel sweep of a backward-Euler time
 /// integrator to advance the rate equations by one (sub-)cycle (dtit).
 inline void step_rate_gauss_seidel(
-  double* dtit, gr_mask_type anydust, double* h2dust, double* rhoH,
-  double* dedot_prev, double* HIdot_prev, gr_mask_type* itmask,
+  double* dtit, IndexRange idx_range, gr_mask_type anydust, double* h2dust,
+  double* rhoH, double* dedot_prev, double* HIdot_prev, gr_mask_type* itmask,
   gr_mask_type* itmask_metal, chemistry_data* my_chemistry,
   grackle_field_data* my_fields, photo_rate_storage my_uvb_rates,
-  IndexRange idx_range,
   grackle::impl::GrainSpeciesCollection grain_growth_rates,
   grackle::impl::SpeciesCollection species_tmpdens,
   grackle::impl::CollisionalRxnRateCollection kcol_buf,