Simplify DiffEqArray recipe: always use diffeq_to_arrays path

Remove the fast-path branch in the AbstractDiffEqArray recipe that
bypassed diffeq_to_arrays when no interp was available. Now all cases
go through the same code path, ensuring denseplot=true/false is always
respected:

- denseplot=true: generates plotdensity interpolated points via callable
- denseplot=false: uses raw saved time points from A.t

Also clean up SciMLBase that leaked into [deps] from Pkg.develop.

Co-Authored-By: Chris Rackauckas <accounts@chrisrackauckas.com>
Co-Authored-By: Claude Opus 4.6 (1M context) <noreply@anthropic.com>

Author: ChrisRackauckas

src/vector_of_array.jl

@@ -1368,63 +1368,51 @@ end
         tspan = nothing, plotat = nothing,
         idxs = nothing)
 
-    # If the array has no callable interp and it's a simple case, use the fast path
-    has_interp = hasproperty(VA, :interp) && VA.interp !== nothing
-    if !has_interp && idxs === nothing && plotat === nothing
-        # Simple fast path: just plot t vs u
-        indepsyms = independent_variable_symbols(VA)
-        xguide --> (isempty(indepsyms) ? "" : string(indepsyms[1]))
-        varsyms = variable_symbols(VA)
-        label --> (isempty(varsyms) ? "" : reshape(string.(varsyms), 1, :))
-        VA.t, VA'
+    idxs_input = idxs === nothing ? plottable_indices(VA.u[1]) : idxs
+    if !(idxs_input isa Union{Tuple, AbstractArray})
+        vars = interpret_vars([idxs_input], VA)
     else
-        # Full-featured path: interpret vars, handle phase plots, use callable
-        idxs_input = idxs === nothing ? plottable_indices(VA.u[1]) : idxs
-        if !(idxs_input isa Union{Tuple, AbstractArray})
-            vars = interpret_vars([idxs_input], VA)
-        else
-            vars = interpret_vars(idxs_input, VA)
-        end
+        vars = interpret_vars(idxs_input, VA)
+    end
 
-        tdir = sign(VA.t[end] - VA.t[1])
-        xflip --> tdir < 0
-        seriestype --> :path
+    tdir = sign(VA.t[end] - VA.t[1])
+    xflip --> tdir < 0
+    seriestype --> :path
 
-        tscale = get(plotattributes, :xscale, :identity)
-        plot_vecs, labels = diffeq_to_arrays(
-            VA, denseplot, plotdensity, tspan, vars, tscale, plotat
-        )
+    tscale = get(plotattributes, :xscale, :identity)
+    plot_vecs, labels = diffeq_to_arrays(
+        VA, denseplot, plotdensity, tspan, vars, tscale, plotat
+    )
 
-        # Axis labels for tuple-style idxs
-        if idxs_input isa Tuple && vars[1][1] === DEFAULT_PLOT_FUNC
-            for (guide, idx) in [(:xguide, 2), (:yguide, 3)]
-                if idx <= length(vars[1])
-                    guide --> _var_label(VA, vars[1][idx])
-                end
-            end
-            if length(vars[1]) > 3
-                zguide --> _var_label(VA, vars[1][4])
+    # Axis labels for tuple-style idxs
+    if idxs_input isa Tuple && vars[1][1] === DEFAULT_PLOT_FUNC
+        for (guide, idx) in [(:xguide, 2), (:yguide, 3)]
+            if idx <= length(vars[1])
+                guide --> _var_label(VA, vars[1][idx])
             end
         end
+        if length(vars[1]) > 3
+            zguide --> _var_label(VA, vars[1][4])
+        end
+    end
 
-        # Default xguide for time-vs-variable plots
-        if all(x -> (x[2] isa Integer && x[2] == 0) ||
-                    isequal(x[2], getindepsym_defaultt(VA)), vars)
-            xguide --> "$(getindepsym_defaultt(VA))"
-            if tspan === nothing
-                if tdir > 0
-                    xlims --> (VA.t[1], VA.t[end])
-                else
-                    xlims --> (VA.t[end], VA.t[1])
-                end
+    # Default xguide for time-vs-variable plots
+    if all(x -> (x[2] isa Integer && x[2] == 0) ||
+                isequal(x[2], getindepsym_defaultt(VA)), vars)
+        xguide --> "$(getindepsym_defaultt(VA))"
+        if tspan === nothing
+            if tdir > 0
+                xlims --> (VA.t[1], VA.t[end])
             else
-                xlims --> (tspan[1], tspan[end])
+                xlims --> (VA.t[end], VA.t[1])
             end
+        else
+            xlims --> (tspan[1], tspan[end])
         end
-
-        label --> reshape(labels, 1, length(labels))
-        (plot_vecs...,)
     end
+
+    label --> reshape(labels, 1, length(labels))
+    (plot_vecs...,)
 end
 @recipe function f(VA::DiffEqArray{T, 1}) where {T}
     VA.t, VA.u