Barry Smith <[email protected]> writes: > You can do MatGetDiagonal() then do a MatMult() followed by > subtracting a VecPointsizeMult().
I would use VecPointsizeMult followed by MatMultAdd. One fewer traversal of a vector.
signature.asc
Description: PGP signature
