def METHOD_NAME(self, freq):
    """Compute the FD delay at an array of frequencies."""
    FD_mapping = self.get_prefix_mapping_component("FD")
    log_freq = np.log(freq.to(u.GHz).value)
    non_finite = np.invert(np.isfinite(log_freq))
    log_freq[non_finite] = 0.0
    FD_coeff = [
        getattr(self, FD_mapping[ii]).value
        for ii in range(self.num_FD_terms, 0, -1)
    ]
    FD_coeff += [0.0]  # Zeroth term of polynomial
    FD_delay = np.polyval(FD_coeff, log_freq)
    return FD_delay * self.FD1.units