38) and additionally generalized for applications with singular matrices ... and VSINGULAR= options to set three singularity criteria for the inversion of the Hessian matrix H. The singularity ...