Waveform Relaxation Time Domain Solver for Subsystem Arrays