use new formula to calculate jacobi rotation constant, refactor comments