Partial Differential Equation Toolbox | ![]() ![]() |
Fast solution of Poisson's equation on a rectangular grid
Syntax
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
Reference
Strang, Gilbert, Introduction to Applied Mathematics, Wellesley-Cambridge Press, Cambridge, MA, 1986, pp. 453-458.
![]() | poimesh | refinemesh | ![]() |