next up previous contents index
Next: Geometric variables, inner and Up: Border(), buildmesh(), polygon() Previous: Border(), buildmesh(), polygon()   Contents   Index

buildmesh and domain decomposition

There is a problem with buildmesh when doing domain decomposition: by default Gfem swap the diagonals at the corners of the domain if the triangle has two boundary edges. This will lead to bad domain decomposition at the sub-domain interfaces.

To solve this, there is a new flag for buildmesh which is optional:

buildmesh(<max_number_of_vertices>, <flag>)
where <flag> =

$\displaystyle \left\{
\begin{array}[l]{l}
= 0 \mbox{ classic way: do diagonal s...
...g}\\
= 1 \mbox{ domain decomposition: no diagonal swaping}
\end{array}\right.$



Christophe Prud'homme 2001-10-20