From 3562f30efaaa6e7c5abc5d46622e218e8d6d95f9 Mon Sep 17 00:00:00 2001
From: Max Luebke <mluebke@gfz-potsdam.de>
Date: Thu, 2 Dec 2021 10:51:47 +0100
Subject: [PATCH] change dx to 1/n-1

---
 src/BTCSDiffusion.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/BTCSDiffusion.cpp b/src/BTCSDiffusion.cpp
index 287fcb7..fcebc92 100644
--- a/src/BTCSDiffusion.cpp
+++ b/src/BTCSDiffusion.cpp
@@ -50,7 +50,7 @@ void BTCSDiffusion::setBoundaryCondition(std::vector<double> input,
 void BTCSDiffusion::simulate(std::vector<double> &c, std::vector<double> &alpha,
                              double timestep) {
   // calculate dx
-  double dx = 1. / this->dim_x;
+  double dx = 1. / (this->dim_x - 1);
 
   // calculate size needed for A matrix and b,x vectors
   int size = this->dim_x + 2;