Why is LinearAlgebra.pinv not type stable? Can it be implemented in a type stable way? Example (thanks to @albheim ) using LinearAlgebra using JET function test(M) A = rand(M, M + 2) pinv(A) end test(10) report = @report_opt test(10) println(report) … ││┌ pinv(A::Matrix{Float64}; atol:…

Methods with keyword arguments are implemented to call hidden helper methods that have the written body, so @code_warntype is really difficult to use with them. Well, I dug out the underlying call and am including the part where things get a little unstable, even though the internal arrays and the …

Pinv not Type Stable?

General Usage Performance

Zentrik September 16, 2023, 8:19pm 14

Yh it would be useful to get this fixed as it seems to cause pinv to not work on the GPU, Support for LinearAlgebra.pinv · Issue #2070 · JuliaGPU/CUDA.jl · GitHub.

2 Likes

Topic		Replies	Views
Type Instability using StaticArrays that I dont get Performance	4	227	January 31, 2024
Bug in LinearAlgebra pinv() General Usage linearalgebra	9	487	February 18, 2022
Type stability with chained parametric types General Usage parametric-types , type-stability	2	498	August 10, 2018
Type stability of vector of matrices General Usage performance	4	223	April 29, 2024
Pinv of matrix transpose? General Usage linearalgebra	5	725	April 6, 2019

Pinv not Type Stable?

Related topics