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 November 3, 2022, 1:45pm 2

Preallocation doesn’t make sense for SMatrix — there is no “allocation” cost in the first place. So what is your motivation?

Topic		Replies	Views
How to preallocate LU factorization Numerics linearalgebra	3	336	November 3, 2022
In place lufact for sparse matrices and prellocation for lufact output Numerics performance , linearalgebra	1	624	September 26, 2017
Replacing normal arrays with static arrays New to Julia linearalgebra	6	1417	August 31, 2019
MArray or SizedArray in StaticArrays package General Usage	5	2503	February 24, 2017
How to reduce large memory allocation by lu! with sparse matrices Performance question , memory-allocation , ode , sparse	8	1501	July 17, 2021

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

Related topics