From 9ec382877ec1879a9fdb24bd61946e6a840eaf50 Mon Sep 17 00:00:00 2001
From: Max Luebke <mluebke@gfz-potsdam.de>
Date: Tue, 1 Mar 2022 11:19:06 +0100
Subject: [PATCH] Fix function parameters.

- Use private data types instead of plain Eigen types
---
 src/BTCSDiffusion.cpp |  6 +++---
 src/BTCSDiffusion.hpp | 10 +++++-----
 2 files changed, 8 insertions(+), 8 deletions(-)

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);