I've also found some independent resources describing this exact
algorithm. In particular, one of them also offers a performance tip: we
can use fract() instead of x - floor() since we don't need the exact
modulus semantics.