Use of the squared dirac operator in variational relativistic calculations