diff --git a/include/cantera/numerics/BandMatrix.h b/include/cantera/numerics/BandMatrix.h
index 80bb60e33..65b0a4909 100644
--- a/include/cantera/numerics/BandMatrix.h
+++ b/include/cantera/numerics/BandMatrix.h
@@ -301,6 +301,15 @@ public:
      */
     virtual size_t checkColumns(doublereal& valueSmall) const;
 
+    //! Change the way the matrix is factored
+    /*!
+     *  @param fAlgorithm   integer
+     *                   0 LU factorization
+     *                   1 QR factorization
+     */
+    virtual void useFactorAlgorithm(int fAlgorithm);
+
+
 protected:
 
     //! Matrix data
diff --git a/src/numerics/BandMatrix.cpp b/src/numerics/BandMatrix.cpp
index 3601d5b6e..95c2a7d49 100644
--- a/src/numerics/BandMatrix.cpp
+++ b/src/numerics/BandMatrix.cpp
@@ -441,4 +441,10 @@ void BandMatrix::copyData(const GeneralMatrix& y)
     (void) memcpy(DATA_PTR(data), yyPtr->ptrColumn(0), n);
 }
 
+void BandMatrix::useFactorAlgorithm(int fAlgorithm)
+{
+//    useQR_ = fAlgorithm;
+}
+
+
 }