def METHOD_NAME(self):
    """Creates an interpolation grid in H_0, omega_m and computes quantities in Dd
    and Ds_Dds.
    :return:
    """
    grid2d = np.dstack(np.meshgrid(self._H0_range, self._omega_m_range)).reshape(
        -1, 2
    )
    H0_grid = grid2d[:, 0]
    omega_m_grid = grid2d[:, 1]
    Dd_grid = np.zeros_like(H0_grid)
    Ds_Dds_grid = np.zeros_like(H0_grid)
    for i in range(len(H0_grid)):
        Dd, Ds_Dds = cosmo2angular_diameter_distances(
            H0_grid[i], omega_m_grid[i], self.z_d, self.z_s
        )
        Dd_grid[i] = Dd
        Ds_Dds_grid[i] = Ds_Dds
    self._f_H0 = interpolate.interp2d(
        Dd_grid,
        Ds_Dds_grid,
        H0_grid,
        kind="linear",
        copy=False,
        bounds_error=False,
        fill_value=-1,
    )
    print("H0 interpolation done")
    self._f_omega_m = interpolate.interp2d(
        Dd_grid,
        Ds_Dds_grid,
        omega_m_grid,
        kind="linear",
        copy=False,
        bounds_error=False,
        fill_value=0,
    )
    print("omega_m interpolation done")