Function richelot 

pub fn richelot(g: &[Quad; 3], p: &BigInt) -> Richelot

The Richelot (2,2)-isogeny from a quadratic splitting G = (G₁, G₂, G₃) over 𝔽_p.