I am struggling to understand the type instability of function foo(x,y,z,r) r_prime = Point3D(x,y,z) result = r - rprime end given as input x = 0.; y = 0.; z = 0.; r = Point3D(1.,1.,1.) @code_warntype foo(x,y,z,r) The output I obtain is MethodInstance for foo(::Float64, ::Float64, ::…

It is very likely that Point3D is somehow an alias for SVector, and the source of type instability is rprime (in result = r - rprime), which is some leftover from the running session (of type Any) and should be r_prime instead.

SVector and Type Stabilty

Specific Domains Modelling & Simulations

Domenico_Lahaye April 23, 2024, 12:56pm 3

const Point3D = SVector{3,Float64}

Avoiding allocations when normalizing a vector

Topic		Replies	Views
Type stability bug in compiler? General Usage question	13	601	January 27, 2023
Type Instability using StaticArrays that I dont get Performance	4	243	January 31, 2024
How can I fix the type instability Performance question	3	383	February 5, 2024
Why is this expression type unstable? General Usage question , type-stability , symbolics	1	918	May 9, 2022
Revisiting type stability General Usage type-stability	5	705	August 23, 2018

SVector and Type Stabilty

Related topics