maincogpolartest.html

有限元学习研究用源代码(老外的),供科研人员参考

<TITLE>COG 2.0: disc</TITLE>

<H1>Intersection with a disc</H1>

 <P>In this variant, instead of the <A HREF="maincogtest.html">3D
cube</A>, the intersection of the 
 <A HREF="cogtest.html">test cogeometry</A> with a 2D disc (shifted in
z-direction) is computed.

<PRE>
#include "<A HREF="coglib.hxx">cog/coglib.hxx</A>"

cogeometry test_cogeometry();

#include "<A HREF="cogtest.cxx">cogtest.cxx</A>"

int main()
{
  cogeometry geom = test_cogeometry();

  cogenOctree base = new <A HREF="cogenoctree.html">CogenOctree</A>();
  base-&gt;setBorder(0.0,1.0, -wzPi,wzPi);
  base-&gt;refinementGlobal(0.1,0.3);
</PRE>

 <P>Now, let's define the coordinates. Let's assume that we want to
compute a non-trivial z-intersection of the circle with the
geometry (z=0.3), and the cogeometry may be a non-trivial 3D
cogeometry:

<PRE>
  base-&gt;addPlaneZ(0.3);
</PRE>

 <P>This requires 3D coordinates!  But polar coordinates are only 2D
coordinates, thus, we have to construct a <B>3D continuation</B>:

<PRE>
  wzcoordinates polar = new wz3DContinuation(new wzPolarCoordinates);
  cogenChart gen = new CogenChart(base-&gt;generator(),polar);
  wzgrid grid = (*gen)(geom);
  grid->write("test.sg");
}
</PRE>