<div dir="ltr">Does PETSc include an efficient implementation for the operation<br>y=(A-diag(A))x or y_i=\sum_{j!=i}A_{ij}x_j on a sparse matrix A?<br><br>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.<br><br>Thanks<br>Ingo<br></div>