There’s an ODE to be solved. It’s a chemical reaction system. There is 863 reactions and 230 species. The system is quite stiff and large. So I ues some algorithms advised by Solving Large Stiff Equations · DifferentialEquations.jl . I tried every algorithms, but the values never change at each step…

This almost ceratinly means you coded your function incorrectly. if I had to guess, it would be that you used the in place form of f(du, u, p, t) but didn’t update du, but it’s hard to say for sure without seeing the code.

Here is the code. It seems du will update. Whatever algorithm I try, the value will not change. function crnn(du,c,k,t) c = reshape(c,(230,1)) k = reshape(k,(863,1)) log_r = w_in * @.log(clamp(c,1e-35,1e15)) log_r = log_r + @.log(clamp(k,1e-35,1e15)) r = @.exp(log_r); du = …

Yeah. My guess was correct. You need du .= w_out' * r; Julia uses “pass by sharing” so when you say du = ... you aren’t changing du you’re just making a new du. Also the return du` isn’t doing anything here since you’re using the in-place formulation.

Could you please help me to solve my problem? Please see the post: Forecasting time series data with neural ordinary differential equations in Julia

Could you please help? [image] Forecasting time series data with neural ordinary differential equations in Julia General Usage Dear all, I am a new user of Julia and following a tutorial on NeuralODE (please see the tutorial: Experiments with Neural ODEs in Julia ). How …

It’s also unclear what is intended with w_out. Is it intentionally global? Might help to do a MWE with a small, stiff, toy system. That can help separate syntax issues from the question of how to solve stiff systems.

Problem with stiff ODE solved by Julia

General Usage

dmjalal90 January 16, 2023, 2:12am 6

Could you please help?

Topic		Replies	Views
Using stiff solvers for a large system of ode Performance question , diffeq , pde	3	480	April 18, 2023
Strange Behaviour of Stiff Solvers in DifferentialEquations.jl Numerics diffeq	2	679	November 23, 2019
Error while using DifferentialEquations.jl stiff methods Numerics diffeq	3	930	March 19, 2019
Julia DifferentialEquations for state-to-state binary mixture kinetics General Usage	1	259	November 26, 2020
ODE system solvable with forward Euler but not general solvers (early stiffness) New to Julia diffeq	9	1498	April 2, 2022

Problem with stiff ODE solved by Julia

Related topics