Test for the new GridGenerator::laplace_transformation function.

git-svn-id: https://svn.dealii.org/trunk@4676 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/tests/deal.II/grid_transform.cc b/tests/deal.II/grid_transform.cc
new file mode 100644 (file)
index 0000000..aa4a396
--- /dev/null
+++ b/tests/deal.II/grid_transform.cc
@@ -0,0 +1,118 @@
+//----------------------------  grid_transform.cc  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 1998, 1999, 2000, 2001 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  grid_transform.cc  ---------------------------
+
+
+
+#include <base/logstream.h>
+#include <grid/tria.h>
+#include <grid/tria_accessor.h>
+#include <grid/tria_boundary_lib.h>
+#include <grid/grid_generator.h>
+#include <grid/grid_out.h>
+#include <fe/mapping_q.h>
+
+#include <fstream>
+
+
+int main ()
+{
+  const unsigned int dim=2;
+  Point<dim> origin;
+  HyperShellBoundary<dim> boundary(origin);
+  MappingQ<dim> mapping(2);
+  Triangulation<dim> tria;
+  tria.set_boundary(0, boundary);
+  const double inner_radius=1.;
+  const double outer_radius=5.;
+  GridGenerator::hyper_shell(tria, origin, inner_radius, outer_radius, 8);
+  tria.refine_global(2);
+  
+  GridOut grid_out;
+  if (false)
+    {
+      ofstream eps_stream1("grid_untransformed.output");
+      grid_out.write_eps(tria, eps_stream1, &mapping);
+    }
+  
+                                  // build up a map of vertex indices
+                                  // of boundary vertices to the new
+                                  // boundary points
+  std::map<unsigned int,Point<dim> > new_points;
+  
+                                  // new center and new radius
+                                  // of the inner circle.
+  const Point<dim> n_center(0,-1);
+  const double n_radius=0.5;
+
+  Triangulation<dim>::cell_iterator cell=tria.begin_active(),
+                                            endc=tria.end();
+  Triangulation<dim>::face_iterator face;
+  for (; cell!=endc; ++cell)
+    {
+      if (cell->at_boundary())
+       for (unsigned int face_no=0; face_no<GeometryInfo<dim>::faces_per_cell; ++face_no)
+         {
+           face=cell->face(face_no);
+           if (face->at_boundary())
+             for (unsigned int vertex_no=0;
+                  vertex_no<GeometryInfo<dim>::vertices_per_face; ++vertex_no)
+               {
+                 const Point<dim> &v=face->vertex(vertex_no);
+                 if (fabs(sqrt(v.square())-outer_radius)<1e-15)
+                   {
+                                                      // leave the
+                                                      // point, where
+                                                      // they are.
+                     new_points.insert(pair<unsigned int, Point<dim> > (
+                       face->vertex_index(vertex_no), v));
+                   }
+                 else if (fabs(sqrt(v.square())-inner_radius)<1e-15)
+                   {
+                                                      // move the
+                                                      // center of
+                                                      // the inner
+                                                      // circle to
+                                                      // (-1,0) and
+                                                      // take half
+                                                      // the radius
+                                                      // of the
+                                                      // circle.
+                     new_points.insert(pair<unsigned int, Point<dim> > (
+                       face->vertex_index(vertex_no), n_radius/inner_radius*v+n_center));
+                     face->set_boundary_indicator(1);
+                   }
+                 else
+                   Assert(false, ExcInternalError());
+               }  
+         }
+    }
+
+                                  // test output
+  if (false)
+    {
+      cout << "New points:" << endl;
+      std::map<unsigned int,Point<dim> >::const_iterator map_iter;
+      for (map_iter=new_points.begin(); map_iter!=new_points.end(); ++map_iter)
+       cout << "vertex " << map_iter->first << ": " << map_iter->second << endl;
+    }
+
+  GridGenerator::laplace_transformation (tria, new_points);
+  HyperBallBoundary<dim> inner_ball(n_center, n_radius);
+  tria.set_boundary(1, inner_ball);
+  ofstream eps_stream2("grid_transform.output");
+  grid_out.write_eps(tria, eps_stream2, &mapping);
+
+  tria.clear();
+  
+  return 0;
+};

diff --git a/tests/deal.II/grid_transform.checked b/tests/deal.II/grid_transform.checked
new file mode 100644 (file)
index 0000000..3922ee5
--- /dev/null
+++ b/tests/deal.II/grid_transform.checked
@@ -0,0 +1,429 @@
+%!PS-Adobe-2.0 EPSF-1.2
+%%Title: deal.II Output
+%%Creator: the deal.II library
+
+%%BoundingBox: 0 0 301 301
+/m {moveto} bind def
+/x {lineto stroke} bind def
+/b {0 0 0 setrgbcolor} def
+/r {1 0 0 setrgbcolor} def
+%%EndProlog
+
+0.5 setlinewidth
+b 288.582 207.403 m 257.034 189.713 x
+b 257.034 189.713 m 224.847 170.233 x
+b 236.58 140.342 m 230.509 154.954 x
+b 230.509 154.954 m 224.847 170.233 x
+b 300 150 m 268.549 145.757 x
+b 268.549 145.757 m 236.58 140.342 x
+b 297.118 179.264 m 263.011 167.491 x
+b 263.011 167.491 m 230.509 154.954 x
+b 268.549 145.757 m 263.011 167.491 x
+b 263.011 167.491 m 257.034 189.713 x
+b 256.066 256.066 m 233.827 229.583 x
+b 233.827 229.583 m 211.221 201.563 x
+b 224.847 170.233 m 217.904 185.384 x
+b 217.904 185.384 m 211.221 201.563 x
+b 274.72 233.336 m 245.625 209.463 x
+b 245.625 209.463 m 217.904 185.384 x
+b 257.034 189.713 m 245.625 209.463 x
+b 245.625 209.463 m 233.827 229.583 x
+b 211.221 201.563 m 187.758 170.568 x
+b 187.758 170.568 m 160.607 130.607 x
+b 224.847 170.233 m 196.457 151.107 x
+b 196.457 151.107 m 163.858 125.74 x
+b 217.904 185.384 m 191.922 160.203 x
+b 191.922 160.203 m 162.472 128.334 x
+b 196.457 151.107 m 191.922 160.203 x
+b 191.922 160.203 m 187.758 170.568 x
+b 236.58 140.342 m 203.398 132.809 x
+b 203.398 132.809 m 165 120 x
+b 230.509 154.954 m 199.624 141.608 x
+b 199.624 141.608 m 164.712 122.926 x
+b 203.398 132.809 m 199.624 141.608 x
+b 199.624 141.608 m 196.457 151.107 x
+b 207.403 288.582 m 194.335 252.412 x
+b 194.335 252.412 m 181.003 214.078 x
+b 211.221 201.563 m 195.952 207.336 x
+b 195.952 207.336 m 181.003 214.078 x
+b 233.336 274.72 m 214.196 240.892 x
+b 214.196 240.892 m 195.952 207.336 x
+b 233.827 229.583 m 214.196 240.892 x
+b 214.196 240.892 m 194.335 252.412 x
+b 150 300 m 150 264.305 x
+b 150 264.305 m 150 226.922 x
+b 181.003 214.078 m 165.522 219.94 x
+b 165.522 219.94 m 150 226.922 x
+b 179.264 297.118 m 172.224 258.278 x
+b 172.224 258.278 m 165.522 219.94 x
+b 194.335 252.412 m 172.224 258.278 x
+b 172.224 258.278 m 150 264.305 x
+b 150 226.922 m 150 186.208 x
+b 150 186.208 m 150 135 x
+b 181.003 214.078 m 169.243 178.321 x
+b 169.243 178.321 m 155.74 133.858 x
+b 165.522 219.94 m 159.662 181.57 x
+b 159.662 181.57 m 152.926 134.712 x
+b 169.243 178.321 m 159.662 181.57 x
+b 159.662 181.57 m 150 186.208 x
+b 195.952 207.336 m 178.257 173.867 x
+b 178.257 173.867 m 158.334 132.472 x
+b 187.758 170.568 m 178.257 173.867 x
+b 178.257 173.867 m 169.243 178.321 x
+b 92.5975 288.582 m 105.665 252.412 x
+b 105.665 252.412 m 118.997 214.078 x
+b 150 226.922 m 134.478 219.94 x
+b 134.478 219.94 m 118.997 214.078 x
+b 120.736 297.118 m 127.776 258.278 x
+b 127.776 258.278 m 134.478 219.94 x
+b 150 264.305 m 127.776 258.278 x
+b 127.776 258.278 m 105.665 252.412 x
+b 43.934 256.066 m 66.1734 229.583 x
+b 66.1734 229.583 m 88.7789 201.563 x
+b 118.997 214.078 m 104.048 207.336 x
+b 104.048 207.336 m 88.7789 201.563 x
+b 66.6645 274.72 m 85.8041 240.892 x
+b 85.8041 240.892 m 104.048 207.336 x
+b 105.665 252.412 m 85.8041 240.892 x
+b 85.8041 240.892 m 66.1734 229.583 x
+b 88.7789 201.563 m 112.242 170.568 x
+b 112.242 170.568 m 139.393 130.607 x
+b 118.997 214.078 m 130.757 178.321 x
+b 130.757 178.321 m 144.26 133.858 x
+b 104.048 207.336 m 121.743 173.867 x
+b 121.743 173.867 m 141.666 132.472 x
+b 130.757 178.321 m 121.743 173.867 x
+b 121.743 173.867 m 112.242 170.568 x
+b 134.478 219.94 m 140.338 181.57 x
+b 140.338 181.57 m 147.074 134.712 x
+b 150 186.208 m 140.338 181.57 x
+b 140.338 181.57 m 130.757 178.321 x
+b 11.4181 207.403 m 42.9656 189.713 x
+b 42.9656 189.713 m 75.1528 170.233 x
+b 88.7789 201.563 m 82.0957 185.384 x
+b 82.0957 185.384 m 75.1528 170.233 x
+b 25.2796 233.336 m 54.3747 209.463 x
+b 54.3747 209.463 m 82.0957 185.384 x
+b 66.1734 229.583 m 54.3747 209.463 x
+b 54.3747 209.463 m 42.9656 189.713 x
+b 0 150 m 31.4513 145.757 x
+b 31.4513 145.757 m 63.4202 140.342 x
+b 75.1528 170.233 m 69.4912 154.954 x
+b 69.4912 154.954 m 63.4202 140.342 x
+b 2.88221 179.264 m 36.9894 167.491 x
+b 36.9894 167.491 m 69.4912 154.954 x
+b 42.9656 189.713 m 36.9894 167.491 x
+b 36.9894 167.491 m 31.4513 145.757 x
+b 63.4202 140.342 m 96.6018 132.809 x
+b 96.6018 132.809 m 135 120 x
+b 75.1528 170.233 m 103.543 151.107 x
+b 103.543 151.107 m 136.142 125.74 x
+b 69.4912 154.954 m 100.376 141.608 x
+b 100.376 141.608 m 135.288 122.926 x
+b 103.543 151.107 m 100.376 141.608 x
+b 100.376 141.608 m 96.6018 132.809 x
+b 82.0957 185.384 m 108.078 160.203 x
+b 108.078 160.203 m 137.528 128.334 x
+b 112.242 170.568 m 108.078 160.203 x
+b 108.078 160.203 m 103.543 151.107 x
+b 11.4181 92.5975 m 42.9656 101.043 x
+b 42.9656 101.043 m 75.1528 108.228 x
+b 63.4202 140.342 m 69.4912 123.909 x
+b 69.4912 123.909 m 75.1528 108.228 x
+b 2.88221 120.736 m 36.9894 123.043 x
+b 36.9894 123.043 m 69.4912 123.909 x
+b 31.4513 145.757 m 36.9894 123.043 x
+b 36.9894 123.043 m 42.9656 101.043 x
+b 43.934 43.934 m 66.1734 61.93 x
+b 66.1734 61.93 m 88.7789 79.1211 x
+b 75.1528 108.228 m 82.0957 93.4791 x
+b 82.0957 93.4791 m 88.7789 79.1211 x
+b 25.2796 66.6645 m 54.3747 81.0713 x
+b 54.3747 81.0713 m 82.0957 93.4791 x
+b 42.9656 101.043 m 54.3747 81.0713 x
+b 54.3747 81.0713 m 66.1734 61.93 x
+b 88.7789 79.1211 m 112.242 95.0513 x
+b 112.242 95.0513 m 139.393 109.393 x
+b 75.1528 108.228 m 103.543 112.621 x
+b 103.543 112.621 m 136.142 114.26 x
+b 82.0957 93.4791 m 108.078 103.689 x
+b 108.078 103.689 m 137.528 111.666 x
+b 103.543 112.621 m 108.078 103.689 x
+b 108.078 103.689 m 112.242 95.0513 x
+b 69.4912 123.909 m 100.376 122.283 x
+b 100.376 122.283 m 135.288 117.074 x
+b 96.6018 132.809 m 100.376 122.283 x
+b 100.376 122.283 m 103.543 112.621 x
+b 92.5975 11.4181 m 105.665 38.3432 x
+b 105.665 38.3432 m 118.997 64.3834 x
+b 88.7789 79.1211 m 104.048 71.5271 x
+b 104.048 71.5271 m 118.997 64.3834 x
+b 66.6645 25.2796 m 85.8041 49.6419 x
+b 85.8041 49.6419 m 104.048 71.5271 x
+b 66.1734 61.93 m 85.8041 49.6419 x
+b 85.8041 49.6419 m 105.665 38.3432 x
+b 150 0 m 150 27.2078 x
+b 150 27.2078 m 150 53.7625 x
+b 118.997 64.3834 m 134.478 58.9226 x
+b 134.478 58.9226 m 150 53.7625 x
+b 120.736 2.88221 m 127.776 32.2565 x
+b 127.776 32.2565 m 134.478 58.9226 x
+b 105.665 38.3432 m 127.776 32.2565 x
+b 127.776 32.2565 m 150 27.2078 x
+b 150 53.7625 m 150 79.4113 x
+b 150 79.4113 m 150 105 x
+b 118.997 64.3834 m 130.757 85.4069 x
+b 130.757 85.4069 m 144.26 106.142 x
+b 134.478 58.9226 m 140.338 82.3218 x
+b 140.338 82.3218 m 147.074 105.288 x
+b 130.757 85.4069 m 140.338 82.3218 x
+b 140.338 82.3218 m 150 79.4113 x
+b 104.048 71.5271 m 121.743 90.0239 x
+b 121.743 90.0239 m 141.666 107.528 x
+b 112.242 95.0513 m 121.743 90.0239 x
+b 121.743 90.0239 m 130.757 85.4069 x
+b 207.403 11.4181 m 194.335 38.3432 x
+b 194.335 38.3432 m 181.003 64.3834 x
+b 150 53.7625 m 165.522 58.9226 x
+b 165.522 58.9226 m 181.003 64.3834 x
+b 179.264 2.88221 m 172.224 32.2565 x
+b 172.224 32.2565 m 165.522 58.9226 x
+b 150 27.2078 m 172.224 32.2565 x
+b 172.224 32.2565 m 194.335 38.3432 x
+b 256.066 43.934 m 233.827 61.93 x
+b 233.827 61.93 m 211.221 79.1211 x
+b 181.003 64.3834 m 195.952 71.5271 x
+b 195.952 71.5271 m 211.221 79.1211 x
+b 233.336 25.2796 m 214.196 49.6419 x
+b 214.196 49.6419 m 195.952 71.5271 x
+b 194.335 38.3432 m 214.196 49.6419 x
+b 214.196 49.6419 m 233.827 61.93 x
+b 211.221 79.1211 m 187.758 95.0513 x
+b 187.758 95.0513 m 160.607 109.393 x
+b 181.003 64.3834 m 169.243 85.4069 x
+b 169.243 85.4069 m 155.74 106.142 x
+b 195.952 71.5271 m 178.257 90.0239 x
+b 178.257 90.0239 m 158.334 107.528 x
+b 169.243 85.4069 m 178.257 90.0239 x
+b 178.257 90.0239 m 187.758 95.0513 x
+b 165.522 58.9226 m 159.662 82.3218 x
+b 159.662 82.3218 m 152.926 105.288 x
+b 150 79.4113 m 159.662 82.3218 x
+b 159.662 82.3218 m 169.243 85.4069 x
+b 288.582 92.5975 m 257.034 101.043 x
+b 257.034 101.043 m 224.847 108.228 x
+b 211.221 79.1211 m 217.904 93.4791 x
+b 217.904 93.4791 m 224.847 108.228 x
+b 274.72 66.6645 m 245.625 81.0713 x
+b 245.625 81.0713 m 217.904 93.4791 x
+b 233.827 61.93 m 245.625 81.0713 x
+b 245.625 81.0713 m 257.034 101.043 x
+b 224.847 108.228 m 230.509 123.909 x
+b 230.509 123.909 m 236.58 140.342 x
+b 297.118 120.736 m 263.011 123.043 x
+b 263.011 123.043 m 230.509 123.909 x
+b 257.034 101.043 m 263.011 123.043 x
+b 263.011 123.043 m 268.549 145.757 x
+b 224.847 108.228 m 196.457 112.621 x
+b 196.457 112.621 m 163.858 114.26 x
+b 230.509 123.909 m 199.624 122.283 x
+b 199.624 122.283 m 164.712 117.074 x
+b 196.457 112.621 m 199.624 122.283 x
+b 199.624 122.283 m 203.398 132.809 x
+b 217.904 93.4791 m 191.922 103.689 x
+b 191.922 103.689 m 162.472 111.666 x
+b 187.758 95.0513 m 191.922 103.689 x
+b 191.922 103.689 m 196.457 112.621 x
+b 300 150 m 299.678 159.817 x
+b 299.678 159.817 m 298.717 169.572 x
+b 298.717 169.572 m 297.118 179.264 x
+b 297.118 179.264 m 294.887 188.83 x
+b 294.887 188.83 m 292.042 198.209 x
+b 292.042 198.209 m 288.582 207.403 x
+b 288.582 207.403 m 284.528 216.35 x
+b 284.528 216.35 m 279.907 224.994 x
+b 279.907 224.994 m 274.72 233.336 x
+b 274.72 233.336 m 268.999 241.32 x
+b 268.999 241.32 m 262.78 248.896 x
+b 262.78 248.896 m 256.066 256.066 x
+b 162.472 128.334 m 161.9 129.132 x
+b 161.9 129.132 m 161.278 129.89 x
+b 161.278 129.89 m 160.607 130.607 x
+b 163.858 125.74 m 163.453 126.635 x
+b 163.453 126.635 m 162.991 127.499 x
+b 162.991 127.499 m 162.472 128.334 x
+b 164.712 122.926 m 164.489 123.883 x
+b 164.489 123.883 m 164.204 124.821 x
+b 164.204 124.821 m 163.858 125.74 x
+b 165 120 m 164.968 120.982 x
+b 164.968 120.982 m 164.872 121.957 x
+b 164.872 121.957 m 164.712 122.926 x
+b 256.066 256.066 m 248.896 262.78 x
+b 248.896 262.78 m 241.32 268.999 x
+b 241.32 268.999 m 233.336 274.72 x
+b 233.336 274.72 m 224.994 279.907 x
+b 224.994 279.907 m 216.35 284.528 x
+b 216.35 284.528 m 207.403 288.582 x
+b 207.403 288.582 m 198.209 292.042 x
+b 198.209 292.042 m 188.83 294.887 x
+b 188.83 294.887 m 179.264 297.118 x
+b 179.264 297.118 m 169.572 298.717 x
+b 169.572 298.717 m 159.817 299.678 x
+b 159.817 299.678 m 150 300 x
+b 152.926 134.712 m 151.957 134.872 x
+b 151.957 134.872 m 150.982 134.968 x
+b 150.982 134.968 m 150 135 x
+b 155.74 133.858 m 154.821 134.204 x
+b 154.821 134.204 m 153.883 134.489 x
+b 153.883 134.489 m 152.926 134.712 x
+b 158.334 132.472 m 157.499 132.991 x
+b 157.499 132.991 m 156.635 133.453 x
+b 156.635 133.453 m 155.74 133.858 x
+b 160.607 130.607 m 159.89 131.278 x
+b 159.89 131.278 m 159.132 131.9 x
+b 159.132 131.9 m 158.334 132.472 x
+b 150 300 m 140.183 299.678 x
+b 140.183 299.678 m 130.428 298.717 x
+b 130.428 298.717 m 120.736 297.118 x
+b 120.736 297.118 m 111.17 294.887 x
+b 111.17 294.887 m 101.791 292.042 x
+b 101.791 292.042 m 92.5975 288.582 x
+b 92.5975 288.582 m 83.6505 284.528 x
+b 83.6505 284.528 m 75.0062 279.907 x
+b 75.0062 279.907 m 66.6645 274.72 x
+b 66.6645 274.72 m 58.6803 268.999 x
+b 58.6803 268.999 m 51.1035 262.78 x
+b 51.1035 262.78 m 43.934 256.066 x
+b 141.666 132.472 m 140.868 131.9 x
+b 140.868 131.9 m 140.11 131.278 x
+b 140.11 131.278 m 139.393 130.607 x
+b 144.26 133.858 m 143.365 133.453 x
+b 143.365 133.453 m 142.501 132.991 x
+b 142.501 132.991 m 141.666 132.472 x
+b 147.074 134.712 m 146.117 134.489 x
+b 146.117 134.489 m 145.179 134.204 x
+b 145.179 134.204 m 144.26 133.858 x
+b 150 135 m 149.018 134.968 x
+b 149.018 134.968 m 148.043 134.872 x
+b 148.043 134.872 m 147.074 134.712 x
+b 43.934 256.066 m 37.2195 248.896 x
+b 37.2195 248.896 m 31.0014 241.32 x
+b 31.0014 241.32 m 25.2796 233.336 x
+b 25.2796 233.336 m 20.0928 224.994 x
+b 20.0928 224.994 m 15.4723 216.35 x
+b 15.4723 216.35 m 11.4181 207.403 x
+b 11.4181 207.403 m 7.95839 198.209 x
+b 7.95839 198.209 m 5.1131 188.83 x
+b 5.1131 188.83 m 2.88221 179.264 x
+b 2.88221 179.264 m 1.28253 169.572 x
+b 1.28253 169.572 m 0.321791 159.817 x
+b 0.321791 159.817 m 0 150 x
+b 135.288 122.926 m 135.128 121.957 x
+b 135.128 121.957 m 135.032 120.982 x
+b 135.032 120.982 m 135 120 x
+b 136.142 125.74 m 135.796 124.821 x
+b 135.796 124.821 m 135.511 123.883 x
+b 135.511 123.883 m 135.288 122.926 x
+b 137.528 128.334 m 137.009 127.499 x
+b 137.009 127.499 m 136.547 126.635 x
+b 136.547 126.635 m 136.142 125.74 x
+b 139.393 130.607 m 138.722 129.89 x
+b 138.722 129.89 m 138.1 129.132 x
+b 138.1 129.132 m 137.528 128.334 x
+b 0 150 m 0.321791 140.183 x
+b 0.321791 140.183 m 1.28253 130.428 x
+b 1.28253 130.428 m 2.88221 120.736 x
+b 2.88221 120.736 m 5.1131 111.17 x
+b 5.1131 111.17 m 7.95839 101.791 x
+b 7.95839 101.791 m 11.4181 92.5975 x
+b 11.4181 92.5975 m 15.4723 83.6505 x
+b 15.4723 83.6505 m 20.0928 75.0062 x
+b 20.0928 75.0062 m 25.2796 66.6645 x
+b 25.2796 66.6645 m 31.0014 58.6803 x
+b 31.0014 58.6803 m 37.2195 51.1035 x
+b 37.2195 51.1035 m 43.934 43.934 x
+b 137.528 111.666 m 138.1 110.868 x
+b 138.1 110.868 m 138.722 110.11 x
+b 138.722 110.11 m 139.393 109.393 x
+b 136.142 114.26 m 136.547 113.365 x
+b 136.547 113.365 m 137.009 112.501 x
+b 137.009 112.501 m 137.528 111.666 x
+b 135.288 117.074 m 135.511 116.117 x
+b 135.511 116.117 m 135.796 115.179 x
+b 135.796 115.179 m 136.142 114.26 x
+b 135 120 m 135.032 119.018 x
+b 135.032 119.018 m 135.128 118.043 x
+b 135.128 118.043 m 135.288 117.074 x
+b 43.934 43.934 m 51.1035 37.2195 x
+b 51.1035 37.2195 m 58.6803 31.0014 x
+b 58.6803 31.0014 m 66.6645 25.2796 x
+b 66.6645 25.2796 m 75.0062 20.0928 x
+b 75.0062 20.0928 m 83.6505 15.4723 x
+b 83.6505 15.4723 m 92.5975 11.4181 x
+b 92.5975 11.4181 m 101.791 7.95839 x
+b 101.791 7.95839 m 111.17 5.1131 x
+b 111.17 5.1131 m 120.736 2.88221 x
+b 120.736 2.88221 m 130.428 1.28253 x
+b 130.428 1.28253 m 140.183 0.321791 x
+b 140.183 0.321791 m 150 0 x
+b 147.074 105.288 m 148.043 105.128 x
+b 148.043 105.128 m 149.018 105.032 x
+b 149.018 105.032 m 150 105 x
+b 144.26 106.142 m 145.179 105.796 x
+b 145.179 105.796 m 146.117 105.511 x
+b 146.117 105.511 m 147.074 105.288 x
+b 141.666 107.528 m 142.501 107.009 x
+b 142.501 107.009 m 143.365 106.547 x
+b 143.365 106.547 m 144.26 106.142 x
+b 139.393 109.393 m 140.11 108.722 x
+b 140.11 108.722 m 140.868 108.1 x
+b 140.868 108.1 m 141.666 107.528 x
+b 150 0 m 159.817 0.321791 x
+b 159.817 0.321791 m 169.572 1.28253 x
+b 169.572 1.28253 m 179.264 2.88221 x
+b 179.264 2.88221 m 188.83 5.1131 x
+b 188.83 5.1131 m 198.209 7.95839 x
+b 198.209 7.95839 m 207.403 11.4181 x
+b 207.403 11.4181 m 216.35 15.4723 x
+b 216.35 15.4723 m 224.994 20.0928 x
+b 224.994 20.0928 m 233.336 25.2796 x
+b 233.336 25.2796 m 241.32 31.0014 x
+b 241.32 31.0014 m 248.896 37.2195 x
+b 248.896 37.2195 m 256.066 43.934 x
+b 158.334 107.528 m 159.132 108.1 x
+b 159.132 108.1 m 159.89 108.722 x
+b 159.89 108.722 m 160.607 109.393 x
+b 155.74 106.142 m 156.635 106.547 x
+b 156.635 106.547 m 157.499 107.009 x
+b 157.499 107.009 m 158.334 107.528 x
+b 152.926 105.288 m 153.883 105.511 x
+b 153.883 105.511 m 154.821 105.796 x
+b 154.821 105.796 m 155.74 106.142 x
+b 150 105 m 150.982 105.032 x
+b 150.982 105.032 m 151.957 105.128 x
+b 151.957 105.128 m 152.926 105.288 x
+b 256.066 43.934 m 262.78 51.1035 x
+b 262.78 51.1035 m 268.999 58.6803 x
+b 268.999 58.6803 m 274.72 66.6645 x
+b 274.72 66.6645 m 279.907 75.0062 x
+b 279.907 75.0062 m 284.528 83.6505 x
+b 284.528 83.6505 m 288.582 92.5975 x
+b 288.582 92.5975 m 292.042 101.791 x
+b 292.042 101.791 m 294.887 111.17 x
+b 294.887 111.17 m 297.118 120.736 x
+b 297.118 120.736 m 298.717 130.428 x
+b 298.717 130.428 m 299.678 140.183 x
+b 299.678 140.183 m 300 150 x
+b 164.712 117.074 m 164.872 118.043 x
+b 164.872 118.043 m 164.968 119.018 x
+b 164.968 119.018 m 165 120 x
+b 163.858 114.26 m 164.204 115.179 x
+b 164.204 115.179 m 164.489 116.117 x
+b 164.489 116.117 m 164.712 117.074 x
+b 162.472 111.666 m 162.991 112.501 x
+b 162.991 112.501 m 163.453 113.365 x
+b 163.453 113.365 m 163.858 114.26 x
+b 160.607 109.393 m 161.278 110.11 x
+b 161.278 110.11 m 161.9 110.868 x
+b 161.9 110.868 m 162.472 111.666 x
+showpage