adding richardson relaxation to monte carlo solver