Regularized Numerical Methods for the Nonlinear Schrödinger Equation with Singular Nonlinearity