This work presents highly parallelized, GPU-accelerated computer software that performs the Bayesian inference of EM (BioEM) formalism efficiently and scales nearly ideally both on pure CPU and on CPU+GPU architectures, thus enabling Bayesian analysis of tens of thousands of images in a reasonable time.