Merge pull request #308 from mschauer/fixsolve

Fix bug in 1x1 solve

Author: c42f

src/solve.jl

@@ -4,7 +4,7 @@
 @inline solve(::Size, ::Size, a, b) = inv(a) * b
 
 @inline function solve(::Size{(1,1)}, ::Size{(1,)}, a::StaticMatrix{<:Any, <:Any, Ta}, b::StaticVector{<:Any, Tb}) where {Ta, Tb}
-    @inbounds return similar_type(b, typeof(b[1] \ a[1]))(b[1] \ a[1])
+    @inbounds return similar_type(b, typeof(a[1] \ b[1]))(a[1] \ b[1])
 end
 
 @inline function solve(::Size{(2,2)}, ::Size{(2,)}, a::StaticMatrix{<:Any, <:Any, Ta}, b::StaticVector{<:Any, Tb}) where {Ta, Tb}

test/solve.jl

@@ -5,8 +5,8 @@ using StaticArrays, Base.Test
             (m, v) in ((SMatrix{n,n}, SVector{n}), (MMatrix{n,n}, MVector{n})),
                 elty in (Float64, Int)
 
-        A = elty.(rand(-99:2:99,n,n))
-        b = A*ones(elty,n)
+        A = elty.(rand(-99:2:99, n, n))
+        b = A * elty.(rand(2:5, n))
         @test m(A)\v(b) ≈ A\b
     end