Code for sampling points on positive dimensional solution sets
in double double and quad double precision added.