4-6 * skip # header * * either abs=4e-11 rel=1e-18 # global 9-13 7 abs=1e-8 # canonical momentum