Merge pull request #275 from avik-pal/patch-1

LM needs to be `_vcat` not `vcat`
SciML · Nov 3, 2023 · b8d43a3 · b8d43a3 · avik-pal · Nov 3, 2023
2 parents a607868 + e634dc0
commit b8d43a3
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/levenberg.jl b/src/levenberg.jl
@@ -356,7 +356,7 @@ function perform_step!(cache::LevenbergMarquardtCache{false, fastls}) where {fas
 
     # Usual Levenberg-Marquardt step ("velocity").
     if fastls
-        cache.mat_tmp = vcat(J, λ * cache.DᵀD)
+        cache.mat_tmp = _vcat(J, λ * cache.DᵀD)
         cache.rhs_tmp[1:length(fu1)] .= -_vec(fu1)
         linres = dolinsolve(alg.precs, linsolve; A = cache.mat_tmp,
             b = cache.rhs_tmp, linu = _vec(cache.v), p = p, reltol = cache.abstol)