From a5e58c80d70978c12d24ff47b2362c9145af325c Mon Sep 17 00:00:00 2001 From: Francois Coppens Date: Thu, 29 Jul 2021 11:48:38 +0200 Subject: [PATCH] Passing argument context instead of creating and passing a locally created context. #25 --- org/qmckl_sherman_morrison_woodbury.org | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/org/qmckl_sherman_morrison_woodbury.org b/org/qmckl_sherman_morrison_woodbury.org index ab7db81..047846c 100644 --- a/org/qmckl_sherman_morrison_woodbury.org +++ b/org/qmckl_sherman_morrison_woodbury.org @@ -5,7 +5,7 @@ Low- and high-level functions that use the Sherman-Morrison and Woodbury matrix inversion formulas to update the inverse of a non-singualr matrix - + * Headers #+begin_src elisp :noexport :results none :exports none (org-babel-lob-ingest "../tools/lib.org") @@ -255,6 +255,7 @@ end function qmckl_sherman_morrison_f #+begin_src c :tangle (eval c) :comments org #include #include "qmckl.h" +#include qmckl_exit_code qmckl_sherman_morrison_c(const qmckl_context context, const uint64_t Dim, @@ -303,6 +304,7 @@ qmckl_exit_code qmckl_sherman_morrison_c(const qmckl_context context, l += 1; } + return QMCKL_SUCCESS; } @@ -911,8 +913,6 @@ qmckl_exit_code qmckl_sherman_morrison_splitting_c(const qmckl_context context, // std::cerr << "Called qmckl_sherman_morrison_splitting with " << N_updates << " updates" << std::endl; // #endif - qmckl_context local_context; - local_context = qmckl_context_create(); qmckl_exit_code rc; double later_updates[Dim * N_updates]; @@ -923,7 +923,7 @@ qmckl_exit_code qmckl_sherman_morrison_splitting_c(const qmckl_context context, breakdown, Slater_inv, later_updates, later_index, &later); if (later > 0) { - rc = qmckl_sherman_morrison_splitting_c(local_context, Dim, later, + rc = qmckl_sherman_morrison_splitting_c(context, Dim, later, later_updates, later_index, breakdown, Slater_inv); }