A full multigrid method and fourth-order compact difference scheme is designed to solve the 3D Helmholtz equation on unequal mesh size. Three dimensional restriction and prolongation operators of the ...