How can I simulate variates from a truncated inverse gaussian distribution?