def METHOD_NAME(typ_nlin, typ_nloc):
    tg = ccl.nl_pt.PTNumberCountsTracer(1., 0., 0., bk2=1.)
    tm = ccl.nl_pt.PTNumberCountsTracer(1., 0., 0.)
    ptc1 = ccl.nl_pt.LagrangianPTCalculator(
        b1_pk_kind=typ_nlin, bk2_pk_kind=typ_nloc, cosmo=COSMO)
    ptc2 = ccl.nl_pt.LagrangianPTCalculator(
        b1_pk_kind=typ_nloc, cosmo=COSMO)
    pkmm = ptc1.get_biased_pk2d(tm, tracer2=tm)
    pkmm2 = ptc2.get_biased_pk2d(tm, tracer2=tm)
    pkgg = ptc1.get_biased_pk2d(tg, tracer2=tg)
    ks = np.geomspace(1E-3, 1E1, 128)
    p1 = pkgg(ks, 1., cosmo=COSMO)
    p2 = pkmm(ks, 1., cosmo=COSMO)+ks**2*pkmm2(ks, 1., cosmo=COSMO)
    assert np.allclose(p1, p2, atol=0, rtol=1E-4)