Text this: Fast numerical integration of relaxation oscillator networks based on singular limit solutions