From 9c632cf1230ae904e89dc06dfba2b63b5c406585 Mon Sep 17 00:00:00 2001
From: Ray Speth <speth@mit.edu>
Date: Sun, 27 Mar 2016 21:45:45 -0400
Subject: [PATCH] [1D] Remove unused / unimplemented methods from StFlow

---
 include/cantera/oneD/StFlow.h | 32 --------------------------------
 src/oneD/StFlow.cpp           |  6 ------
 2 files changed, 38 deletions(-)

diff --git a/include/cantera/oneD/StFlow.h b/include/cantera/oneD/StFlow.h
index 5ebfbdcd9..0922b45fa 100644
--- a/include/cantera/oneD/StFlow.h
+++ b/include/cantera/oneD/StFlow.h
@@ -65,9 +65,6 @@ public:
         return *m_kin;
     }
 
-    virtual void init() {
-    }
-
     /**
      * Set the thermo manager. Note that the flow equations assume
      * the ideal gas equation.
@@ -242,15 +239,11 @@ public:
         return m_do_energy[j];
     }
 
-    void integrateChem(doublereal* x,doublereal dt);
-
     //! Change the grid size. Called after grid refinement.
     void resize(size_t components, size_t points);
 
     virtual void setFixedPoint(int j0, doublereal t0) {}
 
-    void setJac(MultiJac* jac);
-
     //! Set the gas object state to be consistent with the solution at point j.
     void setGas(const doublereal* x, size_t j);
 
@@ -299,18 +292,6 @@ public:
     }
 
 protected:
-    doublereal component(const doublereal* x, size_t i, size_t j) const {
-        return x[index(i,j)];
-    }
-
-    doublereal conc(const doublereal* x, size_t k,size_t j) const {
-        return Y(x,k,j)*density(j)/m_wt[k];
-    }
-
-    doublereal cbar(const doublereal* x, size_t k, size_t j) const {
-        return std::sqrt(8.0*GasConstant * T(x,j) / (Pi * m_wt[k]));
-    }
-
     doublereal wdot(size_t k, size_t j) const {
         return m_wdot(k,j);
     }
@@ -334,17 +315,6 @@ protected:
         }
     }
 
-    //--------------------------------
-    // central-differenced derivatives
-    //--------------------------------
-
-    doublereal cdif2(const doublereal* x, size_t n, size_t j,
-                     const doublereal* f) const {
-        doublereal c1 = (f[j] + f[j-1])*(x[index(n,j)] - x[index(n,j-1)]);
-        doublereal c2 = (f[j+1] + f[j])*(x[index(n,j+1)] - x[index(n,j)]);
-        return (c2/(z(j+1) - z(j)) - c1/(z(j) - z(j-1)))/(z(j+1) - z(j-1));
-    }
-
     //! @name Solution components
     //! @{
 
@@ -470,8 +440,6 @@ protected:
     Kinetics* m_kin;
     Transport* m_trans;
 
-    MultiJac* m_jac;
-
     // boundary emissivities for the radiation calculations
     doublereal m_epsilon_left;
     doublereal m_epsilon_right;
diff --git a/src/oneD/StFlow.cpp b/src/oneD/StFlow.cpp
index face1523d..e5257eed8 100644
--- a/src/oneD/StFlow.cpp
+++ b/src/oneD/StFlow.cpp
@@ -19,7 +19,6 @@ StFlow::StFlow(IdealGasPhase* ph, size_t nsp, size_t points) :
     m_thermo(0),
     m_kin(0),
     m_trans(0),
-    m_jac(0),
     m_epsilon_left(0.0),
     m_epsilon_right(0.0),
     m_do_soret(false),
@@ -820,11 +819,6 @@ XML_Node& StFlow::save(XML_Node& o, const doublereal* const sol)
     return flow;
 }
 
-void StFlow::setJac(MultiJac* jac)
-{
-    m_jac = jac;
-}
-
 void AxiStagnFlow::evalRightBoundary(doublereal* x, doublereal* rsd,
                                      integer* diag, doublereal rdt)
 {