Fast Preconditioner Computation for BICGSTAB-FFT Method of Moments with NURBS in Large Multilayer Structures

Fast computation of the coefficients of the reduced impedance matrix of the method of moment (MM) is proposed by expanding the basis functions (BFs) in pulses and solving an equivalent periodic problem (EPP) for analyzing large multilayer structures with non-uniform rational basis spline (NURBS) modeling of the embedded layout. These coefficients are required by the computation of sparse approximate inverse (SAI) preconditioner, which leads an efficient iterative version of the MM. This reduced coefficient matrix only considers the near field part of the MM matrix. Discrete functions of small sizes are required to implement the pulse expansion and EPP. These discrete functions of small size lead to discrete cyclic convolutions that are computed in a very fast way by fast Fourier transform (FFT)-accelerated matrix-vector multiplication. Results obtained using a conventional laptop show an analysis of very large multilayer structures with resonant layouts, as whole reflectarrays of electrical size 40 times the vacuum wavelengths, where the iterative MM with a SAI preconditioner can be 22.7 times faster than the pure iterative MM without any preconditioner.