  An algorithm for calculating three gauge-invariant helicities (self-, mutual-
and Berger relative helicity) for a magnetic field specified in a rectangular
box is described. The algorithm is tested on a well-known force-free model (Low
and Lou, 1990) presented in vector-potential form.
