Hi,

I am trying to figure out how to use an SDE solver for matrix-valued Stochastic Differential Equations. It looks like the packages GeometricIntegrators.jl, and StochasticIntegrators.jl might have the functionality but I don’t understand how to use them. I have been using DifferentialEquations.jl for univariate SDEs with good success but this fails:

```
using DifferentialEquations, GeometricIntegratorsDiffEq
prob = SDEProblem(drift, diff, u0, tspan, p=p)
sol = solve(prob,GIEuler())
```

throws this error:

ERROR: Incompatible problem+solver pairing.

For example, this can occur if an ODE solver is passed with an SDEProblem.

Solvers are only capable of handling specific problem types. Please double

check that the chosen pairing is capable for handling the given problems.

Problem type: SDEProblem

Solver type: GIEuler

Problem types compatible with the chosen solver: ODEProblem

So is there a way to use Geometric Integrators on SDEs in Julia? If so, could someone share a toy example with the correct syntax?

Thanks,

Simone.