Function sqrt_m1 

pub fn sqrt_m1() -> Fp

sqrt(-1), i.e. 2^((p-1)/4).

2 is a non-square modulo p, so 2^((p-1)/2) = -1 and its square root 2^((p-1)/4) is a square root of -1.