Function mulmod_stable 

pub fn mulmod_stable<T: Float>(a: T, b: T, m: T) -> CoreResult<T>

Stable computation of (a*b) % m avoiding overflow