func PolarDiscriminator(a, b complex128) float64 {
	return cmplx.Phase(a * cmplx.Conj(b))
}
PolarDiscriminator returns the phase angle between two complex vectors
equivalent to arg(a * conj(b)). The returned angle is in the range [-Pi, Pi].
package dsp
import "math/cmplx"