The celerity corresponding to speed \(u\) is defined as
  \(u\gamma\) and the rapidity is
  \(c\cdot\mathrm{atanh}(u/c)\).
Functions celerity_ur() and rapidity_ur() are used for
  the ultrarelativistic case where speeds are very close to the speed of
  light.  Its argument d is the deficit, that is, \(d=c-v\)
  where \(v\) is the speed of the transformation.  Algebraically,
  celerity_ur(c-v) == celerity(v), but if \(d=1-v/c\) is small
  the result of celerity_ur() is more accurate than that of
  celerity().
Things get a bit sticky for celerity and rapidity if \(c\neq
  1\).  The guiding principle in the package is to give the
  celerity and rapidity the same units as \(c\), so if \(u\ll
  c\) we have that all three of celerity(u),
  rapidity(u) and u are approximately equal.  Note
  carefully that, in contrast, \(\gamma\) is dimensionless.  Also
  observe that d in functions celerity_ur() and
  rapidity_ur() has the same units as \(c\).