CUBLAS: Incorrect inversion for matrix with zero pivot
There seems to be a bug in the current CUBLAS library implementation of cublas<t>getrfBatched for matrices of dimension (n) such that 3<=n<=16, when there is a “zero pivot” as you say. A possible workaround is to “identity-extend” your A matrix to be inverted, when n<17, to a size of 17×17 (using matlab nomenclature): LU = … Read more