Evaluating multi-loop Feynman diagrams with infrared and threshold singularities numerically