“…Building on the decoupling between linear electromagnetic and lumped nonlinear blocks [13], a natural approach is to apply a waveform relaxation (WR) algorithm [17], [18], that starting from an initial estimate of the transient port signals to be computed over a given time interval, iteratively refines such estimates until convergence. The procedure stabilizes if the iteration operator is a contraction, regardless of the selected formulation [18], [19], [20], [21]; the number of iterations can be reduced by optimizing the WR decoupling parameters [22], [23], [24]. A general procedure that allows such optimizations is the definition of some reference impedance level, used to build the decoupling network so that exchange of signals between iterations is performed through equivalent scattering waves referenced to this optimal impedance.…”