Hi, Is sparsity of Vec exploited in MatMult? I'm implementing some domain decomposition methods. The residual is non-zero or significant only at very few nodes. Does MatMult knows that most entries of the Vec are zeros? Is there a threshold to regard an value as zero? Thanks, Hui