let c_rem u v =
    sub u (mul v (c_div u v))