Functions or packages to implement Belief propagation?

I look for an efficient way to implement a Belief Propagation (more specific, parallel implementation). Does anyone know about any useful functions or packages?