Fix the order in which recursion occurs so that a) we don't end up dividing by zero if |z| < epsilon, and b) don't lose precision when z is close to epsilon.
Fix test tolerances for new tgamma code.