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

time dependent a/drift and b/noise terms #11

Open
sdementen opened this issue Oct 13, 2023 · 3 comments
Open

time dependent a/drift and b/noise terms #11

sdementen opened this issue Oct 13, 2023 · 3 comments

Comments

@sdementen
Copy link

If we want to have $dX = sigma(t) * dWt$, should we just add an extra variable $dt = 1 *dt$ and then use this variable as time ?
For instance,

# x[0] = time, x[1] = X, dX = (cos(t)+1.1) * dWt
pychastic.sde_problem.SDEProblem(
   lambda x: jnp.array([1., 0.]),
   lambda x: jnp.array([[0, 0], [0, (jnp.cos(x[0]) + 1.1)]]),
   x0=[0, 1],
   tmax=10,
)

It works as expected but would it be better to have for drift and noise, a signature like lambda t,x: ... ?

@RadostW
Copy link
Owner

RadostW commented Oct 13, 2023

Implementing higher order schemes for nonautonomous equations are much more involved. This is a deliberate design choice to ease development.

@sdementen
Copy link
Author

and does the trick of adding the time dimension work ? or is it highly unstable ?
and is the complexity for nonautonomous equations independent of the noise/drift expressions (for instance if they are linear in the other variables or even constant in the other variables, i.e. only time dependent) ?

@sdementen
Copy link
Author

without understanding everything ;-)
https://chat.openai.com/share/2394db4e-d330-4c70-88cb-5038bed14a99

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