“…There are still many technical problems to be solved before BOR-FDTD can realize such a complex electromagnetic field problem [19], [20]. In order to resolve the above problems, numerous of BOR-FDTD methods are proposed, for example, ADI (alternating-direction implicit)-BOR-FDTD [21], [22], LOD (locally one-dimensional)-BOR-FDTD [23], [24], [25], WCS (weakly conditional stability)-BOR-FDTD [26], [27], HIE (hybrid implicit-explicit)-BOR-FDTD [28], [29] and WLP (weighted Laguerre polynomial)-BOR-FDTD [30], [31], [32] and so on. They can well remove the CFL (courant-friedrichslewy) stability condition of the BOR-FDTD method and make the simulation toward BOR structure more efficient, however, the ADI, LOD, WCS and HIE schemes will lead a large numerical dispersion error when the time step is large, and the WLP one will lead to a large amount of memory.…”