Hello, I've run into an issue where the opRestriction and in turn the opExtension constructors don't have an S keyword argument. While the operator on its own works with GPU arrays, this causes an issue when trying to build composite ops for a GPU:
julia> using LinearOperators, CUDA
julia> x = CUDA.rand(256);
julia> op = opRestriction(1:128, 256);
julia> op * x # works
julia> (op * opDiagonal(x)) * x # fails