Update BLAS wrappers

Author: amontoison

src/blas.jl

@@ -86,6 +86,16 @@ export
 
 using ..LinearAlgebra: libblastrampoline, BlasReal, BlasComplex, BlasFloat, BlasInt, DimensionMismatch, checksquare, chkstride1
 
+if USE_BLAS64
+    macro blasfunc(x)
+        return Expr(:quote, Symbol(x, "64_"))
+    end
+else
+    macro blasfunc(x)
+        return Expr(:quote, x)
+    end
+end
+
 include("libblas.jl")
 include("lbt.jl")
 
@@ -108,16 +118,6 @@ Return an object representing the current `libblastrampoline` configuration.
 """
 get_config() = lbt_get_config()
 
-if USE_BLAS64
-    macro blasfunc(x)
-        return Expr(:quote, Symbol(x, "64_"))
-    end
-else
-    macro blasfunc(x)
-        return Expr(:quote, x)
-    end
-end
-
 _tryparse_env_int(key) = tryparse(Int, get(ENV, key, ""))
 
 

src/libblas.jl

@@ -44,6 +44,12 @@ function caxpy(n, ca, cx, incx, cy, incy)
                   Ref{BlasInt}), n, ca, cx, incx, cy, incy)
 end
 
+function caxpby(n, ca, cx, incx, cb, cy, incy)
+    return ccall((@blasfunc(caxpby_), libblastrampoline), Cvoid,
+                 (Ref{BlasInt}, Ref{ComplexF32}, Ptr{ComplexF32}, Ref{BlasInt}, Ref{ComplexF32},
+                  Ptr{ComplexF32}, Ref{BlasInt}), n, ca, cx, incx, cb, cy, incy)
+end
+
 function ccopy(n, cx, incx, cy, incy)
     return ccall((@blasfunc(ccopy_), libblastrampoline), Cvoid,
                  (Ref{BlasInt}, Ptr{ComplexF32}, Ref{BlasInt}, Ptr{ComplexF32}, Ref{BlasInt}), n, cx,
@@ -293,6 +299,12 @@ function daxpy(n, da, dx, incx, dy, incy)
                   Ref{BlasInt}), n, da, dx, incx, dy, incy)
 end
 
+function daxpby(n, da, dx, incx, db, dy, incy)
+    return ccall((@blasfunc(daxpby_), libblastrampoline), Cvoid,
+                 (Ref{BlasInt}, Ref{Float64}, Ptr{Float64}, Ref{BlasInt}, Ref{Float64},
+                  Ptr{Float64}, Ref{BlasInt}), n, da, dx, incx, db, dy, incy)
+end
+
 function dcabs1(z)
     return ccall((@blasfunc(dcabs1_), libblastrampoline), Float64, (Ref{ComplexF64}), z)
 end
@@ -546,6 +558,12 @@ function saxpy(n, sa, sx, incx, sy, incy)
                   Ref{BlasInt}), n, sa, sx, incx, sy, incy)
 end
 
+function saxpby(n, sa, sx, incx, sb, sy, incy)
+    return ccall((@blasfunc(saxpby_), libblastrampoline), Cvoid,
+                 (Ref{BlasInt}, Ref{Float32}, Ptr{Float32}, Ref{BlasInt}, Ref{Float32},
+                  Ptr{Float32}, Ref{BlasInt}), n, sa, sx, incx, sb, sy, incy)
+end
+
 function scabs1(z)
     return ccall((@blasfunc(scabs1_), libblastrampoline), Float32, (Ref{ComplexF32}), z)
 end
@@ -774,6 +792,12 @@ function zaxpy(n, za, zx, incx, zy, incy)
                   Ref{BlasInt}), n, za, zx, incx, zy, incy)
 end
 
+function zaxpby(n, za, zx, incx, zb, zy, incy)
+    return ccall((@blasfunc(zaxpby_), libblastrampoline), Cvoid,
+                 (Ref{BlasInt}, Ref{ComplexF64}, Ptr{ComplexF64}, Ref{BlasInt}, Ref{ComplexF64},
+                  Ptr{ComplexF64}, Ref{BlasInt}), n, za, zx, incx, zb, zy, incy)
+end
+
 function zcopy(n, zx, incx, zy, incy)
     return ccall((@blasfunc(zcopy_), libblastrampoline), Cvoid,
                  (Ref{BlasInt}, Ptr{ComplexF64}, Ref{BlasInt}, Ptr{ComplexF64}, Ref{BlasInt}), n, zx,