Merge pull request #397 from SciML/missing_vec

Missing vec in Descent cache

Author: avik-pal

Project.toml

@@ -1,7 +1,7 @@
 name = "NonlinearSolve"
 uuid = "8913a72c-1f9b-4ce2-8d82-65094dcecaec"
 authors = ["SciML"]
-version = "3.8.3"
+version = "3.8.4"
 
 [deps]
 ADTypes = "47edcb42-4c32-4615-8424-f2b9edc5f35b"

src/descent/newton.jl

@@ -100,7 +100,7 @@ function __internal_solve!(
     if idx === Val(1)
         @bb cache.JᵀJ_cache = transpose(J) × J
     end
-    @bb cache.Jᵀfu_cache = transpose(J) × fu
+    @bb cache.Jᵀfu_cache = transpose(J) × vec(fu)
     @static_timeit cache.timer "linear solve" begin
         δu = cache.lincache(; A = __maybe_symmetric(cache.JᵀJ_cache), b = cache.Jᵀfu_cache,
             kwargs..., linu = _vec(δu), du = _vec(δu),