Compressed hierarchical Schur algorithm for frequency-domain analysis of photonic structures. (English) Zbl 1412.78008

Summary: Three-dimensional finite-difference frequency-domain analyses of partially periodic photonic structures result in large-scale ill-conditioned linear systems. Due to the lack of efficient preconditioner and reordering scheme, existed general-purpose iterative and direct solvers are inadequate to solve these linear systems in time or memory. We propose an efficient direct solver to tackle this problem. By exploring the physical properties, the coefficient matrix structure, and hardware computing efficiency, we extend the concepts of grid geometry manipulation and multi-level Schur method to propose the Compressed Hierarchical Schur algorithm (CHiS). The proposed CHiS algorithm can use less memory and remove redundant computational workloads due to the homogeneity and periodicity of photonic structures. Moreover, CHiS relies on dense BLAS3 operations of sub-matrices that can be computed efficiently with strong scalability by the latest multicore processors or accelerators. The implementation and benchmarks of CHiS demonstrate promising memory usage, timing, and scalability results. The feasibility of future hardware acceleration for CHiS is also addressed using computational data. This high-performance analysis tool can improve the design and modeling capability for various photonic structures.


78M20 Finite difference methods applied to problems in optics and electromagnetic theory
65F05 Direct numerical methods for linear systems and matrix inversion
65Y05 Parallel numerical computation
Full Text: DOI Euclid