Add basic tests for some linear algebra operations (#24)

* Add basic tests for some linear algebra operations

* fix for `Float32` and `Float64`

As suggested by Giordano.

xref #25

---------

Co-authored-by: Neven Sajko <[email protected]>

Author: giordano

src/FixedSizeArrays.jl

@@ -114,4 +114,8 @@ Base.@propagate_inbounds Base.copyto!(dst::FixedSizeArray, src::FixedSizeArray)
 Base.@propagate_inbounds Base.copyto!(dst::FixedSizeArray, src::Memory        ) = copyto2!(dst, src)
 Base.@propagate_inbounds Base.copyto!(dst::Memory        , src::FixedSizeArray) = copyto2!(dst, src)
 
+# unsafe: the native address of the array's storage
+
+Base.unsafe_convert(::Type{Ptr{T}}, a::FixedSizeArray{T}) where {T} = Base.unsafe_convert(Ptr{T}, a.mem)
+
 end # module FixedSizeArrays

test/runtests.jl

@@ -149,4 +149,26 @@ import Aqua
             end
         end
     end
+
+    @testset "LinearAlgebra" begin
+        @testset "$T" for T in (Float16, Float32, Float64)
+            v = randn(T, 8)
+            v_fixed = FixedSizeVector(v)
+            v_sum = @inferred v_fixed + v_fixed
+            @test v_sum isa FixedSizeVector{T}
+            @test v_sum ≈ v + v
+            v_mul = @inferred v_fixed * v_fixed'
+            @test v_mul isa FixedSizeMatrix{T}
+            @test v_mul ≈ v * v'
+
+            M = randn(T, 8, 8)
+            M_fixed = FixedSizeMatrix(M)
+            M_sum = @inferred M_fixed + M_fixed
+            @test M_sum isa FixedSizeMatrix{T}
+            @test M_sum ≈ M + M
+            M_mul = @inferred M_fixed * M_fixed
+            @test M_mul isa FixedSizeMatrix{T}
+            @test M_mul ≈ M * M
+        end
+    end
 end