From 38b4bd0fb2684dc0964a63fb3c6d0f181dbdc802 Mon Sep 17 00:00:00 2001
From: Max Luebke <mluebke@gfz-potsdam.de>
Date: Thu, 27 Jan 2022 09:57:33 +0100
Subject: [PATCH] Update resizing of boundary condition vector

---
 src/BTCSDiffusion.cpp | 17 ++++++-----------
 1 file changed, 6 insertions(+), 11 deletions(-)

diff --git a/src/BTCSDiffusion.cpp b/src/BTCSDiffusion.cpp
index 37f8563..af4c6f7 100644
--- a/src/BTCSDiffusion.cpp
+++ b/src/BTCSDiffusion.cpp
@@ -44,18 +44,13 @@ void BTCSDiffusion::updateInternals() {
     deltas[i] = (double)domain_size[i] / grid_cells[i];
   }
 
-  switch (grid_dim) {
-  case 1:
-    bc.resize(grid_cells[0] + 2, {BTCSDiffusion::BC_CLOSED, 0});
-    break;
-  case 2:
-    bc.resize((grid_cells[0] + 2) * (grid_cells[1] + 2),
-              {BTCSDiffusion::BC_CLOSED, 0});
-    break;
-  case 3:
-    // TODO
-    break;
+  int cells = 1;
+
+  for (int i = 0; i < grid_dim; i++) {
+    cells *= (grid_cells[i] + 2);
   }
+
+  bc.resize(cells, {BTCSDiffusion::BC_CLOSED,0});
 }
 // BTCSDiffusion::BTCSDiffusion(int x) : n_x(x) {
 //   this->grid_dim = 1;