A General Method for the Numerical Computation of Manipulator Singularity Sets