Developer Reference for Intel® oneAPI Math Kernel Library for Fortran

ID 766686
Date 11/07/2023
Public

A newer version of this document is available. Customers should click here to go to the newest version.

Document Table of Contents

Fast Poisson Solver Routines

In addition to the Real Discrete Trigonometric Transforms (TT) interface (refer to Trigonometric Transform Routines), Intel® oneAPI Math Kernel Library (oneMKL) supports thethe Poisson Solver interface. This interface implements a group of routines (Poisson Solver routines) used to compute a solution of Laplace, Poisson, and Helmholtz problems of a special kind using discrete Fourier transforms. Laplace and Poisson problems are special cases of a more general Helmholtz problem. The problems that are solved by the Poisson Solver interface are defined more exactly in Poisson Solver Implementation. The Poisson Solver interface provides much flexibility of use: you can call routines with the default parameter values or adjust routines to your particular needs by manually tuning routine parameters. You can adjust the style of error and warning messages to a Fortrannotation by setting up a dedicated parameter. This adds convenience to debugging, because you can read information in the way that is natural for your code. The Intel® oneAPI Math Kernel Library (oneMKL) Poisson Solver interface currently contains only routines that implement the following solvers:

  • Fast Laplace, Poisson and Helmholtz solvers in a Cartesian coordinate system
  • Fast Poisson and Helmholtz solvers in a spherical coordinate system.

To describe the Intel® oneAPI Math Kernel Library (oneMKL) Poisson Solver interface, the C convention is used. Fortran usage specifics can be found inCalling PDE Support Routines from Fortran.

NOTE:

Fortran users should keep in mind that array indices in Fortran start at 1 instead of 0, as they do in C.