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?
The first step to using the solver
drqsis to form an eigendecomposition of the matrix, using the lapack routinedsyev. If the matrix has infinities in it, thendsyevhangs.This happens if the problem
CERI65IEis run withallow_fallback_method=.true.and the Hybrid option: the M-S method works, butdsyevgoes into an infinite loop prior to sending todtrsCheck for infinities in the Hessian? Is there a better solution?