mainfacetest.html

Delaunay三角形的网格剖分程序

<TITLE>COG 2.1: Boundary conditions</TITLE>

<A HREF="mainreftest.html">next</A>

<H1>COG 2.1: Setting boundary conditions</H1>

 <P>Regions have boundaries between them and outer boundaries.  For
these boundaries, different numbers are useful to distinguish
different boundary conditions.  That's why in COG every boundary face
has a face number.  This face number is part of the
 <A HREF="simplexgrid.html">simplex grid format</A> and may be used by
the application to define the boundary conditions.

 <P>Now, starting with our rectangular geometry, we want to define
different face numbers on different parts on the boundary, like here:

 <P><IMG SRC="facetest.gif">

 <P>Let's consider now how to define these <B>face numbers</B> using
the class <B><A HREF="cogenoctree.html">CogenOctree</A></B>.  We start
with slightly modifying the <A HREF="mainoctreetest.html"> previous
example code</A> - we remember the return value of <B>addBox</B>:

<PRE>
#include "<A HREF="coglib.hxx">cog/coglib.hxx</A>"
int main()
{
  cogenOctree gen = new CogenOctree();
  wzIndex i1 = gen-&gt;addBox(wzRegion(2), 0.0,0.4, 0.0,0.5);
  wzIndex i2 = gen-&gt;addBox(wzRegion(3), 0.2,0.8, 0.3,1.0);
  wzIndex i3 = gen-&gt;addBox(wzRegion(4), 0.4,1.0, 0.7,1.0);
</PRE>

 <P>This return value may be used to set a face number with
<B>setFaceOfBox</B> for the whole rectangle:

<PRE>
  gen-&gt;setFaceOfBox(wzFace(2),i1);
  gen-&gt;setFaceOfBox(wzFace(4),i3);
  gen-&gt;setFaceOfBox(wzFace(3),i2);
</PRE>

 <P>The boundary conditions may be defined for each of the four sides
separately with <B>setFaceOfBoxSide</B>.  The sides are named
<I>wzSideLeft, wzSideRight, wzSideUp, wzSideDown, wzSideNear,
wzSideFar</I>.

<PRE>
  gen-&gt;setFaceOfBoxSide(wzFace(7),i1,wzSideRight);
  gen-&gt;setFaceOfBoxSide(wzFace(6),i2,wzSideLeft);
</PRE>

 <P>The boundary face number may be re defined also inside a box with
<B>setFaceInBox</B> or, as a variant, on some part of a plane with
<B>setFaceOnPlaneX, setFaceOnPlaneY, setFaceOnPlaneZ</B>.

<PRE>
  gen-&gt;setFaceOnPlaneX(wzFace(5),0.2, 0.6,0.7);
  gen-&gt;setFaceOnPlaneY(wzFace(5),0.0, 0.6,0.7);
  gen-&gt;setFaceInBox(wzFace(8), 0.0,0.2, 0.8,1.0);
</PRE>

 <P>Note that the priority depends on the order - the latest call has
priority.  For example, the part of i3 inside i2 gets the face number
of the region i2 instead of i3.

 <P>If you want to write out only some of the faces (probably you
don't want to write out the default face 1) this is possible by
writing out
 <A HREF="mainmaterialtest.html">face type numbers</A> instead of
face numbers and using the "invalid" face type.

<H3>Usage of face numbers to protect edges</H3>

 <P>Note that two of the three regions have the same region number
3. This defines a non-convex region 3. Without any boundary
conditions, the inner edge may be cut off - the grid generator does
not know that this point is something special.  But now the edge is
the place where the boundary condition changes, thus, a special point
which has to be computed exactly.

 <P>In this situation, it is probably also useful to write out
 <A HREF="mainmaterialtest.html">face type numbers</A> instead of face
numbers.

 <HR>

<PRE>
  gen-&gt;refinementGlobal(0.1,0.05);
  wzgrid grid = (*gen)();
  grid->write("test.sg");
}
</PRE>