mainref2test.html

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

<TITLE>COG 2.0: refinement after curved coords</TITLE>

<PRE>
#include "<A HREF="coglib.hxx">cog/coglib.hxx</A>"
int main()
{
  cogenOctree base = new CogenOctree();
  wzIndex i1 = base-&gt;addBox(wzRegion(2), 0.000001,0.4,-0.3,1.5);
  wzIndex i2 = base-&gt;addBox(wzRegion(3), 0.2,0.3);
  wzIndex i3 = base-&gt;addBox(wzRegion(4), 0.2,0.3, 0.3,1.0);
  base-&gt;setFaceOfBox(wzFace(2),i1);
  base-&gt;setFaceOfBox(wzFace(3),i2);
  base-&gt;setFaceOfBox(wzFace(4),i3);
  base-&gt;refinementGlobal(0.05,0.3);
  base-&gt;addPlaneY(wzPi);
  base-&gt;addPlaneY(-wzPi);
  wzcoordinates polar = new wzPolarCoordinates;
  cogenChart chart = new CogenChart(base-&gt;generator(),polar);

  cogenRefine gen = new CogenRefine(chart-&gt;generator());
  gen-&gt;refinementGlobal(0.21,0.21,0.1);
  gen-&gt;faceRefinementGlobal(0.026,1.0);
  gen-&gt;refinementIsotropicalInRegion(0.02,wzRegion(4));
  gen-&gt;refinementIsotropicalNearPoint(0.02,0.04,0.6,0.8,0.0);
  gen-&gt;refinementNearPoint(0.02,0.1,0,  0.2,0.1,0,  0.6,0.4,0);

//  wzgrid grid = (*gen)();
  wzgrid grid = (*chart)();
//  ibgOutput(grid);
  grid->write("test.sg");
//  grid->write("test.sg",wzDefaultSegmentDescription);
}
</PRE>