# \#differentialequation

**URL:** https://discourse.julialang.org/tag/differentialequation/751.md

[Latest](https://discourse.julialang.org/latest.md) · [Categories](https://discourse.julialang.org/categories.md) · [Tags](https://discourse.julialang.org/tags.md)

---

## ["Fixed iteration" DiffEq solvers on the GPU?](https://discourse.julialang.org/t/fixed-iteration-diffeq-solvers-on-the-gpu/139357)

<div class="topic-metadata">

**Author:** [@jl\_enthusiast](https://discourse.julialang.org/u/jl_enthusiast)\
**Replies:** 2\
**Last updated:** [September 11, 2026, 2:37pm UTC](https://discourse.julialang.org/t/fixed-iteration-diffeq-solvers-on-the-gpu/139357 "2026-09-11T14:37:32Z")

</div>

Hey, I am currently trying to (re)write some old hobby project into Julia (related to raytracing) and want to rely on the DiffEq solver interface, ideally inside a KernelAbstractions.jl kernel (as I am performing a lot …

---

## [Does MTK nominal metadata for variables set abstol when solving ODE](https://discourse.julialang.org/t/does-mtk-nominal-metadata-for-variables-set-abstol-when-solving-ode/139123)

<div class="topic-metadata">

**Author:** [@mark.garnett](https://discourse.julialang.org/u/mark.garnett)\
**Replies:** 4\
**Last updated:** [September 4, 2026, 1:33pm UTC](https://discourse.julialang.org/t/does-mtk-nominal-metadata-for-variables-set-abstol-when-solving-ode/139123 "2026-09-04T13:33:06Z")

</div>

when you define a variables @variables begin x(t), \[nominal=1e6\] y(t), \[nominal=1e-3\] end will it essentially do the equivalent of setting abstol to a vector of absolute tolerances rather than a scalar, by multiplying …

---

## [Compiler mis-recognizes contiguous reinterpreted matrices?](https://discourse.julialang.org/t/compiler-mis-recognizes-contiguous-reinterpreted-matrices/138979)

<div class="topic-metadata">

**Author:** [@bremez](https://discourse.julialang.org/u/bremez)\
**Replies:** 11\
**Last updated:** [August 25, 2026, 2:02pm UTC](https://discourse.julialang.org/t/compiler-mis-recognizes-contiguous-reinterpreted-matrices/138979 "2026-08-25T14:02:33Z")

</div>

I am implementing an ODE to use with OrdinaryDiffEq. I am simultaneously evolving in time multiple quantities, which are combinations of vectors and matrices, some real and some complex. The approach I took has been to g…

---

## [Error while installing Juila packages on windows](https://discourse.julialang.org/t/error-while-installing-juila-packages-on-windows/138758)

<div class="topic-metadata">

**Author:** [@Aditya747S](https://discourse.julialang.org/u/Aditya747S)\
**Replies:** 14\
**Last updated:** [August 23, 2026, 1:31pm UTC](https://discourse.julialang.org/t/error-while-installing-juila-packages-on-windows/138758 "2026-08-23T13:31:34Z")

</div>

I am installing julia packages on my windows machine and i am getting errros that go like this → \*.dll is either not designed to run on windows or contains an error. Try installing again using original installation med…

---

## [MTK SDDE ArgumentError when Connecting two Components](https://discourse.julialang.org/t/mtk-sdde-argumenterror-when-connecting-two-components/138861)

<div class="topic-metadata">

**Author:** [@bebrunn](https://discourse.julialang.org/u/bebrunn)\
**Replies:** 7\
**Last updated:** [August 18, 2026, 11:38am UTC](https://discourse.julialang.org/t/mtk-sdde-argumenterror-when-connecting-two-components/138861 "2026-08-18T11:38:45Z")

</div>

Hi, I am trying to connect two MTK components and then solve the resulting SDDE system. However, I am getting an ArgumentsError that I don’t understand. using ModelingToolkit using DifferentialEquations using ModelingT…

---

## ["no method matching" error when running DiffEqGPU example from docs](https://discourse.julialang.org/t/no-method-matching-error-when-running-diffeqgpu-example-from-docs/138471)

<div class="topic-metadata">

**Author:** [@CourtA96](https://discourse.julialang.org/u/CourtA96)\
**Replies:** 3\
**Last updated:** [August 5, 2026, 2:52pm UTC](https://discourse.julialang.org/t/no-method-matching-error-when-running-diffeqgpu-example-from-docs/138471 "2026-08-05T14:52:51Z")

</div>

When I run the following example code from the DiffEqGPU documentation, using DiffEqGPU, OrdinaryDiffEq, StaticArrays, CUDA function lorenz(u, p, t) σ = p\[1\] ρ = p\[2\] β = p\[3\] du1 = σ \* (u\[2\] - u\[1\]) …

---

## [Interpolating solution to constrained BVP](https://discourse.julialang.org/t/interpolating-solution-to-constrained-bvp/137989)

<div class="topic-metadata">

**Author:** [@kar15](https://discourse.julialang.org/u/kar15)\
**Replies:** 3\
**Last updated:** [July 9, 2026, 7:41am UTC](https://discourse.julialang.org/t/interpolating-solution-to-constrained-bvp/137989 "2026-07-09T07:41:02Z")

</div>

Hi all! When trying to find interpolated values of a solution to a constrained BVP using SciMLs collocation method, I get an error. Here is a MWE: using BoundaryValueDiffEq, OptimizationIpopt function f!(du, u, p, t) …

---

## [Issues solving DifferentialEquations EnsembleProblem](https://discourse.julialang.org/t/issues-solving-differentialequations-ensembleproblem/137413)

<div class="topic-metadata">

**Author:** [@Theodore\_Rieger](https://discourse.julialang.org/u/Theodore_Rieger)\
**Replies:** 2\
**Last updated:** [June 3, 2026, 9:50am UTC](https://discourse.julialang.org/t/issues-solving-differentialequations-ensembleproblem/137413 "2026-06-03T09:50:05Z")

</div>

Hello, I was trying to experiment with ensemble problems from differentialequations package. I’m very quickly running into a roadblock in that even the examples from documentation are producing errors. Below is an exampl…

---

## [Tsit5() and DP5() stability](https://discourse.julialang.org/t/tsit5-and-dp5-stability/137137)

<div class="topic-metadata">

**Author:** [@ligeston](https://discourse.julialang.org/u/ligeston)\
**Replies:** 13\
**Last updated:** [May 19, 2026, 2:32pm UTC](https://discourse.julialang.org/t/tsit5-and-dp5-stability/137137 "2026-05-19T14:32:02Z")

</div>

Hi guys, I am wrapping up a script to solve the TOV equation, similar to the TOVsolver\_Julia, and I used ‘CommonSolve’ and ‘DifferentialEquations’ packages. In the solver, unlike the case in TOVsolver\_Julia, I used cal…

---

## [\[BoundaryValueDiffEq.jl\] Reducing allocations](https://discourse.julialang.org/t/boundaryvaluediffeq-jl-reducing-allocations/136255)

<div class="topic-metadata">

**Author:** [@MaxSchelski](https://discourse.julialang.org/u/MaxSchelski)\
**Replies:** 9\
**Last updated:** [March 23, 2026, 2:48pm UTC](https://discourse.julialang.org/t/boundaryvaluediffeq-jl-reducing-allocations/136255 "2026-03-23T14:48:23Z")

</div>

Hi, I started to implement my Python pipeline for solving boundary value problems in Julia using BoundaryValueDiffEq.jl. I tried to optimize the functions but I am wondering if it is possible to increase performance fur…

---

## [Mutation attempt in \`ContinuousCallback\` when solving \`SDEProblem{StaticArraysCore.SVector{2, Float64}}\`](https://discourse.julialang.org/t/mutation-attempt-in-continuouscallback-when-solving-sdeproblem-staticarrayscore-svector-2-float64/136168)

<div class="topic-metadata">

**Author:** [@NoFishLikeIan](https://discourse.julialang.org/u/NoFishLikeIan)\
**Replies:** 1\
**Last updated:** [March 13, 2026, 11:37am UTC](https://discourse.julialang.org/t/mutation-attempt-in-continuouscallback-when-solving-sdeproblem-staticarrayscore-svector-2-float64/136168 "2026-03-13T11:37:29Z")

</div>

I have an SDE where the drift function f is defined out-of-place and the state is a StaticArrays.SVector. I need to terminate! the integration using a ContinuousCallback. This raises an error, as the StochasticDiffEq pac…

---

## [My simulations are suddenly 5-10 times faster. Why?](https://discourse.julialang.org/t/my-simulations-are-suddenly-5-10-times-faster-why/136156)

<div class="topic-metadata">

**Author:** [@ufechner7](https://discourse.julialang.org/u/ufechner7)\
**Replies:** 3\
**Last updated:** [March 11, 2026, 7:36pm UTC](https://discourse.julialang.org/t/my-simulations-are-suddenly-5-10-times-faster-why/136156 "2026-03-11T19:36:40Z")

</div>

I am running a simulation with about 70 states that uses the DFBDF implicit solver. See: KiteModels.jl/examples/bench\_4p.jl at main · OpenSourceAWE/KiteModels.jl · GitHub Recently, I noticed a 5-10x increase in performa…

---

## [Several Compilation errors and continuous recompilation issue](https://discourse.julialang.org/t/several-compilation-errors-and-continuous-recompilation-issue/135894)

<div class="topic-metadata">

**Author:** [@simon79](https://discourse.julialang.org/u/simon79)\
**Replies:** 6\
**Last updated:** [March 6, 2026, 2:18pm UTC](https://discourse.julialang.org/t/several-compilation-errors-and-continuous-recompilation-issue/135894 "2026-03-06T14:18:59Z")

</div>

Hi all, I have been using (Ordinary)DifferentialEquations.jl for a few years without issues. As of today, after an update, a very long list of errors is generated which there seems to be no way to resolve. I paste only …

---

## [A non-classical ODE problem](https://discourse.julialang.org/t/a-non-classical-ode-problem/135868)

<div class="topic-metadata">

**Author:** [@Veenty](https://discourse.julialang.org/u/Veenty)\
**Replies:** 2\
**Last updated:** [February 26, 2026, 5:03pm UTC](https://discourse.julialang.org/t/a-non-classical-ode-problem/135868 "2026-02-26T17:03:03Z")

</div>

I have an evolution problem that I cannot cast in a classical form. \\frac{d}{dt} A u = u + F(u,v) \\\\ G(u) = 0 Where A is a bad matrix to invert, however (A - \\Delta t) is ok to invert. When I’m coding this with an imp…

---

## [DAE‑safe bounding for CoolProp property calls in Julia](https://discourse.julialang.org/t/dae-safe-bounding-for-coolprop-property-calls-in-julia/135416)

<div class="topic-metadata">

**Author:** [@TimGyx](https://discourse.julialang.org/u/TimGyx)\
**Replies:** 7\
**Last updated:** [February 18, 2026, 3:38am UTC](https://discourse.julialang.org/t/dae-safe-bounding-for-coolprop-property-calls-in-julia/135416 "2026-02-18T03:38:45Z")

</div>

I am building a large DAE-based vapor‑compression refrigeration loop in Julia using DifferentialEquations.jl and CoolProp for thermodynamic properties, and I am running into repeated CoolProp crashes during DAE initializ…

---

## [ODE Parameter Estimation with Priors](https://discourse.julialang.org/t/ode-parameter-estimation-with-priors/135691)

<div class="topic-metadata">

**Author:** [@donkeysaddle](https://discourse.julialang.org/u/donkeysaddle)\
**Replies:** 2\
**Last updated:** [February 17, 2026, 5:54pm UTC](https://discourse.julialang.org/t/ode-parameter-estimation-with-priors/135691 "2026-02-17T17:54:42Z")

</div>

I have never done parameter estimation with priors(bayesian inference) before and I am a little unsure if I am going about it the right way. I created a simple example, but will ultimately be doing it for a more complica…

---

## [Unitful.jl and DifferentialEquations.jl: compatibility?](https://discourse.julialang.org/t/unitful-jl-and-differentialequations-jl-compatibility/68594)

<div class="topic-metadata">

**Author:** [@Ickaser](https://discourse.julialang.org/u/Ickaser)\
**Replies:** 16\
**Last updated:** [February 13, 2026, 11:56pm UTC](https://discourse.julialang.org/t/unitful-jl-and-differentialequations-jl-compatibility/68594 "2026-02-13T23:56:04Z")

</div>

I’ve been using Julia for a little while, long enough to have successfully solved ODEs before. I’m trying to use DifferentialEquations together with Unitful, which I thought I had read should work, but it doesn’t seem to…

---

## [Parameter sensitivity of ODE with parameter-dependent event](https://discourse.julialang.org/t/parameter-sensitivity-of-ode-with-parameter-dependent-event/103602)

<div class="topic-metadata">

**Author:** [@Neodym](https://discourse.julialang.org/u/Neodym)\
**Replies:** 4\
**Last updated:** [February 2, 2026, 11:28am UTC](https://discourse.julialang.org/t/parameter-sensitivity-of-ode-with-parameter-dependent-event/103602 "2026-02-02T11:28:17Z")

</div>

Hey, I’m trying to get parameter sensitivities of an ODE problem with a parameter-dependent time-event. Just throwing AD at the problem doesn’t seem to work: using OrdinaryDiffEq using SciMLSensitivity using ForwardDi…

---

## [Singular BVP with integrable boundaries](https://discourse.julialang.org/t/singular-bvp-with-integrable-boundaries/135224)

<div class="topic-metadata">

**Author:** [@NoFishLikeIan](https://discourse.julialang.org/u/NoFishLikeIan)\
**Replies:** 2\
**Last updated:** [January 26, 2026, 7:13am UTC](https://discourse.julialang.org/t/singular-bvp-with-integrable-boundaries/135224 "2026-01-26T07:13:46Z")

</div>

I’m solving a singular first-order, two dimensional ODE arising in reputation games (Faingold & Sannikov, 2011). On \\phi \\in \[0, 1\], I have u' = \\frac{v}{\\phi(1-\\phi)} \\\\ v' = \\frac{1}{\\phi(1-\\phi)} \\left(v + 2\\rho \\f…

---

## [Convergence of SDE solvers](https://discourse.julialang.org/t/convergence-of-sde-solvers/135179)

<div class="topic-metadata">

**Author:** [@diffeqslvr](https://discourse.julialang.org/u/diffeqslvr)\
**Replies:** 6\
**Last updated:** [January 22, 2026, 11:22am UTC](https://discourse.julialang.org/t/convergence-of-sde-solvers/135179 "2026-01-22T11:22:52Z")

</div>

I am wondering how to make sure that an SDE solver has converged on a trajectory level. In other words, I want to make sure that I have chosen abstol and reltol small enough. What I would do in the deterministic case is…

---

## [Difference between tstops and d\_discontinuities](https://discourse.julialang.org/t/difference-between-tstops-and-d-discontinuities/126847)

<div class="topic-metadata">

**Author:** [@cstjean](https://discourse.julialang.org/u/cstjean)\
**Replies:** 3\
**Last updated:** [January 13, 2026, 2:30pm UTC](https://discourse.julialang.org/t/difference-between-tstops-and-d-discontinuities/126847 "2026-01-13T14:30:38Z")

</div>

The docstring says: • tstops: Denotes extra times that the timestepping algorithm must step to. This should be used to help the solver deal with discontinuities and singularities, since stepping exactly at t…

---

## [GPU support for Turing modeling with system of ODEs](https://discourse.julialang.org/t/gpu-support-for-turing-modeling-with-system-of-odes/77410)

<div class="topic-metadata">

**Author:** [@vembha](https://discourse.julialang.org/u/vembha)\
**Replies:** 8\
**Last updated:** [January 5, 2026, 7:38am UTC](https://discourse.julialang.org/t/gpu-support-for-turing-modeling-with-system-of-odes/77410 "2026-01-05T07:38:43Z")

</div>

Hello, I am trying to perform Bayesian analysis on a Turing model with an ecological ODE system embedded; consider it to be a complex variant of a Lotka-Volterra model. As of now, the code is running perfectly fine with…

---

## [Stochastic filtration problem](https://discourse.julialang.org/t/stochastic-filtration-problem/134713)

<div class="topic-metadata">

**Author:** [@ymuv1](https://discourse.julialang.org/u/ymuv1)\
**Replies:** 1\
**Last updated:** [December 23, 2025, 7:13pm UTC](https://discourse.julialang.org/t/stochastic-filtration-problem/134713 "2025-12-23T19:13:46Z")

</div>

Hi all, I’m working with a physical system whose state u(t) follows an SDE du=f(u,t)dt+g(u,t)dW\_t The system produces a measured signal Y\_t given by dY\_t = h(u,t)dt+dW\_t where dW\_t is the same Wiener increment in bo…

---

## [Are there any examples of DiscreteProblem usage?](https://discourse.julialang.org/t/are-there-any-examples-of-discreteproblem-usage/134618)

<div class="topic-metadata">

**Author:** [@teacup775](https://discourse.julialang.org/u/teacup775)\
**Replies:** 3\
**Last updated:** [December 22, 2025, 4:58pm UTC](https://discourse.julialang.org/t/are-there-any-examples-of-discreteproblem-usage/134618 "2025-12-22T16:58:45Z")

</div>

I have searched extensively online and cannot find anywhere a simple example of a set up for a discrete problem. A simply want an example of usage and set up.

---

## [Current state of reaction-diffusion solvers with ModelingToolkit](https://discourse.julialang.org/t/current-state-of-reaction-diffusion-solvers-with-modelingtoolkit/118492)

<div class="topic-metadata">

**Author:** [@Alicanter](https://discourse.julialang.org/u/Alicanter)\
**Replies:** 10\
**Last updated:** [December 18, 2025, 9:49pm UTC](https://discourse.julialang.org/t/current-state-of-reaction-diffusion-solvers-with-modelingtoolkit/118492 "2025-12-18T21:49:48Z")

</div>

I am trying to solve a reaction-diffusion in Julia. I’ve worked with ModelingToolkit and the SciML ecosystem before, but it seems that with PDEs there is not a single best way to solve them. The simplest way I found to…

---

## [Very slow training of Neural ODE with exogeneous input compared to jax using dense interpolated solution](https://discourse.julialang.org/t/very-slow-training-of-neural-ode-with-exogeneous-input-compared-to-jax-using-dense-interpolated-solution/134557)

<div class="topic-metadata">

**Author:** [@johtok](https://discourse.julialang.org/u/johtok)\
**Replies:** 4\
**Last updated:** [December 15, 2025, 11:09pm UTC](https://discourse.julialang.org/t/very-slow-training-of-neural-ode-with-exogeneous-input-compared-to-jax-using-dense-interpolated-solution/134557 "2025-12-15T23:09:39Z")

</div>

Hi guys! Ive been breaking my neck trying to get a julia alternative to this jax neural ode with exogeneous input script working but without luck! The julia version is painfully slow whereas the jax version runs in a c…

---

## [Poincare section for Henon Helies system](https://discourse.julialang.org/t/poincare-section-for-henon-helies-system/124059)

<div class="topic-metadata">

**Author:** [@nn-heckfy](https://discourse.julialang.org/u/nn-heckfy)\
**Replies:** 9\
**Last updated:** [November 10, 2025, 3:12pm UTC](https://discourse.julialang.org/t/poincare-section-for-henon-helies-system/124059 "2025-11-10T15:12:48Z")

</div>

Greetings. I’m trying to solve a problem of plotting a Poincare map for the Henon Heiles’s system. But the figure is completely painted over, but should have cavities. Maybe I’m using the Julia commands incorrectly someh…

---

## [Issue with SS.solve using DynamicSS method in SteadyStateDiffEq](https://discourse.julialang.org/t/issue-with-ss-solve-using-dynamicss-method-in-steadystatediffeq/133712)

<div class="topic-metadata">

**Author:** [@Bill\_Holmes](https://discourse.julialang.org/u/Bill_Holmes)\
**Replies:** 3\
**Last updated:** [November 7, 2025, 5:13pm UTC](https://discourse.julialang.org/t/issue-with-ss-solve-using-dynamicss-method-in-steadystatediffeq/133712 "2025-11-07T17:13:21Z")

</div>

I have a quick question about using SteadyStateDiffEq.jl. I am getting error behavior when trying to use the DynamicSS method (following Steady State Solvers · DifferentialEquations.jl). Below is a MWE. I 1) define a sta…

---

## [Large allocations with MTK event](https://discourse.julialang.org/t/large-allocations-with-mtk-event/133680)

<div class="topic-metadata">

**Author:** [@klinders](https://discourse.julialang.org/u/klinders)\
**Replies:** 1\
**Last updated:** [November 5, 2025, 11:35pm UTC](https://discourse.julialang.org/t/large-allocations-with-mtk-event/133680 "2025-11-05T23:35:22Z")

</div>

Dear all, Recently, you helped me get events working in MTK. Thank you for all the help! Continuing on the topic, I ran into a performance issue with events. I have the following code to measure the time and allocation…

---

## [Solving SemilinearODE without using automatic differentiation](https://discourse.julialang.org/t/solving-semilinearode-without-using-automatic-differentiation/133434)

<div class="topic-metadata">

**Author:** [@sanjeev](https://discourse.julialang.org/u/sanjeev)\
**Replies:** 4\
**Last updated:** [October 27, 2025, 2:10am UTC](https://discourse.julialang.org/t/solving-semilinearode-without-using-automatic-differentiation/133434 "2025-10-27T02:10:13Z")

</div>

I have semilinear ODE with the time-dependent linear implicit part as function f1! So, instead of solving (I - γΔt \* J₁) \* u\_{n+1} = rhs, I want to use Krylov method to solve (I - γΔt \* J1) \* u\_{n+1} = rhs, where J1 \* u…

[Next page](https://discourse.julialang.org/tag/differentialequation/751.md?match_all_tags=true&page=1&tags%5B%5D=differentialequation)