I want to solve a linear system where the coefficient matrix is pentadiagonal block matrix with many of blocks are zero matrices. My matrix is stored in the device memory in row major format with the help of cudaMallocPitch. Any help ll be appreciated.

Thank You