Massively Parallel Hybrid TLM-PEEC Solver and Model Order Reduction for 3D Nonlinear Electromagnetic Transient Analysis