Closes #305. Add ability to set GSI solver convergence tolerance from mixture object

src/gsi/GasSurfaceInteraction.cpp

@@ -210,6 +210,20 @@ void GasSurfaceInteraction::setIterationsSurfaceBalance(const int& iter)
 
 //==============================================================================
 
+void GasSurfaceInteraction::setWriteConvergenceSurfaceBalance(const bool flag)
+{
+    mp_surf->setWriteConvergenceSurfaceBalance(flag);
+}
+
+//==============================================================================
+
+void GasSurfaceInteraction::setConvergenceToleranceSurfaceBalance(const double& tol)
+{
+    mp_surf->setConvergenceToleranceSurfaceBalance(tol);
+}
+
+//==============================================================================
+
 void GasSurfaceInteraction::getMassBlowingRate(double& mdot){
     mdot = mp_surf->massBlowingRate();
 }

src/gsi/GasSurfaceInteraction.h

@@ -172,6 +172,18 @@ class GasSurfaceInteraction
      */
     void setIterationsSurfaceBalance(const int& iter);
 
+    /**
+     * Function which allows to report convergence history for Newton solves.
+     * Default is false.
+     */
+    void setWriteConvergenceSurfaceBalance(const bool flag);
+
+    /**
+     * Function which allows to set convergence tolerance for Newton solves.
+     * Default value is 1e-8.
+     */
+    void setConvergenceToleranceSurfaceBalance(const double& tol);
+
     /**
      * Function which return the total mass blowing flux.
      *

src/gsi/Surface.h

@@ -185,6 +185,32 @@ class Surface
         << "the surface energy balance!";
     }
 
+//==============================================================================
+
+    /**
+     * Purely virtual function to be called to set the flag to output
+     * convergence history for the Newton solver.
+     */
+    virtual void setWriteConvergenceSurfaceBalance(const bool flag)
+    {
+        throw LogicError()
+        << "setWriteConvergenceSurfaceBalance can be called only when solving "
+        << "the surface energy balance!";
+    }
+
+//==============================================================================
+
+    /**
+     * Purely virtual function to be called to set the tolerance
+     * parameter for the Newton solver.
+     */
+    virtual void setConvergenceToleranceSurfaceBalance(const double& tol)
+    {
+        throw LogicError()
+        << "setConvergenceToleranceSurfaceBalance can be called only when solving "
+        << "the surface energy balance!";
+    }
+
 //==============================================================================
 
     /**

src/gsi/SurfaceBalanceSolverMass.cpp

@@ -187,6 +187,20 @@ class SurfaceBalanceSolverMass :
 
     void setIterationsSurfaceBalance(const int& iter){ setMaxIterations(iter); }
 
+//==============================================================================
+
+    void setWriteConvergenceSurfaceBalance(const bool flag)
+    {
+        setWriteConvergenceHistory(flag);
+    }
+
+//==============================================================================
+
+    void setConvergenceToleranceSurfaceBalance(const double& tol)
+    {
+        setEpsilon(tol);
+    }
+
 //==============================================================================
 
     double massBlowingRate() {

src/gsi/SurfaceBalanceSolverMassEnergy.cpp

@@ -230,6 +230,20 @@ class SurfaceBalanceSolverMassEnergy :
 
     void setIterationsSurfaceBalance(const int& iter){ setMaxIterations(iter); }
 
+//==============================================================================
+
+    void setWriteConvergenceSurfaceBalance(const bool flag)
+    {
+        setWriteConvergenceHistory(flag);
+    }
+
+//==============================================================================
+
+    void setConvergenceToleranceSurfaceBalance(const double& tol)
+    {
+        setEpsilon(tol);
+    }
+
 //==============================================================================
 
     double massBlowingRate()