Completely Positive Trace Preserving Numerical Methods for Long-Term Generalized Maxwell-Bloch Simulations