def METHOD_NAME(ctx, lb_method, params):
    config = generate_config(ctx, params)
    fields = generate_fields(config, lb_method.stencil)
    initial_rho = sp.Symbol('rho_0')
    pdfs_setter = lbmpy.macroscopic_value_kernels.macroscopic_values_setter(
        lb_method,
        initial_rho,
        fields['velocity'].center_vector,
        fields['pdfs'].center_vector)
    return pdfs_setter