Structure-Preserving Integration for Magnetic Gaussian Wave Packet Dynamics