This paper introduces an extension of the time-splitting sine-spectral (TSSP)
method for solving damped focusing nonlinear Schr\"{o}dinger equations (NLS).
The method is explicit, unconditionally stable and time transversal invariant.
Moreover, it preserves the exact decay rate for the normalization of the wave
function if linear damping terms are added to the NLS. Extensive numerical
tests are presented for cubic focusing nonlinear Schr\"{o}dinger equations in
2d with a linear, cubic or a quintic damping term. Our numerical results show
that quintic or cubic damping always arrests blowup, while linear damping can
arrest blowup only when the damping parameter \dt is larger than a threshold
value \dt_{\rm th}. We note that our method can also be applied to solve the
3d Gross-Pitaevskii equation with a quintic damping term to model the dynamics
of a collapsing and exploding Bose-Einstein condensate (BEC).Comment: SIAM Journal on Numerical Analysis, to appea