MatZeroRows(), MatZeroRowsIS() and MatZeroRowsLocal(), MatZeroRowsLocalIS() now take optional x and b vectors. The b[index[]] is set to diag*x[index[]] if the b and x are provided. Barry