本文提出了一种多重网格预处理的高阶浸入有限差分求解器，能够精确高效地求解复杂的二维和三维域上的Poisson方程。该求解器采用低阶Shortley-Weller多重网格方法对高阶无矩阵Krylov子空间求解器进行预处理。无矩阵方法能够完全兼容高阶浸入界面法（IIM）离散化的边界和界面条件，以及高阶小波适应的多分辨率网格。通过对二维域的验证和分析表明，该算法能够为具有Dirichlet、Neumann和/或界面跃变边界条件的Laplace和Poisson问题提供高阶精度的结果，并且通过多重网格方法有效预处理。进一步研究表明，所提出的算法能够在复杂三维域上高效求解Laplace和Poisson问题的高阶离散化，并可利用数千个计算核心在多分辨率网格上实现。据我们所知，这项工作解决了迄今为止应用高阶浸入方法解决椭圆偏微分方程的最大问题规模，并首次实现了三维多分辨率自适应网格上的高阶结果。综合来看，这项研究为利用高阶浸入方法求解具有边界或界面条件的各类三维偏微分方程铺平了道路，包括线性和非线性弹性问题、不可压缩Navier-Stokes方程以及流固耦合问题。

This work presents a multigrid preconditioned high order immersed finite difference solver to accurately and efficiently solve the Poisson equation on complex 2D and 3D domains. The solver employs a low order Shortley-Weller multigrid method to precondition a high order matrix-free Krylov subspace solver. The matrix-free approach enables full compatibility with high order IIM discretizations of boundary and interface conditions, as well as high order wavelet-adapted multiresolution grids. Through verification and analysis on 2D domains, we demonstrate the ability of the algorithm to provide high order accurate results to Laplace and Poisson problems with Dirichlet, Neumann, and/or interface jump boundary conditions, all effectively preconditioned using the multigrid method. We further show that the proposed method is able to efficiently solve high order discretizations of Laplace and Poisson problems on complex 3D domains using thousands of compute cores and on multiresolution grids. To our knowledge, this work presents the largest problem sizes tackled with high order immersed methods applied to elliptic partial differential equations, and the first high order results on 3D multiresolution adaptive grids. Together, this work paves the way for employing high order immersed methods to a variety of 3D partial differential equations with boundary or inter-face conditions, including linear and non-linear elasticity problems, the incompressible Navier-Stokes equations, and fluid-structure interactions.