An alternative to StaticArrays is AbstractTensors.jl combined with Grassmann.jl . It’s an experiment where I am trying to build a new foundation for mathematixs in Julia language.
Poblem: StaticArrays.jl is too large of a dependency and has too many extra features.
To simplify my dependencies, I need a simplified variant of SVector, MVector, and SizedVector.
Solution: In AbstractTensors.jl they are replaced with Values, Variables, and FixedVector.
julia> using AbstractTensors # loads 10x faster than StaticArrays
julia> a = Values(1,2,3)
3-element Values{3,Int64} with indices SOneTo(3):
1
2
3
My goal is to have an implementation that is robust as in StaticArrays b…
This was because of a special explicit case for dimensions 1 and 2 and 3, which are now accounted for in Grassmann.jl also.
Also, support has been added for Moore-Penrose inverses for underdetermined and overdetermined linear systems. For underdetermined cases, the exterior product algorithm works ~20x faster than the SMatrix algorithm, and it is numerically stable. For overdetermined equations, the method used is based on the traditional normal equations, and this is prone to more numerical i…
The reason StaticArrays is not in Base is so that its version can be iterated seperately.
5 Likes