WIP AbstractTimeseriesSolution plot recipe

modified:   Project.toml
modified:   src/MechGlueDiffEqBase.jl
new file:   test/Manifest.toml
new file:   test/Project.toml
modified:   test/runtests.jl

Author: hustf

Project.toml

@@ -1,11 +1,12 @@
 name = "MechGlueDiffEqBase"
 uuid = "2532746b-52b5-4539-9431-8bb183ab067f"
 authors = ["hustf <[email protected]> and contributors"]
-version = "0.1.2"
+version = "0.1.3"
 
 [deps]
 DiffEqBase = "2b5f629d-d688-5b77-993f-72d75c75574e"
 MechanicalUnits = "e6be9192-89dc-11e9-36e6-5dbcb28f419e"
+SciMLBase = "0bca4576-84f4-4d90-8ffe-ffa030f20462"
 Unitfu = "5ee08b94-2369-4f4a-b8c7-99333ba35fb0"
 
 [compat]

src/MechGlueDiffEqBase.jl

@@ -1,6 +1,10 @@
 module MechGlueDiffEqBase
-import MechanicalUnits.Unitfu.AbstractQuantity
+import Unitfu: AbstractQuantity, Quantity
 import DiffEqBase: value, ODE_DEFAULT_NORM, UNITLESS_ABS2
+import SciMLBase
+import SciMLBase: AbstractTimeseriesSolution, RecipesBase
+import SciMLBase: AbstractDiscreteProblem, AbstractRODESolution, SensitivityInterpolation
+import RecipesBase.@recipe
 export value, ODE_DEFAULT_NORM, UNITLESS_ABS2, Unitfu
 
 function value(x::Type{AbstractQuantity{T,D,U}}) where {T,D,U}
@@ -23,4 +27,27 @@ end
 end
 
 
+@recipe function f(sol::AbstractTimeseriesSolution{T, N, A};
+    plot_analytic=false,
+    denseplot = (sol.dense ||
+                 typeof(sol.prob) <: AbstractDiscreteProblem) &&
+                 !(typeof(sol) <: AbstractRODESolution) &&
+                 !(hasfield(typeof(sol),:interp) &&
+                   typeof(sol.interp) <: SensitivityInterpolation),
+    plotdensity = min(Int(1e5),sol.tslocation==0 ?
+                 (typeof(sol.prob) <: AbstractDiscreteProblem ?
+                 max(1000,100*length(sol)) :
+                 max(1000,10*length(sol))) :
+                 1000*sol.tslocation),
+    tspan = nothing, axis_safety = 0.1,
+    vars=nothing) where {T<:Quantity, N, A<:Array{<:Quantity}}
+
+    print("----------------------------------MechGlueDiffEqBase")
+    return Unifu.ustrip(sol.u), Unifu.ustrip(sol.t)
+end
+        # T = Quantity{Float64,  ᴸ∙ ᴹ∙ ᵀ⁻², Unitfu.FreeUnits{(N,),  ᴸ∙ ᴹ∙ ᵀ⁻², nothing}}
+    # N = 1
+    # A = Vector{Quantity{Float64,  ᴸ∙ ᴹ∙ ᵀ⁻², Unitfu.FreeUnits{(N,),  ᴸ∙ ᴹ∙ ᵀ⁻², nothing}}}
+
+
 end