Inconsistency when solving system of SDE with diagonal noise

The higher order methods for diagonal noise are much better than those for non-diagonal noise. If you use non-diagonal noise, it needs to switch to an algorithm that has lower stability, worse error estimators, and worse performance. If you have a diagonal noise SDE you really want to make sure you define it as diagonal to get all of the benefits of the much better solvers: non-diagonal SDEs are mathematically a hundred times harder than diagonal SDEs!

1 Like