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

matrix SDE problem with non-diagonal noise #991

Open
maol17 opened this issue Oct 21, 2023 · 1 comment
Open

matrix SDE problem with non-diagonal noise #991

maol17 opened this issue Oct 21, 2023 · 1 comment

Comments

@maol17
Copy link

maol17 commented Oct 21, 2023

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

@ChrisRackauckas
Copy link
Member

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.

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

2 participants