Analysis of Millimeter Wave Scattering by an Electrically Large Metallic Grating Using Wavelet-Based Algebraic Multigrid Preconditioned CG Method

An effective wavelet based multigrid preconditioned conjugate gradient method is developed to solve electromagnetic large matrix problem for millimeter wave scattering application. By using wavelet transformation we restrict the large matrix equation to a relative smaller matrix and which can be solved rapidly. The solution is prolonged as the new improvement for the conjugate gradient (CG) method. Numerical results show that our developed wavelet based multigrid preconditioned CG method can reach large improvement of computational complexity. Due to the automaticity of wavelet transformation, this method is potential to be a block box solver without physical background.     

An Effective Multigrid Preconditioned CG Algorithm for Millimeter Wave Scattering by an Infinite Plane Metallic Grating

An effective multigrid based preconditioned conjugate gradient method is developed to solve electromagnetic large matrix problem for millimeter wave scattering application. By using multigrid technique we restrict the large matrix equation to a relative smaller matrix and which can be solved rapidly. The solution is prolonged as the initial guess for the conjugate gradient (CG) method. Numerical results show that our developed method can reach five times improvement of computational complexity.     

间断问题的并行多重网格预条件共轭梯度法求解

首先比较了四种串行多重网格预条件共轭梯度法求解四类典型间断问题时,数值收敛因子与间断类型、间断程度、网格步长的关系;然后,在算法中引入Schwarz并行,探讨了并行算法收敛因子与子区域个数的关系,以及Schwarz并行的获益。给出了详细的性能比较结果。     

Research of the Algebraic Multigrid Method for Electron Optical Simulator

At present, electron optical simulator (EOS) takes a long time to solve linear FEM systems. The algebraic multigrid preconditioned conjugate gradient (AMGPCG) method can improve the efficiency of solving systems. This paper is focused on the implementation of the AMGPCG method in EOS. The aggregation-based scheme, which uses two passes of a pairwise matching algorithm and the K-cyle scheme, is adopted in the aggregation-based algebraic multigrid method. Numerical experiments show the advantages and disadvantages of the AMG algorithm in peak memory and solving efficiency. The AMGPCG is more efficient than the iterative methods used in the past and only needs one coarsening when EOS computes the particle motion trajectory.     

Object plane detection and phase-amplitude imaging based on transport of intensity equation

By using transport of intensity equation (TIE), phase distribution of an object is retrieved from through-focus intensity images. This technique allows simple and robust phase imaging compared with an interferometric approach. However, it is hard to measure phase distribution when a dynamic object moves in the direction of an optical axis. To clear this problem, autofocusing TIE which is based on local statistics is proposed. The proposed technique achieves the detection of the object plane and the retrieval of a focused object phase distribution simultaneously. In this approach, an object plane is determined by the focusing techniques based on local statistics such as variance, gradient, and Laplacian of amplitude distribution, after phase distribution in an image sensor plane is retrieved by the TIE. The performance of these three statistics is evaluated in numerical and optical experiments, and a suitable focus value is determined for precise phase imaging.     

An h-adaptive finite element solver for the calculations of the electronic structures

In this paper, a framework of using h-adaptive finite element method for the Kohn–Sham equation on the tetrahedron mesh is presented. The Kohn–Sham equation is discretized by the finite element method, and the h-adaptive technique is adopted to optimize the accuracy and the efficiency of the algorithm. The locally optimal block preconditioned conjugate gradient method is employed for solving the generalized eigenvalue problem, and an algebraic multigrid preconditioner is used to accelerate the solver. A variety of numerical experiments demonstrate the effectiveness of our algorithm for both the all-electron and the pseudo-potential calculations.     

Application of Robust Incomplete Cholesky Factorizations Preconditioned CG Algorithms for FEM Analysis of Millimeter Wave Filters

The Incomplete Cholesky factorizations preconditioning scheme is applied to the conjugate gradient (CG) method for solving a large system of linear equations resulting from finite element method (FEM) analysis of millimeter wave filters. As is well known, the convergence of CG method deteriorates with increasing EM wave number and in millimeter wave band the eigen-values of A are more and more scattered between both the right and the left half-plane. The efficient implementation of this preconditioned CG (PCG) algorithm is described in details for Complex coefficient matrix. With incomplete factorization preconditioning scheme in the conjugate gradient algorithm, this PCG approach can reach convergence in 20 times CPU time shorter than CG for several typical millimeter wave structures.     

Restricted Additive Schwarz Preconditioner for Elliptic Equations with Jump Coefficients

This paper provides a proof of robustness of the restricted additive Schwarz preconditioner with harmonic overlap (RASHO) for the second order elliptic problems with jump coefficients. By analyzing the eigenvalue distribution of the RASHO preconditioner, we prove that the convergence rate of preconditioned conjugate gradient method with RASHO preconditioner is uniform with respect to the large jump and mesh size.     

Millimeter Wave Scattering of 2-D Frequency Selective Surface by Efficient Method of Lines with Preconditioned CG Technique

In this paper, both banded and symmetric successive overrelaxation (SSOR) preconditioned conjugate gradient (PCG) techniques are combined with method of lines (MOL) to further enhance the computational efficiency of this semi-analytic method. The electromagnetic wave scattering of 2-D frequency-selective surface is used as the examples to describe its implementation, whose analysis usually needs fast algorithms because of electrically large dimension. For arbitrary incident wave, helmholtz equation and boundary condition are used to calculate the impedance matrix and then to obtain reduced current-voltage linear matrix equation in spatial domain. Both banded and effective symmetric successive overrelaxation preconditioned conjugate gradient iterative method are chosen to solve this matrix equation. Our numerical results show that PCG methods can converge to accurate solution in much fewer iteration steps for analysis of the electromagnetic wave scattering from 2-D frequency-selective surface.     

Analysis of Millimeter Wave Scattering by the Infinite Plane Metallic Grating Using PCG-FFT Technique

In this paper, both fast Fourier transformation (FFT) and preconditioned CG technique are introduced into method of lines (MOL) to further enhance the computational efficiency of this semi-analytic method. Electromagnetic wave scattering by an infinite plane metallic grating is used as the examples to describe its implementation. For arbitrary incident wave, Helmholz equation and boundary condition are first transformed into new ones so that the impedance matrix elements are calculated by FFT technique. As a result, this Topelitz impedance matrix only requires O(N) memory storage for the conjugate gradient FFT method to solve the current distribution with the computational complexity O(N log N) . Our numerical results show that circulate matrix preconditioner can speed up CG-FFT method to converge in much smaller CPU time than the banded matrix preconditioner.     

改进的约束共轭梯度闪光照相图像重建算法

针对闪光照相图像信噪比低的特点,提出了一种改进的约束共轭梯度闪光照相图像重建算法。该算法在约束共轭梯度迭代重建的基础上,提出了新的预优矩阵选取方案,减小了重建图像的轴线噪声,利用松弛迭代步长代替共轭梯度法中的最优迭代步长,并采用了新的收敛准则,在保证算法收敛的同时,减少了重建算法的计算量。数值试验表明,与传统约束共轭梯度重建算法相比,改进算法稳定收敛,迭代速度更快,并能有效提高重建质量。     

Application of SSOR Preconditioning Technique to Method of Lines for Millimeter Wave Scattering

In this paper, symmetric successive overrelaxation (SSOR) preconditioned CG technique are introduced into method of lines (MOL) to further enhance the computational efficiency of this semi-analytic method. Millimeter wave scattering by an infinite plane metallic grating is used as the examples to describe its implementation, whose analysis usually needs fast algorithms because of electrically large dimension. For arbitrary incident wave, Helmholz equation and boundary condition are used to calculate the impedance matrix and then to obtain reduced current-voltage linear matrix equation in spatial domain. An effective symmetric successive overrelaxation preconditioned conjugate gradient iterative method, SSOR-PCG, is chosen to solve this matrix equation. With SSOR as the preconditioner as well as its efficient implementing in CG algorithm, PCG method can converge to accurate solution in much fewer iteration steps.     

改进的约束共轭梯度闪光照相图像重建算法

 针对闪光照相图像信噪比低的特点,提出了一种改进的约束共轭梯度闪光照相图像重建算法。该算法在约束共轭梯度迭代重建的基础上,提出了新的预优矩阵选取方案,减小了重建图像的轴线噪声,利用松弛迭代步长代替共轭梯度法中的最优迭代步长,并采用了新的收敛准则,在保证算法收敛的同时,减少了重建算法的计算量。数值试验表明,与传统约束共轭梯度重建算法相比,改进算法稳定收敛,迭代速度更快,并能有效提高重建质量。     

采用高阶近似逆矩阵的块PCG法

本文应用矩阵元素阶、阶矩阵及消去法的影响域等概念,给出了强主元多对角阵高阶近似求逆的一种快速算法。在强主元条件下,该法可应用于非对称阵和非正定阵。本文将该法与块预处理共轭梯度法相结合,应用于椭圆型方程数值解及类似问题的计算。数值结果表明,该法不仅适用范围较广,也具有较高的计算效率。     

Multigrid algorithms for high-order discontinuous Galerkin discretizations of the compressible Navier–Stokes equations

Multigrid algorithms are developed for systems arising from high-order discontinuous Galerkin discretizations of the compressible Navier–Stokes equations on unstructured meshes. The algorithms are based on coupling both p- and h-multigrid (ph-multigrid) methods which are used in nonlinear or linear forms, and either directly as solvers or as preconditioners to a Newton–Krylov method.The performance of the algorithms are examined in solving the laminar flow over an airfoil configuration. It is shown that the choice of the cycling strategy is crucial in achieving efficient and scalable solvers. For the multigrid solvers, while the order-independent convergence rate is obtained with a proper cycle type, the mesh-independent performance is achieved only if the coarsest problem is solved to a sufficient accuracy. On the other hand, the multigrid preconditioned Newton–GMRES solver appears to be insensitive to this condition and mesh-independent convergence is achieved under the desirable condition that the coarsest problem is solved using a fixed number of multigrid cycles regardless of the size of the problem.It is concluded that the Newton–GMRES solver with the multigrid preconditioning yields the most efficient and robust algorithm among those studied.     

Correction to: Generation of amplitude- and phase-modulated signal beam with a phase-only spatial light modulator and its detection by the transport of intensity equation method for holographic data storage

This study evaluates a novel holographic data storage (HDS) that uses a phase-only spatial light modulator (SLM) for the multilevel complex amplitude modulation of a signal beam and the transport of intensity equation (TIE) method to detect the signal beam without interferometry, to increase the capacity of the HDS, simplifying its optical system, and improving the stability of the signal beam modulation and detection. Both the amplitude and phase of the signal beam were modulated by a computer-generated hologram displayed in a phase-only SLM, a 4-f optical system, and a pinhole placed in the Fourier plane. The complex amplitude-modulated signal beam generated by this scheme does not always perfectly match the target complex amplitude, and deviations from the amplitude and phase of the target complex amplitude may exist. It is unclear whether the TIE method, which is sensitive to the state of the beam intensity and the phase distributions to be detected (such as zero-intensity points and phase discontinuities), can accurately detect a signal beam whose complex amplitude is modulated by the modulation scheme with a phase-only SLM. Here, we demonstrate via numerical simulations and experiments that several methods of complex amplitude generation using a phase-only SLM can achieve multilevel modulation of the amplitude and phase of a signal beam and are suitable for detection by the TIE method in HDS.

     

非对称线性代数方程组的并行算法

对五、七和九对角矩阵在并行-向量处理机上建立了解非对称线性代数方程组的PGCR算法,并分析了它的收敛性。本文给出的算法与多处理机的PGCR算法具有相同的选代次数。利用向量机模拟并行向量计算机并做数值实验,数值结果显示此算法很有效。     

A new symmetrized solution for phase retrieval using the transport of intensity equation

We propose a novel symmetrization method for solving the transport of intensity equation (TIE) using fast Fourier transforms for situations where the input images may or may not exhibit spatial periodicity. The method is derived from the analysis of intensity conservation law and the internal symmetry of the TIE, and is illustrated for both a computational and an experimental data set.     

三维预处理技术初探

本文应用高阶近似LU分解法对三维椭圆型方程的正规七点差分矩阵给出了预处理阵,并对典型问题给出了预处理共轭梯度法的计算结果。文中所采用的分阶方论以及阶矩阵等概念,适用于一般的非奇异对角优势稀疏矩阵。对不同网格点数的计算结果表明,三维预处理共轭梯度法仍具有超线性的收敛速率,在高阶情形下超线性的特征尤为突出,由于非零对角线数增长过快,高阶的三维预处理方法不宜采用,零阶方法和一阶方法是值得推荐的。     

A dispersion minimizing finite difference scheme and preconditioned solver for the 3D Helmholtz equation

In this paper, a new 27-point finite difference method is presented for solving the 3D Helmholtz equation with perfectly matched layer (PML), which is a second order scheme and pointwise consistent with the equation. An error analysis is made between the numerical wavenumber and the exact wavenumber, and a refined choice strategy based on minimizing the numerical dispersion is proposed for choosing weight parameters. A full-coarsening multigrid-based preconditioned Bi-CGSTAB method is developed for solving the linear system stemming from the Helmholtz equation with PML by the finite difference scheme. The shifted-Laplacian is extended to precondition the 3D Helmholtz equation, and a spectral analysis is given. The discrete preconditioned system is solved by the Bi-CGSTAB method, with a multigrid method used to invert the preconditioner approximately. Full-coarsening multigrid is employed, and a new matrix-based prolongation operator is constructed accordingly. Numerical results are presented to demonstrate the efficiency of both the new 27-point finite difference scheme with refined parameters, and the preconditioned Bi-CGSTAB method with the 3D full-coarsening multigrid.