Adjust implementation of laplace_transform to what had been documented.

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

diff --git a/tests/deal.II/grid_transform_02.cc b/tests/deal.II/grid_transform_02.cc
new file mode 100644 (file)
index 0000000..9ceb05b
--- /dev/null
+++ b/tests/deal.II/grid_transform_02.cc
@@ -0,0 +1,89 @@
+// ---------------------------------------------------------------------
+// $Id$
+//
+// Copyright (C) 1998 - 2014 by the deal.II authors
+//
+// This file is part of the deal.II library.
+//
+// The deal.II library is free software; you can use it, redistribute
+// it, and/or modify it under the terms of the GNU Lesser General
+// Public License as published by the Free Software Foundation; either
+// version 2.1 of the License, or (at your option) any later version.
+// The full text of the license can be found in the file LICENSE at
+// the top level of the deal.II distribution.
+//
+// ---------------------------------------------------------------------
+
+
+// Verify that GridTools::laplace_transform can deal with interior
+// nodes being pinned to a new location as well. The test itself
+// doesn't make much sense since it leads to a few inverted cells, but
+// it allows for easy visual inspection that the desired result
+// happens.
+//
+// (Testcase adapted from one by Denis Davydov.)
+
+
+#include "../tests.h"
+#include <deal.II/base/logstream.h>
+#include <deal.II/grid/tria.h>
+#include <deal.II/grid/tria_accessor.h>
+#include <deal.II/grid/tria_boundary_lib.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/grid_out.h>
+#include <deal.II/grid/grid_tools.h>
+#include <deal.II/fe/mapping_q.h>
+
+#include <fstream>
+#include <iomanip>
+
+
+int main ()
+{
+  const int dim = 2;
+
+  Triangulation< dim > tria;
+  std::map< unsigned int, Point< dim > > new_points;
+  const unsigned int N = 16;
+  GridGenerator::subdivided_hyper_cube (tria, N, -5, 5);
+
+  // find the vertex at the origin
+  Triangulation<dim>::active_cell_iterator
+    cell = GridTools::find_active_cell_around_point    (tria,Point<dim>());
+
+  unsigned int best_vertex   = cell->vertex_index(0);//vertex number on local triangulation
+  Point<dim>   best_pos      = cell->vertex(0);
+  double       best_dist     = Point<dim>().distance(best_pos);
+
+  for (unsigned int vertex_no = 1; vertex_no < GeometryInfo<dim>::vertices_per_cell; vertex_no++) {
+    const double dist = Point<dim>().distance(cell->vertex(vertex_no));
+    if (dist < best_dist)
+      {
+       best_pos    = cell->vertex(vertex_no);
+       best_vertex = cell->vertex_index(vertex_no);
+       best_dist   = dist;
+      }
+  }
+  // move the point at the origin by 2 units to the right
+  new_points[best_vertex] = Point<dim>();
+  new_points[best_vertex][0] += 2;
+
+  // now pin all of the points on the boundary
+  cell = tria.begin_active();
+  Triangulation<dim>::active_cell_iterator endc = tria.end();
+
+  for ( ; cell != endc; ++cell)
+    if (cell->at_boundary() == true)
+      for (unsigned int face = 0; face < GeometryInfo<dim>::faces_per_cell;++face)
+       if (cell->face(face)->at_boundary() == true)
+         for (unsigned int v=0; v < GeometryInfo<dim>::vertices_per_face;++v) {
+           unsigned int vertex_number = cell->face(face)->vertex_index(v);
+           new_points[vertex_number] =  cell->face(face)->vertex(v);
+         }
+
+  // then compute new point locations and output the result
+  GridTools::laplace_transform (new_points,tria);
+  std::ofstream out ("output");
+  GridOut grid_out;
+  grid_out.write_eps (tria, out);
+}

diff --git a/tests/deal.II/grid_transform_02.output b/tests/deal.II/grid_transform_02.output
new file mode 100644 (file)
index 0000000..86305e7
--- /dev/null
+++ b/tests/deal.II/grid_transform_02.output
@@ -0,0 +1,653 @@
+%!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 300 150 m 268.547 145.754 x
+b 297.118 179.264 m 264.131 167.89 x
+b 268.547 145.754 m 264.131 167.89 x
+b 297.118 179.264 m 264.131 167.89 x
+b 288.582 207.403 m 257.033 189.709 x
+b 264.131 167.89 m 257.033 189.709 x
+b 268.547 145.754 m 236.579 140.339 x
+b 264.131 167.89 m 230.508 154.954 x
+b 268.547 145.754 m 264.131 167.89 x
+b 236.579 140.339 m 230.508 154.954 x
+b 264.131 167.89 m 230.508 154.954 x
+b 257.033 189.709 m 224.848 170.233 x
+b 264.131 167.89 m 257.033 189.709 x
+b 230.508 154.954 m 224.848 170.233 x
+b 288.582 207.403 m 257.033 189.709 x
+b 274.72 233.336 m 246.575 210.274 x
+b 257.033 189.709 m 246.575 210.274 x
+b 274.72 233.336 m 246.575 210.274 x
+b 256.066 256.066 m 233.826 229.58 x
+b 246.575 210.274 m 233.826 229.58 x
+b 257.033 189.709 m 224.848 170.233 x
+b 246.575 210.274 m 217.903 185.384 x
+b 257.033 189.709 m 246.575 210.274 x
+b 224.848 170.233 m 217.903 185.384 x
+b 246.575 210.274 m 217.903 185.384 x
+b 233.826 229.58 m 211.221 201.559 x
+b 246.575 210.274 m 233.826 229.58 x
+b 217.903 185.384 m 211.221 201.559 x
+b 236.579 140.339 m 203.395 132.808 x
+b 230.508 154.954 m 199.876 141.732 x
+b 236.579 140.339 m 230.508 154.954 x
+b 203.395 132.808 m 199.876 141.732 x
+b 230.508 154.954 m 199.876 141.732 x
+b 224.848 170.233 m 196.452 151.101 x
+b 230.508 154.954 m 224.848 170.233 x
+b 199.876 141.732 m 196.452 151.101 x
+b 203.395 132.808 m 165 120 x
+b 199.876 141.732 m 164.712 122.926 x
+b 203.395 132.808 m 199.876 141.732 x
+b 199.876 141.732 m 164.712 122.926 x
+b 196.452 151.101 m 163.858 125.74 x
+b 199.876 141.732 m 196.452 151.101 x
+b 224.848 170.233 m 196.452 151.101 x
+b 217.903 185.384 m 192.136 160.42 x
+b 224.848 170.233 m 217.903 185.384 x
+b 196.452 151.101 m 192.136 160.42 x
+b 217.903 185.384 m 192.136 160.42 x
+b 211.221 201.559 m 187.756 170.564 x
+b 217.903 185.384 m 211.221 201.559 x
+b 192.136 160.42 m 187.756 170.564 x
+b 196.452 151.101 m 163.858 125.74 x
+b 192.136 160.42 m 162.472 128.334 x
+b 196.452 151.101 m 192.136 160.42 x
+b 192.136 160.42 m 162.472 128.334 x
+b 187.756 170.564 m 160.607 130.607 x
+b 192.136 160.42 m 187.756 170.564 x
+b 256.066 256.066 m 233.826 229.58 x
+b 233.336 274.72 m 214.831 242.018 x
+b 233.826 229.58 m 214.831 242.018 x
+b 233.336 274.72 m 214.831 242.018 x
+b 207.403 288.582 m 194.334 252.407 x
+b 214.831 242.018 m 194.334 252.407 x
+b 233.826 229.58 m 211.221 201.559 x
+b 214.831 242.018 m 195.952 207.335 x
+b 233.826 229.58 m 214.831 242.018 x
+b 211.221 201.559 m 195.952 207.335 x
+b 214.831 242.018 m 195.952 207.335 x
+b 194.334 252.407 m 181.003 214.078 x
+b 214.831 242.018 m 194.334 252.407 x
+b 195.952 207.335 m 181.003 214.078 x
+b 207.403 288.582 m 194.334 252.407 x
+b 179.264 297.118 m 172.447 259.574 x
+b 194.334 252.407 m 172.447 259.574 x
+b 179.264 297.118 m 172.447 259.574 x
+b 150 300 m 150 264.302 x
+b 172.447 259.574 m 150 264.302 x
+b 194.334 252.407 m 181.003 214.078 x
+b 172.447 259.574 m 165.522 219.94 x
+b 194.334 252.407 m 172.447 259.574 x
+b 181.003 214.078 m 165.522 219.94 x
+b 172.447 259.574 m 165.522 219.94 x
+b 150 264.302 m 150 226.918 x
+b 172.447 259.574 m 150 264.302 x
+b 165.522 219.94 m 150 226.918 x
+b 211.221 201.559 m 187.756 170.564 x
+b 195.952 207.335 m 178.4 174.155 x
+b 211.221 201.559 m 195.952 207.335 x
+b 187.756 170.564 m 178.4 174.155 x
+b 195.952 207.335 m 178.4 174.155 x
+b 181.003 214.078 m 169.241 178.312 x
+b 195.952 207.335 m 181.003 214.078 x
+b 178.4 174.155 m 169.241 178.312 x
+b 187.756 170.564 m 160.607 130.607 x
+b 178.4 174.155 m 158.334 132.472 x
+b 187.756 170.564 m 178.4 174.155 x
+b 178.4 174.155 m 158.334 132.472 x
+b 169.241 178.312 m 155.74 133.858 x
+b 178.4 174.155 m 169.241 178.312 x
+b 181.003 214.078 m 169.241 178.312 x
+b 165.522 219.94 m 159.712 181.896 x
+b 181.003 214.078 m 165.522 219.94 x
+b 169.241 178.312 m 159.712 181.896 x
+b 165.522 219.94 m 159.712 181.896 x
+b 150 226.918 m 150 186.203 x
+b 165.522 219.94 m 150 226.918 x
+b 159.712 181.896 m 150 186.203 x
+b 169.241 178.312 m 155.74 133.858 x
+b 159.712 181.896 m 152.926 134.712 x
+b 169.241 178.312 m 159.712 181.896 x
+b 159.712 181.896 m 152.926 134.712 x
+b 150 186.203 m 150 135 x
+b 159.712 181.896 m 150 186.203 x
+b 150 300 m 150 264.302 x
+b 120.736 297.118 m 127.553 259.574 x
+b 150 264.302 m 127.553 259.574 x
+b 120.736 297.118 m 127.553 259.574 x
+b 92.5975 288.582 m 105.666 252.407 x
+b 127.553 259.574 m 105.666 252.407 x
+b 150 264.302 m 150 226.918 x
+b 127.553 259.574 m 134.478 219.94 x
+b 150 264.302 m 127.553 259.574 x
+b 150 226.918 m 134.478 219.94 x
+b 127.553 259.574 m 134.478 219.94 x
+b 105.666 252.407 m 118.997 214.078 x
+b 127.553 259.574 m 105.666 252.407 x
+b 134.478 219.94 m 118.997 214.078 x
+b 92.5975 288.582 m 105.666 252.407 x
+b 66.6645 274.72 m 85.1692 242.018 x
+b 105.666 252.407 m 85.1692 242.018 x
+b 66.6645 274.72 m 85.1692 242.018 x
+b 43.934 256.066 m 66.1743 229.58 x
+b 85.1692 242.018 m 66.1743 229.58 x
+b 105.666 252.407 m 118.997 214.078 x
+b 85.1692 242.018 m 104.048 207.335 x
+b 105.666 252.407 m 85.1692 242.018 x
+b 118.997 214.078 m 104.048 207.335 x
+b 85.1692 242.018 m 104.048 207.335 x
+b 66.1743 229.58 m 88.7794 201.559 x
+b 85.1692 242.018 m 66.1743 229.58 x
+b 104.048 207.335 m 88.7794 201.559 x
+b 150 226.918 m 150 186.203 x
+b 134.478 219.94 m 140.288 181.896 x
+b 150 226.918 m 134.478 219.94 x
+b 150 186.203 m 140.288 181.896 x
+b 134.478 219.94 m 140.288 181.896 x
+b 118.997 214.078 m 130.759 178.312 x
+b 134.478 219.94 m 118.997 214.078 x
+b 140.288 181.896 m 130.759 178.312 x
+b 150 186.203 m 150 135 x
+b 140.288 181.896 m 147.074 134.712 x
+b 150 186.203 m 140.288 181.896 x
+b 140.288 181.896 m 147.074 134.712 x
+b 130.759 178.312 m 144.26 133.858 x
+b 140.288 181.896 m 130.759 178.312 x
+b 118.997 214.078 m 130.759 178.312 x
+b 104.048 207.335 m 121.6 174.155 x
+b 118.997 214.078 m 104.048 207.335 x
+b 130.759 178.312 m 121.6 174.155 x
+b 104.048 207.335 m 121.6 174.155 x
+b 88.7794 201.559 m 112.244 170.564 x
+b 104.048 207.335 m 88.7794 201.559 x
+b 121.6 174.155 m 112.244 170.564 x
+b 130.759 178.312 m 144.26 133.858 x
+b 121.6 174.155 m 141.666 132.472 x
+b 130.759 178.312 m 121.6 174.155 x
+b 121.6 174.155 m 141.666 132.472 x
+b 112.244 170.564 m 139.393 130.607 x
+b 121.6 174.155 m 112.244 170.564 x
+b 43.934 256.066 m 66.1743 229.58 x
+b 25.2796 233.336 m 53.4245 210.274 x
+b 66.1743 229.58 m 53.4245 210.274 x
+b 25.2796 233.336 m 53.4245 210.274 x
+b 11.4181 207.403 m 42.9675 189.709 x
+b 53.4245 210.274 m 42.9675 189.709 x
+b 66.1743 229.58 m 88.7794 201.559 x
+b 53.4245 210.274 m 82.0968 185.384 x
+b 66.1743 229.58 m 53.4245 210.274 x
+b 88.7794 201.559 m 82.0968 185.384 x
+b 53.4245 210.274 m 82.0968 185.384 x
+b 42.9675 189.709 m 75.1517 170.233 x
+b 53.4245 210.274 m 42.9675 189.709 x
+b 82.0968 185.384 m 75.1517 170.233 x
+b 11.4181 207.403 m 42.9675 189.709 x
+b 2.88221 179.264 m 35.8686 167.89 x
+b 42.9675 189.709 m 35.8686 167.89 x
+b 2.88221 179.264 m 35.8686 167.89 x
+b 0 150 m 31.4526 145.754 x
+b 35.8686 167.89 m 31.4526 145.754 x
+b 42.9675 189.709 m 75.1517 170.233 x
+b 35.8686 167.89 m 69.4924 154.954 x
+b 42.9675 189.709 m 35.8686 167.89 x
+b 75.1517 170.233 m 69.4924 154.954 x
+b 35.8686 167.89 m 69.4924 154.954 x
+b 31.4526 145.754 m 63.421 140.339 x
+b 35.8686 167.89 m 31.4526 145.754 x
+b 69.4924 154.954 m 63.421 140.339 x
+b 88.7794 201.559 m 112.244 170.564 x
+b 82.0968 185.384 m 107.864 160.42 x
+b 88.7794 201.559 m 82.0968 185.384 x
+b 112.244 170.564 m 107.864 160.42 x
+b 82.0968 185.384 m 107.864 160.42 x
+b 75.1517 170.233 m 103.548 151.101 x
+b 82.0968 185.384 m 75.1517 170.233 x
+b 107.864 160.42 m 103.548 151.101 x
+b 112.244 170.564 m 139.393 130.607 x
+b 107.864 160.42 m 137.528 128.334 x
+b 112.244 170.564 m 107.864 160.42 x
+b 107.864 160.42 m 137.528 128.334 x
+b 103.548 151.101 m 136.142 125.74 x
+b 107.864 160.42 m 103.548 151.101 x
+b 75.1517 170.233 m 103.548 151.101 x
+b 69.4924 154.954 m 100.124 141.732 x
+b 75.1517 170.233 m 69.4924 154.954 x
+b 103.548 151.101 m 100.124 141.732 x
+b 69.4924 154.954 m 100.124 141.732 x
+b 63.421 140.339 m 96.6048 132.808 x
+b 69.4924 154.954 m 63.421 140.339 x
+b 100.124 141.732 m 96.6048 132.808 x
+b 103.548 151.101 m 136.142 125.74 x
+b 100.124 141.732 m 135.288 122.926 x
+b 103.548 151.101 m 100.124 141.732 x
+b 100.124 141.732 m 135.288 122.926 x
+b 96.6048 132.808 m 135 120 x
+b 100.124 141.732 m 96.6048 132.808 x
+b 0 150 m 31.4526 145.754 x
+b 2.88221 120.736 m 35.8686 122.996 x
+b 31.4526 145.754 m 35.8686 122.996 x
+b 2.88221 120.736 m 35.8686 122.996 x
+b 11.4181 92.5975 m 42.9675 101.04 x
+b 35.8686 122.996 m 42.9675 101.04 x
+b 31.4526 145.754 m 63.421 140.339 x
+b 35.8686 122.996 m 69.4924 123.91 x
+b 31.4526 145.754 m 35.8686 122.996 x
+b 63.421 140.339 m 69.4924 123.91 x
+b 35.8686 122.996 m 69.4924 123.91 x
+b 42.9675 101.04 m 75.1517 108.227 x
+b 35.8686 122.996 m 42.9675 101.04 x
+b 69.4924 123.91 m 75.1517 108.227 x
+b 11.4181 92.5975 m 42.9675 101.04 x
+b 25.2796 66.6645 m 53.4245 80.6122 x
+b 42.9675 101.04 m 53.4245 80.6122 x
+b 25.2796 66.6645 m 53.4245 80.6122 x
+b 43.934 43.934 m 66.1743 61.9287 x
+b 53.4245 80.6122 m 66.1743 61.9287 x
+b 42.9675 101.04 m 75.1517 108.227 x
+b 53.4245 80.6122 m 82.0968 93.4804 x
+b 42.9675 101.04 m 53.4245 80.6122 x
+b 75.1517 108.227 m 82.0968 93.4804 x
+b 53.4245 80.6122 m 82.0968 93.4804 x
+b 66.1743 61.9287 m 88.7794 79.1183 x
+b 53.4245 80.6122 m 66.1743 61.9287 x
+b 82.0968 93.4804 m 88.7794 79.1183 x
+b 63.421 140.339 m 96.6048 132.808 x
+b 69.4924 123.91 m 100.124 122.308 x
+b 63.421 140.339 m 69.4924 123.91 x
+b 96.6048 132.808 m 100.124 122.308 x
+b 69.4924 123.91 m 100.124 122.308 x
+b 75.1517 108.227 m 103.548 112.619 x
+b 69.4924 123.91 m 75.1517 108.227 x
+b 100.124 122.308 m 103.548 112.619 x
+b 96.6048 132.808 m 135 120 x
+b 100.124 122.308 m 135.288 117.074 x
+b 96.6048 132.808 m 100.124 122.308 x
+b 100.124 122.308 m 135.288 117.074 x
+b 103.548 112.619 m 136.142 114.26 x
+b 100.124 122.308 m 103.548 112.619 x
+b 75.1517 108.227 m 103.548 112.619 x
+b 82.0968 93.4804 m 107.864 103.619 x
+b 75.1517 108.227 m 82.0968 93.4804 x
+b 103.548 112.619 m 107.864 103.619 x
+b 82.0968 93.4804 m 107.864 103.619 x
+b 88.7794 79.1183 m 112.244 95.052 x
+b 82.0968 93.4804 m 88.7794 79.1183 x
+b 107.864 103.619 m 112.244 95.052 x
+b 103.548 112.619 m 136.142 114.26 x
+b 107.864 103.619 m 137.528 111.666 x
+b 103.548 112.619 m 107.864 103.619 x
+b 107.864 103.619 m 137.528 111.666 x
+b 112.244 95.052 m 139.393 109.393 x
+b 107.864 103.619 m 112.244 95.052 x
+b 43.934 43.934 m 66.1743 61.9287 x
+b 66.6645 25.2796 m 85.1692 48.8674 x
+b 66.1743 61.9287 m 85.1692 48.8674 x
+b 66.6645 25.2796 m 85.1692 48.8674 x
+b 92.5975 11.4181 m 105.666 38.3417 x
+b 85.1692 48.8674 m 105.666 38.3417 x
+b 66.1743 61.9287 m 88.7794 79.1183 x
+b 85.1692 48.8674 m 104.048 71.5288 x
+b 66.1743 61.9287 m 85.1692 48.8674 x
+b 88.7794 79.1183 m 104.048 71.5288 x
+b 85.1692 48.8674 m 104.048 71.5288 x
+b 105.666 38.3417 m 118.997 64.3816 x
+b 85.1692 48.8674 m 105.666 38.3417 x
+b 104.048 71.5288 m 118.997 64.3816 x
+b 92.5975 11.4181 m 105.666 38.3417 x
+b 120.736 2.88221 m 127.553 31.3115 x
+b 105.666 38.3417 m 127.553 31.3115 x
+b 120.736 2.88221 m 127.553 31.3115 x
+b 150 0 m 150 27.207 x
+b 127.553 31.3115 m 150 27.207 x
+b 105.666 38.3417 m 118.997 64.3816 x
+b 127.553 31.3115 m 134.478 58.9245 x
+b 105.666 38.3417 m 127.553 31.3115 x
+b 118.997 64.3816 m 134.478 58.9245 x
+b 127.553 31.3115 m 134.478 58.9245 x
+b 150 27.207 m 150 53.7599 x
+b 127.553 31.3115 m 150 27.207 x
+b 134.478 58.9245 m 150 53.7599 x
+b 88.7794 79.1183 m 112.244 95.052 x
+b 104.048 71.5288 m 121.6 89.8844 x
+b 88.7794 79.1183 m 104.048 71.5288 x
+b 112.244 95.052 m 121.6 89.8844 x
+b 104.048 71.5288 m 121.6 89.8844 x
+b 118.997 64.3816 m 130.759 85.4079 x
+b 104.048 71.5288 m 118.997 64.3816 x
+b 121.6 89.8844 m 130.759 85.4079 x
+b 112.244 95.052 m 139.393 109.393 x
+b 121.6 89.8844 m 141.666 107.528 x
+b 112.244 95.052 m 121.6 89.8844 x
+b 121.6 89.8844 m 141.666 107.528 x
+b 130.759 85.4079 m 144.26 106.142 x
+b 121.6 89.8844 m 130.759 85.4079 x
+b 118.997 64.3816 m 130.759 85.4079 x
+b 134.478 58.9245 m 140.288 82.1435 x
+b 118.997 64.3816 m 134.478 58.9245 x
+b 130.759 85.4079 m 140.288 82.1435 x
+b 134.478 58.9245 m 140.288 82.1435 x
+b 150 53.7599 m 150 79.4129 x
+b 134.478 58.9245 m 150 53.7599 x
+b 140.288 82.1435 m 150 79.4129 x
+b 130.759 85.4079 m 144.26 106.142 x
+b 140.288 82.1435 m 147.074 105.288 x
+b 130.759 85.4079 m 140.288 82.1435 x
+b 140.288 82.1435 m 147.074 105.288 x
+b 150 79.4129 m 150 105 x
+b 140.288 82.1435 m 150 79.4129 x
+b 150 0 m 150 27.207 x
+b 179.264 2.88221 m 172.447 31.3115 x
+b 150 27.207 m 172.447 31.3115 x
+b 179.264 2.88221 m 172.447 31.3115 x
+b 207.403 11.4181 m 194.334 38.3417 x
+b 172.447 31.3115 m 194.334 38.3417 x
+b 150 27.207 m 150 53.7599 x
+b 172.447 31.3115 m 165.522 58.9245 x
+b 150 27.207 m 172.447 31.3115 x
+b 150 53.7599 m 165.522 58.9245 x
+b 172.447 31.3115 m 165.522 58.9245 x
+b 194.334 38.3417 m 181.003 64.3816 x
+b 172.447 31.3115 m 194.334 38.3417 x
+b 165.522 58.9245 m 181.003 64.3816 x
+b 207.403 11.4181 m 194.334 38.3417 x
+b 233.336 25.2796 m 214.831 48.8674 x
+b 194.334 38.3417 m 214.831 48.8674 x
+b 233.336 25.2796 m 214.831 48.8674 x
+b 256.066 43.934 m 233.826 61.9287 x
+b 214.831 48.8674 m 233.826 61.9287 x
+b 194.334 38.3417 m 181.003 64.3816 x
+b 214.831 48.8674 m 195.952 71.5288 x
+b 194.334 38.3417 m 214.831 48.8674 x
+b 181.003 64.3816 m 195.952 71.5288 x
+b 214.831 48.8674 m 195.952 71.5288 x
+b 233.826 61.9287 m 211.221 79.1183 x
+b 214.831 48.8674 m 233.826 61.9287 x
+b 195.952 71.5288 m 211.221 79.1183 x
+b 150 53.7599 m 150 79.4129 x
+b 165.522 58.9245 m 159.712 82.1435 x
+b 150 53.7599 m 165.522 58.9245 x
+b 150 79.4129 m 159.712 82.1435 x
+b 165.522 58.9245 m 159.712 82.1435 x
+b 181.003 64.3816 m 169.241 85.4079 x
+b 165.522 58.9245 m 181.003 64.3816 x
+b 159.712 82.1435 m 169.241 85.4079 x
+b 150 79.4129 m 150 105 x
+b 159.712 82.1435 m 152.926 105.288 x
+b 150 79.4129 m 159.712 82.1435 x
+b 159.712 82.1435 m 152.926 105.288 x
+b 169.241 85.4079 m 155.74 106.142 x
+b 159.712 82.1435 m 169.241 85.4079 x
+b 181.003 64.3816 m 169.241 85.4079 x
+b 195.952 71.5288 m 178.4 89.8844 x
+b 181.003 64.3816 m 195.952 71.5288 x
+b 169.241 85.4079 m 178.4 89.8844 x
+b 195.952 71.5288 m 178.4 89.8844 x
+b 211.221 79.1183 m 187.756 95.052 x
+b 195.952 71.5288 m 211.221 79.1183 x
+b 178.4 89.8844 m 187.756 95.052 x
+b 169.241 85.4079 m 155.74 106.142 x
+b 178.4 89.8844 m 158.334 107.528 x
+b 169.241 85.4079 m 178.4 89.8844 x
+b 178.4 89.8844 m 158.334 107.528 x
+b 187.756 95.052 m 160.607 109.393 x
+b 178.4 89.8844 m 187.756 95.052 x
+b 256.066 43.934 m 233.826 61.9287 x
+b 274.72 66.6645 m 246.575 80.6122 x
+b 233.826 61.9287 m 246.575 80.6122 x
+b 274.72 66.6645 m 246.575 80.6122 x
+b 288.582 92.5975 m 257.033 101.04 x
+b 246.575 80.6122 m 257.033 101.04 x
+b 233.826 61.9287 m 211.221 79.1183 x
+b 246.575 80.6122 m 217.903 93.4804 x
+b 233.826 61.9287 m 246.575 80.6122 x
+b 211.221 79.1183 m 217.903 93.4804 x
+b 246.575 80.6122 m 217.903 93.4804 x
+b 257.033 101.04 m 224.848 108.227 x
+b 246.575 80.6122 m 257.033 101.04 x
+b 217.903 93.4804 m 224.848 108.227 x
+b 288.582 92.5975 m 257.033 101.04 x
+b 297.118 120.736 m 264.131 122.996 x
+b 257.033 101.04 m 264.131 122.996 x
+b 297.118 120.736 m 264.131 122.996 x
+b 300 150 m 268.547 145.754 x
+b 264.131 122.996 m 268.547 145.754 x
+b 257.033 101.04 m 224.848 108.227 x
+b 264.131 122.996 m 230.508 123.91 x
+b 257.033 101.04 m 264.131 122.996 x
+b 224.848 108.227 m 230.508 123.91 x
+b 264.131 122.996 m 230.508 123.91 x
+b 268.547 145.754 m 236.579 140.339 x
+b 264.131 122.996 m 268.547 145.754 x
+b 230.508 123.91 m 236.579 140.339 x
+b 211.221 79.1183 m 187.756 95.052 x
+b 217.903 93.4804 m 192.136 103.619 x
+b 211.221 79.1183 m 217.903 93.4804 x
+b 187.756 95.052 m 192.136 103.619 x
+b 217.903 93.4804 m 192.136 103.619 x
+b 224.848 108.227 m 196.452 112.619 x
+b 217.903 93.4804 m 224.848 108.227 x
+b 192.136 103.619 m 196.452 112.619 x
+b 187.756 95.052 m 160.607 109.393 x
+b 192.136 103.619 m 162.472 111.666 x
+b 187.756 95.052 m 192.136 103.619 x
+b 192.136 103.619 m 162.472 111.666 x
+b 196.452 112.619 m 163.858 114.26 x
+b 192.136 103.619 m 196.452 112.619 x
+b 224.848 108.227 m 196.452 112.619 x
+b 230.508 123.91 m 199.876 122.308 x
+b 224.848 108.227 m 230.508 123.91 x
+b 196.452 112.619 m 199.876 122.308 x
+b 230.508 123.91 m 199.876 122.308 x
+b 236.579 140.339 m 203.395 132.808 x
+b 230.508 123.91 m 236.579 140.339 x
+b 199.876 122.308 m 203.395 132.808 x
+b 196.452 112.619 m 163.858 114.26 x
+b 199.876 122.308 m 164.712 117.074 x
+b 196.452 112.619 m 199.876 122.308 x
+b 199.876 122.308 m 164.712 117.074 x
+b 203.395 132.808 m 165 120 x
+b 199.876 122.308 m 203.395 132.808 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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
+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 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 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
+showpage