Plotting ODEProblem solutions with Makie

airpmb · March 27, 2021, 9:29pm

Whatever (plot recipes?) allows for the solution output from solving an ODEProblem in ModelingToolkit.jl (or maybe just OrdinaryDiffEq.jl?) to be just plotted using Plots.jl is wonderful and magical. I love how it’s possible to just do sol = solve(problem); plot(sol) as well as plotting individual states by name (plot(sol, vars=[(x, y), (w, z)]) etc).

However overall I’m finding Makie to be my favorite tool for plotting/visualization. Is there anything out there yet that provides similar power and expressivity but using Makie instead of Plots?

ChrisRackauckas · March 28, 2021, 5:01pm

I want to move everything to Makie. Please help us. The recipes are in here:

github.com

SciML/SciMLBase.jl/blob/master/src/solutions/solution_interface.jl#L151


      
                  end
              elseif i isa Base.Integer || i isa AbstractRange || i isa AbstractVector{<:Base.Integer}
                  A[i, args...]
              else
                  error("Invalid indexing of solution")
              end
          end
          
          function observed(A::AbstractTimeseriesSolution, sym, i::Int)
              getobserved(A)(sym, A[i], A.prob.p, A.t[i])
          end
          
          function observed(A::AbstractTimeseriesSolution, sym, i::AbstractArray{Int})
              getobserved(A).((sym,), A.u[i], (A.prob.p,), A.t[i])
          end
          
          function observed(A::AbstractTimeseriesSolution, sym, i::Colon)
              getobserved(A).((sym,), A.u, (A.prob.p,), A.t)
          end
          
          Base.@propagate_inbounds function Base.getindex(A::AbstractNoTimeSolution, sym)

I think @asinghvi17 had a working version calling the recipe from Makie, but I forget how he did it. The main issue is that type recipes in Makie don’t allow custom keyword arguments.

That might be the last blocker.

rveltz · March 31, 2021, 8:19am

I am interested in this as well. Any progress has been made?

EvoArt · April 2, 2021, 11:46pm

This is pretty big!

I’ve been wondering which plotting package(s) to invest time in learning. If SciML is backing Makie, then it seems like a no-brainer.

kool7d · April 4, 2021, 8:18am

This looks interesting. I have been working on a few things that would benefit immensely from incorporating SciML. If anyone wants to chat about moving this stuff to Makie I think that could be beneficial.

rveltz · February 21, 2024, 7:48pm

is it closed as an “issue”?

asinghvi17 · February 21, 2024, 8:39pm

We have basic support but it’s a bit barebones at the moment. Work is ongoing on the Makie end to make this more possible, though!