Partial Differential Equation Toolbox |

poimesh

Make regular mesh on a rectangular geometry

Syntax

-     [p,e,t]=poimesh(g,nx,ny)
      [p,e,t]=poimesh(g,n)
      [p,e,t]=poimesh(g)

Description

`[p,e,t]=poimesh(g,nx,ny)` constructs a regular mesh on the rectangular geometry specified by g, by dividing the "x edge" into `nx` pieces and the "y edge" into `ny` pieces, and placing `(nx+1)*(ny+1)` points at the intersections.

The "x edge" is the one that makes the smallest angle with the x-axis.

`[p,e,t]=poimesh(g,n)` uses `nx=ny=n`, and `[p,e,t]=poimesh(g)` uses `nx=ny=1`.

The triangular mesh is described by the mesh data `p`, `e`, and `t`. For details on the mesh data representation, see `initmesh`.

For best performance with `poisolv`, the larger of `nx` and `ny` should be a power of 2.

If `g` does not seem to describe a rectangle, `p` is zero on return.

Examples

Try the demo command `pdedemo8`. The solution of Poisson's equation over a rectangular grid with boundary condition given by the file `squareb4` is returned. The solution time is compared to the usual Finite Element Method (FEM) approach.

See Also

`initmesh`, `poisolv`

| poiindex |  | poisolv |