Parametrization Jun 1, 2018 fitting penalty=0.1 transition functions penalty=10 (+local Gaussian distance weighting) optimization took 27’11’’ varying scale