From 42e679f530180b02914d9b3d9e7b1437e352659c Mon Sep 17 00:00:00 2001 From: Ralf Hartmann Date: Fri, 18 May 2001 17:29:12 +0000 Subject: [PATCH] Test for the new GridGenerator::laplace_transformation function. git-svn-id: https://svn.dealii.org/trunk@4676 0785d39b-7218-0410-832d-ea1e28bc413d --- tests/deal.II/grid_transform.cc | 118 ++++++++ tests/deal.II/grid_transform.checked | 429 +++++++++++++++++++++++++++ 2 files changed, 547 insertions(+) create mode 100644 tests/deal.II/grid_transform.cc create mode 100644 tests/deal.II/grid_transform.checked diff --git a/tests/deal.II/grid_transform.cc b/tests/deal.II/grid_transform.cc new file mode 100644 index 0000000000..aa4a396fe6 --- /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 +#include +#include +#include +#include +#include +#include + +#include + + +int main () +{ + const unsigned int dim=2; + Point origin; + HyperShellBoundary boundary(origin); + MappingQ mapping(2); + Triangulation 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 > new_points; + + // new center and new radius + // of the inner circle. + const Point n_center(0,-1); + const double n_radius=0.5; + + Triangulation::cell_iterator cell=tria.begin_active(), + endc=tria.end(); + Triangulation::face_iterator face; + for (; cell!=endc; ++cell) + { + if (cell->at_boundary()) + for (unsigned int face_no=0; face_no::faces_per_cell; ++face_no) + { + face=cell->face(face_no); + if (face->at_boundary()) + for (unsigned int vertex_no=0; + vertex_no::vertices_per_face; ++vertex_no) + { + const Point &v=face->vertex(vertex_no); + if (fabs(sqrt(v.square())-outer_radius)<1e-15) + { + // leave the + // point, where + // they are. + new_points.insert(pair > ( + 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 > ( + 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 >::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 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 index 0000000000..3922ee5411 --- /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 -- 2.39.5