Normalises radians to (−pi, pi].
radians
(−pi, pi]
static double normaliseShort(double radians) { final w = radians % tau; // always in [0, 2*pi) return w > pi ? w - tau : w; // remap (pi, 2*pi) to (−pi, 0) }