- replace the OpenMP based multi-core loop with a single-threaded numerical
  integration approximation of the sum in the posterior.  my test case
  yields essentially identical marginalized posterior results in a bit less
  time using only one core of a four-core machine that the OpenMP based
  loop saturates.