put sumTol(drBarPars, drBarPars) into physpars.{h,cpp}

This avoids a circular dependence between utils.{h,cpp} and
physpars.{h,cpp}.  Now physpars depends on utils (via linalg),
but utils does not depend on physpars.