[petsc-users] Off-diagonal matrix-vector product y=(A-diag(A))x

Ingo Gaertner ingogaertner.tus at gmail.com
Fri Apr 14 03:55:55 CDT 2017


Does PETSc include an efficient implementation for the operation
y=(A-diag(A))x or y_i=\sum_{j!=i}A_{ij}x_j on a sparse matrix A?

In words, I need a matrix-vector product after the matrix diagonal has been
set to zero. For efficiency reasons I can't copy and modify the matrix or
first calculate the full product and then subtract the diagonal
contribution.

Thanks
Ingo
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20170414/d4d60aae/attachment.html>


More information about the petsc-users mailing list