def METHOD_NAME():
    # Check P(k) for linear tracer with b1=1 is the same
    # as matter P(k)
    ptc = ccl.nl_pt.LagrangianPTCalculator(cosmo=COSMO,
                                           b1_pk_kind='pt')
    tg = ccl.nl_pt.PTNumberCountsTracer(b1=1.0)
    tm = ccl.nl_pt.PTMatterTracer()
    ks = np.geomspace(1E-3, 10, 64)
    pk1 = ptc.get_biased_pk2d(tm)(ks, 1.0, cosmo=COSMO)
    pk2 = ptc.get_biased_pk2d(tg)(ks, 1.0, cosmo=COSMO)
    assert np.allclose(pk1, pk2, atol=0, rtol=1E-3)