From fbd57782484f463db00e3cd63754dda4596ae5e7 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Mon, 24 Mar 2014 03:14:05 +0000 Subject: [PATCH] Adjust implementation of laplace_transform to what had been documented. git-svn-id: https://svn.dealii.org/trunk@32685 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/deal.II/grid_transform_02.cc | 89 ++++ tests/deal.II/grid_transform_02.output | 653 +++++++++++++++++++++++++ 2 files changed, 742 insertions(+) create mode 100644 tests/deal.II/grid_transform_02.cc create mode 100644 tests/deal.II/grid_transform_02.output diff --git a/tests/deal.II/grid_transform_02.cc b/tests/deal.II/grid_transform_02.cc new file mode 100644 index 0000000000..9ceb05b33c --- /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 +#include +#include +#include +#include +#include +#include +#include + +#include +#include + + +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::active_cell_iterator + cell = GridTools::find_active_cell_around_point (tria,Point()); + + unsigned int best_vertex = cell->vertex_index(0);//vertex number on local triangulation + Point best_pos = cell->vertex(0); + double best_dist = Point().distance(best_pos); + + for (unsigned int vertex_no = 1; vertex_no < GeometryInfo::vertices_per_cell; vertex_no++) { + const double dist = Point().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(); + new_points[best_vertex][0] += 2; + + // now pin all of the points on the boundary + cell = tria.begin_active(); + Triangulation::active_cell_iterator endc = tria.end(); + + for ( ; cell != endc; ++cell) + if (cell->at_boundary() == true) + for (unsigned int face = 0; face < GeometryInfo::faces_per_cell;++face) + if (cell->face(face)->at_boundary() == true) + for (unsigned int v=0; v < GeometryInfo::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 index 0000000000..86305e75be --- /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 -- 2.39.5