defined SetOperatorMat.

This function sets the matrix of the linear system to solve. we only
have to set it up only once.