Unexpected allocations in loops

[lmiq]

Array{F}

Is an abstract type. Probably you want just Vector{F}.

Edit: actually since you want to define the dimension at runtime, you may need to define:

@t-bltg

struct Species{A}
  q::F
  ρ::A
end

Such that you get a concrete species.