“…In this framework, a constant demand exists for reliable and efficient numerical tools to characterize such media over large frequency bandwidth (from direct current to tens of GHz), including but not restricted to applications such as electromagnetic interference/compatibility (EMI/EMC), 5G networks, more electrical aircrafts/automobiles. The state-of-the-art regarding fullwave electromagnetic (EM) modeling shows that various methods were successfully used to simulate composite fiber-reinforced media including models in frequency domain (finite element method [1], integral equation [2]), and time domain (finite difference [3], discontinuous Galerkin [4], finite integral technique (FIT) [5]). However, due to the intrinsic variability of composite materials (e.g., conductivity, locations or shapes of inclusions, .…”