Function sinmx_accurate 

fn sinmx_accurate(
    y: DoubleDouble,
    sin_k: DoubleDouble,
    cos_k: DoubleDouble,
    x: f64,
) -> f64