Function phase 

pub(crate) fn phase(vector: [Complex; 3], point: [f64; 3]) -> Complex

exp(i k·r) for a complex wavevector k and a real point r, from the real and imaginary parts of k·r separately.