static sf_double_complex dsdot(int n, const sf_complex *cx, const sf_complex *cy) /* Hermitian dot product */ { ... }