algorithm should work

add openmp and run on cluster

then add importance sampling