def METHOD_NAME(sampler, machine, parameters, state):
    pdf = jnp.absolute(
        to_array(sampler.hilbert, machine.apply, parameters) ** sampler.machine_pow
    )
    pdf = pdf / pdf.sum()
    return state.replace(pdf=pdf)