EnsembleSummary does not work on solution of complex SDE

[HumpyBlumpy]

Running
summ = EnsembleSummary(sol,0.0:0.1:10.0)
where sol is a solution of EnsembleProblem of SDEProblem returns
MethodError: no method matching isless(::Complex{Float64}, ::Complex{Float64})

[ChrisRackauckas-Claude]

Still reproduces verbatim:

MethodError: no method matching isless(::ComplexF64, ::ComplexF64)
  [1] lt(o::Base.Order.ForwardOrdering, a::ComplexF64, b::ComplexF64)

An identical SDE with real-valued output summarizes fine (u[end]=15388.95, qlow[end]=1.394), which isolates it to the complex values rather than the ensemble machinery.

Worth noting where it fails, though: the quantile/ordering path. Quantiles have no natural definition for complex numbers, so this probably needs a designed semantic (component-wise? magnitude? refuse with a clear error?) rather than a straightforward fix.

Verified on DifferentialEquations v8.0.2 / OrdinaryDiffEq v7.1.3 / Julia 1.11.9.