Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Ensemble plot, single component #976

Open
knappa opened this issue Aug 16, 2023 · 3 comments
Open

Ensemble plot, single component #976

knappa opened this issue Aug 16, 2023 · 3 comments

Comments

@knappa
Copy link

knappa commented Aug 16, 2023

It's conceivable that I'm just confused, but I've been trying to plot a single component of an EnsembleSummary for 18 dimensional SDE following the documentation here:
https://docs.sciml.ai/DiffEqDocs/stable/features/ensemble/#Example-4:-Using-the-Analysis-Tools
where it says:

Since there are 8 components to the differential equation, this can get messy, so let's only plot the 3rd component:

using Plots;
plot(summ; idxs = 3);

However, both the example plot in the documentation and the result in my own code display all of the components. (The components of my SDE are at different scales, so a co-mingled plot isn't very useful.)

@ChrisRackauckas
Copy link
Member

I think there is currently a bug in the plot recipe that needs to be handled.

@knappa
Copy link
Author

knappa commented Aug 23, 2023

I'm a julia neophyte, but I'd be happy to look into fixing this bug myself if given a little direction.

@ChrisRackauckas
Copy link
Member

Ref SciML/DiffEqDocs.jl#683.

idxs is the thing in the individual plot recipe https://github.com/SciML/SciMLBase.jl/blob/master/src/solutions/solution_interface.jl#L247

The ensemble plot recipe is here: https://github.com/SciML/SciMLBase.jl/blob/master/src/ensemble/ensemble_solutions.jl#L157.

I believe this should do it: SciML/SciMLBase.jl#493

If you or @ErikQQY can test it out that would be great.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants