Fix Aqua error

Author: ErikQQY

lib/BoundaryValueDiffEqCore/src/utils.jl

@@ -623,6 +623,12 @@ end
 @inline __default_cost(f) = f
 @inline __default_cost(fun::BVPFunction) = __default_cost(fun.cost)
 
+"""
+    __construct_internal_problem
+
+Constructs the internal problem based on the type of the boundary value problem and the
+algorithm used. It returns either a `NonlinearProblem` or an `OptimizationProblem`.
+"""
 function __construct_internal_problem(
         prob::BVProblem, alg, loss, jac, jac_prototype, resid_prototype, y, p, M, N)
     T = eltype(y)

lib/BoundaryValueDiffEqMIRK/Project.toml

@@ -42,10 +42,12 @@ FastClosures = "0.3.2"
 ForwardDiff = "0.10.38, 1"
 Hwloc = "3"
 InteractiveUtils = "<0.0.1, 1"
+Ipopt = "1.10.6"
 JET = "0.9.18"
 LinearAlgebra = "1.10"
 LinearSolve = "2.36.2, 3"
 Mooncake = "0.4.108"
+OptimizationMOI = ""
 OrdinaryDiffEqRosenbrock = "1"
 PreallocationTools = "0.4.24"
 PrecompileTools = "1.2"

lib/BoundaryValueDiffEqMIRK/src/adaptivity.jl

@@ -498,8 +498,7 @@ end
     new_prob = remake(prob, u0 = high_sol)
     high_cache = SciMLBase.__init(new_prob, alg, adaptive = false)
 
-    high_nlprob = __construct_nlproblem(
-        high_cache, vec(high_sol), VectorOfArray(high_sol.u))
+    high_nlprob = __construct_problem(high_cache, vec(high_sol), VectorOfArray(high_sol.u))
     high_sol_original = __solve(
         high_nlprob, nlsolve_alg; cache.nlsolve_kwargs..., alias_u0 = true)
     recursive_unflatten!(high_sol, high_sol_original.u)
@@ -518,7 +517,7 @@ end
     new_prob = remake(prob, u0 = high_sol)
     high_cache = SciMLBase.__init(new_prob, __high_order_method(alg), adaptive = false)
 
-    high_nlprob = __construct_nlproblem(high_cache, sol.u, high_sol)
+    high_nlprob = __construct_problem(high_cache, sol.u, high_sol)
     high_sol_nlprob = __solve(
         high_nlprob, nlsolve_alg; cache.nlsolve_kwargs..., alias_u0 = true)
     recursive_unflatten!(high_sol, high_sol_nlprob)