diff --git a/src/BTCSDiffusion.cpp b/src/BTCSDiffusion.cpp
index 71ff04b..f18916b 100644
--- a/src/BTCSDiffusion.cpp
+++ b/src/BTCSDiffusion.cpp
@@ -251,8 +251,8 @@ void Diffusion::BTCSDiffusion::simulate2D(
 }
 
 inline void Diffusion::BTCSDiffusion::fillMatrixFromRow(
-    const Eigen::VectorXd &alpha,
-    const Eigen::Vector<Diffusion::boundary_condition, Eigen::Dynamic> &bc,
+    const DVectorRowMajor &alpha,
+    const BCVectorRowMajor &bc,
     int size, double dx, double time_step) {
 
   Diffusion::boundary_condition left = bc[0];
@@ -289,7 +289,7 @@ inline void Diffusion::BTCSDiffusion::fillMatrixFromRow(
 }
 
 inline void Diffusion::BTCSDiffusion::fillVectorFromRowADI(
-    const DVectorRowMajor &c, const Eigen::VectorXd alpha,
+    const DVectorRowMajor &c, const DVectorRowMajor alpha,
     const BCVectorRowMajor &bc, const DVectorRowMajor &t0_c, int size,
     double dx, double time_step) {
 
diff --git a/src/BTCSDiffusion.hpp b/src/BTCSDiffusion.hpp
index ff64189..fcbab40 100644
--- a/src/BTCSDiffusion.hpp
+++ b/src/BTCSDiffusion.hpp
@@ -128,12 +128,12 @@ private:
                   Eigen::Map<const DMatrixRowMajor> &alpha,
                   Eigen::Map<const BCMatrixRowMajor> &bc);
 
-  inline void fillMatrixFromRow(
-      const Eigen::VectorXd &alpha,
-      const Eigen::Vector<Diffusion::boundary_condition, Eigen::Dynamic> &bc,
-      int size, double dx, double time_step);
+  inline void fillMatrixFromRow(const DVectorRowMajor &alpha,
+                                const BCVectorRowMajor &bc, int size, double dx,
+                                double time_step);
+
   inline void fillVectorFromRowADI(const DVectorRowMajor &c,
-                                   const Eigen::VectorXd alpha,
+                                   const DVectorRowMajor alpha,
                                    const BCVectorRowMajor &bc,
                                    const DVectorRowMajor &t0_c, int size,
                                    double dx, double time_step);