We analyzed the largest eigenvalues λ max of the A n matrix for the standard CHARMM36 cholesterol model as well as the HMR and VIS ones using our script. This shows that also for atomistic systems in which all bonds are constrained, the analysis of the eigenvalues of the A n matrix is a valuable diagnosis tool. Our script allows detecting potential convergence issues of the LINCS algorithm and estimating the required LINCS settings. We emphasize that the limitation of the Gromacs MD engine to evaluate the eigenvalues is not a conceptual but a practical matter that stems from its highly optimized nature.