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

Investigate the possibility of StaticArrays tableau methods #114

Closed
ChrisRackauckas opened this issue Aug 6, 2017 · 1 comment
Closed

Comments

@ChrisRackauckas
Copy link
Member

When I measured before, tableau-based methods were really really bad compared to hand-unrolled methods. But this was with Arrays and FixedSizeArrays. But now there's StaticArrays.jl, could it work?

The issue is that there's tons of specialization for a lot of the methods. Each of them "skip true zeros" in the tableaus (which, by the assumptions used to do the derivation, can be a lot in higher order methods). They have different error estimation schemes (DP8 and BS5 are good examples there). They have wildly different dense output schemes. Etc. So I don't know what we'll really get out of it because we'll still have to specialize a lot, but maybe it could be useful.

@ChrisRackauckas
Copy link
Member Author

SDIRKs and Rosenbrocks can probably use it too. The issue is that doing this will make it so that .= doesn't put as much in a single call, which will hurt immensely hurt GPU kernels.

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

1 participant