Merge pull request #366 from martinholters/mh/mul_0.7

Fixes for mul! on Julia 0.7

Author: andyferris

src/matrix_multiply.jl

@@ -28,7 +28,7 @@ const StaticVecOrMat{T} = Union{StaticVector{<:Any, T}, StaticMatrix{<:Any, <:An
     @inline At_mul_Bt!(dest::StaticVecOrMat, A::StaticVecOrMat, B::StaticVecOrMat) = mul!(dest, transpose(A), transpose(B))
     @inline At_mul_B!(dest::StaticVecOrMat, A::StaticVecOrMat, B::StaticVecOrMat) = mul!(dest, transpose(A), B)
 else
-    import LinearAlgebra: BlasFloat, matprod
+    import LinearAlgebra: BlasFloat, matprod, mul!
 end
 
 
@@ -342,6 +342,32 @@ end
             gemm = :cgemm_
         end
 
+        if VERSION < v"0.7-"
+            blascall = quote
+                ccall((Base.BLAS.@blasfunc($gemm), Base.BLAS.libblas), Nothing,
+                    (Ptr{UInt8}, Ptr{UInt8}, Ptr{Base.BLAS.BlasInt}, Ptr{Base.BLAS.BlasInt},
+                     Ptr{Base.BLAS.BlasInt}, Ptr{$T}, Ptr{$T}, Ptr{Base.BLAS.BlasInt},
+                     Ptr{$T}, Ptr{Base.BLAS.BlasInt}, Ptr{$T}, Ptr{$T},
+                     Ptr{Base.BLAS.BlasInt}),
+                     &transA, &transB, &m, &n,
+                     &ka, &alpha, a, &strideA,
+                     b, &strideB, &beta, c,
+                     &strideC)
+             end
+         else
+             blascall = quote
+                 ccall((LinearAlgebra.BLAS.@blasfunc($gemm), LinearAlgebra.BLAS.libblas), Nothing,
+                     (Ref{UInt8}, Ref{UInt8}, Ref{LinearAlgebra.BLAS.BlasInt}, Ref{LinearAlgebra.BLAS.BlasInt},
+                      Ref{LinearAlgebra.BLAS.BlasInt}, Ref{$T}, Ptr{$T}, Ref{LinearAlgebra.BLAS.BlasInt},
+                      Ptr{$T}, Ref{LinearAlgebra.BLAS.BlasInt}, Ref{$T}, Ptr{$T},
+                      Ref{LinearAlgebra.BLAS.BlasInt}),
+                      transA, transB, m, n,
+                      ka, alpha, a, strideA,
+                      b, strideB, beta, c,
+                      strideC)
+            end
+        end
+
         return quote
             alpha = one(T)
             beta = zero(T)
@@ -355,15 +381,8 @@ end
             strideB = $(sb[1])
             strideC = $(s[1])
 
-            ccall((Base.BLAS.@blasfunc($gemm), Base.BLAS.libblas), Nothing,
-                (Ptr{UInt8}, Ptr{UInt8}, Ptr{Base.BLAS.BlasInt}, Ptr{Base.BLAS.BlasInt},
-                 Ptr{Base.BLAS.BlasInt}, Ptr{$T}, Ptr{$T}, Ptr{Base.BLAS.BlasInt},
-                 Ptr{$T}, Ptr{Base.BLAS.BlasInt}, Ptr{$T}, Ptr{$T},
-                 Ptr{Base.BLAS.BlasInt}),
-                 &transA, &transB, &m, &n,
-                 &ka, &alpha, a, &strideA,
-                 b, &strideB, &beta, c,
-                 &strideC)
+            $blascall
+
             return c
         end
     else