Random field Ising model in two dimensions: Bethe approximation, cluster variational method and message passing algorithms