Once the iterative diagonalization algorithm has converged, the matrix of displacement vectors S and the corresponding potential energy Hessian-vector products Y E are projected into the full redundant internal coordinate basis and used to update the approximate potential energy Hessian B. This is done with the multi-secant TS-BFGS Hessian update formula, 1,28-31