Partial Differential Equation Toolbox |

poisolv

Fast solution of Poisson's equation on a rectangular grid

Syntax

-     u=poisolv(b,p,e,t,f)

Description

`u=poisolv(b,p,e,t,f)` solves Poisson's equation with Dirichlet boundary conditions on a regular rectangular grid. A combination of sine transforms and tridiagonal solutions is used for increased performance.

The boundary conditions `b` must specify Dirichlet conditions for all boundary points.

The mesh `p`, `e`, and `t` must be a regular rectangular grid. For details on the mesh data representation, see `initmesh`.

`f` gives the right-hand side of Poisson's equation.

Apart from roundoff errors, the result should be the same as `u=assempde(b,p,e,t,1,0,f)`.

See Also

`poicalc`, `poimesh`

Reference

Strang, Gilbert, Introduction to Applied Mathematics, Wellesley-Cambridge Press, Cambridge, MA, 1986, pp. 453-458.

| poimesh |  | refinemesh |