A bit more clean up

jchristopherson · Feb 28, 2024 · 5dfab86 · 5dfab86
1 parent 7031790
commit 5dfab86
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/src/levenberg_marquardt.f90 b/src/levenberg_marquardt.f90
@@ -1,7 +1,8 @@
 submodule (fstats) levenberg_marquardt
 ! REFERENCES:
 ! 1. https://people.duke.edu/~hpgavin/ExperimentalSystems/lm.pdf
-    use linalg, only : mtx_inverse, lu_factor, mtx_mult
+    use linalg, only : mtx_inverse, lu_factor, mtx_mult, solve_lu, &
+        diag_mtx_mult, rank1_update
     use fstats_errors
     implicit none
 contains
@@ -453,7 +454,7 @@ subroutine broyden_update(pOld, yOld, jac, p, y, dp, dy)
         dp = p - pOld
         h2 = dot_product(dp, dp)
         dy = y - yOld - matmul(jac, dp)
-        call recip_mult_array(h2, dy)   ! compute dy / h2
+        dy = dy / h2
         call rank1_update(1.0d0, dy, dp, jac)
     end subroutine