I have a type-stability question. In my 2D axisymmetric ray tracing code, I define an array of surfaces that has the type Union{Cone, Plane}, like surfaces = Matrix{Union{Cone, Plane}}(undef, 100, 2) All of them should in principle be cones, but due to small errors from a discretization earlier on…

Cone and Plane are abstract types. Only unions of isbits types are fast. If you are using 64 bits floats then you probably want the type to be Union{Cone{Float64},Plane{Float64}}.

Array of union causing 10x slowdown

jling December 7, 2022, 2:56am 2

Array of union causing 10x slowdown