fix: handle non-finite SED flux

                wave = self.lamb_qso * (1.0 + obj.z)

        else:

            raise ValueError("Object type not known")

        # check SED table

        bad = ~(np.isfinite(wave) & np.isfinite(flux))

        if np.any(bad):

            idx = int(np.argmax(bad))

            raise ValueError(

                f"[WARN] Skip {obj.type} (id={obj.id}): SED has {int(bad.sum())} "

                f"non-finite flux (first at wave={wave[idx]:g}, flux={flux[idx]:g})"

            )

        speci = interpolate.interp1d(wave, flux, fill_value=0.0, bounds_error=False)

        lamb = np.arange(2000, 11001 + 0.5, 0.5)

        y = speci(lamb)