[CL] Dropped time step size setter.

ChemistryLib/ChemicalSolverInterface.h

@@ -20,8 +20,6 @@ public:
     virtual void doWaterChemistryCalculation(
         std::vector<GlobalVector*>& process_solutions, double const dt) = 0;
 
-    virtual void setTimeStepSize(double const dt) = 0;
-
     virtual ~ChemicalSolverInterface() = default;
 };
 }  // namespace ChemistryLib

ChemistryLib/PhreeqcIO.cpp

@@ -87,9 +87,8 @@ void PhreeqcIO::doWaterChemistryCalculation(
 {
     setAqueousSolutionsOrUpdateProcessSolutions(
         process_solutions, Status::SettingAqueousSolutions);
-    setTimeStepSize(dt);
 
-    writeInputsToFile();
+    writeInputsToFile(dt);
 
     execute();
 
@@ -175,7 +174,7 @@ void PhreeqcIO::setAqueousSolutionsOrUpdateProcessSolutions(
     }
 }
 
-void PhreeqcIO::writeInputsToFile()
+void PhreeqcIO::writeInputsToFile(double const dt)
 {
     DBUG("Writing phreeqc inputs into file '%s'.", _phreeqc_input_file.c_str());
     std::ofstream out(_phreeqc_input_file, std::ofstream::out);
@@ -186,7 +185,7 @@ void PhreeqcIO::writeInputsToFile()
                   _phreeqc_input_file.c_str());
     }
 
-    out << *this;
+    out << (*this << dt);
 
     if (!out)
     {

ChemistryLib/PhreeqcIO.h

@@ -51,9 +51,7 @@ public:
         std::vector<GlobalVector*> const& process_solutions,
         Status const status);
 
-    void setTimeStepSize(double const dt) override { _dt = dt; }
-
-    void writeInputsToFile();
+    void writeInputsToFile(double const dt);
 
     void execute();
 
@@ -67,6 +65,12 @@ public:
     std::string const _phreeqc_input_file;
 
 private:
+    PhreeqcIO& operator<<(int dt)
+    {
+        _dt = dt;
+        return *this;
+    }
+
     std::string const _database;
     std::vector<AqueousSolution> _aqueous_solutions;
     std::vector<EquilibriumPhase> _equilibrium_phases;

ChemistryLib/PhreeqcKernel.h

@@ -47,8 +47,6 @@ public:
     void setAqueousSolutions(
         std::vector<GlobalVector*> const& process_solutions);
 
-    void setTimeStepSize(double const dt) override;
-
     void execute(std::vector<GlobalVector*>& process_solutions);
 
     void updateNodalProcessSolutions(
@@ -70,6 +68,8 @@ private:
         return strcmp(element, "H") == 0 ? true : false;
     }
 
+    void setTimeStepSize(double const dt);
+
     std::map<int, struct master*> _process_id_to_master_map;
     cxxISolution _initial_aqueous_solution;
     std::vector<ReactionRate> const _reaction_rates;