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
Weak Order 2.0 SRK methods #182
Comments
* added constant cache part of DRI1 method * changed sampling of random variables and added RI1 * changed sampling of random variables and added RI1 * n point distributions from weak_utils * order conditions based on tableau for roessler RI * type fixes and scalar test from the paper * added the srk_weak_final file to runtests and # the tests in weak_convergence * Update runtests.jl * add missing test dependency Co-authored-by: Frank Schaefer <frank.schaefer@unibas.ch> Co-authored-by: Christopher Rackauckas <accounts@chrisrackauckas.com>
Note that for implementing the adaptive schemes, all that needs to be done is the addition of the error estimator: https://github.com/SciML/StochasticDiffEq.jl/blob/master/src/perform_step/sri.jl#L111-L118 . If you notice the implementation of EnsembleGPUArray and EnsembleCPUArray will then automatically use the normed pooled error, so that recovers the scheme of Rossler (the Brownian bridge is handled automatically by the machinery of the package. |
Probably, we should also add
and
|
GG @frankschae |
The text was updated successfully, but these errors were encountered: