Can Intel's Poisson solver routines be applied to the following general Poisson equation?

greiner08 — Fri, 19 Oct 2018 12:34:38 GMT

Hi,

I see in the documentation that Intels' Poisson solver routines can be used to solve equations of the form

∇.∇u = f

Can they also (perhaps via some additional steps) be applied to an equation of the form

∇.ε∇u = f

Thanks

The possible algorithm

mecej4 — Mon, 22 Oct 2018 00:20:52 GMT

The possible algorithm choices depend on whether ε is (i) a constant, (ii) position dependent, or (iii) a known function of u. In the first two cases, the finite-difference equations are still linear, and can be solved using Pardiso as in the case when ε = 1. In the third case, you may use Kirchoff's transformation, φ = \int ε du, which leads to ∇.∇φ = f; you solve for φ, and recover u from φ at the end.

topic The possible algorithm in Intel® oneAPI Math Kernel Library

Can Intel's Poisson solver routines be applied to the following general Poisson equation?

The possible algorithm