Hi everybody, I am trying to solve the following system of ODE with DifferentialEquations.jl function newton!(t, y, dy, μ) r = norm(y[1:3]) dy[1:3] = y[4:6] dy[4:6] = -μ * y[1:3] / r^3 end Without units this works perfectly: function test() r0 = [1131.340, -2282.343, 6672.423] …

DifferentialEquations.jl and systems with heterogenous units

ChrisRackauckas July 15, 2020, 1:16pm 15

Probably the broadcast overload + constant prop is helping

Topic		Replies	Views
Unitful + DiffEq + different physical quantities + derivatives? General Usage question , package , diffeq , performance	1	396	June 9, 2021
Unitful and automatic differentiation? General Usage forwarddiff , unitful	18	1282	November 9, 2021
Unitful.jl and DifferentialEquations.jl: compatibility? New to Julia unitful , differentialequation	15	2356	January 28, 2022
[ANN] UnitfulTensors.jl: Efficient arrays with units Package Announcements linearalgebra , arrays , unitful	17	1341	March 26, 2025
Still looking for a way to use units meaningfully with SciML ecosystem General Usage forwarddiff , unitful	17	1386	January 29, 2025