Add generalized hitandrun() interface

 - handles both equality and inequality constraints transparently
 - also add simplexConstraints() to make sampling from the simplex easier
 - documentation updated to reflect changes