On the computation of solution spaces in high dimensions