choose Float64 as default type for trust region adaptation parameters…

… for now
SciML · Sep 26, 2023 · 0190764 · 0190764
1 parent 177def2
commit 0190764
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/trustRegion.jl b/src/trustRegion.jl
@@ -242,7 +242,7 @@ function SciMLBase.__init(prob::NonlinearProblem{uType, iip}, alg::TrustRegion,
     radius_update_scheme = alg.radius_update_scheme
 
     # set default type for all trust region parameters
-    trustType = eltype(u) #typeof(alg.initial_trust_radius)
+    trustType = Float64 #typeof(alg.initial_trust_radius)
     max_trust_radius = convert(trustType, alg.max_trust_radius)
     if iszero(max_trust_radius)
         max_trust_radius = convert(trustType, max(norm(fu1), maximum(u) - minimum(u)))