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 Instability using StaticArrays that I dont get Performance	4	225	January 31, 2024
Type instability with tuples of functions and parameters Performance question , code_warntype , type-stability	13	627	May 25, 2022
Type unstable in NonlinearSolve.jl Optimization (Mathematical) code_warntype , type-stability	3	75	September 18, 2024
Strange type instability General Usage	2	330	December 22, 2023
Ominous type instability Performance bug , code_warntype , type-stability	10	1170	August 20, 2022

SVector and Type Stabilty

Related topics