fix: partial revert of previous fix

Author: avik-pal

src/internal/helpers.jl

@@ -1,9 +1,8 @@
 # Evaluate the residual function at a given point
 function evaluate_f(prob::AbstractNonlinearProblem{uType, iip}, u) where {uType, iip}
-    (; f, u0, p) = prob
+    (; f, p) = prob
     if iip
-        fu = f.resid_prototype === nothing ? zero(u) :
-             promote_type(eltype(u), eltype(f.resid_prototype)).(f.resid_prototype)
+        fu = f.resid_prototype === nothing ? zero(u) : similar(f.resid_prototype)
         f(fu, u, p)
     else
         fu = f(u, p)

src/internal/jacobian.jl

@@ -92,11 +92,10 @@ function JacobianCache(prob, alg, f::F, fu_, u, p; stats, autodiff = nothing,
                 end
             end
         else
-            # NOTE: don't use similar because that might lead to an unwrapped array
             if eltype(f.jac_prototype) <: Bool
-                promote_type(eltype(fu), eltype(u)).(f.jac_prototype)
+                similar(f.jac_prototype, promote_type(eltype(fu), eltype(u)))
             else
-                copy(f.jac_prototype)
+                similar(f.jac_prototype)
             end
         end
     end