Comparison of numerical methods for modeling laser mode locking with saturable gain