DifferentialEquations.jl icon indicating copy to clipboard operation
DifferentialEquations.jl copied to clipboard

matrix SDE problem with non-diagonal noise

Open maol17 opened this issue 2 years ago • 1 comments

Hi,

I'm trying to solve the matrix SDE problem $dA = AD_1 dW_1+D_2AdW_2$, where $A$, $D_1$ and $D_2$ are $n\times n$ matrices. I'm wondering what noise prototype I should use? I tried noise_rate_prototype=zeros(n,n,2) but it doesn't work

maol17 avatar Oct 21 '23 04:10 maol17

For the noise prototype you need to consider it in the vec(u) form as a linear operator, i.e. in the column-order definition of u as a vector. Not the best but it will work.

ChrisRackauckas avatar Oct 27 '23 04:10 ChrisRackauckas