Galahad solver if the Jacobian has infinities

[tyronerees]

The first step to using the solver drqs is to form an eigendecomposition of the matrix, using the lapack routine dsyev. If the matrix has infinities in it, then dsyev hangs.

This happens if the problem CERI65IE is run with allow_fallback_method=.true. and the Hybrid option: the M-S method works, but dsyev goes into an infinite loop prior to sending to dtrs

Check for infinities in the Hessian? Is there a better solution?