def is_real_pole(s):
    """Apply a threshold to determine if a pole should be treated as purely
    real, in which case its conjugate should not be used in models"""
    return abs(s.imag) < 0.001 * abs(s.real)
