From af64b91ca6d78d50b4c23994e95dd6954e2ff8b8 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Fran=C3=A7ois=20Coppens?= <mail@fmgjcoppens.nl>
Date: Fri, 16 Apr 2021 17:33:15 +0200
Subject: [PATCH] Started debugging Slagel-splitting in MaponiA3.

---
 src/SM_MaponiA3S.cpp | 9 +++++----
 1 file changed, 5 insertions(+), 4 deletions(-)

diff --git a/src/SM_MaponiA3S.cpp b/src/SM_MaponiA3S.cpp
index bb02872..634ba3a 100644
--- a/src/SM_MaponiA3S.cpp
+++ b/src/SM_MaponiA3S.cpp
@@ -129,10 +129,6 @@ void MaponiA3S(double *Slater_inv, unsigned int Dim, unsigned int N_updates,
   }
   memcpy(Slater_inv, last, Dim * Dim * sizeof(double));
 
-  if (later > 0) {
-    MaponiA3S(Slater_inv, Dim, later, later_updates, later_index);
-  }
-
   /*
   CLEANUP MEMORY
   */
@@ -144,6 +140,11 @@ void MaponiA3S(double *Slater_inv, unsigned int Dim, unsigned int N_updates,
     delete[] ylk[l];
   }
   delete[] Al, next, p;
+
+  if (later > 0) {
+    std::cout << "Entering recursive loop with " << l << " updates" << std::endl;
+    MaponiA3S(Slater_inv, Dim, later, later_updates, later_index);
+  }
 }
 
 extern "C" {