Fixed reducedim operations for CuQRPackedQ

Author: Ralph Kube

lib/cusolver/linalg.jl

@@ -27,6 +27,14 @@ CUDA.CuMatrix(A::CuQRPackedQ) = orgqr!(copy(A.factors), A.τ)
 CUDA.CuArray(A::CuQRPackedQ) = CuMatrix(A)
 Base.Matrix(A::CuQRPackedQ) = Matrix(CuMatrix(A))
 
+# CuQRPackedQ isn't an AbstractGPUArray, so we miss out on lots of GPU-compatible methods.
+# add a couple common methods here that avoid GPU-incompatible fallbacks.
+Base.similar(A::CuQRPackedQ{S,T}) where {S,T} = CuArray{S, ndims(a)}(undef, size(a))
+Base.sum(xs::CuQRPackedQ) = sum(CuArray(xs))
+Base.prod(xs::CuQRPackedQ) = prod(CuArray(xs))
+Base.maximum(xs::CuQRPackedQ) = maximum(CuArray(xs))
+Base.minimum(xs::CuQRPackedQ) = minimum(CuArray(xs))
+
 function Base.getproperty(A::CuQR, d::Symbol)
     m, n = size(getfield(A, :factors))
     if d == :R