In this article, hierarchical finite element method (FEM) based on curvilinear elements is used to study three-dimensional (3D) Figure 5 Measured results of IMD3 for the fixed uneven DPA and the proposed optimum uneven DPA. [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com] Figure 6 Measured results of ACLR for the fixed uneven DPA and the proposed optimum uneven DPA. [Color figure can be viewed in the online issue, which is available at wileyonlinelibrary.com]electromagnetic problems. The incomplete Cholesky preconditioned loose generalized minimal residual solver (LGMRES) based on decomposition algorithm (DA) is applied to solve the FEM equations. The efficiency of the proposed approach is studied on several numerical problems. Numerical results demonstrate that the DA-based LGMRES is especially effective when the hierarchical FEM is used to solve 3D electromagnetic problems.