fix a small numerical bug
try to parallelize, but fail due to qhull
re-compute pairwise IBS only when interacting
commend out debugbox