I was wondering if anyone had ever tried using cuBlas as a substitute for something like MKL with PETSc. I've been wondering if it would give better performance than MKL for my direct solves with cholmod even though the block sizes are small for cholmod i.e. 32x32 is the default I believe. If so, were there any tricky aspects to using cuBlas in this way? Thanks, Dave