Function double 

pub fn double<Fe>(_: (Fe, Fe, Fe), curve_equation_b: Fe) -> (Fe, Fe, Fe)
where
    Fe: Field,

Implements the exception-free point doubling formula from Renes-Costello-Batina 2015 (Algorithm 6).