From dda44a719288ef5fcc931e03355885179797d3e1 Mon Sep 17 00:00:00 2001 From: Guido Kanschat Date: Thu, 31 Jul 2014 15:27:03 +0200 Subject: [PATCH] all existing multigrid tests verified with fixed MGConstrainedDoFs new test for Neumann data --- .../deal.II/multigrid/mg_constrained_dofs.h | 29 +- tests/multigrid/mg_output_dirichlet.cc | 299 ++++++++++++++++++ tests/multigrid/mg_output_dirichlet.output | 37 +++ .../{mg_output.cc => mg_output_neumann.cc} | 39 +-- ...output.output => mg_output_neumann.output} | 0 tests/multigrid/step-16-bdry1.output | 42 +-- tests/multigrid/step-39-02a.cc | 8 +- 7 files changed, 396 insertions(+), 58 deletions(-) create mode 100644 tests/multigrid/mg_output_dirichlet.cc create mode 100644 tests/multigrid/mg_output_dirichlet.output rename tests/multigrid/{mg_output.cc => mg_output_neumann.cc} (89%) rename tests/multigrid/{mg_output.output => mg_output_neumann.output} (100%) diff --git a/include/deal.II/multigrid/mg_constrained_dofs.h b/include/deal.II/multigrid/mg_constrained_dofs.h index 843737b690..b9795fc0f5 100644 --- a/include/deal.II/multigrid/mg_constrained_dofs.h +++ b/include/deal.II/multigrid/mg_constrained_dofs.h @@ -42,6 +42,20 @@ class MGConstrainedDoFs : public Subscriptor public: typedef std::vector >::size_type size_dof; + /** + * Fill the internal data structures with hanging node constraints + * extracted from the dof handler object. Works with natural + * boundary conditions only. There exists a sister function setting + * up boundary constraints as well. + * + * This function ensures that on every level, degrees of freedom at interior + * edges of a refinement level are treated corrected but leaves degrees of + * freedom at the boundary of the domain untouched assuming that no + * Dirichlet boundary conditions for them exist. + */ + template + void initialize (const DoFHandler &dof); + /** * Fill the internal data structures with values extracted from the dof * handler object and apply the boundary values provided. @@ -137,21 +151,6 @@ public: private: - /** - * @warning This function generates an inconsistent object if not - * called from the other initialize() in this class. - * - * Fill the internal data structures with values extracted from the dof - * handler object. - * - * This function ensures that on every level, degrees of freedom at interior - * edges of a refinement level are treated corrected but leaves degrees of - * freedom at the boundary of the domain untouched assuming that no - * Dirichlet boundary conditions for them exist. - */ - template - void initialize (const DoFHandler &dof); - /** * The indices of boundary dofs for each level. */ diff --git a/tests/multigrid/mg_output_dirichlet.cc b/tests/multigrid/mg_output_dirichlet.cc new file mode 100644 index 0000000000..bf44280a83 --- /dev/null +++ b/tests/multigrid/mg_output_dirichlet.cc @@ -0,0 +1,299 @@ +// --------------------------------------------------------------------- +// $Id$ +// +// Copyright (C) 2000 - 2013 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. +// +// --------------------------------------------------------------------- + + +// like _01 but on adaptively refined grid + +//TODO:[GK] Add checks for RT again! + +#include "../tests.h" +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include + +#include +#include +#include +#include + +using namespace std; + +template +void +reinit_vector (const dealii::DoFHandler &mg_dof, + MGLevelObject > &v) +{ + for (unsigned int level=v.min_level(); + level<=v.max_level(); ++level) + { + unsigned int n = mg_dof.n_dofs (level); + v[level].reinit(n); + } +} + + +template +void refine_mesh (Triangulation &triangulation) +{ + bool cell_refined = false; + for (typename Triangulation::active_cell_iterator + cell = triangulation.begin_active(); + cell != triangulation.end(); ++cell) + { + for (unsigned int vertex=0; + vertex < GeometryInfo::vertices_per_cell; + ++vertex) + { + const Point p = cell->vertex(vertex); + const Point origin = (dim == 2 ? + Point(0,0) : + Point(0,0,0)); + const double dist = p.distance(origin); + if (dist<0.25/numbers::PI) + { + cell->set_refine_flag (); + cell_refined = true; + break; + } + } + } + if (!cell_refined) //if no cell was selected for refinement, refine global + for (typename Triangulation::active_cell_iterator + cell = triangulation.begin_active(); + cell != triangulation.end(); ++cell) + cell->set_refine_flag(); + triangulation.execute_coarsening_and_refinement (); +} + +template +void initialize (const DoFHandler &dof, + Vector &u) +{ + unsigned int counter=0; + const unsigned int dofs_per_cell = dof.get_fe().dofs_per_cell; + std::vector dof_indices(dofs_per_cell); + for (typename DoFHandler::active_cell_iterator + cell = dof.begin_active(); + cell != dof.end(); ++cell) + { + cell->get_dof_indices(dof_indices); + for (unsigned int i=0; i +void initialize (const DoFHandler &dof, + MGLevelObject > &u) +{ + const unsigned int dofs_per_cell = dof.get_fe().dofs_per_cell; + const unsigned int dofs_per_face = dof.get_fe().dofs_per_face; + std::vector dof_indices(dofs_per_cell); + std::vector face_indices(dofs_per_face); + + for (unsigned int l=0; l::cell_iterator + cell = dof.begin_mg(l); + cell != dof.end_mg(l); ++cell) + { + cell->get_mg_dof_indices(dof_indices); + for (unsigned int f=0; f::faces_per_cell; ++f) + { + cell->face(f)->get_mg_dof_indices(cell->level(),face_indices); + if (cell->face(f)->at_boundary()) + { + for (unsigned int i=0; ilevel()](face_indices[i]) = 1; + } + else + { + for (unsigned int i=0; ilevel()](face_indices[i]) != 1) + u[cell->level()](face_indices[i]) = 0; + } + } + } + } +} + +template +void print (const DoFHandler &dof, + MGLevelObject > &u) +{ + const unsigned int dofs_per_cell = dof.get_fe().dofs_per_cell; + std::vector dof_indices(dofs_per_cell); + for (unsigned int l=0; l::cell_iterator + cell = dof.begin_mg(l); + cell != dof.end_mg(l); ++cell) + { + cell->get_mg_dof_indices(dof_indices); + for (unsigned int i=0; i +void print_diff (const DoFHandler &dof_1, const DoFHandler &dof_2, + const Vector &u, const Vector &v) +{ + Vector diff; + diff.reinit (u); + const unsigned int dofs_per_cell = dof_1.get_fe().dofs_per_cell; + std::vector dof_indices_1(dofs_per_cell); + std::vector dof_indices_2(dofs_per_cell); + for (typename DoFHandler::active_cell_iterator + cell_1 = dof_1.begin_active(), cell_2 = dof_2.begin_active(); + cell_1 != dof_1.end(); ++cell_1, ++cell_2) + { + cell_1->get_dof_indices(dof_indices_1); + cell_2->get_dof_indices(dof_indices_2); + for (unsigned int i=0; i +void check_simple(const FiniteElement &fe) +{ + deallog << fe.get_name() << std::endl; + + Triangulation tr(Triangulation::limit_level_difference_at_vertices); + GridGenerator::hyper_cube(tr,-1,1); + typename FunctionMap::type dirichlet_boundary_functions; + ZeroFunction homogeneous_dirichlet_bc (1); + dirichlet_boundary_functions[0] = &homogeneous_dirichlet_bc; + + tr.refine_global (1); + refine_mesh(tr); + refine_mesh(tr); + refine_mesh(tr); + refine_mesh(tr); + refine_mesh(tr); + refine_mesh(tr); + + DoFHandler mgdof(tr); + mgdof.distribute_dofs(fe); + mgdof.distribute_mg_dofs(fe); + + DoFHandler mgdof_renumbered(tr); + mgdof_renumbered.distribute_dofs(fe); + mgdof_renumbered.distribute_mg_dofs(fe); + + ConstraintMatrix hnc, hnc_renumbered; + hnc.clear (); + hnc_renumbered.clear (); + DoFTools::make_hanging_node_constraints ( + mgdof, + hnc); + + DoFTools::make_hanging_node_constraints ( + mgdof_renumbered, + hnc_renumbered); + hnc.close (); + hnc_renumbered.close (); + + + std::vector block_component (4,0); + block_component[2] = 1; + block_component[3] = 1; + + DoFRenumbering::component_wise (mgdof_renumbered, block_component); + for (unsigned int level=0; level > transfer(hnc, mg_constrained_dofs); + transfer.build_matrices(mgdof); + + MGTransferPrebuilt > transfer_renumbered(hnc_renumbered, mg_constrained_dofs_renumbered); + transfer_renumbered.build_matrices(mgdof_renumbered); + + Vector u(mgdof.n_dofs()); + initialize(mgdof,u); + + MGLevelObject > v(0, tr.n_levels()-1); + reinit_vector(mgdof, v); + + transfer.copy_to_mg(mgdof, v, u); + print(mgdof, v); + + u=0; + initialize(mgdof, v); + for (unsigned int l=0; l diff = u; + + initialize(mgdof_renumbered,u); + reinit_vector(mgdof_renumbered, v); + transfer_renumbered.copy_to_mg(mgdof_renumbered, v, u); + deallog << "copy_to_mg" << std::endl; + print(mgdof_renumbered, v); + + u=0; + initialize(mgdof_renumbered, v); + for (unsigned int l=0; l(FE_Q<2>(1), 2)); + check_simple (FESystem<2>(FE_Q<2>(1), 4)); +} diff --git a/tests/multigrid/mg_output_dirichlet.output b/tests/multigrid/mg_output_dirichlet.output new file mode 100644 index 0000000000..3959b73cd6 --- /dev/null +++ b/tests/multigrid/mg_output_dirichlet.output @@ -0,0 +1,37 @@ + +DEAL::FESystem<2>[FE_Q<2>(1)^4] +DEAL:: +DEAL::Level 0 +DEAL:: 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::Level 1 +DEAL:: 0 0 0 0 0 0 0 0 0 0 0 0 589199 589384 589568 589753 0 0 0 0 0 0 0 0 589199 589384 589568 589753 0 0 0 0 0 0 0 0 589199 589384 589568 589753 0 0 0 0 0 0 0 0 589199 589384 589568 589753 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::Level 2 +DEAL:: 0 0 0 0 0 0 0 0 0 0 0 0 8378 8384 8391 8398 0 0 0 0 0 0 0 0 8378 8384 8391 8398 35911 35928 35946 35964 0 0 0 0 8378 8384 8391 8398 0 0 0 0 38331 38348 38366 38383 8378 8384 8391 8398 35911 35928 35946 35964 38331 38348 38366 38383 499066 499208 499351 499493 0 0 0 0 0 0 0 0 35911 35928 35946 35964 8998 9004 9011 9018 0 0 0 0 0 0 0 0 8998 9004 9011 9018 0 0 0 0 35911 35928 35946 35964 8998 9004 9011 9018 499066 499208 499351 499493 42677 42695 42712 42730 8998 9004 9011 9018 0 0 0 0 42677 42695 42712 42730 0 0 0 0 0 0 0 0 38331 38348 38366 38383 0 0 0 0 9603 9610 9616 9623 38331 38348 38366 38383 499066 499208 499351 499493 9603 9610 9616 9623 44745 44762 44780 44797 0 0 0 0 9603 9610 9616 9623 0 0 0 0 0 0 0 0 9603 9610 9616 9623 44745 44762 44780 44797 0 0 0 0 0 0 0 0 499066 499208 499351 499493 42677 42695 42712 42730 44745 44762 44780 44797 10225 10231 10238 10245 42677 42695 42712 42730 0 0 0 0 10225 10231 10238 10245 0 0 0 0 44745 44762 44780 44797 10225 10231 10238 10245 0 0 0 0 0 0 0 0 10225 10231 10238 10245 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::Level 3 +DEAL:: 1613 1614 1615 1617 2504 2506 2508 2510 2542 2544 2546 2548 13976 13983 13990 13997 2504 2506 2508 2510 1802 1804 1806 1808 13976 13983 13990 13997 48352 48370 48388 48406 2542 2544 2546 2548 13976 13983 13990 13997 2034 2036 2038 2040 51881 51899 51917 51935 13976 13983 13990 13997 48352 48370 48388 48406 51881 51899 51917 51935 373589 373688 373787 373887 1802 1804 1806 1808 2618 2620 2622 2624 48352 48370 48388 48406 14964 14971 14978 14985 2618 2620 2622 2624 1713 1714 1715 1717 14964 14971 14978 14985 2706 2708 2710 2712 48352 48370 48388 48406 14964 14971 14978 14985 373589 373688 373787 373887 58446 58465 58483 58501 14964 14971 14978 14985 2706 2708 2710 2712 58446 58465 58483 58501 2310 2312 2314 2316 2034 2036 2038 2040 51881 51899 51917 51935 2754 2756 2758 2760 15893 15900 15907 15914 51881 51899 51917 51935 373589 373688 373787 373887 15893 15900 15907 15914 61411 61429 61447 61465 2754 2756 2758 2760 15893 15900 15907 15914 1813 1814 1815 1817 2848 2850 2852 2854 15893 15900 15907 15914 61411 61429 61447 61465 2848 2850 2852 2854 2482 2484 2486 2488 373589 373688 373787 373887 58446 58465 58483 58501 61411 61429 61447 61465 16888 16895 16902 16909 58446 58465 58483 58501 2310 2312 2314 2316 16888 16895 16902 16909 2918 2920 2922 2924 61411 61429 61447 61465 16888 16895 16902 16909 2482 2484 2486 2488 2962 2964 2966 2968 16888 16895 16902 16909 2918 2920 2922 2924 2962 2964 2966 2968 1913 1914 1915 1917 533 534 535 536 17 18 19 20 565 566 567 568 49 50 51 52 17 18 19 20 129 130 131 132 49 50 51 52 161 162 163 164 565 566 567 568 49 50 51 52 1613 1614 1615 1617 2504 2506 2508 2510 49 50 51 52 161 162 163 164 2504 2506 2508 2510 1802 1804 1806 1808 553 554 555 556 569 570 571 572 97 98 99 100 113 114 115 116 569 570 571 572 1613 1614 1615 1617 113 114 115 116 2542 2544 2546 2548 97 98 99 100 113 114 115 116 257 258 259 260 273 274 275 276 113 114 115 116 2542 2544 2546 2548 273 274 275 276 2034 2036 2038 2040 129 130 131 132 145 146 147 148 161 162 163 164 177 178 179 180 145 146 147 148 577 578 579 580 177 178 179 180 609 610 611 612 161 162 163 164 177 178 179 180 1802 1804 1806 1808 2618 2620 2622 2624 177 178 179 180 609 610 611 612 2618 2620 2622 2624 1713 1714 1715 1717 1713 1714 1715 1717 633 634 635 636 2706 2708 2710 2712 241 242 243 244 633 634 635 636 637 638 639 640 241 242 243 244 245 246 247 248 2706 2708 2710 2712 241 242 243 244 2310 2312 2314 2316 401 402 403 404 241 242 243 244 245 246 247 248 401 402 403 404 405 406 407 408 257 258 259 260 273 274 275 276 289 290 291 292 305 306 307 308 273 274 275 276 2034 2036 2038 2040 305 306 307 308 2754 2756 2758 2760 289 290 291 292 305 306 307 308 641 642 643 644 657 658 659 660 305 306 307 308 2754 2756 2758 2760 657 658 659 660 1813 1814 1815 1817 1813 1814 1815 1817 2848 2850 2852 2854 693 694 695 696 369 370 371 372 2848 2850 2852 2854 2482 2484 2486 2488 369 370 371 372 481 482 483 484 693 694 695 696 369 370 371 372 701 702 703 704 377 378 379 380 369 370 371 372 481 482 483 484 377 378 379 380 489 490 491 492 2310 2312 2314 2316 401 402 403 404 2918 2920 2922 2924 433 434 435 436 401 402 403 404 405 406 407 408 433 434 435 436 437 438 439 440 2918 2920 2922 2924 433 434 435 436 1913 1914 1915 1917 721 722 723 724 433 434 435 436 437 438 439 440 721 722 723 724 725 726 727 728 2482 2484 2486 2488 2962 2964 2966 2968 481 482 483 484 497 498 499 500 2962 2964 2966 2968 1913 1914 1915 1917 497 498 499 500 737 738 739 740 481 482 483 484 497 498 499 500 489 490 491 492 505 506 507 508 497 498 499 500 737 738 739 740 505 506 507 508 745 746 747 748 513 514 515 516 529 530 531 532 545 546 547 548 561 562 563 564 529 530 531 532 533 534 535 536 561 562 563 564 565 566 567 568 545 546 547 548 561 562 563 564 553 554 555 556 569 570 571 572 561 562 563 564 565 566 567 568 569 570 571 572 1613 1614 1615 1617 577 578 579 580 593 594 595 596 609 610 611 612 625 626 627 628 593 594 595 596 597 598 599 600 625 626 627 628 629 630 631 632 609 610 611 612 625 626 627 628 1713 1714 1715 1717 633 634 635 636 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 657 658 659 660 673 674 675 676 689 690 691 692 657 658 659 660 1813 1814 1815 1817 689 690 691 692 693 694 695 696 673 674 675 676 689 690 691 692 681 682 683 684 697 698 699 700 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 1913 1914 1915 1917 721 722 723 724 737 738 739 740 753 754 755 756 721 722 723 724 725 726 727 728 753 754 755 756 757 758 759 760 737 738 739 740 753 754 755 756 745 746 747 748 761 762 763 764 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 +DEAL::Level 4 +DEAL:: 2573 2574 2575 2577 4040 4042 4044 4046 4078 4080 4082 4084 21010 21018 21026 21034 4040 4042 4044 4046 3338 3340 3342 3344 21010 21018 21026 21034 57784 57803 57821 57840 4078 4080 4082 4084 21010 21018 21026 21034 3570 3572 3574 3576 62884 62902 62921 62939 21010 21018 21026 21034 57784 57803 57821 57840 62884 62902 62921 62939 214204 214259 214313 214367 3338 3340 3342 3344 4154 4156 4158 4160 57784 57803 57821 57840 23469 23477 23485 23493 4154 4156 4158 4160 2673 2674 2675 2677 23469 23477 23485 23493 4242 4244 4246 4248 57784 57803 57821 57840 23469 23477 23485 23493 214204 214259 214313 214367 72406 72425 72443 72462 23469 23477 23485 23493 4242 4244 4246 4248 72406 72425 72443 72462 3846 3848 3850 3852 3570 3572 3574 3576 62884 62902 62921 62939 4290 4292 4294 4296 25693 25701 25709 25717 62884 62902 62921 62939 214204 214259 214313 214367 25693 25701 25709 25717 76516 76535 76553 76572 4290 4292 4294 4296 25693 25701 25709 25717 2773 2774 2775 2777 4384 4386 4388 4390 25693 25701 25709 25717 76516 76535 76553 76572 4384 4386 4388 4390 4018 4020 4022 4024 214204 214259 214313 214367 72406 72425 72443 72462 76516 76535 76553 76572 28181 28189 28197 28205 72406 72425 72443 72462 3846 3848 3850 3852 28181 28189 28197 28205 4454 4456 4458 4460 76516 76535 76553 76572 28181 28189 28197 28205 4018 4020 4022 4024 4498 4500 4502 4504 28181 28189 28197 28205 4454 4456 4458 4460 4498 4500 4502 4504 2873 2874 2875 2877 0 0 0 0 0 0 0 0 1333 1334 1335 1336 817 818 819 820 0 0 0 0 0 0 0 0 817 818 819 820 929 930 931 932 1333 1334 1335 1336 817 818 819 820 2573 2574 2575 2577 4040 4042 4044 4046 817 818 819 820 929 930 931 932 4040 4042 4044 4046 3338 3340 3342 3344 0 0 0 0 1337 1338 1339 1340 0 0 0 0 881 882 883 884 1337 1338 1339 1340 2573 2574 2575 2577 881 882 883 884 4078 4080 4082 4084 0 0 0 0 881 882 883 884 0 0 0 0 1041 1042 1043 1044 881 882 883 884 4078 4080 4082 4084 1041 1042 1043 1044 3570 3572 3574 3576 0 0 0 0 0 0 0 0 929 930 931 932 945 946 947 948 0 0 0 0 0 0 0 0 945 946 947 948 1377 1378 1379 1380 929 930 931 932 945 946 947 948 3338 3340 3342 3344 4154 4156 4158 4160 945 946 947 948 1377 1378 1379 1380 4154 4156 4158 4160 2673 2674 2675 2677 2673 2674 2675 2677 1401 1402 1403 1404 4242 4244 4246 4248 1009 1010 1011 1012 1401 1402 1403 1404 0 0 0 0 1009 1010 1011 1012 0 0 0 0 4242 4244 4246 4248 1009 1010 1011 1012 3846 3848 3850 3852 1169 1170 1171 1172 1009 1010 1011 1012 0 0 0 0 1169 1170 1171 1172 0 0 0 0 0 0 0 0 1041 1042 1043 1044 0 0 0 0 1073 1074 1075 1076 1041 1042 1043 1044 3570 3572 3574 3576 1073 1074 1075 1076 4290 4292 4294 4296 0 0 0 0 1073 1074 1075 1076 0 0 0 0 1425 1426 1427 1428 1073 1074 1075 1076 4290 4292 4294 4296 1425 1426 1427 1428 2773 2774 2775 2777 2773 2774 2775 2777 4384 4386 4388 4390 1461 1462 1463 1464 1137 1138 1139 1140 4384 4386 4388 4390 4018 4020 4022 4024 1137 1138 1139 1140 1249 1250 1251 1252 1461 1462 1463 1464 1137 1138 1139 1140 0 0 0 0 0 0 0 0 1137 1138 1139 1140 1249 1250 1251 1252 0 0 0 0 0 0 0 0 3846 3848 3850 3852 1169 1170 1171 1172 4454 4456 4458 4460 1201 1202 1203 1204 1169 1170 1171 1172 0 0 0 0 1201 1202 1203 1204 0 0 0 0 4454 4456 4458 4460 1201 1202 1203 1204 2873 2874 2875 2877 1489 1490 1491 1492 1201 1202 1203 1204 0 0 0 0 1489 1490 1491 1492 0 0 0 0 4018 4020 4022 4024 4498 4500 4502 4504 1249 1250 1251 1252 1265 1266 1267 1268 4498 4500 4502 4504 2873 2874 2875 2877 1265 1266 1267 1268 1505 1506 1507 1508 1249 1250 1251 1252 1265 1266 1267 1268 0 0 0 0 0 0 0 0 1265 1266 1267 1268 1505 1506 1507 1508 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1329 1330 1331 1332 0 0 0 0 0 0 0 0 1329 1330 1331 1332 1333 1334 1335 1336 0 0 0 0 1329 1330 1331 1332 0 0 0 0 1337 1338 1339 1340 1329 1330 1331 1332 1333 1334 1335 1336 1337 1338 1339 1340 2573 2574 2575 2577 0 0 0 0 0 0 0 0 1377 1378 1379 1380 1393 1394 1395 1396 0 0 0 0 0 0 0 0 1393 1394 1395 1396 0 0 0 0 1377 1378 1379 1380 1393 1394 1395 1396 2673 2674 2675 2677 1401 1402 1403 1404 1393 1394 1395 1396 0 0 0 0 1401 1402 1403 1404 0 0 0 0 0 0 0 0 1425 1426 1427 1428 0 0 0 0 1457 1458 1459 1460 1425 1426 1427 1428 2773 2774 2775 2777 1457 1458 1459 1460 1461 1462 1463 1464 0 0 0 0 1457 1458 1459 1460 0 0 0 0 0 0 0 0 1457 1458 1459 1460 1461 1462 1463 1464 0 0 0 0 0 0 0 0 2873 2874 2875 2877 1489 1490 1491 1492 1505 1506 1507 1508 1521 1522 1523 1524 1489 1490 1491 1492 0 0 0 0 1521 1522 1523 1524 0 0 0 0 1505 1506 1507 1508 1521 1522 1523 1524 0 0 0 0 0 0 0 0 1521 1522 1523 1524 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::Level 5 +DEAL:: 3213 3214 3215 3217 6446 6448 6451 6453 6608 6610 6613 6615 31336 31346 31356 31367 6446 6448 6451 6453 8811 8814 8817 8820 31336 31346 31356 31367 46902 46916 46930 46944 6608 6610 6613 6615 31336 31346 31356 31367 9651 9654 9657 9660 51742 51756 51770 51784 31336 31346 31356 31367 46902 46916 46930 46944 51742 51756 51770 51784 63888 63904 63920 63936 8811 8814 8817 8820 7132 7134 7137 7139 46902 46916 46930 46944 36905 36915 36925 36936 7132 7134 7137 7139 3313 3314 3315 3317 36905 36915 36925 36936 7587 7589 7592 7594 46902 46916 46930 46944 36905 36915 36925 36936 63888 63904 63920 63936 59698 59712 59726 59740 36905 36915 36925 36936 7587 7589 7592 7594 59698 59712 59726 59740 11275 11278 11281 11284 9651 9654 9657 9660 51742 51756 51770 51784 7796 7798 7801 7803 42194 42204 42214 42225 51742 51756 51770 51784 63888 63904 63920 63936 42194 42204 42214 42225 63134 63148 63162 63176 7796 7798 7801 7803 42194 42204 42214 42225 3413 3414 3415 3417 8260 8262 8265 8267 42194 42204 42214 42225 63134 63148 63162 63176 8260 8262 8265 8267 11835 11838 11841 11844 63888 63904 63920 63936 59698 59712 59726 59740 63134 63148 63162 63176 47879 47889 47899 47910 59698 59712 59726 59740 11275 11278 11281 11284 47879 47889 47899 47910 8775 8777 8780 8782 63134 63148 63162 63176 47879 47889 47899 47910 11835 11838 11841 11844 8946 8948 8951 8953 47879 47889 47899 47910 8775 8777 8780 8782 8946 8948 8951 8953 3513 3514 3515 3517 0 0 0 0 0 0 0 0 2101 2102 2103 2104 1585 1586 1587 1588 0 0 0 0 0 0 0 0 1585 1586 1587 1588 1697 1698 1699 1700 2101 2102 2103 2104 1585 1586 1587 1588 3213 3214 3215 3217 6446 6448 6451 6453 1585 1586 1587 1588 1697 1698 1699 1700 6446 6448 6451 6453 8811 8814 8817 8820 0 0 0 0 2105 2106 2107 2108 0 0 0 0 1649 1650 1651 1652 2105 2106 2107 2108 3213 3214 3215 3217 1649 1650 1651 1652 6608 6610 6613 6615 0 0 0 0 1649 1650 1651 1652 0 0 0 0 1809 1810 1811 1812 1649 1650 1651 1652 6608 6610 6613 6615 1809 1810 1811 1812 9651 9654 9657 9660 0 0 0 0 0 0 0 0 1697 1698 1699 1700 1713 1714 1715 1716 0 0 0 0 0 0 0 0 1713 1714 1715 1716 2145 2146 2147 2148 1697 1698 1699 1700 1713 1714 1715 1716 8811 8814 8817 8820 7132 7134 7137 7139 1713 1714 1715 1716 2145 2146 2147 2148 7132 7134 7137 7139 3313 3314 3315 3317 3313 3314 3315 3317 2169 2170 2171 2172 7587 7589 7592 7594 1777 1778 1779 1780 2169 2170 2171 2172 0 0 0 0 1777 1778 1779 1780 0 0 0 0 7587 7589 7592 7594 1777 1778 1779 1780 11275 11278 11281 11284 1937 1938 1939 1940 1777 1778 1779 1780 0 0 0 0 1937 1938 1939 1940 0 0 0 0 0 0 0 0 1809 1810 1811 1812 0 0 0 0 1841 1842 1843 1844 1809 1810 1811 1812 9651 9654 9657 9660 1841 1842 1843 1844 7796 7798 7801 7803 0 0 0 0 1841 1842 1843 1844 0 0 0 0 2193 2194 2195 2196 1841 1842 1843 1844 7796 7798 7801 7803 2193 2194 2195 2196 3413 3414 3415 3417 3413 3414 3415 3417 8260 8262 8265 8267 2229 2230 2231 2232 1905 1906 1907 1908 8260 8262 8265 8267 11835 11838 11841 11844 1905 1906 1907 1908 2017 2018 2019 2020 2229 2230 2231 2232 1905 1906 1907 1908 0 0 0 0 0 0 0 0 1905 1906 1907 1908 2017 2018 2019 2020 0 0 0 0 0 0 0 0 11275 11278 11281 11284 1937 1938 1939 1940 8775 8777 8780 8782 1969 1970 1971 1972 1937 1938 1939 1940 0 0 0 0 1969 1970 1971 1972 0 0 0 0 8775 8777 8780 8782 1969 1970 1971 1972 3513 3514 3515 3517 2257 2258 2259 2260 1969 1970 1971 1972 0 0 0 0 2257 2258 2259 2260 0 0 0 0 11835 11838 11841 11844 8946 8948 8951 8953 2017 2018 2019 2020 2033 2034 2035 2036 8946 8948 8951 8953 3513 3514 3515 3517 2033 2034 2035 2036 2273 2274 2275 2276 2017 2018 2019 2020 2033 2034 2035 2036 0 0 0 0 0 0 0 0 2033 2034 2035 2036 2273 2274 2275 2276 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2097 2098 2099 2100 0 0 0 0 0 0 0 0 2097 2098 2099 2100 2101 2102 2103 2104 0 0 0 0 2097 2098 2099 2100 0 0 0 0 2105 2106 2107 2108 2097 2098 2099 2100 2101 2102 2103 2104 2105 2106 2107 2108 3213 3214 3215 3217 0 0 0 0 0 0 0 0 2145 2146 2147 2148 2161 2162 2163 2164 0 0 0 0 0 0 0 0 2161 2162 2163 2164 0 0 0 0 2145 2146 2147 2148 2161 2162 2163 2164 3313 3314 3315 3317 2169 2170 2171 2172 2161 2162 2163 2164 0 0 0 0 2169 2170 2171 2172 0 0 0 0 0 0 0 0 2193 2194 2195 2196 0 0 0 0 2225 2226 2227 2228 2193 2194 2195 2196 3413 3414 3415 3417 2225 2226 2227 2228 2229 2230 2231 2232 0 0 0 0 2225 2226 2227 2228 0 0 0 0 0 0 0 0 2225 2226 2227 2228 2229 2230 2231 2232 0 0 0 0 0 0 0 0 3513 3514 3515 3517 2257 2258 2259 2260 2273 2274 2275 2276 2289 2290 2291 2292 2257 2258 2259 2260 0 0 0 0 2289 2290 2291 2292 0 0 0 0 2273 2274 2275 2276 2289 2290 2291 2292 0 0 0 0 0 0 0 0 2289 2290 2291 2292 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::Level 6 +DEAL:: 0 0 0 0 0 0 0 0 3533 3534 3535 3537 5786 5788 5790 5792 0 0 0 0 0 0 0 0 5786 5788 5790 5792 6586 6588 6590 6592 3533 3534 3535 3537 5786 5788 5790 5792 8414 8417 8420 8423 12324 12328 12332 12336 5786 5788 5790 5792 6586 6588 6590 6592 12324 12328 12332 12336 13444 13448 13452 13456 0 0 0 0 3693 3694 3695 3697 0 0 0 0 6034 6036 6038 6040 3693 3694 3695 3697 8414 8417 8420 8423 6034 6036 6038 6040 12596 12600 12604 12608 0 0 0 0 6034 6036 6038 6040 0 0 0 0 7506 7508 7510 7512 6034 6036 6038 6040 12596 12600 12604 12608 7506 7508 7510 7512 15060 15064 15068 15072 8414 8417 8420 8423 12324 12328 12332 12336 12596 12600 12604 12608 12964 12968 12972 12976 12324 12328 12332 12336 13444 13448 13452 13456 12964 12968 12972 12976 13988 13992 13996 14000 12596 12600 12604 12608 12964 12968 12972 12976 15060 15064 15068 15072 15396 15400 15404 15408 12964 12968 12972 12976 13988 13992 13996 14000 15396 15400 15404 15408 16420 16424 16428 16432 0 0 0 0 0 0 0 0 6586 6588 6590 6592 6810 6812 6814 6816 0 0 0 0 0 0 0 0 6810 6812 6814 6816 4273 4274 4275 4277 6586 6588 6590 6592 6810 6812 6814 6816 13444 13448 13452 13456 13988 13992 13996 14000 6810 6812 6814 6816 4273 4274 4275 4277 13988 13992 13996 14000 10046 10049 10052 10055 13444 13448 13452 13456 13988 13992 13996 14000 13988 13992 13996 14000 14500 14504 14508 14512 13988 13992 13996 14000 10046 10049 10052 10055 14500 14504 14508 14512 14996 15000 15004 15008 13988 13992 13996 14000 14500 14504 14508 14512 16420 16424 16428 16432 16932 16936 16940 16944 14500 14504 14508 14512 14996 15000 15004 15008 16932 16936 16940 16944 17524 17528 17532 17536 10046 10049 10052 10055 4673 4674 4675 4677 14996 15000 15004 15008 7590 7592 7594 7596 4673 4674 4675 4677 0 0 0 0 7590 7592 7594 7596 0 0 0 0 14996 15000 15004 15008 7590 7592 7594 7596 17524 17528 17532 17536 9062 9064 9066 9068 7590 7592 7594 7596 0 0 0 0 9062 9064 9066 9068 0 0 0 0 0 0 0 0 7506 7508 7510 7512 0 0 0 0 7826 7828 7830 7832 7506 7508 7510 7512 15060 15064 15068 15072 7826 7828 7830 7832 16180 16184 16188 16192 0 0 0 0 7826 7828 7830 7832 0 0 0 0 4933 4934 4935 4937 7826 7828 7830 7832 16180 16184 16188 16192 4933 4934 4935 4937 11502 11505 11508 11511 15060 15064 15068 15072 15396 15400 15404 15408 16180 16184 16188 16192 16548 16552 16556 16560 15396 15400 15404 15408 16420 16424 16428 16432 16548 16552 16556 16560 17572 17576 17580 17584 16180 16184 16188 16192 16548 16552 16556 16560 11502 11505 11508 11511 17060 17064 17068 17072 16548 16552 16556 16560 17572 17576 17580 17584 17060 17064 17068 17072 18372 18376 18380 18384 11502 11505 11508 11511 17060 17064 17068 17072 5333 5334 5335 5337 8642 8644 8646 8648 17060 17064 17068 17072 18372 18376 18380 18384 8642 8644 8646 8648 9442 9444 9446 9448 5333 5334 5335 5337 8642 8644 8646 8648 0 0 0 0 0 0 0 0 8642 8644 8646 8648 9442 9444 9446 9448 0 0 0 0 0 0 0 0 16420 16424 16428 16432 16932 16936 16940 16944 17572 17576 17580 17584 18084 18088 18092 18096 16932 16936 16940 16944 17524 17528 17532 17536 18084 18088 18092 18096 18580 18584 18588 18592 17572 17576 17580 17584 18084 18088 18092 18096 18372 18376 18380 18384 18980 18984 18988 18992 18084 18088 18092 18096 18580 18584 18588 18592 18980 18984 18988 18992 13122 13125 13128 13131 17524 17528 17532 17536 9062 9064 9066 9068 18580 18584 18588 18592 9382 9384 9386 9388 9062 9064 9066 9068 0 0 0 0 9382 9384 9386 9388 0 0 0 0 18580 18584 18588 18592 9382 9384 9386 9388 13122 13125 13128 13131 5913 5914 5915 5917 9382 9384 9386 9388 0 0 0 0 5913 5914 5915 5917 0 0 0 0 18372 18376 18380 18384 18980 18984 18988 18992 9442 9444 9446 9448 9666 9668 9670 9672 18980 18984 18988 18992 13122 13125 13128 13131 9666 9668 9670 9672 6073 6074 6075 6077 9442 9444 9446 9448 9666 9668 9670 9672 0 0 0 0 0 0 0 0 9666 9668 9670 9672 6073 6074 6075 6077 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2609 2610 2611 2612 0 0 0 0 0 0 0 0 2609 2610 2611 2612 3533 3534 3535 3537 0 0 0 0 2609 2610 2611 2612 0 0 0 0 3693 3694 3695 3697 2609 2610 2611 2612 3533 3534 3535 3537 3693 3694 3695 3697 8414 8417 8420 8423 0 0 0 0 0 0 0 0 4273 4274 4275 4277 2673 2674 2675 2676 0 0 0 0 0 0 0 0 2673 2674 2675 2676 0 0 0 0 4273 4274 4275 4277 2673 2674 2675 2676 10046 10049 10052 10055 4673 4674 4675 4677 2673 2674 2675 2676 0 0 0 0 4673 4674 4675 4677 0 0 0 0 0 0 0 0 4933 4934 4935 4937 0 0 0 0 2737 2738 2739 2740 4933 4934 4935 4937 11502 11505 11508 11511 2737 2738 2739 2740 5333 5334 5335 5337 0 0 0 0 2737 2738 2739 2740 0 0 0 0 0 0 0 0 2737 2738 2739 2740 5333 5334 5335 5337 0 0 0 0 0 0 0 0 13122 13125 13128 13131 5913 5914 5915 5917 6073 6074 6075 6077 2801 2802 2803 2804 5913 5914 5915 5917 0 0 0 0 2801 2802 2803 2804 0 0 0 0 6073 6074 6075 6077 2801 2802 2803 2804 0 0 0 0 0 0 0 0 2801 2802 2803 2804 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::Level 7 +DEAL:: 0 0 0 0 0 0 0 0 0 0 0 0 2865 2866 2867 2868 0 0 0 0 0 0 0 0 2865 2866 2867 2868 2913 2914 2915 2916 0 0 0 0 2865 2866 2867 2868 0 0 0 0 3089 3090 3091 3092 2865 2866 2867 2868 2913 2914 2915 2916 3089 3090 3091 3092 3137 3138 3139 3140 0 0 0 0 0 0 0 0 2913 2914 2915 2916 2929 2930 2931 2932 0 0 0 0 0 0 0 0 2929 2930 2931 2932 3361 3362 3363 3364 2913 2914 2915 2916 2929 2930 2931 2932 3137 3138 3139 3140 3153 3154 3155 3156 2929 2930 2931 2932 3361 3362 3363 3364 3153 3154 3155 3156 3457 3458 3459 3460 0 0 0 0 0 0 0 0 0 0 0 0 2993 2994 2995 2996 0 0 0 0 0 0 0 0 2993 2994 2995 2996 3105 3106 3107 3108 0 0 0 0 2993 2994 2995 2996 0 0 0 0 3025 3026 3027 3028 2993 2994 2995 2996 3105 3106 3107 3108 3025 3026 3027 3028 3201 3202 3203 3204 0 0 0 0 3025 3026 3027 3028 0 0 0 0 3057 3058 3059 3060 3025 3026 3027 3028 3201 3202 3203 3204 3057 3058 3059 3060 3233 3234 3235 3236 0 0 0 0 3057 3058 3059 3060 0 0 0 0 3857 3858 3859 3860 3057 3058 3059 3060 3233 3234 3235 3236 3857 3858 3859 3860 3969 3970 3971 3972 0 0 0 0 3089 3090 3091 3092 3105 3106 3107 3108 3121 3122 3123 3124 3089 3090 3091 3092 3137 3138 3139 3140 3121 3122 3123 3124 3169 3170 3171 3172 3105 3106 3107 3108 3121 3122 3123 3124 3201 3202 3203 3204 3217 3218 3219 3220 3121 3122 3123 3124 3169 3170 3171 3172 3217 3218 3219 3220 3265 3266 3267 3268 3137 3138 3139 3140 3153 3154 3155 3156 3169 3170 3171 3172 3185 3186 3187 3188 3153 3154 3155 3156 3457 3458 3459 3460 3185 3186 3187 3188 3489 3490 3491 3492 3169 3170 3171 3172 3185 3186 3187 3188 3265 3266 3267 3268 3281 3282 3283 3284 3185 3186 3187 3188 3489 3490 3491 3492 3281 3282 3283 3284 3585 3586 3587 3588 3201 3202 3203 3204 3217 3218 3219 3220 3233 3234 3235 3236 3249 3250 3251 3252 3217 3218 3219 3220 3265 3266 3267 3268 3249 3250 3251 3252 3297 3298 3299 3300 3233 3234 3235 3236 3249 3250 3251 3252 3969 3970 3971 3972 3985 3986 3987 3988 3249 3250 3251 3252 3297 3298 3299 3300 3985 3986 3987 3988 4033 4034 4035 4036 3265 3266 3267 3268 3281 3282 3283 3284 3297 3298 3299 3300 3313 3314 3315 3316 3281 3282 3283 3284 3585 3586 3587 3588 3313 3314 3315 3316 3617 3618 3619 3620 3297 3298 3299 3300 3313 3314 3315 3316 4033 4034 4035 4036 4049 4050 4051 4052 3313 3314 3315 3316 3617 3618 3619 3620 4049 4050 4051 4052 4353 4354 4355 4356 0 0 0 0 0 0 0 0 3361 3362 3363 3364 3377 3378 3379 3380 0 0 0 0 0 0 0 0 3377 3378 3379 3380 3425 3426 3427 3428 3361 3362 3363 3364 3377 3378 3379 3380 3457 3458 3459 3460 3473 3474 3475 3476 3377 3378 3379 3380 3425 3426 3427 3428 3473 3474 3475 3476 3521 3522 3523 3524 0 0 0 0 0 0 0 0 3425 3426 3427 3428 3441 3442 3443 3444 0 0 0 0 0 0 0 0 3441 3442 3443 3444 0 0 0 0 3425 3426 3427 3428 3441 3442 3443 3444 3521 3522 3523 3524 3537 3538 3539 3540 3441 3442 3443 3444 0 0 0 0 3537 3538 3539 3540 0 0 0 0 3457 3458 3459 3460 3473 3474 3475 3476 3489 3490 3491 3492 3505 3506 3507 3508 3473 3474 3475 3476 3521 3522 3523 3524 3505 3506 3507 3508 3553 3554 3555 3556 3489 3490 3491 3492 3505 3506 3507 3508 3585 3586 3587 3588 3601 3602 3603 3604 3505 3506 3507 3508 3553 3554 3555 3556 3601 3602 3603 3604 3649 3650 3651 3652 3521 3522 3523 3524 3537 3538 3539 3540 3553 3554 3555 3556 3569 3570 3571 3572 3537 3538 3539 3540 0 0 0 0 3569 3570 3571 3572 3745 3746 3747 3748 3553 3554 3555 3556 3569 3570 3571 3572 3649 3650 3651 3652 3665 3666 3667 3668 3569 3570 3571 3572 3745 3746 3747 3748 3665 3666 3667 3668 3777 3778 3779 3780 3585 3586 3587 3588 3601 3602 3603 3604 3617 3618 3619 3620 3633 3634 3635 3636 3601 3602 3603 3604 3649 3650 3651 3652 3633 3634 3635 3636 3681 3682 3683 3684 3617 3618 3619 3620 3633 3634 3635 3636 4353 4354 4355 4356 4369 4370 4371 4372 3633 3634 3635 3636 3681 3682 3683 3684 4369 4370 4371 4372 4417 4418 4419 4420 3649 3650 3651 3652 3665 3666 3667 3668 3681 3682 3683 3684 3697 3698 3699 3700 3665 3666 3667 3668 3777 3778 3779 3780 3697 3698 3699 3700 3809 3810 3811 3812 3681 3682 3683 3684 3697 3698 3699 3700 4417 4418 4419 4420 4433 4434 4435 4436 3697 3698 3699 3700 3809 3810 3811 3812 4433 4434 4435 4436 4609 4610 4611 4612 0 0 0 0 0 0 0 0 3745 3746 3747 3748 3761 3762 3763 3764 0 0 0 0 0 0 0 0 3761 3762 3763 3764 0 0 0 0 3745 3746 3747 3748 3761 3762 3763 3764 3777 3778 3779 3780 3793 3794 3795 3796 3761 3762 3763 3764 0 0 0 0 3793 3794 3795 3796 0 0 0 0 3777 3778 3779 3780 3793 3794 3795 3796 3809 3810 3811 3812 3825 3826 3827 3828 3793 3794 3795 3796 0 0 0 0 3825 3826 3827 3828 0 0 0 0 3809 3810 3811 3812 3825 3826 3827 3828 4609 4610 4611 4612 4625 4626 4627 4628 3825 3826 3827 3828 0 0 0 0 4625 4626 4627 4628 0 0 0 0 0 0 0 0 3857 3858 3859 3860 0 0 0 0 3889 3890 3891 3892 3857 3858 3859 3860 3969 3970 3971 3972 3889 3890 3891 3892 4001 4002 4003 4004 0 0 0 0 3889 3890 3891 3892 0 0 0 0 3921 3922 3923 3924 3889 3890 3891 3892 4001 4002 4003 4004 3921 3922 3923 3924 4097 4098 4099 4100 0 0 0 0 3921 3922 3923 3924 0 0 0 0 3953 3954 3955 3956 3921 3922 3923 3924 4097 4098 4099 4100 3953 3954 3955 3956 4129 4130 4131 4132 0 0 0 0 3953 3954 3955 3956 0 0 0 0 0 0 0 0 3953 3954 3955 3956 4129 4130 4131 4132 0 0 0 0 0 0 0 0 3969 3970 3971 3972 3985 3986 3987 3988 4001 4002 4003 4004 4017 4018 4019 4020 3985 3986 3987 3988 4033 4034 4035 4036 4017 4018 4019 4020 4065 4066 4067 4068 4001 4002 4003 4004 4017 4018 4019 4020 4097 4098 4099 4100 4113 4114 4115 4116 4017 4018 4019 4020 4065 4066 4067 4068 4113 4114 4115 4116 4161 4162 4163 4164 4033 4034 4035 4036 4049 4050 4051 4052 4065 4066 4067 4068 4081 4082 4083 4084 4049 4050 4051 4052 4353 4354 4355 4356 4081 4082 4083 4084 4385 4386 4387 4388 4065 4066 4067 4068 4081 4082 4083 4084 4161 4162 4163 4164 4177 4178 4179 4180 4081 4082 4083 4084 4385 4386 4387 4388 4177 4178 4179 4180 4481 4482 4483 4484 4097 4098 4099 4100 4113 4114 4115 4116 4129 4130 4131 4132 4145 4146 4147 4148 4113 4114 4115 4116 4161 4162 4163 4164 4145 4146 4147 4148 4193 4194 4195 4196 4129 4130 4131 4132 4145 4146 4147 4148 0 0 0 0 4241 4242 4243 4244 4145 4146 4147 4148 4193 4194 4195 4196 4241 4242 4243 4244 4289 4290 4291 4292 4161 4162 4163 4164 4177 4178 4179 4180 4193 4194 4195 4196 4209 4210 4211 4212 4177 4178 4179 4180 4481 4482 4483 4484 4209 4210 4211 4212 4513 4514 4515 4516 4193 4194 4195 4196 4209 4210 4211 4212 4289 4290 4291 4292 4305 4306 4307 4308 4209 4210 4211 4212 4513 4514 4515 4516 4305 4306 4307 4308 4737 4738 4739 4740 0 0 0 0 4241 4242 4243 4244 0 0 0 0 4273 4274 4275 4276 4241 4242 4243 4244 4289 4290 4291 4292 4273 4274 4275 4276 4321 4322 4323 4324 0 0 0 0 4273 4274 4275 4276 0 0 0 0 0 0 0 0 4273 4274 4275 4276 4321 4322 4323 4324 0 0 0 0 0 0 0 0 4289 4290 4291 4292 4305 4306 4307 4308 4321 4322 4323 4324 4337 4338 4339 4340 4305 4306 4307 4308 4737 4738 4739 4740 4337 4338 4339 4340 4769 4770 4771 4772 4321 4322 4323 4324 4337 4338 4339 4340 0 0 0 0 0 0 0 0 4337 4338 4339 4340 4769 4770 4771 4772 0 0 0 0 0 0 0 0 4353 4354 4355 4356 4369 4370 4371 4372 4385 4386 4387 4388 4401 4402 4403 4404 4369 4370 4371 4372 4417 4418 4419 4420 4401 4402 4403 4404 4449 4450 4451 4452 4385 4386 4387 4388 4401 4402 4403 4404 4481 4482 4483 4484 4497 4498 4499 4500 4401 4402 4403 4404 4449 4450 4451 4452 4497 4498 4499 4500 4545 4546 4547 4548 4417 4418 4419 4420 4433 4434 4435 4436 4449 4450 4451 4452 4465 4466 4467 4468 4433 4434 4435 4436 4609 4610 4611 4612 4465 4466 4467 4468 4641 4642 4643 4644 4449 4450 4451 4452 4465 4466 4467 4468 4545 4546 4547 4548 4561 4562 4563 4564 4465 4466 4467 4468 4641 4642 4643 4644 4561 4562 4563 4564 4673 4674 4675 4676 4481 4482 4483 4484 4497 4498 4499 4500 4513 4514 4515 4516 4529 4530 4531 4532 4497 4498 4499 4500 4545 4546 4547 4548 4529 4530 4531 4532 4577 4578 4579 4580 4513 4514 4515 4516 4529 4530 4531 4532 4737 4738 4739 4740 4753 4754 4755 4756 4529 4530 4531 4532 4577 4578 4579 4580 4753 4754 4755 4756 4801 4802 4803 4804 4545 4546 4547 4548 4561 4562 4563 4564 4577 4578 4579 4580 4593 4594 4595 4596 4561 4562 4563 4564 4673 4674 4675 4676 4593 4594 4595 4596 4705 4706 4707 4708 4577 4578 4579 4580 4593 4594 4595 4596 4801 4802 4803 4804 4817 4818 4819 4820 4593 4594 4595 4596 4705 4706 4707 4708 4817 4818 4819 4820 0 0 0 0 4609 4610 4611 4612 4625 4626 4627 4628 4641 4642 4643 4644 4657 4658 4659 4660 4625 4626 4627 4628 0 0 0 0 4657 4658 4659 4660 0 0 0 0 4641 4642 4643 4644 4657 4658 4659 4660 4673 4674 4675 4676 4689 4690 4691 4692 4657 4658 4659 4660 0 0 0 0 4689 4690 4691 4692 0 0 0 0 4673 4674 4675 4676 4689 4690 4691 4692 4705 4706 4707 4708 4721 4722 4723 4724 4689 4690 4691 4692 0 0 0 0 4721 4722 4723 4724 0 0 0 0 4705 4706 4707 4708 4721 4722 4723 4724 0 0 0 0 0 0 0 0 4721 4722 4723 4724 0 0 0 0 0 0 0 0 0 0 0 0 4737 4738 4739 4740 4753 4754 4755 4756 4769 4770 4771 4772 4785 4786 4787 4788 4753 4754 4755 4756 4801 4802 4803 4804 4785 4786 4787 4788 4833 4834 4835 4836 4769 4770 4771 4772 4785 4786 4787 4788 0 0 0 0 0 0 0 0 4785 4786 4787 4788 4833 4834 4835 4836 0 0 0 0 0 0 0 0 4801 4802 4803 4804 4817 4818 4819 4820 4833 4834 4835 4836 4849 4850 4851 4852 4817 4818 4819 4820 0 0 0 0 4849 4850 4851 4852 0 0 0 0 4833 4834 4835 4836 4849 4850 4851 4852 0 0 0 0 0 0 0 0 4849 4850 4851 4852 0 0 0 0 0 0 0 0 0 0 0 0copy_to_mg +DEAL:: +DEAL::Level 0 +DEAL:: 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::Level 1 +DEAL:: 0 0 0 0 0 0 0 0 0 0 0 0 589199 589384 589568 589753 0 0 0 0 0 0 0 0 589199 589384 589568 589753 0 0 0 0 0 0 0 0 589199 589384 589568 589753 0 0 0 0 0 0 0 0 589199 589384 589568 589753 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::Level 2 +DEAL:: 0 0 0 0 0 0 0 0 0 0 0 0 8378 8384 8391 8398 0 0 0 0 0 0 0 0 8378 8384 8391 8398 35911 35928 35946 35964 0 0 0 0 8378 8384 8391 8398 0 0 0 0 38331 38348 38366 38383 8378 8384 8391 8398 35911 35928 35946 35964 38331 38348 38366 38383 499066 499208 499351 499493 0 0 0 0 0 0 0 0 35911 35928 35946 35964 8998 9004 9011 9018 0 0 0 0 0 0 0 0 8998 9004 9011 9018 0 0 0 0 35911 35928 35946 35964 8998 9004 9011 9018 499066 499208 499351 499493 42677 42695 42712 42730 8998 9004 9011 9018 0 0 0 0 42677 42695 42712 42730 0 0 0 0 0 0 0 0 38331 38348 38366 38383 0 0 0 0 9603 9610 9616 9623 38331 38348 38366 38383 499066 499208 499351 499493 9603 9610 9616 9623 44745 44762 44780 44797 0 0 0 0 9603 9610 9616 9623 0 0 0 0 0 0 0 0 9603 9610 9616 9623 44745 44762 44780 44797 0 0 0 0 0 0 0 0 499066 499208 499351 499493 42677 42695 42712 42730 44745 44762 44780 44797 10225 10231 10238 10245 42677 42695 42712 42730 0 0 0 0 10225 10231 10238 10245 0 0 0 0 44745 44762 44780 44797 10225 10231 10238 10245 0 0 0 0 0 0 0 0 10225 10231 10238 10245 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::Level 3 +DEAL:: 1613 1614 1615 1617 2504 2506 2508 2510 2542 2544 2546 2548 13976 13983 13990 13997 2504 2506 2508 2510 1802 1804 1806 1808 13976 13983 13990 13997 48352 48370 48388 48406 2542 2544 2546 2548 13976 13983 13990 13997 2034 2036 2038 2040 51881 51899 51917 51935 13976 13983 13990 13997 48352 48370 48388 48406 51881 51899 51917 51935 373589 373688 373787 373887 1802 1804 1806 1808 2618 2620 2622 2624 48352 48370 48388 48406 14964 14971 14978 14985 2618 2620 2622 2624 1713 1714 1715 1717 14964 14971 14978 14985 2706 2708 2710 2712 48352 48370 48388 48406 14964 14971 14978 14985 373589 373688 373787 373887 58446 58465 58483 58501 14964 14971 14978 14985 2706 2708 2710 2712 58446 58465 58483 58501 2310 2312 2314 2316 2034 2036 2038 2040 51881 51899 51917 51935 2754 2756 2758 2760 15893 15900 15907 15914 51881 51899 51917 51935 373589 373688 373787 373887 15893 15900 15907 15914 61411 61429 61447 61465 2754 2756 2758 2760 15893 15900 15907 15914 1813 1814 1815 1817 2848 2850 2852 2854 15893 15900 15907 15914 61411 61429 61447 61465 2848 2850 2852 2854 2482 2484 2486 2488 373589 373688 373787 373887 58446 58465 58483 58501 61411 61429 61447 61465 16888 16895 16902 16909 58446 58465 58483 58501 2310 2312 2314 2316 16888 16895 16902 16909 2918 2920 2922 2924 61411 61429 61447 61465 16888 16895 16902 16909 2482 2484 2486 2488 2962 2964 2966 2968 16888 16895 16902 16909 2918 2920 2922 2924 2962 2964 2966 2968 1913 1914 1915 1917 533 534 535 536 17 18 19 20 565 566 567 568 49 50 51 52 17 18 19 20 129 130 131 132 49 50 51 52 161 162 163 164 565 566 567 568 49 50 51 52 1613 1614 1615 1617 2504 2506 2508 2510 49 50 51 52 161 162 163 164 2504 2506 2508 2510 1802 1804 1806 1808 553 554 555 556 569 570 571 572 97 98 99 100 113 114 115 116 569 570 571 572 1613 1614 1615 1617 113 114 115 116 2542 2544 2546 2548 97 98 99 100 113 114 115 116 257 258 259 260 273 274 275 276 113 114 115 116 2542 2544 2546 2548 273 274 275 276 2034 2036 2038 2040 129 130 131 132 145 146 147 148 161 162 163 164 177 178 179 180 145 146 147 148 577 578 579 580 177 178 179 180 609 610 611 612 161 162 163 164 177 178 179 180 1802 1804 1806 1808 2618 2620 2622 2624 177 178 179 180 609 610 611 612 2618 2620 2622 2624 1713 1714 1715 1717 1713 1714 1715 1717 633 634 635 636 2706 2708 2710 2712 241 242 243 244 633 634 635 636 637 638 639 640 241 242 243 244 245 246 247 248 2706 2708 2710 2712 241 242 243 244 2310 2312 2314 2316 401 402 403 404 241 242 243 244 245 246 247 248 401 402 403 404 405 406 407 408 257 258 259 260 273 274 275 276 289 290 291 292 305 306 307 308 273 274 275 276 2034 2036 2038 2040 305 306 307 308 2754 2756 2758 2760 289 290 291 292 305 306 307 308 641 642 643 644 657 658 659 660 305 306 307 308 2754 2756 2758 2760 657 658 659 660 1813 1814 1815 1817 1813 1814 1815 1817 2848 2850 2852 2854 693 694 695 696 369 370 371 372 2848 2850 2852 2854 2482 2484 2486 2488 369 370 371 372 481 482 483 484 693 694 695 696 369 370 371 372 701 702 703 704 377 378 379 380 369 370 371 372 481 482 483 484 377 378 379 380 489 490 491 492 2310 2312 2314 2316 401 402 403 404 2918 2920 2922 2924 433 434 435 436 401 402 403 404 405 406 407 408 433 434 435 436 437 438 439 440 2918 2920 2922 2924 433 434 435 436 1913 1914 1915 1917 721 722 723 724 433 434 435 436 437 438 439 440 721 722 723 724 725 726 727 728 2482 2484 2486 2488 2962 2964 2966 2968 481 482 483 484 497 498 499 500 2962 2964 2966 2968 1913 1914 1915 1917 497 498 499 500 737 738 739 740 481 482 483 484 497 498 499 500 489 490 491 492 505 506 507 508 497 498 499 500 737 738 739 740 505 506 507 508 745 746 747 748 513 514 515 516 529 530 531 532 545 546 547 548 561 562 563 564 529 530 531 532 533 534 535 536 561 562 563 564 565 566 567 568 545 546 547 548 561 562 563 564 553 554 555 556 569 570 571 572 561 562 563 564 565 566 567 568 569 570 571 572 1613 1614 1615 1617 577 578 579 580 593 594 595 596 609 610 611 612 625 626 627 628 593 594 595 596 597 598 599 600 625 626 627 628 629 630 631 632 609 610 611 612 625 626 627 628 1713 1714 1715 1717 633 634 635 636 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 657 658 659 660 673 674 675 676 689 690 691 692 657 658 659 660 1813 1814 1815 1817 689 690 691 692 693 694 695 696 673 674 675 676 689 690 691 692 681 682 683 684 697 698 699 700 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 1913 1914 1915 1917 721 722 723 724 737 738 739 740 753 754 755 756 721 722 723 724 725 726 727 728 753 754 755 756 757 758 759 760 737 738 739 740 753 754 755 756 745 746 747 748 761 762 763 764 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 +DEAL::Level 4 +DEAL:: 2573 2574 2575 2577 4040 4042 4044 4046 4078 4080 4082 4084 21010 21018 21026 21034 4040 4042 4044 4046 3338 3340 3342 3344 21010 21018 21026 21034 57784 57803 57821 57840 4078 4080 4082 4084 21010 21018 21026 21034 3570 3572 3574 3576 62884 62902 62921 62939 21010 21018 21026 21034 57784 57803 57821 57840 62884 62902 62921 62939 214204 214259 214313 214367 3338 3340 3342 3344 4154 4156 4158 4160 57784 57803 57821 57840 23469 23477 23485 23493 4154 4156 4158 4160 2673 2674 2675 2677 23469 23477 23485 23493 4242 4244 4246 4248 57784 57803 57821 57840 23469 23477 23485 23493 214204 214259 214313 214367 72406 72425 72443 72462 23469 23477 23485 23493 4242 4244 4246 4248 72406 72425 72443 72462 3846 3848 3850 3852 3570 3572 3574 3576 62884 62902 62921 62939 4290 4292 4294 4296 25693 25701 25709 25717 62884 62902 62921 62939 214204 214259 214313 214367 25693 25701 25709 25717 76516 76535 76553 76572 4290 4292 4294 4296 25693 25701 25709 25717 2773 2774 2775 2777 4384 4386 4388 4390 25693 25701 25709 25717 76516 76535 76553 76572 4384 4386 4388 4390 4018 4020 4022 4024 214204 214259 214313 214367 72406 72425 72443 72462 76516 76535 76553 76572 28181 28189 28197 28205 72406 72425 72443 72462 3846 3848 3850 3852 28181 28189 28197 28205 4454 4456 4458 4460 76516 76535 76553 76572 28181 28189 28197 28205 4018 4020 4022 4024 4498 4500 4502 4504 28181 28189 28197 28205 4454 4456 4458 4460 4498 4500 4502 4504 2873 2874 2875 2877 0 0 0 0 0 0 0 0 1333 1334 1335 1336 817 818 819 820 0 0 0 0 0 0 0 0 817 818 819 820 929 930 931 932 1333 1334 1335 1336 817 818 819 820 2573 2574 2575 2577 4040 4042 4044 4046 817 818 819 820 929 930 931 932 4040 4042 4044 4046 3338 3340 3342 3344 0 0 0 0 1337 1338 1339 1340 0 0 0 0 881 882 883 884 1337 1338 1339 1340 2573 2574 2575 2577 881 882 883 884 4078 4080 4082 4084 0 0 0 0 881 882 883 884 0 0 0 0 1041 1042 1043 1044 881 882 883 884 4078 4080 4082 4084 1041 1042 1043 1044 3570 3572 3574 3576 0 0 0 0 0 0 0 0 929 930 931 932 945 946 947 948 0 0 0 0 0 0 0 0 945 946 947 948 1377 1378 1379 1380 929 930 931 932 945 946 947 948 3338 3340 3342 3344 4154 4156 4158 4160 945 946 947 948 1377 1378 1379 1380 4154 4156 4158 4160 2673 2674 2675 2677 2673 2674 2675 2677 1401 1402 1403 1404 4242 4244 4246 4248 1009 1010 1011 1012 1401 1402 1403 1404 0 0 0 0 1009 1010 1011 1012 0 0 0 0 4242 4244 4246 4248 1009 1010 1011 1012 3846 3848 3850 3852 1169 1170 1171 1172 1009 1010 1011 1012 0 0 0 0 1169 1170 1171 1172 0 0 0 0 0 0 0 0 1041 1042 1043 1044 0 0 0 0 1073 1074 1075 1076 1041 1042 1043 1044 3570 3572 3574 3576 1073 1074 1075 1076 4290 4292 4294 4296 0 0 0 0 1073 1074 1075 1076 0 0 0 0 1425 1426 1427 1428 1073 1074 1075 1076 4290 4292 4294 4296 1425 1426 1427 1428 2773 2774 2775 2777 2773 2774 2775 2777 4384 4386 4388 4390 1461 1462 1463 1464 1137 1138 1139 1140 4384 4386 4388 4390 4018 4020 4022 4024 1137 1138 1139 1140 1249 1250 1251 1252 1461 1462 1463 1464 1137 1138 1139 1140 0 0 0 0 0 0 0 0 1137 1138 1139 1140 1249 1250 1251 1252 0 0 0 0 0 0 0 0 3846 3848 3850 3852 1169 1170 1171 1172 4454 4456 4458 4460 1201 1202 1203 1204 1169 1170 1171 1172 0 0 0 0 1201 1202 1203 1204 0 0 0 0 4454 4456 4458 4460 1201 1202 1203 1204 2873 2874 2875 2877 1489 1490 1491 1492 1201 1202 1203 1204 0 0 0 0 1489 1490 1491 1492 0 0 0 0 4018 4020 4022 4024 4498 4500 4502 4504 1249 1250 1251 1252 1265 1266 1267 1268 4498 4500 4502 4504 2873 2874 2875 2877 1265 1266 1267 1268 1505 1506 1507 1508 1249 1250 1251 1252 1265 1266 1267 1268 0 0 0 0 0 0 0 0 1265 1266 1267 1268 1505 1506 1507 1508 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1329 1330 1331 1332 0 0 0 0 0 0 0 0 1329 1330 1331 1332 1333 1334 1335 1336 0 0 0 0 1329 1330 1331 1332 0 0 0 0 1337 1338 1339 1340 1329 1330 1331 1332 1333 1334 1335 1336 1337 1338 1339 1340 2573 2574 2575 2577 0 0 0 0 0 0 0 0 1377 1378 1379 1380 1393 1394 1395 1396 0 0 0 0 0 0 0 0 1393 1394 1395 1396 0 0 0 0 1377 1378 1379 1380 1393 1394 1395 1396 2673 2674 2675 2677 1401 1402 1403 1404 1393 1394 1395 1396 0 0 0 0 1401 1402 1403 1404 0 0 0 0 0 0 0 0 1425 1426 1427 1428 0 0 0 0 1457 1458 1459 1460 1425 1426 1427 1428 2773 2774 2775 2777 1457 1458 1459 1460 1461 1462 1463 1464 0 0 0 0 1457 1458 1459 1460 0 0 0 0 0 0 0 0 1457 1458 1459 1460 1461 1462 1463 1464 0 0 0 0 0 0 0 0 2873 2874 2875 2877 1489 1490 1491 1492 1505 1506 1507 1508 1521 1522 1523 1524 1489 1490 1491 1492 0 0 0 0 1521 1522 1523 1524 0 0 0 0 1505 1506 1507 1508 1521 1522 1523 1524 0 0 0 0 0 0 0 0 1521 1522 1523 1524 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::Level 5 +DEAL:: 3213 3214 3215 3217 6446 6448 6451 6453 6608 6610 6613 6615 31336 31346 31356 31367 6446 6448 6451 6453 8811 8814 8817 8820 31336 31346 31356 31367 46902 46916 46930 46944 6608 6610 6613 6615 31336 31346 31356 31367 9651 9654 9657 9660 51742 51756 51770 51784 31336 31346 31356 31367 46902 46916 46930 46944 51742 51756 51770 51784 63888 63904 63920 63936 8811 8814 8817 8820 7132 7134 7137 7139 46902 46916 46930 46944 36905 36915 36925 36936 7132 7134 7137 7139 3313 3314 3315 3317 36905 36915 36925 36936 7587 7589 7592 7594 46902 46916 46930 46944 36905 36915 36925 36936 63888 63904 63920 63936 59698 59712 59726 59740 36905 36915 36925 36936 7587 7589 7592 7594 59698 59712 59726 59740 11275 11278 11281 11284 9651 9654 9657 9660 51742 51756 51770 51784 7796 7798 7801 7803 42194 42204 42214 42225 51742 51756 51770 51784 63888 63904 63920 63936 42194 42204 42214 42225 63134 63148 63162 63176 7796 7798 7801 7803 42194 42204 42214 42225 3413 3414 3415 3417 8260 8262 8265 8267 42194 42204 42214 42225 63134 63148 63162 63176 8260 8262 8265 8267 11835 11838 11841 11844 63888 63904 63920 63936 59698 59712 59726 59740 63134 63148 63162 63176 47879 47889 47899 47910 59698 59712 59726 59740 11275 11278 11281 11284 47879 47889 47899 47910 8775 8777 8780 8782 63134 63148 63162 63176 47879 47889 47899 47910 11835 11838 11841 11844 8946 8948 8951 8953 47879 47889 47899 47910 8775 8777 8780 8782 8946 8948 8951 8953 3513 3514 3515 3517 0 0 0 0 0 0 0 0 2101 2102 2103 2104 1585 1586 1587 1588 0 0 0 0 0 0 0 0 1585 1586 1587 1588 1697 1698 1699 1700 2101 2102 2103 2104 1585 1586 1587 1588 3213 3214 3215 3217 6446 6448 6451 6453 1585 1586 1587 1588 1697 1698 1699 1700 6446 6448 6451 6453 8811 8814 8817 8820 0 0 0 0 2105 2106 2107 2108 0 0 0 0 1649 1650 1651 1652 2105 2106 2107 2108 3213 3214 3215 3217 1649 1650 1651 1652 6608 6610 6613 6615 0 0 0 0 1649 1650 1651 1652 0 0 0 0 1809 1810 1811 1812 1649 1650 1651 1652 6608 6610 6613 6615 1809 1810 1811 1812 9651 9654 9657 9660 0 0 0 0 0 0 0 0 1697 1698 1699 1700 1713 1714 1715 1716 0 0 0 0 0 0 0 0 1713 1714 1715 1716 2145 2146 2147 2148 1697 1698 1699 1700 1713 1714 1715 1716 8811 8814 8817 8820 7132 7134 7137 7139 1713 1714 1715 1716 2145 2146 2147 2148 7132 7134 7137 7139 3313 3314 3315 3317 3313 3314 3315 3317 2169 2170 2171 2172 7587 7589 7592 7594 1777 1778 1779 1780 2169 2170 2171 2172 0 0 0 0 1777 1778 1779 1780 0 0 0 0 7587 7589 7592 7594 1777 1778 1779 1780 11275 11278 11281 11284 1937 1938 1939 1940 1777 1778 1779 1780 0 0 0 0 1937 1938 1939 1940 0 0 0 0 0 0 0 0 1809 1810 1811 1812 0 0 0 0 1841 1842 1843 1844 1809 1810 1811 1812 9651 9654 9657 9660 1841 1842 1843 1844 7796 7798 7801 7803 0 0 0 0 1841 1842 1843 1844 0 0 0 0 2193 2194 2195 2196 1841 1842 1843 1844 7796 7798 7801 7803 2193 2194 2195 2196 3413 3414 3415 3417 3413 3414 3415 3417 8260 8262 8265 8267 2229 2230 2231 2232 1905 1906 1907 1908 8260 8262 8265 8267 11835 11838 11841 11844 1905 1906 1907 1908 2017 2018 2019 2020 2229 2230 2231 2232 1905 1906 1907 1908 0 0 0 0 0 0 0 0 1905 1906 1907 1908 2017 2018 2019 2020 0 0 0 0 0 0 0 0 11275 11278 11281 11284 1937 1938 1939 1940 8775 8777 8780 8782 1969 1970 1971 1972 1937 1938 1939 1940 0 0 0 0 1969 1970 1971 1972 0 0 0 0 8775 8777 8780 8782 1969 1970 1971 1972 3513 3514 3515 3517 2257 2258 2259 2260 1969 1970 1971 1972 0 0 0 0 2257 2258 2259 2260 0 0 0 0 11835 11838 11841 11844 8946 8948 8951 8953 2017 2018 2019 2020 2033 2034 2035 2036 8946 8948 8951 8953 3513 3514 3515 3517 2033 2034 2035 2036 2273 2274 2275 2276 2017 2018 2019 2020 2033 2034 2035 2036 0 0 0 0 0 0 0 0 2033 2034 2035 2036 2273 2274 2275 2276 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2097 2098 2099 2100 0 0 0 0 0 0 0 0 2097 2098 2099 2100 2101 2102 2103 2104 0 0 0 0 2097 2098 2099 2100 0 0 0 0 2105 2106 2107 2108 2097 2098 2099 2100 2101 2102 2103 2104 2105 2106 2107 2108 3213 3214 3215 3217 0 0 0 0 0 0 0 0 2145 2146 2147 2148 2161 2162 2163 2164 0 0 0 0 0 0 0 0 2161 2162 2163 2164 0 0 0 0 2145 2146 2147 2148 2161 2162 2163 2164 3313 3314 3315 3317 2169 2170 2171 2172 2161 2162 2163 2164 0 0 0 0 2169 2170 2171 2172 0 0 0 0 0 0 0 0 2193 2194 2195 2196 0 0 0 0 2225 2226 2227 2228 2193 2194 2195 2196 3413 3414 3415 3417 2225 2226 2227 2228 2229 2230 2231 2232 0 0 0 0 2225 2226 2227 2228 0 0 0 0 0 0 0 0 2225 2226 2227 2228 2229 2230 2231 2232 0 0 0 0 0 0 0 0 3513 3514 3515 3517 2257 2258 2259 2260 2273 2274 2275 2276 2289 2290 2291 2292 2257 2258 2259 2260 0 0 0 0 2289 2290 2291 2292 0 0 0 0 2273 2274 2275 2276 2289 2290 2291 2292 0 0 0 0 0 0 0 0 2289 2290 2291 2292 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::Level 6 +DEAL:: 0 0 0 0 0 0 0 0 3533 3534 3535 3537 5786 5788 5790 5792 0 0 0 0 0 0 0 0 5786 5788 5790 5792 6586 6588 6590 6592 3533 3534 3535 3537 5786 5788 5790 5792 8414 8417 8420 8423 12324 12328 12332 12336 5786 5788 5790 5792 6586 6588 6590 6592 12324 12328 12332 12336 13444 13448 13452 13456 0 0 0 0 3693 3694 3695 3697 0 0 0 0 6034 6036 6038 6040 3693 3694 3695 3697 8414 8417 8420 8423 6034 6036 6038 6040 12596 12600 12604 12608 0 0 0 0 6034 6036 6038 6040 0 0 0 0 7506 7508 7510 7512 6034 6036 6038 6040 12596 12600 12604 12608 7506 7508 7510 7512 15060 15064 15068 15072 8414 8417 8420 8423 12324 12328 12332 12336 12596 12600 12604 12608 12964 12968 12972 12976 12324 12328 12332 12336 13444 13448 13452 13456 12964 12968 12972 12976 13988 13992 13996 14000 12596 12600 12604 12608 12964 12968 12972 12976 15060 15064 15068 15072 15396 15400 15404 15408 12964 12968 12972 12976 13988 13992 13996 14000 15396 15400 15404 15408 16420 16424 16428 16432 0 0 0 0 0 0 0 0 6586 6588 6590 6592 6810 6812 6814 6816 0 0 0 0 0 0 0 0 6810 6812 6814 6816 4273 4274 4275 4277 6586 6588 6590 6592 6810 6812 6814 6816 13444 13448 13452 13456 13988 13992 13996 14000 6810 6812 6814 6816 4273 4274 4275 4277 13988 13992 13996 14000 10046 10049 10052 10055 13444 13448 13452 13456 13988 13992 13996 14000 13988 13992 13996 14000 14500 14504 14508 14512 13988 13992 13996 14000 10046 10049 10052 10055 14500 14504 14508 14512 14996 15000 15004 15008 13988 13992 13996 14000 14500 14504 14508 14512 16420 16424 16428 16432 16932 16936 16940 16944 14500 14504 14508 14512 14996 15000 15004 15008 16932 16936 16940 16944 17524 17528 17532 17536 10046 10049 10052 10055 4673 4674 4675 4677 14996 15000 15004 15008 7590 7592 7594 7596 4673 4674 4675 4677 0 0 0 0 7590 7592 7594 7596 0 0 0 0 14996 15000 15004 15008 7590 7592 7594 7596 17524 17528 17532 17536 9062 9064 9066 9068 7590 7592 7594 7596 0 0 0 0 9062 9064 9066 9068 0 0 0 0 0 0 0 0 7506 7508 7510 7512 0 0 0 0 7826 7828 7830 7832 7506 7508 7510 7512 15060 15064 15068 15072 7826 7828 7830 7832 16180 16184 16188 16192 0 0 0 0 7826 7828 7830 7832 0 0 0 0 4933 4934 4935 4937 7826 7828 7830 7832 16180 16184 16188 16192 4933 4934 4935 4937 11502 11505 11508 11511 15060 15064 15068 15072 15396 15400 15404 15408 16180 16184 16188 16192 16548 16552 16556 16560 15396 15400 15404 15408 16420 16424 16428 16432 16548 16552 16556 16560 17572 17576 17580 17584 16180 16184 16188 16192 16548 16552 16556 16560 11502 11505 11508 11511 17060 17064 17068 17072 16548 16552 16556 16560 17572 17576 17580 17584 17060 17064 17068 17072 18372 18376 18380 18384 11502 11505 11508 11511 17060 17064 17068 17072 5333 5334 5335 5337 8642 8644 8646 8648 17060 17064 17068 17072 18372 18376 18380 18384 8642 8644 8646 8648 9442 9444 9446 9448 5333 5334 5335 5337 8642 8644 8646 8648 0 0 0 0 0 0 0 0 8642 8644 8646 8648 9442 9444 9446 9448 0 0 0 0 0 0 0 0 16420 16424 16428 16432 16932 16936 16940 16944 17572 17576 17580 17584 18084 18088 18092 18096 16932 16936 16940 16944 17524 17528 17532 17536 18084 18088 18092 18096 18580 18584 18588 18592 17572 17576 17580 17584 18084 18088 18092 18096 18372 18376 18380 18384 18980 18984 18988 18992 18084 18088 18092 18096 18580 18584 18588 18592 18980 18984 18988 18992 13122 13125 13128 13131 17524 17528 17532 17536 9062 9064 9066 9068 18580 18584 18588 18592 9382 9384 9386 9388 9062 9064 9066 9068 0 0 0 0 9382 9384 9386 9388 0 0 0 0 18580 18584 18588 18592 9382 9384 9386 9388 13122 13125 13128 13131 5913 5914 5915 5917 9382 9384 9386 9388 0 0 0 0 5913 5914 5915 5917 0 0 0 0 18372 18376 18380 18384 18980 18984 18988 18992 9442 9444 9446 9448 9666 9668 9670 9672 18980 18984 18988 18992 13122 13125 13128 13131 9666 9668 9670 9672 6073 6074 6075 6077 9442 9444 9446 9448 9666 9668 9670 9672 0 0 0 0 0 0 0 0 9666 9668 9670 9672 6073 6074 6075 6077 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2609 2610 2611 2612 0 0 0 0 0 0 0 0 2609 2610 2611 2612 3533 3534 3535 3537 0 0 0 0 2609 2610 2611 2612 0 0 0 0 3693 3694 3695 3697 2609 2610 2611 2612 3533 3534 3535 3537 3693 3694 3695 3697 8414 8417 8420 8423 0 0 0 0 0 0 0 0 4273 4274 4275 4277 2673 2674 2675 2676 0 0 0 0 0 0 0 0 2673 2674 2675 2676 0 0 0 0 4273 4274 4275 4277 2673 2674 2675 2676 10046 10049 10052 10055 4673 4674 4675 4677 2673 2674 2675 2676 0 0 0 0 4673 4674 4675 4677 0 0 0 0 0 0 0 0 4933 4934 4935 4937 0 0 0 0 2737 2738 2739 2740 4933 4934 4935 4937 11502 11505 11508 11511 2737 2738 2739 2740 5333 5334 5335 5337 0 0 0 0 2737 2738 2739 2740 0 0 0 0 0 0 0 0 2737 2738 2739 2740 5333 5334 5335 5337 0 0 0 0 0 0 0 0 13122 13125 13128 13131 5913 5914 5915 5917 6073 6074 6075 6077 2801 2802 2803 2804 5913 5914 5915 5917 0 0 0 0 2801 2802 2803 2804 0 0 0 0 6073 6074 6075 6077 2801 2802 2803 2804 0 0 0 0 0 0 0 0 2801 2802 2803 2804 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::Level 7 +DEAL:: 0 0 0 0 0 0 0 0 0 0 0 0 2865 2866 2867 2868 0 0 0 0 0 0 0 0 2865 2866 2867 2868 2913 2914 2915 2916 0 0 0 0 2865 2866 2867 2868 0 0 0 0 3089 3090 3091 3092 2865 2866 2867 2868 2913 2914 2915 2916 3089 3090 3091 3092 3137 3138 3139 3140 0 0 0 0 0 0 0 0 2913 2914 2915 2916 2929 2930 2931 2932 0 0 0 0 0 0 0 0 2929 2930 2931 2932 3361 3362 3363 3364 2913 2914 2915 2916 2929 2930 2931 2932 3137 3138 3139 3140 3153 3154 3155 3156 2929 2930 2931 2932 3361 3362 3363 3364 3153 3154 3155 3156 3457 3458 3459 3460 0 0 0 0 0 0 0 0 0 0 0 0 2993 2994 2995 2996 0 0 0 0 0 0 0 0 2993 2994 2995 2996 3105 3106 3107 3108 0 0 0 0 2993 2994 2995 2996 0 0 0 0 3025 3026 3027 3028 2993 2994 2995 2996 3105 3106 3107 3108 3025 3026 3027 3028 3201 3202 3203 3204 0 0 0 0 3025 3026 3027 3028 0 0 0 0 3057 3058 3059 3060 3025 3026 3027 3028 3201 3202 3203 3204 3057 3058 3059 3060 3233 3234 3235 3236 0 0 0 0 3057 3058 3059 3060 0 0 0 0 3857 3858 3859 3860 3057 3058 3059 3060 3233 3234 3235 3236 3857 3858 3859 3860 3969 3970 3971 3972 0 0 0 0 3089 3090 3091 3092 3105 3106 3107 3108 3121 3122 3123 3124 3089 3090 3091 3092 3137 3138 3139 3140 3121 3122 3123 3124 3169 3170 3171 3172 3105 3106 3107 3108 3121 3122 3123 3124 3201 3202 3203 3204 3217 3218 3219 3220 3121 3122 3123 3124 3169 3170 3171 3172 3217 3218 3219 3220 3265 3266 3267 3268 3137 3138 3139 3140 3153 3154 3155 3156 3169 3170 3171 3172 3185 3186 3187 3188 3153 3154 3155 3156 3457 3458 3459 3460 3185 3186 3187 3188 3489 3490 3491 3492 3169 3170 3171 3172 3185 3186 3187 3188 3265 3266 3267 3268 3281 3282 3283 3284 3185 3186 3187 3188 3489 3490 3491 3492 3281 3282 3283 3284 3585 3586 3587 3588 3201 3202 3203 3204 3217 3218 3219 3220 3233 3234 3235 3236 3249 3250 3251 3252 3217 3218 3219 3220 3265 3266 3267 3268 3249 3250 3251 3252 3297 3298 3299 3300 3233 3234 3235 3236 3249 3250 3251 3252 3969 3970 3971 3972 3985 3986 3987 3988 3249 3250 3251 3252 3297 3298 3299 3300 3985 3986 3987 3988 4033 4034 4035 4036 3265 3266 3267 3268 3281 3282 3283 3284 3297 3298 3299 3300 3313 3314 3315 3316 3281 3282 3283 3284 3585 3586 3587 3588 3313 3314 3315 3316 3617 3618 3619 3620 3297 3298 3299 3300 3313 3314 3315 3316 4033 4034 4035 4036 4049 4050 4051 4052 3313 3314 3315 3316 3617 3618 3619 3620 4049 4050 4051 4052 4353 4354 4355 4356 0 0 0 0 0 0 0 0 3361 3362 3363 3364 3377 3378 3379 3380 0 0 0 0 0 0 0 0 3377 3378 3379 3380 3425 3426 3427 3428 3361 3362 3363 3364 3377 3378 3379 3380 3457 3458 3459 3460 3473 3474 3475 3476 3377 3378 3379 3380 3425 3426 3427 3428 3473 3474 3475 3476 3521 3522 3523 3524 0 0 0 0 0 0 0 0 3425 3426 3427 3428 3441 3442 3443 3444 0 0 0 0 0 0 0 0 3441 3442 3443 3444 0 0 0 0 3425 3426 3427 3428 3441 3442 3443 3444 3521 3522 3523 3524 3537 3538 3539 3540 3441 3442 3443 3444 0 0 0 0 3537 3538 3539 3540 0 0 0 0 3457 3458 3459 3460 3473 3474 3475 3476 3489 3490 3491 3492 3505 3506 3507 3508 3473 3474 3475 3476 3521 3522 3523 3524 3505 3506 3507 3508 3553 3554 3555 3556 3489 3490 3491 3492 3505 3506 3507 3508 3585 3586 3587 3588 3601 3602 3603 3604 3505 3506 3507 3508 3553 3554 3555 3556 3601 3602 3603 3604 3649 3650 3651 3652 3521 3522 3523 3524 3537 3538 3539 3540 3553 3554 3555 3556 3569 3570 3571 3572 3537 3538 3539 3540 0 0 0 0 3569 3570 3571 3572 3745 3746 3747 3748 3553 3554 3555 3556 3569 3570 3571 3572 3649 3650 3651 3652 3665 3666 3667 3668 3569 3570 3571 3572 3745 3746 3747 3748 3665 3666 3667 3668 3777 3778 3779 3780 3585 3586 3587 3588 3601 3602 3603 3604 3617 3618 3619 3620 3633 3634 3635 3636 3601 3602 3603 3604 3649 3650 3651 3652 3633 3634 3635 3636 3681 3682 3683 3684 3617 3618 3619 3620 3633 3634 3635 3636 4353 4354 4355 4356 4369 4370 4371 4372 3633 3634 3635 3636 3681 3682 3683 3684 4369 4370 4371 4372 4417 4418 4419 4420 3649 3650 3651 3652 3665 3666 3667 3668 3681 3682 3683 3684 3697 3698 3699 3700 3665 3666 3667 3668 3777 3778 3779 3780 3697 3698 3699 3700 3809 3810 3811 3812 3681 3682 3683 3684 3697 3698 3699 3700 4417 4418 4419 4420 4433 4434 4435 4436 3697 3698 3699 3700 3809 3810 3811 3812 4433 4434 4435 4436 4609 4610 4611 4612 0 0 0 0 0 0 0 0 3745 3746 3747 3748 3761 3762 3763 3764 0 0 0 0 0 0 0 0 3761 3762 3763 3764 0 0 0 0 3745 3746 3747 3748 3761 3762 3763 3764 3777 3778 3779 3780 3793 3794 3795 3796 3761 3762 3763 3764 0 0 0 0 3793 3794 3795 3796 0 0 0 0 3777 3778 3779 3780 3793 3794 3795 3796 3809 3810 3811 3812 3825 3826 3827 3828 3793 3794 3795 3796 0 0 0 0 3825 3826 3827 3828 0 0 0 0 3809 3810 3811 3812 3825 3826 3827 3828 4609 4610 4611 4612 4625 4626 4627 4628 3825 3826 3827 3828 0 0 0 0 4625 4626 4627 4628 0 0 0 0 0 0 0 0 3857 3858 3859 3860 0 0 0 0 3889 3890 3891 3892 3857 3858 3859 3860 3969 3970 3971 3972 3889 3890 3891 3892 4001 4002 4003 4004 0 0 0 0 3889 3890 3891 3892 0 0 0 0 3921 3922 3923 3924 3889 3890 3891 3892 4001 4002 4003 4004 3921 3922 3923 3924 4097 4098 4099 4100 0 0 0 0 3921 3922 3923 3924 0 0 0 0 3953 3954 3955 3956 3921 3922 3923 3924 4097 4098 4099 4100 3953 3954 3955 3956 4129 4130 4131 4132 0 0 0 0 3953 3954 3955 3956 0 0 0 0 0 0 0 0 3953 3954 3955 3956 4129 4130 4131 4132 0 0 0 0 0 0 0 0 3969 3970 3971 3972 3985 3986 3987 3988 4001 4002 4003 4004 4017 4018 4019 4020 3985 3986 3987 3988 4033 4034 4035 4036 4017 4018 4019 4020 4065 4066 4067 4068 4001 4002 4003 4004 4017 4018 4019 4020 4097 4098 4099 4100 4113 4114 4115 4116 4017 4018 4019 4020 4065 4066 4067 4068 4113 4114 4115 4116 4161 4162 4163 4164 4033 4034 4035 4036 4049 4050 4051 4052 4065 4066 4067 4068 4081 4082 4083 4084 4049 4050 4051 4052 4353 4354 4355 4356 4081 4082 4083 4084 4385 4386 4387 4388 4065 4066 4067 4068 4081 4082 4083 4084 4161 4162 4163 4164 4177 4178 4179 4180 4081 4082 4083 4084 4385 4386 4387 4388 4177 4178 4179 4180 4481 4482 4483 4484 4097 4098 4099 4100 4113 4114 4115 4116 4129 4130 4131 4132 4145 4146 4147 4148 4113 4114 4115 4116 4161 4162 4163 4164 4145 4146 4147 4148 4193 4194 4195 4196 4129 4130 4131 4132 4145 4146 4147 4148 0 0 0 0 4241 4242 4243 4244 4145 4146 4147 4148 4193 4194 4195 4196 4241 4242 4243 4244 4289 4290 4291 4292 4161 4162 4163 4164 4177 4178 4179 4180 4193 4194 4195 4196 4209 4210 4211 4212 4177 4178 4179 4180 4481 4482 4483 4484 4209 4210 4211 4212 4513 4514 4515 4516 4193 4194 4195 4196 4209 4210 4211 4212 4289 4290 4291 4292 4305 4306 4307 4308 4209 4210 4211 4212 4513 4514 4515 4516 4305 4306 4307 4308 4737 4738 4739 4740 0 0 0 0 4241 4242 4243 4244 0 0 0 0 4273 4274 4275 4276 4241 4242 4243 4244 4289 4290 4291 4292 4273 4274 4275 4276 4321 4322 4323 4324 0 0 0 0 4273 4274 4275 4276 0 0 0 0 0 0 0 0 4273 4274 4275 4276 4321 4322 4323 4324 0 0 0 0 0 0 0 0 4289 4290 4291 4292 4305 4306 4307 4308 4321 4322 4323 4324 4337 4338 4339 4340 4305 4306 4307 4308 4737 4738 4739 4740 4337 4338 4339 4340 4769 4770 4771 4772 4321 4322 4323 4324 4337 4338 4339 4340 0 0 0 0 0 0 0 0 4337 4338 4339 4340 4769 4770 4771 4772 0 0 0 0 0 0 0 0 4353 4354 4355 4356 4369 4370 4371 4372 4385 4386 4387 4388 4401 4402 4403 4404 4369 4370 4371 4372 4417 4418 4419 4420 4401 4402 4403 4404 4449 4450 4451 4452 4385 4386 4387 4388 4401 4402 4403 4404 4481 4482 4483 4484 4497 4498 4499 4500 4401 4402 4403 4404 4449 4450 4451 4452 4497 4498 4499 4500 4545 4546 4547 4548 4417 4418 4419 4420 4433 4434 4435 4436 4449 4450 4451 4452 4465 4466 4467 4468 4433 4434 4435 4436 4609 4610 4611 4612 4465 4466 4467 4468 4641 4642 4643 4644 4449 4450 4451 4452 4465 4466 4467 4468 4545 4546 4547 4548 4561 4562 4563 4564 4465 4466 4467 4468 4641 4642 4643 4644 4561 4562 4563 4564 4673 4674 4675 4676 4481 4482 4483 4484 4497 4498 4499 4500 4513 4514 4515 4516 4529 4530 4531 4532 4497 4498 4499 4500 4545 4546 4547 4548 4529 4530 4531 4532 4577 4578 4579 4580 4513 4514 4515 4516 4529 4530 4531 4532 4737 4738 4739 4740 4753 4754 4755 4756 4529 4530 4531 4532 4577 4578 4579 4580 4753 4754 4755 4756 4801 4802 4803 4804 4545 4546 4547 4548 4561 4562 4563 4564 4577 4578 4579 4580 4593 4594 4595 4596 4561 4562 4563 4564 4673 4674 4675 4676 4593 4594 4595 4596 4705 4706 4707 4708 4577 4578 4579 4580 4593 4594 4595 4596 4801 4802 4803 4804 4817 4818 4819 4820 4593 4594 4595 4596 4705 4706 4707 4708 4817 4818 4819 4820 0 0 0 0 4609 4610 4611 4612 4625 4626 4627 4628 4641 4642 4643 4644 4657 4658 4659 4660 4625 4626 4627 4628 0 0 0 0 4657 4658 4659 4660 0 0 0 0 4641 4642 4643 4644 4657 4658 4659 4660 4673 4674 4675 4676 4689 4690 4691 4692 4657 4658 4659 4660 0 0 0 0 4689 4690 4691 4692 0 0 0 0 4673 4674 4675 4676 4689 4690 4691 4692 4705 4706 4707 4708 4721 4722 4723 4724 4689 4690 4691 4692 0 0 0 0 4721 4722 4723 4724 0 0 0 0 4705 4706 4707 4708 4721 4722 4723 4724 0 0 0 0 0 0 0 0 4721 4722 4723 4724 0 0 0 0 0 0 0 0 0 0 0 0 4737 4738 4739 4740 4753 4754 4755 4756 4769 4770 4771 4772 4785 4786 4787 4788 4753 4754 4755 4756 4801 4802 4803 4804 4785 4786 4787 4788 4833 4834 4835 4836 4769 4770 4771 4772 4785 4786 4787 4788 0 0 0 0 0 0 0 0 4785 4786 4787 4788 4833 4834 4835 4836 0 0 0 0 0 0 0 0 4801 4802 4803 4804 4817 4818 4819 4820 4833 4834 4835 4836 4849 4850 4851 4852 4817 4818 4819 4820 0 0 0 0 4849 4850 4851 4852 0 0 0 0 4833 4834 4835 4836 4849 4850 4851 4852 0 0 0 0 0 0 0 0 4849 4850 4851 4852 0 0 0 0 0 0 0 0 0 0 0 0 +DEAL::diff 0 diff --git a/tests/multigrid/mg_output.cc b/tests/multigrid/mg_output_neumann.cc similarity index 89% rename from tests/multigrid/mg_output.cc rename to tests/multigrid/mg_output_neumann.cc index c82d8560e2..ef92fedf2f 100644 --- a/tests/multigrid/mg_output.cc +++ b/tests/multigrid/mg_output_neumann.cc @@ -38,7 +38,6 @@ #include #include #include -#include #include #include #include @@ -52,7 +51,7 @@ using namespace std; template void -reinit_vector (const dealii::MGDoFHandler &mg_dof, +reinit_vector (const dealii::DoFHandler &mg_dof, MGLevelObject > &v) { for (unsigned int level=v.min_level(); @@ -98,13 +97,13 @@ void refine_mesh (Triangulation &triangulation) } template -void initialize (const MGDoFHandler &dof, +void initialize (const DoFHandler &dof, Vector &u) { unsigned int counter=0; const unsigned int dofs_per_cell = dof.get_fe().dofs_per_cell; std::vector dof_indices(dofs_per_cell); - for (typename MGDoFHandler::active_cell_iterator + for (typename DoFHandler::active_cell_iterator cell = dof.begin_active(); cell != dof.end(); ++cell) { @@ -115,7 +114,7 @@ void initialize (const MGDoFHandler &dof, } template -void initialize (const MGDoFHandler &dof, +void initialize (const DoFHandler &dof, MGLevelObject > &u) { const unsigned int dofs_per_cell = dof.get_fe().dofs_per_cell; @@ -125,9 +124,9 @@ void initialize (const MGDoFHandler &dof, for (unsigned int l=0; l::cell_iterator - cell = dof.begin(l); - cell != dof.end(l); ++cell) + for (typename DoFHandler::cell_iterator + cell = dof.begin_mg(l); + cell != dof.end_mg(l); ++cell) { cell->get_mg_dof_indices(dof_indices); for (unsigned int f=0; f::faces_per_cell; ++f) @@ -150,7 +149,7 @@ void initialize (const MGDoFHandler &dof, } template -void print (const MGDoFHandler &dof, +void print (const DoFHandler &dof, MGLevelObject > &u) { const unsigned int dofs_per_cell = dof.get_fe().dofs_per_cell; @@ -159,9 +158,9 @@ void print (const MGDoFHandler &dof, { deallog << std::endl; deallog << "Level " << l << std::endl; - for (typename MGDoFHandler::cell_iterator - cell = dof.begin(l); - cell != dof.end(l); ++cell) + for (typename DoFHandler::cell_iterator + cell = dof.begin_mg(l); + cell != dof.end_mg(l); ++cell) { cell->get_mg_dof_indices(dof_indices); for (unsigned int i=0; i &dof, } template -void print_diff (const MGDoFHandler &dof_1, const MGDoFHandler &dof_2, +void print_diff (const DoFHandler &dof_1, const DoFHandler &dof_2, const Vector &u, const Vector &v) { Vector diff; @@ -179,7 +178,7 @@ void print_diff (const MGDoFHandler &dof_1, const MGDoFHandler &dof_2, const unsigned int dofs_per_cell = dof_1.get_fe().dofs_per_cell; std::vector dof_indices_1(dofs_per_cell); std::vector dof_indices_2(dofs_per_cell); - for (typename MGDoFHandler::active_cell_iterator + for (typename DoFHandler::active_cell_iterator cell_1 = dof_1.begin_active(), cell_2 = dof_2.begin_active(); cell_1 != dof_1.end(); ++cell_1, ++cell_2) { @@ -199,6 +198,8 @@ void check_simple(const FiniteElement &fe) Triangulation tr(Triangulation::limit_level_difference_at_vertices); GridGenerator::hyper_cube(tr,-1,1); + typename FunctionMap::type dirichlet_boundary_functions; + tr.refine_global (1); refine_mesh(tr); refine_mesh(tr); @@ -207,11 +208,13 @@ void check_simple(const FiniteElement &fe) refine_mesh(tr); refine_mesh(tr); - MGDoFHandler mgdof(tr); + DoFHandler mgdof(tr); mgdof.distribute_dofs(fe); + mgdof.distribute_mg_dofs(fe); - MGDoFHandler mgdof_renumbered(tr); + DoFHandler mgdof_renumbered(tr); mgdof_renumbered.distribute_dofs(fe); + mgdof_renumbered.distribute_mg_dofs(fe); ConstraintMatrix hnc, hnc_renumbered; hnc.clear (); @@ -237,10 +240,10 @@ void check_simple(const FiniteElement &fe) MGConstrainedDoFs mg_constrained_dofs; - mg_constrained_dofs.initialize(mgdof); + mg_constrained_dofs.initialize(mgdof, dirichlet_boundary_functions); MGConstrainedDoFs mg_constrained_dofs_renumbered; - mg_constrained_dofs_renumbered.initialize(mgdof_renumbered); + mg_constrained_dofs_renumbered.initialize(mgdof_renumbered, dirichlet_boundary_functions); MGTransferPrebuilt > transfer(hnc, mg_constrained_dofs); transfer.build_matrices(mgdof); diff --git a/tests/multigrid/mg_output.output b/tests/multigrid/mg_output_neumann.output similarity index 100% rename from tests/multigrid/mg_output.output rename to tests/multigrid/mg_output_neumann.output diff --git a/tests/multigrid/step-16-bdry1.output b/tests/multigrid/step-16-bdry1.output index a0be674368..17d2f3bb22 100644 --- a/tests/multigrid/step-16-bdry1.output +++ b/tests/multigrid/step-16-bdry1.output @@ -3,55 +3,55 @@ DEAL::Cycle 0: DEAL:: Number of active cells: 4 DEAL::Number of degrees of freedom: 9 L0: 4 L1: 9 DEAL:: Number of degrees of freedom: 9 (by level: 4, 9) -DEAL:cg::Starting value 1.00000 -DEAL:cg::Convergence step 1 value 0.00000 -DEAL:: 1 CG iterations needed to obtain convergence. +DEAL:cg::Starting value 1.36931 +DEAL:cg::Convergence step 4 value 8.55703e-14 +DEAL:: 4 CG iterations needed to obtain convergence. DEAL::Cycle 1: DEAL:: Number of active cells: 16 DEAL::Number of degrees of freedom: 25 L0: 4 L1: 9 L2: 25 DEAL:: Number of degrees of freedom: 25 (by level: 4, 9, 25) -DEAL:cg::Starting value 0.750000 -DEAL:cg::Convergence step 5 value 1.17215e-16 -DEAL:: 5 CG iterations needed to obtain convergence. +DEAL:cg::Starting value 0.843171 +DEAL:cg::Convergence step 6 value 9.78112e-15 +DEAL:: 6 CG iterations needed to obtain convergence. DEAL::Cycle 2: DEAL:: Number of active cells: 28 DEAL::Number of degrees of freedom: 41 L0: 4 L1: 9 L2: 25 L3: 25 DEAL:: Number of degrees of freedom: 41 (by level: 4, 9, 25, 25) -DEAL:cg::Starting value 0.628894 -DEAL:cg::Convergence step 7 value 2.44070e-14 -DEAL:: 7 CG iterations needed to obtain convergence. +DEAL:cg::Starting value 0.737526 +DEAL:cg::Convergence step 8 value 8.81541e-15 +DEAL:: 8 CG iterations needed to obtain convergence. DEAL::Cycle 3: DEAL:: Number of active cells: 40 DEAL::Number of degrees of freedom: 57 L0: 4 L1: 9 L2: 25 L3: 25 L4: 25 DEAL:: Number of degrees of freedom: 57 (by level: 4, 9, 25, 25, 25) -DEAL:cg::Starting value 0.620541 -DEAL:cg::Convergence step 7 value 5.62398e-13 -DEAL:: 7 CG iterations needed to obtain convergence. +DEAL:cg::Starting value 0.730417 +DEAL:cg::Convergence step 8 value 9.36194e-14 +DEAL:: 8 CG iterations needed to obtain convergence. DEAL::Cycle 4: DEAL:: Number of active cells: 52 DEAL::Number of degrees of freedom: 73 L0: 4 L1: 9 L2: 25 L3: 25 L4: 25 L5: 25 DEAL:: Number of degrees of freedom: 73 (by level: 4, 9, 25, 25, 25, 25) -DEAL:cg::Starting value 0.620015 -DEAL:cg::Convergence step 8 value 4.43441e-14 +DEAL:cg::Starting value 0.729970 +DEAL:cg::Convergence step 8 value 2.81558e-13 DEAL:: 8 CG iterations needed to obtain convergence. DEAL::Cycle 5: DEAL:: Number of active cells: 160 DEAL::Number of degrees of freedom: 201 L0: 4 L1: 9 L2: 25 L3: 65 L4: 65 L5: 65 L6: 65 DEAL:: Number of degrees of freedom: 201 (by level: 4, 9, 25, 65, 65, 65, 65) -DEAL:cg::Starting value 0.415999 -DEAL:cg::Convergence step 9 value 4.46942e-14 +DEAL:cg::Starting value 0.492333 +DEAL:cg::Convergence step 9 value 5.07798e-13 DEAL:: 9 CG iterations needed to obtain convergence. DEAL::Cycle 6: DEAL:: Number of active cells: 304 DEAL::Number of degrees of freedom: 357 L0: 4 L1: 9 L2: 25 L3: 81 L4: 81 L5: 81 L6: 81 L7: 153 DEAL:: Number of degrees of freedom: 357 (by level: 4, 9, 25, 81, 81, 81, 81, 153) -DEAL:cg::Starting value 0.377840 -DEAL:cg::Convergence step 8 value 1.00219e-13 +DEAL:cg::Starting value 0.404672 +DEAL:cg::Convergence step 8 value 5.14673e-13 DEAL:: 8 CG iterations needed to obtain convergence. DEAL::Cycle 7: DEAL:: Number of active cells: 676 DEAL::Number of degrees of freedom: 761 L0: 4 L1: 9 L2: 25 L3: 81 L4: 81 L5: 81 L6: 121 L7: 209 L8: 465 DEAL:: Number of degrees of freedom: 761 (by level: 4, 9, 25, 81, 81, 81, 121, 209, 465) -DEAL:cg::Starting value 0.377713 -DEAL:cg::Convergence step 8 value 5.29691e-13 -DEAL:: 8 CG iterations needed to obtain convergence. +DEAL:cg::Starting value 0.404553 +DEAL:cg::Convergence step 9 value 6.34913e-14 +DEAL:: 9 CG iterations needed to obtain convergence. diff --git a/tests/multigrid/step-39-02a.cc b/tests/multigrid/step-39-02a.cc index 6eb5dafa25..ae740381c3 100644 --- a/tests/multigrid/step-39-02a.cc +++ b/tests/multigrid/step-39-02a.cc @@ -550,10 +550,10 @@ namespace Step39 mg_smoother.set_symmetric(true); mg_smoother.set_variable(false); - MGMatrix, Vector > mgmatrix(&mg_matrix); - MGMatrix, Vector > mgdown(&mg_matrix_dg_down); - MGMatrix, Vector > mgup(&mg_matrix_dg_up); - MGMatrix, Vector > mgedge(&mg_matrix_in_out); + mg::Matrix > mgmatrix(mg_matrix); + mg::Matrix > mgdown(mg_matrix_dg_down); + mg::Matrix > mgup(mg_matrix_dg_up); + mg::Matrix > mgedge(mg_matrix_in_out); Multigrid > mg(mg_dof_handler, mgmatrix, mg_coarse, mg_transfer, -- 2.39.5