Merge 04145c4 into 83d2440

src/pdmat.jl

@@ -34,7 +34,7 @@ function pdadd!(r::Matrix, a::Matrix, b::PDMat, c)
     _addscal!(r, a, b.mat, c)
 end
 
-*{T<:Real}(a::PDMat{T}, c::T) = PDMat(a.mat * c)
+*{S<:Real, T<:Real}(a::PDMat{S}, c::T) = PDMat(a.mat * c)
 *(a::PDMat, x::StridedVecOrMat) = a.mat * x
 \(a::PDMat, x::StridedVecOrMat) = a.chol \ x
 

test/generics.jl

@@ -0,0 +1,24 @@
+
+# test operators with pd matrix types
+using PDMats
+using Base.Test
+
+# test scalar multiplication 
+print_with_color(:blue, "Testing scalar multiplication\n")
+pm1 = PDMat(eye(3))
+pm2 = PDiagMat(ones(3))
+pm3 = ScalMat(3,1)
+
+pm1a = PDMat(3.0 .* eye(3))
+pm2a = PDiagMat(3.0 .* ones(3))
+pm3a = ScalMat(3, 3)
+
+pmats = Any[pm1, pm2, pm3]
+pmatsa= Any[pm1a,pm2a,pm3a]
+
+for i in 1:length(pmats)
+    @test full(3.0 * pmats[i]) == full(pmatsa[i])
+    @test full(pmats[i] * 3.0) == full(pmatsa[i])
+    @test full(3 * pmats[i])   == full(pmatsa[i])
+    @test full(pmats[i] * 3)   == full(pmatsa[i])
+end

test/runtests.jl

@@ -1,4 +1,4 @@
-tests = ["pdmtypes", "addition"]
+tests = ["pdmtypes", "addition", "generics"]
 println("Running tests ...")
 
 for t in tests