Fix MIRK QA public API checks

Co-Authored-By: Chris Rackauckas <accounts@chrisrackauckas.com>

Author: ChrisRackauckas

lib/BoundaryValueDiffEqMIRK/src/BoundaryValueDiffEqMIRK.jl

@@ -12,8 +12,8 @@ using BoundaryValueDiffEqCore: BoundaryValueDiffEqCore,
     __concrete_solve_algorithm, diff!, EvalSol,
     concrete_jacobian_algorithm, eval_bc_residual,
     eval_bc_residual!, __maybe_matmul!, __resize!,
-    __extract_problem_details, interval,
-    __needs_diffcache, __maybe_allocate_diffcache,
+    __extract_problem_details, interval, nodual_value,
+    __maybe_allocate_diffcache,
     __restructure_sol, __cache_trait, __get_bcresid_prototype,
     safe_similar, __vec, __vec_f, __vec_f!, __vec_bc, __vec_bc!,
     recursive_flatten_twopoint!,

lib/BoundaryValueDiffEqMIRK/src/interpolation.jl

@@ -218,8 +218,10 @@ function __stage_weighted_zero(prototype, k, weights)
     return fill!(similar(prototype, T), zero(T))
 end
 
-@inline __primal_value(x) = x
-@inline __primal_value(x::ForwardDiff.Dual) = __primal_value(ForwardDiff.value(x))
+@inline function __primal_value(x)
+    y = nodual_value(x)
+    return y === x ? y : __primal_value(y)
+end
 
 """
     EvalSol