[petsc-users] MatMult doesn't give the correct result