Getting an openmp version working.  As expected, roughly factor of 2 increase on 2-4 cores.
Blocking out the algorithm for an MPI version