If you call lu on an SMatrix, you get a static LU factorization—cool: using StaticArrays using LinearAlgebra A = SMatrix{3, 3}(rand(3, 3)) lu(A) # StaticArrays.LU # L factor: # 3×3 LowerTriangular{Float64, SMatrix{3, 3, Float64, 9}}: # 1.0 ⋅ ⋅ # 0.706739 1.0 ⋅ # 0.64896…

[image] maxkapur: But sometimes, you want to preallocate space for the L and U factors. Preallocation doesn’t make sense for SMatrix — there is no “allocation” cost in the first place. So what is your motivation?

How to preallocate a static array for use with `lu!`

General Usage

stevengj Split this topic November 3, 2022, 2:28pm 5

2 posts were split to a new topic: How to preallocate LU factorization

Topic		Replies	Views
How to preallocate LU factorization Numerics linearalgebra	3	277	November 3, 2022
Allocations while using static arrays in a struct General Usage memory-allocation	3	1569	December 12, 2018
Allocations when calling eigen for hermitian 3x3 SMatrix General Usage linearalgebra , staticarrays	5	507	May 15, 2022
Operations on `StaticArrays` without intermediate allocation General Usage	9	1226	June 13, 2017
Allocation during construction of compile-time sized structure Performance staticarrays	1	333	September 18, 2022

How to preallocate a static array for use with `lu!`

Related topics