Here is a prototype of my problem: using OrdinaryDiffEq n = 2^18; u0 = rand(ComplexF64, n); function f(du, u, p, t) for i in 1:length(u) du[i] = -1.0im * u[i] end return end mutable struct A u0::Vector{ComplexF64} integrator function A(u0) a = new() …

Confusing memory allocations when using the integrator of DifferentialEquations.jl

Oscar_Smith January 11, 2022, 9:42pm 2

Your struct isn’t type stable. Specifically, the compiler can’t figure out what the type of a.integrator is.

Topic		Replies	Views
A lot of memory allocations which I do not understand Performance	2	561	March 23, 2020
Unexpected Allocations in Julia 0.5 General Usage	4	573	December 5, 2016
Speed up solution system of ODEs New to Julia diffeq	23	4430	June 2, 2019
Help reduce large gc time Performance diffeq	32	3790	November 18, 2018
Controlling memory allocation using external struct as placeholder General Usage question	5	895	March 11, 2017