diff --git a/demo/CMakeLists.txt b/demo/CMakeLists.txt index decc784d3f5f148adf6ced04986810b056103f1a..2d9c54973d4a507099ac284099a75e473e7d1eda 100644 --- a/demo/CMakeLists.txt +++ b/demo/CMakeLists.txt @@ -105,6 +105,16 @@ if (TOOLS_DIR) add_executable("drivenCavity" ${drivenCavity}) target_link_libraries("drivenCavity" ${BASIS_LIBS}) + include_directories(/home/spraetor/projects/src/common) + set(navierStokes ${TOOLS_DIR}/diffuseDomain/POperators.cc +# ${TOOLS_DIR}/baseProblems/NavierStokes_TaylorHood.cc + ${TOOLS_DIR}/baseProblems/NavierStokes_TaylorHood_RB.cc + ${TOOLS_DIR}/baseProblems/time/ExtendedRosenbrockStationary.cc + /home/spraetor/projects/src/common/GeometryTools.cc + src/navierStokes.cc) + add_executable("navierStokes" ${navierStokes}) + target_link_libraries("navierStokes" ${BASIS_LIBS}) + endif() else() message(WARNING "No tools directory specified! Some demos will not be build.") diff --git a/demo/init/elliptImplicit.dat.2d b/demo/init/elliptImplicit.dat.2d index eb6b42ee80811a00f4d1cdb8255fd9597a8a62a2..11f9b6950e78c5e438b7217b80aca3e6884f9617 100644 --- a/demo/init/elliptImplicit.dat.2d +++ b/demo/init/elliptImplicit.dat.2d @@ -12,8 +12,8 @@ mesh->refinement->interface width: 0.01 ellipt->mesh: elliptMesh ellipt->dim: 2 ellipt->components: 1 -ellipt->polynomial degree[0]: 1 - +ellipt->polynomial degree[0]: 2 + ellipt->solver: cg ellipt->solver->ell: 1 ellipt->solver->max iteration: 1000 @@ -33,7 +33,7 @@ ellipt->marker[0]->MSGamma: 0.5 ellipt->adapt[0]->tolerance: 1e-6 ellipt->adapt[0]->refine bisections: 2 -ellipt->adapt->max iteration: 5 +ellipt->adapt->max iteration: 10 ellipt->output->filename: output/elliptImplicit.2d ellipt->output->ParaView format: 1 diff --git a/demo/init/navierStokes.dat.2d b/demo/init/navierStokes.dat.2d new file mode 100644 index 0000000000000000000000000000000000000000..ab21693d0d08455af80f31daaa28fd7a4bc42c52 --- /dev/null +++ b/demo/init/navierStokes.dat.2d @@ -0,0 +1,50 @@ +dimension of world: 2 + +% ====================== VARIABLES ======================== +output_folder: output +output_postfix: _ns +mesh_name: mesh +polynomial-degree: 1 + +% ====================== MESH ============================= +mesh->H: 4.1 + +% ====================== MESH ============================= +${mesh_name}->macro file name: macro/kanal_msh.2d +${mesh_name}->global refinements: 0 +${mesh_name}->check: 0 + +% ====================== INCLUDES ========================= +#include "init/navierStokes_TaylorHood.inc.2d" + +% ====================== USER_PARAMETER - NS ============== + +ns->viscosity: 1/100 +ns->theta: 0.5 +%ns->force: [0.0, 9.81] % gravitational force [m/s^2] + +ns->force dirichlet bc: 0 + +ns->Um: 1.5 + +ns->initial velocity: 0 +ns->initial velocity value: 0.0 +ns->laplace operator: 0 % 0... div(nu*grad(u)), 1... div(0.5*nu*(grad(u)+grad(u)^T)) [sehr langsam] +ns->non-linear term: 2 % 0... u^old*grad(u_i^old), 1... u'*grad(u_i^old), 2... u^old*grad(u'_i) + +% ====================== TIMESTEPS ======================== +adapt->max iteration: 1 +adapt->max timestep iteration: 1 +adapt->max time iteration: 1 + +adapt->timestep: 5.e-4 +adapt->max timestep: 1e+10 +adapt->min timestep: 1e-6 +adapt->start time: 0.0 +adapt->end time: 10 + +% ====================== ESTIMATORS ======================= +adapt->strategy: 0 % 0=explicit, 1=implicit + +WAIT: 1 + diff --git a/demo/init/navierStokes_Chorin.inc.2d b/demo/init/navierStokes_Chorin.inc.2d new file mode 100644 index 0000000000000000000000000000000000000000..9c5c23b7b26c1d044a0b0d942362e5a1cf8d0651 --- /dev/null +++ b/demo/init/navierStokes_Chorin.inc.2d @@ -0,0 +1,111 @@ +${ns}->predictor->mesh: ${mesh_name} +${ns}->pressure->mesh: ${mesh_name} +${ns}->corrector->mesh: ${mesh_name} + +% ============ USER_PARAMETERS - NS ================================= +${ns}->viscosity: 1 +${ns}->beta: 1 +${ns}->sigma: 0.072 +${ns}->exponent: 2 +%${ns}->force: [0.0, -9.81] % gravitational force + +${ns}->viscosity1: 1 +${ns}->viscosity2: ${${ns}->viscosity1}/1000 + +${ns}->force dirichlet bc: 0 +${ns}->poisson problem pertubation: 0 % 1..applySingularPertubation, 2..applySingularDBC, 0..nothing +${ns}->simple algorithm: 1 +${ns}->calculate pressure: 0 +${ns}->initial velocity: 3 +${ns}->initial velocity value: 0.1 +${ns}->laplace operator: 1 % 0... div(nu*grad(u)), 1... div(0.5*nu*(grad(u)+grad(u)^T)) [sehr langsam] +${ns}->non-linear term: 2 % 0... u^old*grad(u_i^old), 1... u'*grad(u_i^old), 2... u^old*grad(u'_i) + +% =========== OUTPUT ============================================== + +${ns}->predictor->output->filename: ${output_folder}/velocity/predictor${output_postfix}_ +${ns}->pressure->output->filename: ${output_folder}/velocity/pressure${output_postfix}_ +${ns}->corrector->output->filename: ${output_folder}/velocity/velocity${output_postfix}_ + +% ============= PROBLEM-SPACES ================================== +${ns}->predictor->components: 2 +${ns}->pressure->components: 1 +${ns}->corrector->components: 2 + +${ns}->predictor->polynomial degree[0]: ${polynomial-degree} +${ns}->predictor->polynomial degree[1]: ${polynomial-degree} +${ns}->pressure->polynomial degree[0]: ${polynomial-degree} +${ns}->corrector->polynomial degree[0]: ${polynomial-degree} +${ns}->corrector->polynomial degree[1]: ${polynomial-degree} + +${ns}->predictor->dim: 2 +${ns}->pressure->dim: 2 +${ns}->corrector->dim: 2 + +${ns}->pressure->name: pressure +${ns}->corrector->name: velocity + +% ================== SOLVER ====================================== +${ns}->predictor->solver: umfpack +${ns}->predictor->solver->symmetric strategy: 0 +${ns}->predictor->solver->ell: 3 +${ns}->predictor->solver->max iteration: 500 +${ns}->predictor->solver->restart: 10 % only used for GMRES +${ns}->predictor->solver->tolerance: 1.e-8 +${ns}->predictor->solver->info: 1 +${ns}->predictor->solver->left precon: ilu + +${ns}->predictor->petsc->ksp_type: preonly +${ns}->predictor->petsc->pc_type: lu +${ns}->predictor->petsc->pc_factor_mat_solver_package: mumps + +${ns}->pressure->solver: umfpack +${ns}->pressure->solver->symmetric strategy: 0 +${ns}->pressure->solver->ell: 3 +${ns}->pressure->solver->max iteration: 500 +${ns}->pressure->solver->tolerance: 1.e-8 +${ns}->pressure->solver->info: 1 +${ns}->pressure->solver->left precon: ilu + +${ns}->pressure->petsc->ksp_type: preonly +${ns}->pressure->petsc->pc_type: lu +${ns}->pressure->petsc->pc_factor_mat_solver_package: mumps + +${ns}->corrector->solver: umfpack +${ns}->corrector->solver->symmetric strategy: 0 +${ns}->corrector->solver->ell: 3 +${ns}->corrector->solver->max iteration: 500 +${ns}->corrector->solver->tolerance: 1.e-8 +${ns}->corrector->solver->info: 1 +${ns}->corrector->solver->left precon: ilu + +${ns}->corrector->petsc->ksp_type: preonly +${ns}->corrector->petsc->pc_type: lu +${ns}->corrector->petsc->pc_factor_mat_solver_package: mumps + +% =================== OUTPUT ========================================= +${ns}->predictor->output->ParaView animation: 0 +${ns}->predictor->output->ParaView format: 0 +${ns}->predictor->output->write every i-th timestep: 1 +%${ns}->predictor->output->compression: gzip +${ns}->predictor->output->append index: 1 +${ns}->predictor->output->index length: 7 +${ns}->predictor->output->index decimals: 5 + +${ns}->pressure->output->ParaView animation: 0 +${ns}->pressure->output->ParaView format: 0 +${ns}->pressure->output->write every i-th timestep: 1 +%${ns}->pressure->output->compression: gzip +${ns}->pressure->output->append index: 1 +${ns}->pressure->output->index length: 7 +${ns}->pressure->output->index decimals: 5 + +${ns}->corrector->output->ParaView animation: 1 +${ns}->corrector->output->ParaView vector format: 1 +${ns}->corrector->output->write vector as 3d vector: 1 +${ns}->corrector->output->write every i-th timestep: 1 +%${ns}->corrector->output->compression: gzip +${ns}->corrector->output->append index: 1 +${ns}->corrector->output->index length: 7 +${ns}->corrector->output->index decimals: 5 + diff --git a/demo/init/navierStokes_rb.dat.2d b/demo/init/navierStokes_rb.dat.2d new file mode 100644 index 0000000000000000000000000000000000000000..6e51f5337ac6f850531642adaa62ce1fdaae1d3c --- /dev/null +++ b/demo/init/navierStokes_rb.dat.2d @@ -0,0 +1,83 @@ +dimension of world: 2 + +% ====================== VARIABLES ======================== +output_folder: output +output_postfix: _refinement3 +mesh_name: mesh +polynomial-degree: 1 + +% ====================== MESH ============================= +mesh->H: 4.1 + +obstacle->num vertices: 7 +obstacle->vertex[0]: [1.5, 2] +obstacle->vertex[1]: [2, 2.5] +% obstacle->vertex[2]: [4.5, 2] +% obstacle->vertex[3]: [2, 1.5] +obstacle->vertex[2]: [2.5, 2.05] +obstacle->vertex[3]: [6.5, 2.05] +obstacle->vertex[4]: [6.5, 1.95] +obstacle->vertex[5]: [2.5, 1.95] +obstacle->vertex[6]: [2, 1.5] + +mesh->refinement->initial level: 1 +mesh->refinement->level on interface: 6 +mesh->refinement->level in inner domain: 1 +mesh->refinement->level in outer domain: 1 +mesh->refinement->interface width: 0.5 +mesh->refinement->fade out width: 1.0 + +% ====================== MESH ============================= +${mesh_name}->macro file name: macro/kanal_square_fin.2d +${mesh_name}->global refinements: 0 +${mesh_name}->check: 0 + +% ====================== INCLUDES ========================= +#include "init/navierStokes_TaylorHood.inc.2d" + +% ====================== USER_PARAMETER - NS ============== + +ns->viscosity: 1/500 +ns->theta: 0.5 +%ns->force: [0.0, 9.81] % gravitational force [m/s^2] + +ns->force dirichlet bc: 0 + +ns->Um: 1.5 + +ns->initial velocity: 0 +ns->initial velocity value: 0.0 +ns->laplace operator: 0 % 0... div(nu*grad(u)), 1... div(0.5*nu*(grad(u)+grad(u)^T)) [sehr langsam] +ns->non-linear term: 2 % 0... u^old*grad(u_i^old), 1... u'*grad(u_i^old), 2... u^old*grad(u'_i) + +% ====================== TIMESTEPS ======================== + +adapt->rosenbrock method: rodasp +adapt->fix first timesteps: 0 +adapt->rosenbrock->timestep study: 0 +adapt->rosenbrock->timestep study steps: 0 +adapt->rosenbrock->error weights: [1,1,1] + +adapt[0]->time tolerance: 1.e-3 + +adapt->timestep: 5.e-4 +adapt->max timestep: 1e+1 +adapt->min timestep: 1e-4 +adapt->start time: 0.0 +adapt->end time: 25.0 + +% ====================== ESTIMATORS ======================= +adapt->strategy: 0 % 0=explicit, 1=implicit + +%ns->space->estimator[0]: residual +%ns->space->estimator[0]->error norm: 1 % 1: H1_NORM, 2: L2_NORM +%ns->space->estimator[0]->C0: 0.1 % constant of element residual +%ns->space->estimator[0]->C1: 0.1 % constant of jump residual + +ns->space->marker[0]->strategy: 0 % 0: no adaption 1: GR 2: MS 3: ES 4:GERS + +adapt[0]->tolerance: 1e-3 +adapt->max iteration: 0 + +WAIT: 1 + diff --git a/demo/init/pfc.dat.2d b/demo/init/pfc.dat.2d index 9af20ab72952348e8a92944e9839d1e5367f4e30..dbbb04a952a4d5e7e3d43c59e5de31890a7dc58f 100644 --- a/demo/init/pfc.dat.2d +++ b/demo/init/pfc.dat.2d @@ -39,7 +39,11 @@ adapt->timestep: 1.e-1 adapt->max timestep: 1e+10 adapt->min timestep: 1e-10 adapt->start time: 0.0 -adapt->end time: 50.0 +adapt->end time: 10000.0 + +adapt->sequence->calc end time: 0 +adapt->sequence->timesteps: [1.e-2,2.e-2,5.e-2,1.e-1] +adapt->sequence->number of timesteps: [500,1500,1500,1500] % =================== OUTPUT ========================================= diff --git a/demo/macro/kanal_square.2d b/demo/macro/kanal_square.2d new file mode 100644 index 0000000000000000000000000000000000000000..8965f8f0d34f04a7cf80cd4040ea3f82d7ab7934 --- /dev/null +++ b/demo/macro/kanal_square.2d @@ -0,0 +1,5648 @@ +DIM: 2 +DIM_OF_WORLD: 2 + +number of elements: 1592 +number of vertices: 858 + +vertex coordinates: +0 0 +0.5 0 +1 0 +1.5 0 +2 0 +2.5 0 +3 0 +3.5 0 +4 0 +4.5 0 +5 0 +5.5 0 +6 0 +6.5 0 +7 0 +7.5 0 +8 0 +8.5 0 +9 0 +9.5 0 +10 0 +10.5 0 +11 0 +11.5 0 +12 0 +12.5 0 +13 0 +13.5 0 +14 0 +14.5 0 +15 0 +15.5 0 +16 0 +16.5 0 +17 0 +17.5 0 +18 0 +18.5 0 +19 0 +19.5 0 +20 0 +20.5 0 +21 0 +21.5 0 +22 0 +22.5 0 +23 0 +23.5 0 +24 0 +24.5 0 +25 0 +0 0.5 +0.5 0.5 +1 0.5 +1.5 0.5 +2 0.5 +2.5 0.5 +3 0.5 +3.5 0.5 +4 0.5 +4.5 0.5 +5 0.5 +5.5 0.5 +6 0.5 +6.5 0.5 +7 0.5 +7.5 0.5 +8 0.5 +8.5 0.5 +9 0.5 +9.5 0.5 +10 0.5 +10.5 0.5 +11 0.5 +11.5 0.5 +12 0.5 +12.5 0.5 +13 0.5 +13.5 0.5 +14 0.5 +14.5 0.5 +15 0.5 +15.5 0.5 +16 0.5 +16.5 0.5 +17 0.5 +17.5 0.5 +18 0.5 +18.5 0.5 +19 0.5 +19.5 0.5 +20 0.5 +20.5 0.5 +21 0.5 +21.5 0.5 +22 0.5 +22.5 0.5 +23 0.5 +23.5 0.5 +24 0.5 +24.5 0.5 +25 0.5 +0 1 +0.5 1 +1 1 +1.5 1 +2 1 +2.5 1 +3 1 +3.5 1 +4 1 +4.5 1 +5 1 +5.5 1 +6 1 +6.5 1 +7 1 +7.5 1 +8 1 +8.5 1 +9 1 +9.5 1 +10 1 +10.5 1 +11 1 +11.5 1 +12 1 +12.5 1 +13 1 +13.5 1 +14 1 +14.5 1 +15 1 +15.5 1 +16 1 +16.5 1 +17 1 +17.5 1 +18 1 +18.5 1 +19 1 +19.5 1 +20 1 +20.5 1 +21 1 +21.5 1 +22 1 +22.5 1 +23 1 +23.5 1 +24 1 +24.5 1 +25 1 +0 1.5 +0.5 1.5 +1 1.5 +1.5 1.5 +2 1.5 +2.5 1.5 +3 1.5 +3.5 1.5 +4 1.5 +4.5 1.5 +5 1.5 +5.5 1.5 +6 1.5 +6.5 1.5 +7 1.5 +7.5 1.5 +8 1.5 +8.5 1.5 +9 1.5 +9.5 1.5 +10 1.5 +10.5 1.5 +11 1.5 +11.5 1.5 +12 1.5 +12.5 1.5 +13 1.5 +13.5 1.5 +14 1.5 +14.5 1.5 +15 1.5 +15.5 1.5 +16 1.5 +16.5 1.5 +17 1.5 +17.5 1.5 +18 1.5 +18.5 1.5 +19 1.5 +19.5 1.5 +20 1.5 +20.5 1.5 +21 1.5 +21.5 1.5 +22 1.5 +22.5 1.5 +23 1.5 +23.5 1.5 +24 1.5 +24.5 1.5 +25 1.5 +0 2 +0.5 2 +1 2 +1.5 2 +2.5 2 +3 2 +3.5 2 +4 2 +4.5 2 +5 2 +5.5 2 +6 2 +6.5 2 +7 2 +7.5 2 +8 2 +8.5 2 +9 2 +9.5 2 +10 2 +10.5 2 +11 2 +11.5 2 +12 2 +12.5 2 +13 2 +13.5 2 +14 2 +14.5 2 +15 2 +15.5 2 +16 2 +16.5 2 +17 2 +17.5 2 +18 2 +18.5 2 +19 2 +19.5 2 +20 2 +20.5 2 +21 2 +21.5 2 +22 2 +22.5 2 +23 2 +23.5 2 +24 2 +24.5 2 +25 2 +0 2.5 +0.5 2.5 +1 2.5 +1.5 2.5 +2.5 2.5 +3 2.5 +3.5 2.5 +4 2.5 +4.5 2.5 +5 2.5 +5.5 2.5 +6 2.5 +6.5 2.5 +7 2.5 +7.5 2.5 +8 2.5 +8.5 2.5 +9 2.5 +9.5 2.5 +10 2.5 +10.5 2.5 +11 2.5 +11.5 2.5 +12 2.5 +12.5 2.5 +13 2.5 +13.5 2.5 +14 2.5 +14.5 2.5 +15 2.5 +15.5 2.5 +16 2.5 +16.5 2.5 +17 2.5 +17.5 2.5 +18 2.5 +18.5 2.5 +19 2.5 +19.5 2.5 +20 2.5 +20.5 2.5 +21 2.5 +21.5 2.5 +22 2.5 +22.5 2.5 +23 2.5 +23.5 2.5 +24 2.5 +24.5 2.5 +25 2.5 +2 2.5 +0 3.033333333 +0.5 3.033333333 +1 3.033333333 +1.5 3.033333333 +2 3.033333333 +2.5 3.033333333 +3 3.033333333 +3.5 3.033333333 +4 3.033333333 +4.5 3.033333333 +5 3.033333333 +5.5 3.033333333 +6 3.033333333 +6.5 3.033333333 +7 3.033333333 +7.5 3.033333333 +8 3.033333333 +8.5 3.033333333 +9 3.033333333 +9.5 3.033333333 +10 3.033333333 +10.5 3.033333333 +11 3.033333333 +11.5 3.033333333 +12 3.033333333 +12.5 3.033333333 +13 3.033333333 +13.5 3.033333333 +14 3.033333333 +14.5 3.033333333 +15 3.033333333 +15.5 3.033333333 +16 3.033333333 +16.5 3.033333333 +17 3.033333333 +17.5 3.033333333 +18 3.033333333 +18.5 3.033333333 +19 3.033333333 +19.5 3.033333333 +20 3.033333333 +20.5 3.033333333 +21 3.033333333 +21.5 3.033333333 +22 3.033333333 +22.5 3.033333333 +23 3.033333333 +23.5 3.033333333 +24 3.033333333 +24.5 3.033333333 +25 3.033333333 +0 3.566666667 +0.5 3.566666667 +1 3.566666667 +1.5 3.566666667 +2 3.566666667 +2.5 3.566666667 +3 3.566666667 +3.5 3.566666667 +4 3.566666667 +4.5 3.566666667 +5 3.566666667 +5.5 3.566666667 +6 3.566666667 +6.5 3.566666667 +7 3.566666667 +7.5 3.566666667 +8 3.566666667 +8.5 3.566666667 +9 3.566666667 +9.5 3.566666667 +10 3.566666667 +10.5 3.566666667 +11 3.566666667 +11.5 3.566666667 +12 3.566666667 +12.5 3.566666667 +13 3.566666667 +13.5 3.566666667 +14 3.566666667 +14.5 3.566666667 +15 3.566666667 +15.5 3.566666667 +16 3.566666667 +16.5 3.566666667 +17 3.566666667 +17.5 3.566666667 +18 3.566666667 +18.5 3.566666667 +19 3.566666667 +19.5 3.566666667 +20 3.566666667 +20.5 3.566666667 +21 3.566666667 +21.5 3.566666667 +22 3.566666667 +22.5 3.566666667 +23 3.566666667 +23.5 3.566666667 +24 3.566666667 +24.5 3.566666667 +25 3.566666667 +0 4.1 +0.5 4.1 +1 4.1 +1.5 4.1 +2 4.1 +2.5 4.1 +3 4.1 +3.5 4.1 +4 4.1 +4.5 4.1 +5 4.1 +5.5 4.1 +6 4.1 +6.5 4.1 +7 4.1 +7.5 4.1 +8 4.1 +8.5 4.1 +9 4.1 +9.5 4.1 +10 4.1 +10.5 4.1 +11 4.1 +11.5 4.1 +12 4.1 +12.5 4.1 +13 4.1 +13.5 4.1 +14 4.1 +14.5 4.1 +15 4.1 +15.5 4.1 +16 4.1 +16.5 4.1 +17 4.1 +17.5 4.1 +18 4.1 +18.5 4.1 +19 4.1 +19.5 4.1 +20 4.1 +20.5 4.1 +21 4.1 +21.5 4.1 +22 4.1 +22.5 4.1 +23 4.1 +23.5 4.1 +24 4.1 +24.5 4.1 +25 4.1 +0.25 0.25 +0.75 0.25 +1.25 0.25 +1.75 0.25 +2.25 0.25 +2.75 0.25 +3.25 0.25 +3.75 0.25 +4.25 0.25 +4.75 0.25 +5.25 0.25 +5.75 0.25 +6.25 0.25 +6.75 0.25 +7.25 0.25 +7.75 0.25 +8.25 0.25 +8.75 0.25 +9.25 0.25 +9.75 0.25 +10.25 0.25 +10.75 0.25 +11.25 0.25 +11.75 0.25 +12.25 0.25 +12.75 0.25 +13.25 0.25 +13.75 0.25 +14.25 0.25 +14.75 0.25 +15.25 0.25 +15.75 0.25 +16.25 0.25 +16.75 0.25 +17.25 0.25 +17.75 0.25 +18.25 0.25 +18.75 0.25 +19.25 0.25 +19.75 0.25 +20.25 0.25 +20.75 0.25 +21.25 0.25 +21.75 0.25 +22.25 0.25 +22.75 0.25 +23.25 0.25 +23.75 0.25 +24.25 0.25 +24.75 0.25 +0.25 0.75 +0.75 0.75 +1.25 0.75 +1.75 0.75 +2.25 0.75 +2.75 0.75 +3.25 0.75 +3.75 0.75 +4.25 0.75 +4.75 0.75 +5.25 0.75 +5.75 0.75 +6.25 0.75 +6.75 0.75 +7.25 0.75 +7.75 0.75 +8.25 0.75 +8.75 0.75 +9.25 0.75 +9.75 0.75 +10.25 0.75 +10.75 0.75 +11.25 0.75 +11.75 0.75 +12.25 0.75 +12.75 0.75 +13.25 0.75 +13.75 0.75 +14.25 0.75 +14.75 0.75 +15.25 0.75 +15.75 0.75 +16.25 0.75 +16.75 0.75 +17.25 0.75 +17.75 0.75 +18.25 0.75 +18.75 0.75 +19.25 0.75 +19.75 0.75 +20.25 0.75 +20.75 0.75 +21.25 0.75 +21.75 0.75 +22.25 0.75 +22.75 0.75 +23.25 0.75 +23.75 0.75 +24.25 0.75 +24.75 0.75 +0.25 1.25 +0.75 1.25 +1.25 1.25 +1.75 1.25 +2.25 1.25 +2.75 1.25 +3.25 1.25 +3.75 1.25 +4.25 1.25 +4.75 1.25 +5.25 1.25 +5.75 1.25 +6.25 1.25 +6.75 1.25 +7.25 1.25 +7.75 1.25 +8.25 1.25 +8.75 1.25 +9.25 1.25 +9.75 1.25 +10.25 1.25 +10.75 1.25 +11.25 1.25 +11.75 1.25 +12.25 1.25 +12.75 1.25 +13.25 1.25 +13.75 1.25 +14.25 1.25 +14.75 1.25 +15.25 1.25 +15.75 1.25 +16.25 1.25 +16.75 1.25 +17.25 1.25 +17.75 1.25 +18.25 1.25 +18.75 1.25 +19.25 1.25 +19.75 1.25 +20.25 1.25 +20.75 1.25 +21.25 1.25 +21.75 1.25 +22.25 1.25 +22.75 1.25 +23.25 1.25 +23.75 1.25 +24.25 1.25 +24.75 1.25 +0.25 1.75 +0.75 1.75 +1.25 1.75 +2.75 1.75 +3.25 1.75 +3.75 1.75 +4.25 1.75 +4.75 1.75 +5.25 1.75 +5.75 1.75 +6.25 1.75 +6.75 1.75 +7.25 1.75 +7.75 1.75 +8.25 1.75 +8.75 1.75 +9.25 1.75 +9.75 1.75 +10.25 1.75 +10.75 1.75 +11.25 1.75 +11.75 1.75 +12.25 1.75 +12.75 1.75 +13.25 1.75 +13.75 1.75 +14.25 1.75 +14.75 1.75 +15.25 1.75 +15.75 1.75 +16.25 1.75 +16.75 1.75 +17.25 1.75 +17.75 1.75 +18.25 1.75 +18.75 1.75 +19.25 1.75 +19.75 1.75 +20.25 1.75 +20.75 1.75 +21.25 1.75 +21.75 1.75 +22.25 1.75 +22.75 1.75 +23.25 1.75 +23.75 1.75 +24.25 1.75 +24.75 1.75 +0.25 2.25 +0.75 2.25 +1.25 2.25 +2.75 2.25 +3.25 2.25 +3.75 2.25 +4.25 2.25 +4.75 2.25 +5.25 2.25 +5.75 2.25 +6.25 2.25 +6.75 2.25 +7.25 2.25 +7.75 2.25 +8.25 2.25 +8.75 2.25 +9.25 2.25 +9.75 2.25 +10.25 2.25 +10.75 2.25 +11.25 2.25 +11.75 2.25 +12.25 2.25 +12.75 2.25 +13.25 2.25 +13.75 2.25 +14.25 2.25 +14.75 2.25 +15.25 2.25 +15.75 2.25 +16.25 2.25 +16.75 2.25 +17.25 2.25 +17.75 2.25 +18.25 2.25 +18.75 2.25 +19.25 2.25 +19.75 2.25 +20.25 2.25 +20.75 2.25 +21.25 2.25 +21.75 2.25 +22.25 2.25 +22.75 2.25 +23.25 2.25 +23.75 2.25 +24.25 2.25 +24.75 2.25 +0.25 2.766666667 +0.75 2.766666667 +1.25 2.766666667 +1.75 2.766666667 +2.25 2.766666667 +2.75 2.766666667 +3.25 2.766666667 +3.75 2.766666667 +4.25 2.766666667 +4.75 2.766666667 +5.25 2.766666667 +5.75 2.766666667 +6.25 2.766666667 +6.75 2.766666667 +7.25 2.766666667 +7.75 2.766666667 +8.25 2.766666667 +8.75 2.766666667 +9.25 2.766666667 +9.75 2.766666667 +10.25 2.766666667 +10.75 2.766666667 +11.25 2.766666667 +11.75 2.766666667 +12.25 2.766666667 +12.75 2.766666667 +13.25 2.766666667 +13.75 2.766666667 +14.25 2.766666667 +14.75 2.766666667 +15.25 2.766666667 +15.75 2.766666667 +16.25 2.766666667 +16.75 2.766666667 +17.25 2.766666667 +17.75 2.766666667 +18.25 2.766666667 +18.75 2.766666667 +19.25 2.766666667 +19.75 2.766666667 +20.25 2.766666667 +20.75 2.766666667 +21.25 2.766666667 +21.75 2.766666667 +22.25 2.766666667 +22.75 2.766666667 +23.25 2.766666667 +23.75 2.766666667 +24.25 2.766666667 +24.75 2.766666667 +0.25 3.3 +0.75 3.3 +1.25 3.3 +1.75 3.3 +2.25 3.3 +2.75 3.3 +3.25 3.3 +3.75 3.3 +4.25 3.3 +4.75 3.3 +5.25 3.3 +5.75 3.3 +6.25 3.3 +6.75 3.3 +7.25 3.3 +7.75 3.3 +8.25 3.3 +8.75 3.3 +9.25 3.3 +9.75 3.3 +10.25 3.3 +10.75 3.3 +11.25 3.3 +11.75 3.3 +12.25 3.3 +12.75 3.3 +13.25 3.3 +13.75 3.3 +14.25 3.3 +14.75 3.3 +15.25 3.3 +15.75 3.3 +16.25 3.3 +16.75 3.3 +17.25 3.3 +17.75 3.3 +18.25 3.3 +18.75 3.3 +19.25 3.3 +19.75 3.3 +20.25 3.3 +20.75 3.3 +21.25 3.3 +21.75 3.3 +22.25 3.3 +22.75 3.3 +23.25 3.3 +23.75 3.3 +24.25 3.3 +24.75 3.3 +0.25 3.833333333 +0.75 3.833333333 +1.25 3.833333333 +1.75 3.833333333 +2.25 3.833333333 +2.75 3.833333333 +3.25 3.833333333 +3.75 3.833333333 +4.25 3.833333333 +4.75 3.833333333 +5.25 3.833333333 +5.75 3.833333333 +6.25 3.833333333 +6.75 3.833333333 +7.25 3.833333333 +7.75 3.833333333 +8.25 3.833333333 +8.75 3.833333333 +9.25 3.833333333 +9.75 3.833333333 +10.25 3.833333333 +10.75 3.833333333 +11.25 3.833333333 +11.75 3.833333333 +12.25 3.833333333 +12.75 3.833333333 +13.25 3.833333333 +13.75 3.833333333 +14.25 3.833333333 +14.75 3.833333333 +15.25 3.833333333 +15.75 3.833333333 +16.25 3.833333333 +16.75 3.833333333 +17.25 3.833333333 +17.75 3.833333333 +18.25 3.833333333 +18.75 3.833333333 +19.25 3.833333333 +19.75 3.833333333 +20.25 3.833333333 +20.75 3.833333333 +21.25 3.833333333 +21.75 3.833333333 +22.25 3.833333333 +22.75 3.833333333 +23.25 3.833333333 +23.75 3.833333333 +24.25 3.833333333 +24.75 3.833333333 +1.75 2.25 +2.25 2.25 +1.75 1.75 +2.25 1.75 + +element vertices: +51 0 458 +52 51 458 +1 52 458 +0 1 458 +52 1 459 +53 52 459 +2 53 459 +1 2 459 +53 2 460 +54 53 460 +3 54 460 +2 3 460 +54 3 461 +55 54 461 +4 55 461 +3 4 461 +55 4 462 +56 55 462 +5 56 462 +4 5 462 +56 5 463 +57 56 463 +6 57 463 +5 6 463 +57 6 464 +58 57 464 +7 58 464 +6 7 464 +58 7 465 +59 58 465 +8 59 465 +7 8 465 +59 8 466 +60 59 466 +9 60 466 +8 9 466 +60 9 467 +61 60 467 +10 61 467 +9 10 467 +61 10 468 +62 61 468 +11 62 468 +10 11 468 +62 11 469 +63 62 469 +12 63 469 +11 12 469 +63 12 470 +64 63 470 +13 64 470 +12 13 470 +64 13 471 +65 64 471 +14 65 471 +13 14 471 +65 14 472 +66 65 472 +15 66 472 +14 15 472 +66 15 473 +67 66 473 +16 67 473 +15 16 473 +67 16 474 +68 67 474 +17 68 474 +16 17 474 +68 17 475 +69 68 475 +18 69 475 +17 18 475 +69 18 476 +70 69 476 +19 70 476 +18 19 476 +70 19 477 +71 70 477 +20 71 477 +19 20 477 +71 20 478 +72 71 478 +21 72 478 +20 21 478 +72 21 479 +73 72 479 +22 73 479 +21 22 479 +73 22 480 +74 73 480 +23 74 480 +22 23 480 +74 23 481 +75 74 481 +24 75 481 +23 24 481 +75 24 482 +76 75 482 +25 76 482 +24 25 482 +76 25 483 +77 76 483 +26 77 483 +25 26 483 +77 26 484 +78 77 484 +27 78 484 +26 27 484 +78 27 485 +79 78 485 +28 79 485 +27 28 485 +79 28 486 +80 79 486 +29 80 486 +28 29 486 +80 29 487 +81 80 487 +30 81 487 +29 30 487 +81 30 488 +82 81 488 +31 82 488 +30 31 488 +82 31 489 +83 82 489 +32 83 489 +31 32 489 +83 32 490 +84 83 490 +33 84 490 +32 33 490 +84 33 491 +85 84 491 +34 85 491 +33 34 491 +85 34 492 +86 85 492 +35 86 492 +34 35 492 +86 35 493 +87 86 493 +36 87 493 +35 36 493 +87 36 494 +88 87 494 +37 88 494 +36 37 494 +88 37 495 +89 88 495 +38 89 495 +37 38 495 +89 38 496 +90 89 496 +39 90 496 +38 39 496 +90 39 497 +91 90 497 +40 91 497 +39 40 497 +91 40 498 +92 91 498 +41 92 498 +40 41 498 +92 41 499 +93 92 499 +42 93 499 +41 42 499 +93 42 500 +94 93 500 +43 94 500 +42 43 500 +94 43 501 +95 94 501 +44 95 501 +43 44 501 +95 44 502 +96 95 502 +45 96 502 +44 45 502 +96 45 503 +97 96 503 +46 97 503 +45 46 503 +97 46 504 +98 97 504 +47 98 504 +46 47 504 +98 47 505 +99 98 505 +48 99 505 +47 48 505 +99 48 506 +100 99 506 +49 100 506 +48 49 506 +100 49 507 +101 100 507 +50 101 507 +49 50 507 +102 51 508 +103 102 508 +52 103 508 +51 52 508 +103 52 509 +104 103 509 +53 104 509 +52 53 509 +104 53 510 +105 104 510 +54 105 510 +53 54 510 +105 54 511 +106 105 511 +55 106 511 +54 55 511 +106 55 512 +107 106 512 +56 107 512 +55 56 512 +107 56 513 +108 107 513 +57 108 513 +56 57 513 +108 57 514 +109 108 514 +58 109 514 +57 58 514 +109 58 515 +110 109 515 +59 110 515 +58 59 515 +110 59 516 +111 110 516 +60 111 516 +59 60 516 +111 60 517 +112 111 517 +61 112 517 +60 61 517 +112 61 518 +113 112 518 +62 113 518 +61 62 518 +113 62 519 +114 113 519 +63 114 519 +62 63 519 +114 63 520 +115 114 520 +64 115 520 +63 64 520 +115 64 521 +116 115 521 +65 116 521 +64 65 521 +116 65 522 +117 116 522 +66 117 522 +65 66 522 +117 66 523 +118 117 523 +67 118 523 +66 67 523 +118 67 524 +119 118 524 +68 119 524 +67 68 524 +119 68 525 +120 119 525 +69 120 525 +68 69 525 +120 69 526 +121 120 526 +70 121 526 +69 70 526 +121 70 527 +122 121 527 +71 122 527 +70 71 527 +122 71 528 +123 122 528 +72 123 528 +71 72 528 +123 72 529 +124 123 529 +73 124 529 +72 73 529 +124 73 530 +125 124 530 +74 125 530 +73 74 530 +125 74 531 +126 125 531 +75 126 531 +74 75 531 +126 75 532 +127 126 532 +76 127 532 +75 76 532 +127 76 533 +128 127 533 +77 128 533 +76 77 533 +128 77 534 +129 128 534 +78 129 534 +77 78 534 +129 78 535 +130 129 535 +79 130 535 +78 79 535 +130 79 536 +131 130 536 +80 131 536 +79 80 536 +131 80 537 +132 131 537 +81 132 537 +80 81 537 +132 81 538 +133 132 538 +82 133 538 +81 82 538 +133 82 539 +134 133 539 +83 134 539 +82 83 539 +134 83 540 +135 134 540 +84 135 540 +83 84 540 +135 84 541 +136 135 541 +85 136 541 +84 85 541 +136 85 542 +137 136 542 +86 137 542 +85 86 542 +137 86 543 +138 137 543 +87 138 543 +86 87 543 +138 87 544 +139 138 544 +88 139 544 +87 88 544 +139 88 545 +140 139 545 +89 140 545 +88 89 545 +140 89 546 +141 140 546 +90 141 546 +89 90 546 +141 90 547 +142 141 547 +91 142 547 +90 91 547 +142 91 548 +143 142 548 +92 143 548 +91 92 548 +143 92 549 +144 143 549 +93 144 549 +92 93 549 +144 93 550 +145 144 550 +94 145 550 +93 94 550 +145 94 551 +146 145 551 +95 146 551 +94 95 551 +146 95 552 +147 146 552 +96 147 552 +95 96 552 +147 96 553 +148 147 553 +97 148 553 +96 97 553 +148 97 554 +149 148 554 +98 149 554 +97 98 554 +149 98 555 +150 149 555 +99 150 555 +98 99 555 +150 99 556 +151 150 556 +100 151 556 +99 100 556 +151 100 557 +152 151 557 +101 152 557 +100 101 557 +153 102 558 +154 153 558 +103 154 558 +102 103 558 +154 103 559 +155 154 559 +104 155 559 +103 104 559 +155 104 560 +156 155 560 +105 156 560 +104 105 560 +156 105 561 +157 156 561 +106 157 561 +105 106 561 +157 106 562 +158 157 562 +107 158 562 +106 107 562 +158 107 563 +159 158 563 +108 159 563 +107 108 563 +159 108 564 +160 159 564 +109 160 564 +108 109 564 +160 109 565 +161 160 565 +110 161 565 +109 110 565 +161 110 566 +162 161 566 +111 162 566 +110 111 566 +162 111 567 +163 162 567 +112 163 567 +111 112 567 +163 112 568 +164 163 568 +113 164 568 +112 113 568 +164 113 569 +165 164 569 +114 165 569 +113 114 569 +165 114 570 +166 165 570 +115 166 570 +114 115 570 +166 115 571 +167 166 571 +116 167 571 +115 116 571 +167 116 572 +168 167 572 +117 168 572 +116 117 572 +168 117 573 +169 168 573 +118 169 573 +117 118 573 +169 118 574 +170 169 574 +119 170 574 +118 119 574 +170 119 575 +171 170 575 +120 171 575 +119 120 575 +171 120 576 +172 171 576 +121 172 576 +120 121 576 +172 121 577 +173 172 577 +122 173 577 +121 122 577 +173 122 578 +174 173 578 +123 174 578 +122 123 578 +174 123 579 +175 174 579 +124 175 579 +123 124 579 +175 124 580 +176 175 580 +125 176 580 +124 125 580 +176 125 581 +177 176 581 +126 177 581 +125 126 581 +177 126 582 +178 177 582 +127 178 582 +126 127 582 +178 127 583 +179 178 583 +128 179 583 +127 128 583 +179 128 584 +180 179 584 +129 180 584 +128 129 584 +180 129 585 +181 180 585 +130 181 585 +129 130 585 +181 130 586 +182 181 586 +131 182 586 +130 131 586 +182 131 587 +183 182 587 +132 183 587 +131 132 587 +183 132 588 +184 183 588 +133 184 588 +132 133 588 +184 133 589 +185 184 589 +134 185 589 +133 134 589 +185 134 590 +186 185 590 +135 186 590 +134 135 590 +186 135 591 +187 186 591 +136 187 591 +135 136 591 +187 136 592 +188 187 592 +137 188 592 +136 137 592 +188 137 593 +189 188 593 +138 189 593 +137 138 593 +189 138 594 +190 189 594 +139 190 594 +138 139 594 +190 139 595 +191 190 595 +140 191 595 +139 140 595 +191 140 596 +192 191 596 +141 192 596 +140 141 596 +192 141 597 +193 192 597 +142 193 597 +141 142 597 +193 142 598 +194 193 598 +143 194 598 +142 143 598 +194 143 599 +195 194 599 +144 195 599 +143 144 599 +195 144 600 +196 195 600 +145 196 600 +144 145 600 +196 145 601 +197 196 601 +146 197 601 +145 146 601 +197 146 602 +198 197 602 +147 198 602 +146 147 602 +198 147 603 +199 198 603 +148 199 603 +147 148 603 +199 148 604 +200 199 604 +149 200 604 +148 149 604 +200 149 605 +201 200 605 +150 201 605 +149 150 605 +201 150 606 +202 201 606 +151 202 606 +150 151 606 +202 151 607 +203 202 607 +152 203 607 +151 152 607 +204 153 608 +205 204 608 +154 205 608 +153 154 608 +205 154 609 +206 205 609 +155 206 609 +154 155 609 +206 155 610 +207 206 610 +156 207 610 +155 156 610 +208 158 611 +209 208 611 +159 209 611 +158 159 611 +209 159 612 +210 209 612 +160 210 612 +159 160 612 +210 160 613 +211 210 613 +161 211 613 +160 161 613 +211 161 614 +212 211 614 +162 212 614 +161 162 614 +212 162 615 +213 212 615 +163 213 615 +162 163 615 +213 163 616 +214 213 616 +164 214 616 +163 164 616 +214 164 617 +215 214 617 +165 215 617 +164 165 617 +215 165 618 +216 215 618 +166 216 618 +165 166 618 +216 166 619 +217 216 619 +167 217 619 +166 167 619 +217 167 620 +218 217 620 +168 218 620 +167 168 620 +218 168 621 +219 218 621 +169 219 621 +168 169 621 +219 169 622 +220 219 622 +170 220 622 +169 170 622 +220 170 623 +221 220 623 +171 221 623 +170 171 623 +221 171 624 +222 221 624 +172 222 624 +171 172 624 +222 172 625 +223 222 625 +173 223 625 +172 173 625 +223 173 626 +224 223 626 +174 224 626 +173 174 626 +224 174 627 +225 224 627 +175 225 627 +174 175 627 +225 175 628 +226 225 628 +176 226 628 +175 176 628 +226 176 629 +227 226 629 +177 227 629 +176 177 629 +227 177 630 +228 227 630 +178 228 630 +177 178 630 +228 178 631 +229 228 631 +179 229 631 +178 179 631 +229 179 632 +230 229 632 +180 230 632 +179 180 632 +230 180 633 +231 230 633 +181 231 633 +180 181 633 +231 181 634 +232 231 634 +182 232 634 +181 182 634 +232 182 635 +233 232 635 +183 233 635 +182 183 635 +233 183 636 +234 233 636 +184 234 636 +183 184 636 +234 184 637 +235 234 637 +185 235 637 +184 185 637 +235 185 638 +236 235 638 +186 236 638 +185 186 638 +236 186 639 +237 236 639 +187 237 639 +186 187 639 +237 187 640 +238 237 640 +188 238 640 +187 188 640 +238 188 641 +239 238 641 +189 239 641 +188 189 641 +239 189 642 +240 239 642 +190 240 642 +189 190 642 +240 190 643 +241 240 643 +191 241 643 +190 191 643 +241 191 644 +242 241 644 +192 242 644 +191 192 644 +242 192 645 +243 242 645 +193 243 645 +192 193 645 +243 193 646 +244 243 646 +194 244 646 +193 194 646 +244 194 647 +245 244 647 +195 245 647 +194 195 647 +245 195 648 +246 245 648 +196 246 648 +195 196 648 +246 196 649 +247 246 649 +197 247 649 +196 197 649 +247 197 650 +248 247 650 +198 248 650 +197 198 650 +248 198 651 +249 248 651 +199 249 651 +198 199 651 +249 199 652 +250 249 652 +200 250 652 +199 200 652 +250 200 653 +251 250 653 +201 251 653 +200 201 653 +251 201 654 +252 251 654 +202 252 654 +201 202 654 +252 202 655 +253 252 655 +203 253 655 +202 203 655 +254 204 656 +255 254 656 +205 255 656 +204 205 656 +255 205 657 +256 255 657 +206 256 657 +205 206 657 +256 206 658 +257 256 658 +207 257 658 +206 207 658 +258 208 659 +259 258 659 +209 259 659 +208 209 659 +259 209 660 +260 259 660 +210 260 660 +209 210 660 +260 210 661 +261 260 661 +211 261 661 +210 211 661 +261 211 662 +262 261 662 +212 262 662 +211 212 662 +262 212 663 +263 262 663 +213 263 663 +212 213 663 +263 213 664 +264 263 664 +214 264 664 +213 214 664 +264 214 665 +265 264 665 +215 265 665 +214 215 665 +265 215 666 +266 265 666 +216 266 666 +215 216 666 +266 216 667 +267 266 667 +217 267 667 +216 217 667 +267 217 668 +268 267 668 +218 268 668 +217 218 668 +268 218 669 +269 268 669 +219 269 669 +218 219 669 +269 219 670 +270 269 670 +220 270 670 +219 220 670 +270 220 671 +271 270 671 +221 271 671 +220 221 671 +271 221 672 +272 271 672 +222 272 672 +221 222 672 +272 222 673 +273 272 673 +223 273 673 +222 223 673 +273 223 674 +274 273 674 +224 274 674 +223 224 674 +274 224 675 +275 274 675 +225 275 675 +224 225 675 +275 225 676 +276 275 676 +226 276 676 +225 226 676 +276 226 677 +277 276 677 +227 277 677 +226 227 677 +277 227 678 +278 277 678 +228 278 678 +227 228 678 +278 228 679 +279 278 679 +229 279 679 +228 229 679 +279 229 680 +280 279 680 +230 280 680 +229 230 680 +280 230 681 +281 280 681 +231 281 681 +230 231 681 +281 231 682 +282 281 682 +232 282 682 +231 232 682 +282 232 683 +283 282 683 +233 283 683 +232 233 683 +283 233 684 +284 283 684 +234 284 684 +233 234 684 +284 234 685 +285 284 685 +235 285 685 +234 235 685 +285 235 686 +286 285 686 +236 286 686 +235 236 686 +286 236 687 +287 286 687 +237 287 687 +236 237 687 +287 237 688 +288 287 688 +238 288 688 +237 238 688 +288 238 689 +289 288 689 +239 289 689 +238 239 689 +289 239 690 +290 289 690 +240 290 690 +239 240 690 +290 240 691 +291 290 691 +241 291 691 +240 241 691 +291 241 692 +292 291 692 +242 292 692 +241 242 692 +292 242 693 +293 292 693 +243 293 693 +242 243 693 +293 243 694 +294 293 694 +244 294 694 +243 244 694 +294 244 695 +295 294 695 +245 295 695 +244 245 695 +295 245 696 +296 295 696 +246 296 696 +245 246 696 +296 246 697 +297 296 697 +247 297 697 +246 247 697 +297 247 698 +298 297 698 +248 298 698 +247 248 698 +298 248 699 +299 298 699 +249 299 699 +248 249 699 +299 249 700 +300 299 700 +250 300 700 +249 250 700 +300 250 701 +301 300 701 +251 301 701 +250 251 701 +301 251 702 +302 301 702 +252 302 702 +251 252 702 +302 252 703 +303 302 703 +253 303 703 +252 253 703 +305 254 704 +306 305 704 +255 306 704 +254 255 704 +306 255 705 +307 306 705 +256 307 705 +255 256 705 +307 256 706 +308 307 706 +257 308 706 +256 257 706 +308 257 707 +309 308 707 +304 309 707 +257 304 707 +309 304 708 +310 309 708 +258 310 708 +304 258 708 +310 258 709 +311 310 709 +259 311 709 +258 259 709 +311 259 710 +312 311 710 +260 312 710 +259 260 710 +312 260 711 +313 312 711 +261 313 711 +260 261 711 +313 261 712 +314 313 712 +262 314 712 +261 262 712 +314 262 713 +315 314 713 +263 315 713 +262 263 713 +315 263 714 +316 315 714 +264 316 714 +263 264 714 +316 264 715 +317 316 715 +265 317 715 +264 265 715 +317 265 716 +318 317 716 +266 318 716 +265 266 716 +318 266 717 +319 318 717 +267 319 717 +266 267 717 +319 267 718 +320 319 718 +268 320 718 +267 268 718 +320 268 719 +321 320 719 +269 321 719 +268 269 719 +321 269 720 +322 321 720 +270 322 720 +269 270 720 +322 270 721 +323 322 721 +271 323 721 +270 271 721 +323 271 722 +324 323 722 +272 324 722 +271 272 722 +324 272 723 +325 324 723 +273 325 723 +272 273 723 +325 273 724 +326 325 724 +274 326 724 +273 274 724 +326 274 725 +327 326 725 +275 327 725 +274 275 725 +327 275 726 +328 327 726 +276 328 726 +275 276 726 +328 276 727 +329 328 727 +277 329 727 +276 277 727 +329 277 728 +330 329 728 +278 330 728 +277 278 728 +330 278 729 +331 330 729 +279 331 729 +278 279 729 +331 279 730 +332 331 730 +280 332 730 +279 280 730 +332 280 731 +333 332 731 +281 333 731 +280 281 731 +333 281 732 +334 333 732 +282 334 732 +281 282 732 +334 282 733 +335 334 733 +283 335 733 +282 283 733 +335 283 734 +336 335 734 +284 336 734 +283 284 734 +336 284 735 +337 336 735 +285 337 735 +284 285 735 +337 285 736 +338 337 736 +286 338 736 +285 286 736 +338 286 737 +339 338 737 +287 339 737 +286 287 737 +339 287 738 +340 339 738 +288 340 738 +287 288 738 +340 288 739 +341 340 739 +289 341 739 +288 289 739 +341 289 740 +342 341 740 +290 342 740 +289 290 740 +342 290 741 +343 342 741 +291 343 741 +290 291 741 +343 291 742 +344 343 742 +292 344 742 +291 292 742 +344 292 743 +345 344 743 +293 345 743 +292 293 743 +345 293 744 +346 345 744 +294 346 744 +293 294 744 +346 294 745 +347 346 745 +295 347 745 +294 295 745 +347 295 746 +348 347 746 +296 348 746 +295 296 746 +348 296 747 +349 348 747 +297 349 747 +296 297 747 +349 297 748 +350 349 748 +298 350 748 +297 298 748 +350 298 749 +351 350 749 +299 351 749 +298 299 749 +351 299 750 +352 351 750 +300 352 750 +299 300 750 +352 300 751 +353 352 751 +301 353 751 +300 301 751 +353 301 752 +354 353 752 +302 354 752 +301 302 752 +354 302 753 +355 354 753 +303 355 753 +302 303 753 +356 305 754 +357 356 754 +306 357 754 +305 306 754 +357 306 755 +358 357 755 +307 358 755 +306 307 755 +358 307 756 +359 358 756 +308 359 756 +307 308 756 +359 308 757 +360 359 757 +309 360 757 +308 309 757 +360 309 758 +361 360 758 +310 361 758 +309 310 758 +361 310 759 +362 361 759 +311 362 759 +310 311 759 +362 311 760 +363 362 760 +312 363 760 +311 312 760 +363 312 761 +364 363 761 +313 364 761 +312 313 761 +364 313 762 +365 364 762 +314 365 762 +313 314 762 +365 314 763 +366 365 763 +315 366 763 +314 315 763 +366 315 764 +367 366 764 +316 367 764 +315 316 764 +367 316 765 +368 367 765 +317 368 765 +316 317 765 +368 317 766 +369 368 766 +318 369 766 +317 318 766 +369 318 767 +370 369 767 +319 370 767 +318 319 767 +370 319 768 +371 370 768 +320 371 768 +319 320 768 +371 320 769 +372 371 769 +321 372 769 +320 321 769 +372 321 770 +373 372 770 +322 373 770 +321 322 770 +373 322 771 +374 373 771 +323 374 771 +322 323 771 +374 323 772 +375 374 772 +324 375 772 +323 324 772 +375 324 773 +376 375 773 +325 376 773 +324 325 773 +376 325 774 +377 376 774 +326 377 774 +325 326 774 +377 326 775 +378 377 775 +327 378 775 +326 327 775 +378 327 776 +379 378 776 +328 379 776 +327 328 776 +379 328 777 +380 379 777 +329 380 777 +328 329 777 +380 329 778 +381 380 778 +330 381 778 +329 330 778 +381 330 779 +382 381 779 +331 382 779 +330 331 779 +382 331 780 +383 382 780 +332 383 780 +331 332 780 +383 332 781 +384 383 781 +333 384 781 +332 333 781 +384 333 782 +385 384 782 +334 385 782 +333 334 782 +385 334 783 +386 385 783 +335 386 783 +334 335 783 +386 335 784 +387 386 784 +336 387 784 +335 336 784 +387 336 785 +388 387 785 +337 388 785 +336 337 785 +388 337 786 +389 388 786 +338 389 786 +337 338 786 +389 338 787 +390 389 787 +339 390 787 +338 339 787 +390 339 788 +391 390 788 +340 391 788 +339 340 788 +391 340 789 +392 391 789 +341 392 789 +340 341 789 +392 341 790 +393 392 790 +342 393 790 +341 342 790 +393 342 791 +394 393 791 +343 394 791 +342 343 791 +394 343 792 +395 394 792 +344 395 792 +343 344 792 +395 344 793 +396 395 793 +345 396 793 +344 345 793 +396 345 794 +397 396 794 +346 397 794 +345 346 794 +397 346 795 +398 397 795 +347 398 795 +346 347 795 +398 347 796 +399 398 796 +348 399 796 +347 348 796 +399 348 797 +400 399 797 +349 400 797 +348 349 797 +400 349 798 +401 400 798 +350 401 798 +349 350 798 +401 350 799 +402 401 799 +351 402 799 +350 351 799 +402 351 800 +403 402 800 +352 403 800 +351 352 800 +403 352 801 +404 403 801 +353 404 801 +352 353 801 +404 353 802 +405 404 802 +354 405 802 +353 354 802 +405 354 803 +406 405 803 +355 406 803 +354 355 803 +407 356 804 +408 407 804 +357 408 804 +356 357 804 +408 357 805 +409 408 805 +358 409 805 +357 358 805 +409 358 806 +410 409 806 +359 410 806 +358 359 806 +410 359 807 +411 410 807 +360 411 807 +359 360 807 +411 360 808 +412 411 808 +361 412 808 +360 361 808 +412 361 809 +413 412 809 +362 413 809 +361 362 809 +413 362 810 +414 413 810 +363 414 810 +362 363 810 +414 363 811 +415 414 811 +364 415 811 +363 364 811 +415 364 812 +416 415 812 +365 416 812 +364 365 812 +416 365 813 +417 416 813 +366 417 813 +365 366 813 +417 366 814 +418 417 814 +367 418 814 +366 367 814 +418 367 815 +419 418 815 +368 419 815 +367 368 815 +419 368 816 +420 419 816 +369 420 816 +368 369 816 +420 369 817 +421 420 817 +370 421 817 +369 370 817 +421 370 818 +422 421 818 +371 422 818 +370 371 818 +422 371 819 +423 422 819 +372 423 819 +371 372 819 +423 372 820 +424 423 820 +373 424 820 +372 373 820 +424 373 821 +425 424 821 +374 425 821 +373 374 821 +425 374 822 +426 425 822 +375 426 822 +374 375 822 +426 375 823 +427 426 823 +376 427 823 +375 376 823 +427 376 824 +428 427 824 +377 428 824 +376 377 824 +428 377 825 +429 428 825 +378 429 825 +377 378 825 +429 378 826 +430 429 826 +379 430 826 +378 379 826 +430 379 827 +431 430 827 +380 431 827 +379 380 827 +431 380 828 +432 431 828 +381 432 828 +380 381 828 +432 381 829 +433 432 829 +382 433 829 +381 382 829 +433 382 830 +434 433 830 +383 434 830 +382 383 830 +434 383 831 +435 434 831 +384 435 831 +383 384 831 +435 384 832 +436 435 832 +385 436 832 +384 385 832 +436 385 833 +437 436 833 +386 437 833 +385 386 833 +437 386 834 +438 437 834 +387 438 834 +386 387 834 +438 387 835 +439 438 835 +388 439 835 +387 388 835 +439 388 836 +440 439 836 +389 440 836 +388 389 836 +440 389 837 +441 440 837 +390 441 837 +389 390 837 +441 390 838 +442 441 838 +391 442 838 +390 391 838 +442 391 839 +443 442 839 +392 443 839 +391 392 839 +443 392 840 +444 443 840 +393 444 840 +392 393 840 +444 393 841 +445 444 841 +394 445 841 +393 394 841 +445 394 842 +446 445 842 +395 446 842 +394 395 842 +446 395 843 +447 446 843 +396 447 843 +395 396 843 +447 396 844 +448 447 844 +397 448 844 +396 397 844 +448 397 845 +449 448 845 +398 449 845 +397 398 845 +449 398 846 +450 449 846 +399 450 846 +398 399 846 +450 399 847 +451 450 847 +400 451 847 +399 400 847 +451 400 848 +452 451 848 +401 452 848 +400 401 848 +452 401 849 +453 452 849 +402 453 849 +401 402 849 +453 402 850 +454 453 850 +403 454 850 +402 403 850 +454 403 851 +455 454 851 +404 455 851 +403 404 851 +455 404 852 +456 455 852 +405 456 852 +404 405 852 +456 405 853 +457 456 853 +406 457 853 +405 406 853 +257 207 854 +304 257 854 +258 304 855 +208 258 855 +158 208 857 +157 158 857 +156 157 856 +207 156 856 + +element boundaries: +0 0 2 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 3 +0 0 4 +0 0 2 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 3 +0 0 0 +0 0 2 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 3 +0 0 0 +0 0 2 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 3 +0 0 0 +0 0 2 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 3 +0 0 0 +0 0 2 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 3 +0 0 0 +0 0 2 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 3 +0 0 0 +0 0 2 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 3 +0 0 0 +1 0 0 +0 1 0 +1 0 0 +0 1 0 +1 0 0 +0 1 0 +1 0 0 +0 1 0 + +element neighbours: +3 1 -1 +0 2 203 +1 3 4 +2 0 -1 +7 5 2 +4 6 207 +5 7 8 +6 4 -1 +11 9 6 +8 10 211 +9 11 12 +10 8 -1 +15 13 10 +12 14 215 +13 15 16 +14 12 -1 +19 17 14 +16 18 219 +17 19 20 +18 16 -1 +23 21 18 +20 22 223 +21 23 24 +22 20 -1 +27 25 22 +24 26 227 +25 27 28 +26 24 -1 +31 29 26 +28 30 231 +29 31 32 +30 28 -1 +35 33 30 +32 34 235 +33 35 36 +34 32 -1 +39 37 34 +36 38 239 +37 39 40 +38 36 -1 +43 41 38 +40 42 243 +41 43 44 +42 40 -1 +47 45 42 +44 46 247 +45 47 48 +46 44 -1 +51 49 46 +48 50 251 +49 51 52 +50 48 -1 +55 53 50 +52 54 255 +53 55 56 +54 52 -1 +59 57 54 +56 58 259 +57 59 60 +58 56 -1 +63 61 58 +60 62 263 +61 63 64 +62 60 -1 +67 65 62 +64 66 267 +65 67 68 +66 64 -1 +71 69 66 +68 70 271 +69 71 72 +70 68 -1 +75 73 70 +72 74 275 +73 75 76 +74 72 -1 +79 77 74 +76 78 279 +77 79 80 +78 76 -1 +83 81 78 +80 82 283 +81 83 84 +82 80 -1 +87 85 82 +84 86 287 +85 87 88 +86 84 -1 +91 89 86 +88 90 291 +89 91 92 +90 88 -1 +95 93 90 +92 94 295 +93 95 96 +94 92 -1 +99 97 94 +96 98 299 +97 99 100 +98 96 -1 +103 101 98 +100 102 303 +101 103 104 +102 100 -1 +107 105 102 +104 106 307 +105 107 108 +106 104 -1 +111 109 106 +108 110 311 +109 111 112 +110 108 -1 +115 113 110 +112 114 315 +113 115 116 +114 112 -1 +119 117 114 +116 118 319 +117 119 120 +118 116 -1 +123 121 118 +120 122 323 +121 123 124 +122 120 -1 +127 125 122 +124 126 327 +125 127 128 +126 124 -1 +131 129 126 +128 130 331 +129 131 132 +130 128 -1 +135 133 130 +132 134 335 +133 135 136 +134 132 -1 +139 137 134 +136 138 339 +137 139 140 +138 136 -1 +143 141 138 +140 142 343 +141 143 144 +142 140 -1 +147 145 142 +144 146 347 +145 147 148 +146 144 -1 +151 149 146 +148 150 351 +149 151 152 +150 148 -1 +155 153 150 +152 154 355 +153 155 156 +154 152 -1 +159 157 154 +156 158 359 +157 159 160 +158 156 -1 +163 161 158 +160 162 363 +161 163 164 +162 160 -1 +167 165 162 +164 166 367 +165 167 168 +166 164 -1 +171 169 166 +168 170 371 +169 171 172 +170 168 -1 +175 173 170 +172 174 375 +173 175 176 +174 172 -1 +179 177 174 +176 178 379 +177 179 180 +178 176 -1 +183 181 178 +180 182 383 +181 183 184 +182 180 -1 +187 185 182 +184 186 387 +185 187 188 +186 184 -1 +191 189 186 +188 190 391 +189 191 192 +190 188 -1 +195 193 190 +192 194 395 +193 195 196 +194 192 -1 +199 197 194 +196 198 399 +197 199 -1 +198 196 -1 +203 201 -1 +200 202 403 +201 203 204 +202 200 1 +207 205 202 +204 206 407 +205 207 208 +206 204 5 +211 209 206 +208 210 411 +209 211 212 +210 208 9 +215 213 210 +212 214 415 +213 215 216 +214 212 13 +219 217 214 +216 218 419 +217 219 220 +218 216 17 +223 221 218 +220 222 423 +221 223 224 +222 220 21 +227 225 222 +224 226 427 +225 227 228 +226 224 25 +231 229 226 +228 230 431 +229 231 232 +230 228 29 +235 233 230 +232 234 435 +233 235 236 +234 232 33 +239 237 234 +236 238 439 +237 239 240 +238 236 37 +243 241 238 +240 242 443 +241 243 244 +242 240 41 +247 245 242 +244 246 447 +245 247 248 +246 244 45 +251 249 246 +248 250 451 +249 251 252 +250 248 49 +255 253 250 +252 254 455 +253 255 256 +254 252 53 +259 257 254 +256 258 459 +257 259 260 +258 256 57 +263 261 258 +260 262 463 +261 263 264 +262 260 61 +267 265 262 +264 266 467 +265 267 268 +266 264 65 +271 269 266 +268 270 471 +269 271 272 +270 268 69 +275 273 270 +272 274 475 +273 275 276 +274 272 73 +279 277 274 +276 278 479 +277 279 280 +278 276 77 +283 281 278 +280 282 483 +281 283 284 +282 280 81 +287 285 282 +284 286 487 +285 287 288 +286 284 85 +291 289 286 +288 290 491 +289 291 292 +290 288 89 +295 293 290 +292 294 495 +293 295 296 +294 292 93 +299 297 294 +296 298 499 +297 299 300 +298 296 97 +303 301 298 +300 302 503 +301 303 304 +302 300 101 +307 305 302 +304 306 507 +305 307 308 +306 304 105 +311 309 306 +308 310 511 +309 311 312 +310 308 109 +315 313 310 +312 314 515 +313 315 316 +314 312 113 +319 317 314 +316 318 519 +317 319 320 +318 316 117 +323 321 318 +320 322 523 +321 323 324 +322 320 121 +327 325 322 +324 326 527 +325 327 328 +326 324 125 +331 329 326 +328 330 531 +329 331 332 +330 328 129 +335 333 330 +332 334 535 +333 335 336 +334 332 133 +339 337 334 +336 338 539 +337 339 340 +338 336 137 +343 341 338 +340 342 543 +341 343 344 +342 340 141 +347 345 342 +344 346 547 +345 347 348 +346 344 145 +351 349 346 +348 350 551 +349 351 352 +350 348 149 +355 353 350 +352 354 555 +353 355 356 +354 352 153 +359 357 354 +356 358 559 +357 359 360 +358 356 157 +363 361 358 +360 362 563 +361 363 364 +362 360 161 +367 365 362 +364 366 567 +365 367 368 +366 364 165 +371 369 366 +368 370 571 +369 371 372 +370 368 169 +375 373 370 +372 374 575 +373 375 376 +374 372 173 +379 377 374 +376 378 579 +377 379 380 +378 376 177 +383 381 378 +380 382 583 +381 383 384 +382 380 181 +387 385 382 +384 386 587 +385 387 388 +386 384 185 +391 389 386 +388 390 591 +389 391 392 +390 388 189 +395 393 390 +392 394 595 +393 395 396 +394 392 193 +399 397 394 +396 398 599 +397 399 -1 +398 396 197 +403 401 -1 +400 402 603 +401 403 404 +402 400 201 +407 405 402 +404 406 607 +405 407 408 +406 404 205 +411 409 406 +408 410 611 +409 411 412 +410 408 209 +415 413 410 +412 414 1590 +413 415 416 +414 412 213 +419 417 414 +416 418 1589 +417 419 420 +418 416 217 +423 421 418 +420 422 615 +421 423 424 +422 420 221 +427 425 422 +424 426 619 +425 427 428 +426 424 225 +431 429 426 +428 430 623 +429 431 432 +430 428 229 +435 433 430 +432 434 627 +433 435 436 +434 432 233 +439 437 434 +436 438 631 +437 439 440 +438 436 237 +443 441 438 +440 442 635 +441 443 444 +442 440 241 +447 445 442 +444 446 639 +445 447 448 +446 444 245 +451 449 446 +448 450 643 +449 451 452 +450 448 249 +455 453 450 +452 454 647 +453 455 456 +454 452 253 +459 457 454 +456 458 651 +457 459 460 +458 456 257 +463 461 458 +460 462 655 +461 463 464 +462 460 261 +467 465 462 +464 466 659 +465 467 468 +466 464 265 +471 469 466 +468 470 663 +469 471 472 +470 468 269 +475 473 470 +472 474 667 +473 475 476 +474 472 273 +479 477 474 +476 478 671 +477 479 480 +478 476 277 +483 481 478 +480 482 675 +481 483 484 +482 480 281 +487 485 482 +484 486 679 +485 487 488 +486 484 285 +491 489 486 +488 490 683 +489 491 492 +490 488 289 +495 493 490 +492 494 687 +493 495 496 +494 492 293 +499 497 494 +496 498 691 +497 499 500 +498 496 297 +503 501 498 +500 502 695 +501 503 504 +502 500 301 +507 505 502 +504 506 699 +505 507 508 +506 504 305 +511 509 506 +508 510 703 +509 511 512 +510 508 309 +515 513 510 +512 514 707 +513 515 516 +514 512 313 +519 517 514 +516 518 711 +517 519 520 +518 516 317 +523 521 518 +520 522 715 +521 523 524 +522 520 321 +527 525 522 +524 526 719 +525 527 528 +526 524 325 +531 529 526 +528 530 723 +529 531 532 +530 528 329 +535 533 530 +532 534 727 +533 535 536 +534 532 333 +539 537 534 +536 538 731 +537 539 540 +538 536 337 +543 541 538 +540 542 735 +541 543 544 +542 540 341 +547 545 542 +544 546 739 +545 547 548 +546 544 345 +551 549 546 +548 550 743 +549 551 552 +550 548 349 +555 553 550 +552 554 747 +553 555 556 +554 552 353 +559 557 554 +556 558 751 +557 559 560 +558 556 357 +563 561 558 +560 562 755 +561 563 564 +562 560 361 +567 565 562 +564 566 759 +565 567 568 +566 564 365 +571 569 566 +568 570 763 +569 571 572 +570 568 369 +575 573 570 +572 574 767 +573 575 576 +574 572 373 +579 577 574 +576 578 771 +577 579 580 +578 576 377 +583 581 578 +580 582 775 +581 583 584 +582 580 381 +587 585 582 +584 586 779 +585 587 588 +586 584 385 +591 589 586 +588 590 783 +589 591 592 +590 588 389 +595 593 590 +592 594 787 +593 595 596 +594 592 393 +599 597 594 +596 598 791 +597 599 -1 +598 596 397 +603 601 -1 +600 602 795 +601 603 604 +602 600 401 +607 605 602 +604 606 799 +605 607 608 +606 604 405 +611 609 606 +608 610 803 +609 611 1591 +610 608 409 +615 613 1588 +612 614 807 +613 615 616 +614 612 421 +619 617 614 +616 618 811 +617 619 620 +618 616 425 +623 621 618 +620 622 815 +621 623 624 +622 620 429 +627 625 622 +624 626 819 +625 627 628 +626 624 433 +631 629 626 +628 630 823 +629 631 632 +630 628 437 +635 633 630 +632 634 827 +633 635 636 +634 632 441 +639 637 634 +636 638 831 +637 639 640 +638 636 445 +643 641 638 +640 642 835 +641 643 644 +642 640 449 +647 645 642 +644 646 839 +645 647 648 +646 644 453 +651 649 646 +648 650 843 +649 651 652 +650 648 457 +655 653 650 +652 654 847 +653 655 656 +654 652 461 +659 657 654 +656 658 851 +657 659 660 +658 656 465 +663 661 658 +660 662 855 +661 663 664 +662 660 469 +667 665 662 +664 666 859 +665 667 668 +666 664 473 +671 669 666 +668 670 863 +669 671 672 +670 668 477 +675 673 670 +672 674 867 +673 675 676 +674 672 481 +679 677 674 +676 678 871 +677 679 680 +678 676 485 +683 681 678 +680 682 875 +681 683 684 +682 680 489 +687 685 682 +684 686 879 +685 687 688 +686 684 493 +691 689 686 +688 690 883 +689 691 692 +690 688 497 +695 693 690 +692 694 887 +693 695 696 +694 692 501 +699 697 694 +696 698 891 +697 699 700 +698 696 505 +703 701 698 +700 702 895 +701 703 704 +702 700 509 +707 705 702 +704 706 899 +705 707 708 +706 704 513 +711 709 706 +708 710 903 +709 711 712 +710 708 517 +715 713 710 +712 714 907 +713 715 716 +714 712 521 +719 717 714 +716 718 911 +717 719 720 +718 716 525 +723 721 718 +720 722 915 +721 723 724 +722 720 529 +727 725 722 +724 726 919 +725 727 728 +726 724 533 +731 729 726 +728 730 923 +729 731 732 +730 728 537 +735 733 730 +732 734 927 +733 735 736 +734 732 541 +739 737 734 +736 738 931 +737 739 740 +738 736 545 +743 741 738 +740 742 935 +741 743 744 +742 740 549 +747 745 742 +744 746 939 +745 747 748 +746 744 553 +751 749 746 +748 750 943 +749 751 752 +750 748 557 +755 753 750 +752 754 947 +753 755 756 +754 752 561 +759 757 754 +756 758 951 +757 759 760 +758 756 565 +763 761 758 +760 762 955 +761 763 764 +762 760 569 +767 765 762 +764 766 959 +765 767 768 +766 764 573 +771 769 766 +768 770 963 +769 771 772 +770 768 577 +775 773 770 +772 774 967 +773 775 776 +774 772 581 +779 777 774 +776 778 971 +777 779 780 +778 776 585 +783 781 778 +780 782 975 +781 783 784 +782 780 589 +787 785 782 +784 786 979 +785 787 788 +786 784 593 +791 789 786 +788 790 983 +789 791 -1 +790 788 597 +795 793 -1 +792 794 987 +793 795 796 +794 792 601 +799 797 794 +796 798 991 +797 799 800 +798 796 605 +803 801 798 +800 802 995 +801 803 1584 +802 800 609 +807 805 1587 +804 806 1007 +805 807 808 +806 804 613 +811 809 806 +808 810 1011 +809 811 812 +810 808 617 +815 813 810 +812 814 1015 +813 815 816 +814 812 621 +819 817 814 +816 818 1019 +817 819 820 +818 816 625 +823 821 818 +820 822 1023 +821 823 824 +822 820 629 +827 825 822 +824 826 1027 +825 827 828 +826 824 633 +831 829 826 +828 830 1031 +829 831 832 +830 828 637 +835 833 830 +832 834 1035 +833 835 836 +834 832 641 +839 837 834 +836 838 1039 +837 839 840 +838 836 645 +843 841 838 +840 842 1043 +841 843 844 +842 840 649 +847 845 842 +844 846 1047 +845 847 848 +846 844 653 +851 849 846 +848 850 1051 +849 851 852 +850 848 657 +855 853 850 +852 854 1055 +853 855 856 +854 852 661 +859 857 854 +856 858 1059 +857 859 860 +858 856 665 +863 861 858 +860 862 1063 +861 863 864 +862 860 669 +867 865 862 +864 866 1067 +865 867 868 +866 864 673 +871 869 866 +868 870 1071 +869 871 872 +870 868 677 +875 873 870 +872 874 1075 +873 875 876 +874 872 681 +879 877 874 +876 878 1079 +877 879 880 +878 876 685 +883 881 878 +880 882 1083 +881 883 884 +882 880 689 +887 885 882 +884 886 1087 +885 887 888 +886 884 693 +891 889 886 +888 890 1091 +889 891 892 +890 888 697 +895 893 890 +892 894 1095 +893 895 896 +894 892 701 +899 897 894 +896 898 1099 +897 899 900 +898 896 705 +903 901 898 +900 902 1103 +901 903 904 +902 900 709 +907 905 902 +904 906 1107 +905 907 908 +906 904 713 +911 909 906 +908 910 1111 +909 911 912 +910 908 717 +915 913 910 +912 914 1115 +913 915 916 +914 912 721 +919 917 914 +916 918 1119 +917 919 920 +918 916 725 +923 921 918 +920 922 1123 +921 923 924 +922 920 729 +927 925 922 +924 926 1127 +925 927 928 +926 924 733 +931 929 926 +928 930 1131 +929 931 932 +930 928 737 +935 933 930 +932 934 1135 +933 935 936 +934 932 741 +939 937 934 +936 938 1139 +937 939 940 +938 936 745 +943 941 938 +940 942 1143 +941 943 944 +942 940 749 +947 945 942 +944 946 1147 +945 947 948 +946 944 753 +951 949 946 +948 950 1151 +949 951 952 +950 948 757 +955 953 950 +952 954 1155 +953 955 956 +954 952 761 +959 957 954 +956 958 1159 +957 959 960 +958 956 765 +963 961 958 +960 962 1163 +961 963 964 +962 960 769 +967 965 962 +964 966 1167 +965 967 968 +966 964 773 +971 969 966 +968 970 1171 +969 971 972 +970 968 777 +975 973 970 +972 974 1175 +973 975 976 +974 972 781 +979 977 974 +976 978 1179 +977 979 980 +978 976 785 +983 981 978 +980 982 1183 +981 983 -1 +982 980 789 +987 985 -1 +984 986 1187 +985 987 988 +986 984 793 +991 989 986 +988 990 1191 +989 991 992 +990 988 797 +995 993 990 +992 994 1195 +993 995 996 +994 992 801 +999 997 994 +996 998 1199 +997 999 1000 +998 996 1585 +1003 1001 998 +1000 1002 1203 +1001 1003 1004 +1002 1000 1586 +1007 1005 1002 +1004 1006 1207 +1005 1007 1008 +1006 1004 805 +1011 1009 1006 +1008 1010 1211 +1009 1011 1012 +1010 1008 809 +1015 1013 1010 +1012 1014 1215 +1013 1015 1016 +1014 1012 813 +1019 1017 1014 +1016 1018 1219 +1017 1019 1020 +1018 1016 817 +1023 1021 1018 +1020 1022 1223 +1021 1023 1024 +1022 1020 821 +1027 1025 1022 +1024 1026 1227 +1025 1027 1028 +1026 1024 825 +1031 1029 1026 +1028 1030 1231 +1029 1031 1032 +1030 1028 829 +1035 1033 1030 +1032 1034 1235 +1033 1035 1036 +1034 1032 833 +1039 1037 1034 +1036 1038 1239 +1037 1039 1040 +1038 1036 837 +1043 1041 1038 +1040 1042 1243 +1041 1043 1044 +1042 1040 841 +1047 1045 1042 +1044 1046 1247 +1045 1047 1048 +1046 1044 845 +1051 1049 1046 +1048 1050 1251 +1049 1051 1052 +1050 1048 849 +1055 1053 1050 +1052 1054 1255 +1053 1055 1056 +1054 1052 853 +1059 1057 1054 +1056 1058 1259 +1057 1059 1060 +1058 1056 857 +1063 1061 1058 +1060 1062 1263 +1061 1063 1064 +1062 1060 861 +1067 1065 1062 +1064 1066 1267 +1065 1067 1068 +1066 1064 865 +1071 1069 1066 +1068 1070 1271 +1069 1071 1072 +1070 1068 869 +1075 1073 1070 +1072 1074 1275 +1073 1075 1076 +1074 1072 873 +1079 1077 1074 +1076 1078 1279 +1077 1079 1080 +1078 1076 877 +1083 1081 1078 +1080 1082 1283 +1081 1083 1084 +1082 1080 881 +1087 1085 1082 +1084 1086 1287 +1085 1087 1088 +1086 1084 885 +1091 1089 1086 +1088 1090 1291 +1089 1091 1092 +1090 1088 889 +1095 1093 1090 +1092 1094 1295 +1093 1095 1096 +1094 1092 893 +1099 1097 1094 +1096 1098 1299 +1097 1099 1100 +1098 1096 897 +1103 1101 1098 +1100 1102 1303 +1101 1103 1104 +1102 1100 901 +1107 1105 1102 +1104 1106 1307 +1105 1107 1108 +1106 1104 905 +1111 1109 1106 +1108 1110 1311 +1109 1111 1112 +1110 1108 909 +1115 1113 1110 +1112 1114 1315 +1113 1115 1116 +1114 1112 913 +1119 1117 1114 +1116 1118 1319 +1117 1119 1120 +1118 1116 917 +1123 1121 1118 +1120 1122 1323 +1121 1123 1124 +1122 1120 921 +1127 1125 1122 +1124 1126 1327 +1125 1127 1128 +1126 1124 925 +1131 1129 1126 +1128 1130 1331 +1129 1131 1132 +1130 1128 929 +1135 1133 1130 +1132 1134 1335 +1133 1135 1136 +1134 1132 933 +1139 1137 1134 +1136 1138 1339 +1137 1139 1140 +1138 1136 937 +1143 1141 1138 +1140 1142 1343 +1141 1143 1144 +1142 1140 941 +1147 1145 1142 +1144 1146 1347 +1145 1147 1148 +1146 1144 945 +1151 1149 1146 +1148 1150 1351 +1149 1151 1152 +1150 1148 949 +1155 1153 1150 +1152 1154 1355 +1153 1155 1156 +1154 1152 953 +1159 1157 1154 +1156 1158 1359 +1157 1159 1160 +1158 1156 957 +1163 1161 1158 +1160 1162 1363 +1161 1163 1164 +1162 1160 961 +1167 1165 1162 +1164 1166 1367 +1165 1167 1168 +1166 1164 965 +1171 1169 1166 +1168 1170 1371 +1169 1171 1172 +1170 1168 969 +1175 1173 1170 +1172 1174 1375 +1173 1175 1176 +1174 1172 973 +1179 1177 1174 +1176 1178 1379 +1177 1179 1180 +1178 1176 977 +1183 1181 1178 +1180 1182 1383 +1181 1183 -1 +1182 1180 981 +1187 1185 -1 +1184 1186 1387 +1185 1187 1188 +1186 1184 985 +1191 1189 1186 +1188 1190 1391 +1189 1191 1192 +1190 1188 989 +1195 1193 1190 +1192 1194 1395 +1193 1195 1196 +1194 1192 993 +1199 1197 1194 +1196 1198 1399 +1197 1199 1200 +1198 1196 997 +1203 1201 1198 +1200 1202 1403 +1201 1203 1204 +1202 1200 1001 +1207 1205 1202 +1204 1206 1407 +1205 1207 1208 +1206 1204 1005 +1211 1209 1206 +1208 1210 1411 +1209 1211 1212 +1210 1208 1009 +1215 1213 1210 +1212 1214 1415 +1213 1215 1216 +1214 1212 1013 +1219 1217 1214 +1216 1218 1419 +1217 1219 1220 +1218 1216 1017 +1223 1221 1218 +1220 1222 1423 +1221 1223 1224 +1222 1220 1021 +1227 1225 1222 +1224 1226 1427 +1225 1227 1228 +1226 1224 1025 +1231 1229 1226 +1228 1230 1431 +1229 1231 1232 +1230 1228 1029 +1235 1233 1230 +1232 1234 1435 +1233 1235 1236 +1234 1232 1033 +1239 1237 1234 +1236 1238 1439 +1237 1239 1240 +1238 1236 1037 +1243 1241 1238 +1240 1242 1443 +1241 1243 1244 +1242 1240 1041 +1247 1245 1242 +1244 1246 1447 +1245 1247 1248 +1246 1244 1045 +1251 1249 1246 +1248 1250 1451 +1249 1251 1252 +1250 1248 1049 +1255 1253 1250 +1252 1254 1455 +1253 1255 1256 +1254 1252 1053 +1259 1257 1254 +1256 1258 1459 +1257 1259 1260 +1258 1256 1057 +1263 1261 1258 +1260 1262 1463 +1261 1263 1264 +1262 1260 1061 +1267 1265 1262 +1264 1266 1467 +1265 1267 1268 +1266 1264 1065 +1271 1269 1266 +1268 1270 1471 +1269 1271 1272 +1270 1268 1069 +1275 1273 1270 +1272 1274 1475 +1273 1275 1276 +1274 1272 1073 +1279 1277 1274 +1276 1278 1479 +1277 1279 1280 +1278 1276 1077 +1283 1281 1278 +1280 1282 1483 +1281 1283 1284 +1282 1280 1081 +1287 1285 1282 +1284 1286 1487 +1285 1287 1288 +1286 1284 1085 +1291 1289 1286 +1288 1290 1491 +1289 1291 1292 +1290 1288 1089 +1295 1293 1290 +1292 1294 1495 +1293 1295 1296 +1294 1292 1093 +1299 1297 1294 +1296 1298 1499 +1297 1299 1300 +1298 1296 1097 +1303 1301 1298 +1300 1302 1503 +1301 1303 1304 +1302 1300 1101 +1307 1305 1302 +1304 1306 1507 +1305 1307 1308 +1306 1304 1105 +1311 1309 1306 +1308 1310 1511 +1309 1311 1312 +1310 1308 1109 +1315 1313 1310 +1312 1314 1515 +1313 1315 1316 +1314 1312 1113 +1319 1317 1314 +1316 1318 1519 +1317 1319 1320 +1318 1316 1117 +1323 1321 1318 +1320 1322 1523 +1321 1323 1324 +1322 1320 1121 +1327 1325 1322 +1324 1326 1527 +1325 1327 1328 +1326 1324 1125 +1331 1329 1326 +1328 1330 1531 +1329 1331 1332 +1330 1328 1129 +1335 1333 1330 +1332 1334 1535 +1333 1335 1336 +1334 1332 1133 +1339 1337 1334 +1336 1338 1539 +1337 1339 1340 +1338 1336 1137 +1343 1341 1338 +1340 1342 1543 +1341 1343 1344 +1342 1340 1141 +1347 1345 1342 +1344 1346 1547 +1345 1347 1348 +1346 1344 1145 +1351 1349 1346 +1348 1350 1551 +1349 1351 1352 +1350 1348 1149 +1355 1353 1350 +1352 1354 1555 +1353 1355 1356 +1354 1352 1153 +1359 1357 1354 +1356 1358 1559 +1357 1359 1360 +1358 1356 1157 +1363 1361 1358 +1360 1362 1563 +1361 1363 1364 +1362 1360 1161 +1367 1365 1362 +1364 1366 1567 +1365 1367 1368 +1366 1364 1165 +1371 1369 1366 +1368 1370 1571 +1369 1371 1372 +1370 1368 1169 +1375 1373 1370 +1372 1374 1575 +1373 1375 1376 +1374 1372 1173 +1379 1377 1374 +1376 1378 1579 +1377 1379 1380 +1378 1376 1177 +1383 1381 1378 +1380 1382 1583 +1381 1383 -1 +1382 1380 1181 +1387 1385 -1 +1384 1386 -1 +1385 1387 1388 +1386 1384 1185 +1391 1389 1386 +1388 1390 -1 +1389 1391 1392 +1390 1388 1189 +1395 1393 1390 +1392 1394 -1 +1393 1395 1396 +1394 1392 1193 +1399 1397 1394 +1396 1398 -1 +1397 1399 1400 +1398 1396 1197 +1403 1401 1398 +1400 1402 -1 +1401 1403 1404 +1402 1400 1201 +1407 1405 1402 +1404 1406 -1 +1405 1407 1408 +1406 1404 1205 +1411 1409 1406 +1408 1410 -1 +1409 1411 1412 +1410 1408 1209 +1415 1413 1410 +1412 1414 -1 +1413 1415 1416 +1414 1412 1213 +1419 1417 1414 +1416 1418 -1 +1417 1419 1420 +1418 1416 1217 +1423 1421 1418 +1420 1422 -1 +1421 1423 1424 +1422 1420 1221 +1427 1425 1422 +1424 1426 -1 +1425 1427 1428 +1426 1424 1225 +1431 1429 1426 +1428 1430 -1 +1429 1431 1432 +1430 1428 1229 +1435 1433 1430 +1432 1434 -1 +1433 1435 1436 +1434 1432 1233 +1439 1437 1434 +1436 1438 -1 +1437 1439 1440 +1438 1436 1237 +1443 1441 1438 +1440 1442 -1 +1441 1443 1444 +1442 1440 1241 +1447 1445 1442 +1444 1446 -1 +1445 1447 1448 +1446 1444 1245 +1451 1449 1446 +1448 1450 -1 +1449 1451 1452 +1450 1448 1249 +1455 1453 1450 +1452 1454 -1 +1453 1455 1456 +1454 1452 1253 +1459 1457 1454 +1456 1458 -1 +1457 1459 1460 +1458 1456 1257 +1463 1461 1458 +1460 1462 -1 +1461 1463 1464 +1462 1460 1261 +1467 1465 1462 +1464 1466 -1 +1465 1467 1468 +1466 1464 1265 +1471 1469 1466 +1468 1470 -1 +1469 1471 1472 +1470 1468 1269 +1475 1473 1470 +1472 1474 -1 +1473 1475 1476 +1474 1472 1273 +1479 1477 1474 +1476 1478 -1 +1477 1479 1480 +1478 1476 1277 +1483 1481 1478 +1480 1482 -1 +1481 1483 1484 +1482 1480 1281 +1487 1485 1482 +1484 1486 -1 +1485 1487 1488 +1486 1484 1285 +1491 1489 1486 +1488 1490 -1 +1489 1491 1492 +1490 1488 1289 +1495 1493 1490 +1492 1494 -1 +1493 1495 1496 +1494 1492 1293 +1499 1497 1494 +1496 1498 -1 +1497 1499 1500 +1498 1496 1297 +1503 1501 1498 +1500 1502 -1 +1501 1503 1504 +1502 1500 1301 +1507 1505 1502 +1504 1506 -1 +1505 1507 1508 +1506 1504 1305 +1511 1509 1506 +1508 1510 -1 +1509 1511 1512 +1510 1508 1309 +1515 1513 1510 +1512 1514 -1 +1513 1515 1516 +1514 1512 1313 +1519 1517 1514 +1516 1518 -1 +1517 1519 1520 +1518 1516 1317 +1523 1521 1518 +1520 1522 -1 +1521 1523 1524 +1522 1520 1321 +1527 1525 1522 +1524 1526 -1 +1525 1527 1528 +1526 1524 1325 +1531 1529 1526 +1528 1530 -1 +1529 1531 1532 +1530 1528 1329 +1535 1533 1530 +1532 1534 -1 +1533 1535 1536 +1534 1532 1333 +1539 1537 1534 +1536 1538 -1 +1537 1539 1540 +1538 1536 1337 +1543 1541 1538 +1540 1542 -1 +1541 1543 1544 +1542 1540 1341 +1547 1545 1542 +1544 1546 -1 +1545 1547 1548 +1546 1544 1345 +1551 1549 1546 +1548 1550 -1 +1549 1551 1552 +1550 1548 1349 +1555 1553 1550 +1552 1554 -1 +1553 1555 1556 +1554 1552 1353 +1559 1557 1554 +1556 1558 -1 +1557 1559 1560 +1558 1556 1357 +1563 1561 1558 +1560 1562 -1 +1561 1563 1564 +1562 1560 1361 +1567 1565 1562 +1564 1566 -1 +1565 1567 1568 +1566 1564 1365 +1571 1569 1566 +1568 1570 -1 +1569 1571 1572 +1570 1568 1369 +1575 1573 1570 +1572 1574 -1 +1573 1575 1576 +1574 1572 1373 +1579 1577 1574 +1576 1578 -1 +1577 1579 1580 +1578 1576 1377 +1583 1581 1578 +1580 1582 -1 +1581 1583 -1 +1582 1580 1381 +-1 1585 802 +1584 -1 999 +-1 1587 1003 +1586 -1 804 +-1 1589 612 +1588 -1 417 +-1 1591 413 +1590 -1 610 + diff --git a/demo/macro/kanal_square_fin.2d b/demo/macro/kanal_square_fin.2d new file mode 100644 index 0000000000000000000000000000000000000000..eb686d80ad770a9d27d408b51e29beadff83a996 --- /dev/null +++ b/demo/macro/kanal_square_fin.2d @@ -0,0 +1,5652 @@ +DIM: 2 +DIM_OF_WORLD: 2 + +number of elements: 1592 +number of vertices: 862 + +vertex coordinates: +0 0 +0.5 0 +1 0 +1.5 0 +2 0 +2.5 0 +3 0 +3.5 0 +4 0 +4.5 0 +5 0 +5.5 0 +6 0 +6.5 0 +7 0 +7.5 0 +8 0 +8.5 0 +9 0 +9.5 0 +10 0 +10.5 0 +11 0 +11.5 0 +12 0 +12.5 0 +13 0 +13.5 0 +14 0 +14.5 0 +15 0 +15.5 0 +16 0 +16.5 0 +17 0 +17.5 0 +18 0 +18.5 0 +19 0 +19.5 0 +20 0 +20.5 0 +21 0 +21.5 0 +22 0 +22.5 0 +23 0 +23.5 0 +24 0 +24.5 0 +25 0 +0 0.5 +0.5 0.5 +1 0.5 +1.5 0.5 +2 0.5 +2.5 0.5 +3 0.5 +3.5 0.5 +4 0.5 +4.5 0.5 +5 0.5 +5.5 0.5 +6 0.5 +6.5 0.5 +7 0.5 +7.5 0.5 +8 0.5 +8.5 0.5 +9 0.5 +9.5 0.5 +10 0.5 +10.5 0.5 +11 0.5 +11.5 0.5 +12 0.5 +12.5 0.5 +13 0.5 +13.5 0.5 +14 0.5 +14.5 0.5 +15 0.5 +15.5 0.5 +16 0.5 +16.5 0.5 +17 0.5 +17.5 0.5 +18 0.5 +18.5 0.5 +19 0.5 +19.5 0.5 +20 0.5 +20.5 0.5 +21 0.5 +21.5 0.5 +22 0.5 +22.5 0.5 +23 0.5 +23.5 0.5 +24 0.5 +24.5 0.5 +25 0.5 +0 1 +0.5 1 +1 1 +1.5 1 +2 1 +2.5 1 +3 1 +3.5 1 +4 1 +4.5 1 +5 1 +5.5 1 +6 1 +6.5 1 +7 1 +7.5 1 +8 1 +8.5 1 +9 1 +9.5 1 +10 1 +10.5 1 +11 1 +11.5 1 +12 1 +12.5 1 +13 1 +13.5 1 +14 1 +14.5 1 +15 1 +15.5 1 +16 1 +16.5 1 +17 1 +17.5 1 +18 1 +18.5 1 +19 1 +19.5 1 +20 1 +20.5 1 +21 1 +21.5 1 +22 1 +22.5 1 +23 1 +23.5 1 +24 1 +24.5 1 +25 1 +0 1.5 +0.5 1.5 +1 1.5 +1.5 1.5 +2 1.5 +2.5 1.5 +3 1.5 +3.5 1.5 +4 1.5 +4.5 1.5 +5 1.5 +5.5 1.5 +6 1.5 +6.5 1.5 +7 1.5 +7.5 1.5 +8 1.5 +8.5 1.5 +9 1.5 +9.5 1.5 +10 1.5 +10.5 1.5 +11 1.5 +11.5 1.5 +12 1.5 +12.5 1.5 +13 1.5 +13.5 1.5 +14 1.5 +14.5 1.5 +15 1.5 +15.5 1.5 +16 1.5 +16.5 1.5 +17 1.5 +17.5 1.5 +18 1.5 +18.5 1.5 +19 1.5 +19.5 1.5 +20 1.5 +20.5 1.5 +21 1.5 +21.5 1.5 +22 1.5 +22.5 1.5 +23 1.5 +23.5 1.5 +24 1.5 +24.5 1.5 +25 1.5 +0 2 +0.5 2 +1 2 +1.5 2 +2.5 1.95 +3 1.95 +3.5 1.95 +4 1.95 +4.5 2 +5 2 +5.5 2 +6 2 +6.5 2 +7 2 +7.5 2 +8 2 +8.5 2 +9 2 +9.5 2 +10 2 +10.5 2 +11 2 +11.5 2 +12 2 +12.5 2 +13 2 +13.5 2 +14 2 +14.5 2 +15 2 +15.5 2 +16 2 +16.5 2 +17 2 +17.5 2 +18 2 +18.5 2 +19 2 +19.5 2 +20 2 +20.5 2 +21 2 +21.5 2 +22 2 +22.5 2 +23 2 +23.5 2 +24 2 +24.5 2 +25 2 +0 2.5 +0.5 2.5 +1 2.5 +1.5 2.5 +2.5 2.5 +3 2.5 +3.5 2.5 +4 2.5 +4.5 2.5 +5 2.5 +5.5 2.5 +6 2.5 +6.5 2.5 +7 2.5 +7.5 2.5 +8 2.5 +8.5 2.5 +9 2.5 +9.5 2.5 +10 2.5 +10.5 2.5 +11 2.5 +11.5 2.5 +12 2.5 +12.5 2.5 +13 2.5 +13.5 2.5 +14 2.5 +14.5 2.5 +15 2.5 +15.5 2.5 +16 2.5 +16.5 2.5 +17 2.5 +17.5 2.5 +18 2.5 +18.5 2.5 +19 2.5 +19.5 2.5 +20 2.5 +20.5 2.5 +21 2.5 +21.5 2.5 +22 2.5 +22.5 2.5 +23 2.5 +23.5 2.5 +24 2.5 +24.5 2.5 +25 2.5 +2 2.5 +0 3.033333333 +0.5 3.033333333 +1 3.033333333 +1.5 3.033333333 +2 3.033333333 +2.5 3.033333333 +3 3.033333333 +3.5 3.033333333 +4 3.033333333 +4.5 3.033333333 +5 3.033333333 +5.5 3.033333333 +6 3.033333333 +6.5 3.033333333 +7 3.033333333 +7.5 3.033333333 +8 3.033333333 +8.5 3.033333333 +9 3.033333333 +9.5 3.033333333 +10 3.033333333 +10.5 3.033333333 +11 3.033333333 +11.5 3.033333333 +12 3.033333333 +12.5 3.033333333 +13 3.033333333 +13.5 3.033333333 +14 3.033333333 +14.5 3.033333333 +15 3.033333333 +15.5 3.033333333 +16 3.033333333 +16.5 3.033333333 +17 3.033333333 +17.5 3.033333333 +18 3.033333333 +18.5 3.033333333 +19 3.033333333 +19.5 3.033333333 +20 3.033333333 +20.5 3.033333333 +21 3.033333333 +21.5 3.033333333 +22 3.033333333 +22.5 3.033333333 +23 3.033333333 +23.5 3.033333333 +24 3.033333333 +24.5 3.033333333 +25 3.033333333 +0 3.566666667 +0.5 3.566666667 +1 3.566666667 +1.5 3.566666667 +2 3.566666667 +2.5 3.566666667 +3 3.566666667 +3.5 3.566666667 +4 3.566666667 +4.5 3.566666667 +5 3.566666667 +5.5 3.566666667 +6 3.566666667 +6.5 3.566666667 +7 3.566666667 +7.5 3.566666667 +8 3.566666667 +8.5 3.566666667 +9 3.566666667 +9.5 3.566666667 +10 3.566666667 +10.5 3.566666667 +11 3.566666667 +11.5 3.566666667 +12 3.566666667 +12.5 3.566666667 +13 3.566666667 +13.5 3.566666667 +14 3.566666667 +14.5 3.566666667 +15 3.566666667 +15.5 3.566666667 +16 3.566666667 +16.5 3.566666667 +17 3.566666667 +17.5 3.566666667 +18 3.566666667 +18.5 3.566666667 +19 3.566666667 +19.5 3.566666667 +20 3.566666667 +20.5 3.566666667 +21 3.566666667 +21.5 3.566666667 +22 3.566666667 +22.5 3.566666667 +23 3.566666667 +23.5 3.566666667 +24 3.566666667 +24.5 3.566666667 +25 3.566666667 +0 4.1 +0.5 4.1 +1 4.1 +1.5 4.1 +2 4.1 +2.5 4.1 +3 4.1 +3.5 4.1 +4 4.1 +4.5 4.1 +5 4.1 +5.5 4.1 +6 4.1 +6.5 4.1 +7 4.1 +7.5 4.1 +8 4.1 +8.5 4.1 +9 4.1 +9.5 4.1 +10 4.1 +10.5 4.1 +11 4.1 +11.5 4.1 +12 4.1 +12.5 4.1 +13 4.1 +13.5 4.1 +14 4.1 +14.5 4.1 +15 4.1 +15.5 4.1 +16 4.1 +16.5 4.1 +17 4.1 +17.5 4.1 +18 4.1 +18.5 4.1 +19 4.1 +19.5 4.1 +20 4.1 +20.5 4.1 +21 4.1 +21.5 4.1 +22 4.1 +22.5 4.1 +23 4.1 +23.5 4.1 +24 4.1 +24.5 4.1 +25 4.1 +0.25 0.25 +0.75 0.25 +1.25 0.25 +1.75 0.25 +2.25 0.25 +2.75 0.25 +3.25 0.25 +3.75 0.25 +4.25 0.25 +4.75 0.25 +5.25 0.25 +5.75 0.25 +6.25 0.25 +6.75 0.25 +7.25 0.25 +7.75 0.25 +8.25 0.25 +8.75 0.25 +9.25 0.25 +9.75 0.25 +10.25 0.25 +10.75 0.25 +11.25 0.25 +11.75 0.25 +12.25 0.25 +12.75 0.25 +13.25 0.25 +13.75 0.25 +14.25 0.25 +14.75 0.25 +15.25 0.25 +15.75 0.25 +16.25 0.25 +16.75 0.25 +17.25 0.25 +17.75 0.25 +18.25 0.25 +18.75 0.25 +19.25 0.25 +19.75 0.25 +20.25 0.25 +20.75 0.25 +21.25 0.25 +21.75 0.25 +22.25 0.25 +22.75 0.25 +23.25 0.25 +23.75 0.25 +24.25 0.25 +24.75 0.25 +0.25 0.75 +0.75 0.75 +1.25 0.75 +1.75 0.75 +2.25 0.75 +2.75 0.75 +3.25 0.75 +3.75 0.75 +4.25 0.75 +4.75 0.75 +5.25 0.75 +5.75 0.75 +6.25 0.75 +6.75 0.75 +7.25 0.75 +7.75 0.75 +8.25 0.75 +8.75 0.75 +9.25 0.75 +9.75 0.75 +10.25 0.75 +10.75 0.75 +11.25 0.75 +11.75 0.75 +12.25 0.75 +12.75 0.75 +13.25 0.75 +13.75 0.75 +14.25 0.75 +14.75 0.75 +15.25 0.75 +15.75 0.75 +16.25 0.75 +16.75 0.75 +17.25 0.75 +17.75 0.75 +18.25 0.75 +18.75 0.75 +19.25 0.75 +19.75 0.75 +20.25 0.75 +20.75 0.75 +21.25 0.75 +21.75 0.75 +22.25 0.75 +22.75 0.75 +23.25 0.75 +23.75 0.75 +24.25 0.75 +24.75 0.75 +0.25 1.25 +0.75 1.25 +1.25 1.25 +1.75 1.25 +2.25 1.25 +2.75 1.25 +3.25 1.25 +3.75 1.25 +4.25 1.25 +4.75 1.25 +5.25 1.25 +5.75 1.25 +6.25 1.25 +6.75 1.25 +7.25 1.25 +7.75 1.25 +8.25 1.25 +8.75 1.25 +9.25 1.25 +9.75 1.25 +10.25 1.25 +10.75 1.25 +11.25 1.25 +11.75 1.25 +12.25 1.25 +12.75 1.25 +13.25 1.25 +13.75 1.25 +14.25 1.25 +14.75 1.25 +15.25 1.25 +15.75 1.25 +16.25 1.25 +16.75 1.25 +17.25 1.25 +17.75 1.25 +18.25 1.25 +18.75 1.25 +19.25 1.25 +19.75 1.25 +20.25 1.25 +20.75 1.25 +21.25 1.25 +21.75 1.25 +22.25 1.25 +22.75 1.25 +23.25 1.25 +23.75 1.25 +24.25 1.25 +24.75 1.25 +0.25 1.75 +0.75 1.75 +1.25 1.75 +2.75 1.75 +3.25 1.75 +3.75 1.75 +4.25 1.75 +4.75 1.75 +5.25 1.75 +5.75 1.75 +6.25 1.75 +6.75 1.75 +7.25 1.75 +7.75 1.75 +8.25 1.75 +8.75 1.75 +9.25 1.75 +9.75 1.75 +10.25 1.75 +10.75 1.75 +11.25 1.75 +11.75 1.75 +12.25 1.75 +12.75 1.75 +13.25 1.75 +13.75 1.75 +14.25 1.75 +14.75 1.75 +15.25 1.75 +15.75 1.75 +16.25 1.75 +16.75 1.75 +17.25 1.75 +17.75 1.75 +18.25 1.75 +18.75 1.75 +19.25 1.75 +19.75 1.75 +20.25 1.75 +20.75 1.75 +21.25 1.75 +21.75 1.75 +22.25 1.75 +22.75 1.75 +23.25 1.75 +23.75 1.75 +24.25 1.75 +24.75 1.75 +0.25 2.25 +0.75 2.25 +1.25 2.25 +2.75 2.25 +3.25 2.25 +3.75 2.25 +4.25 2.25 +4.75 2.25 +5.25 2.25 +5.75 2.25 +6.25 2.25 +6.75 2.25 +7.25 2.25 +7.75 2.25 +8.25 2.25 +8.75 2.25 +9.25 2.25 +9.75 2.25 +10.25 2.25 +10.75 2.25 +11.25 2.25 +11.75 2.25 +12.25 2.25 +12.75 2.25 +13.25 2.25 +13.75 2.25 +14.25 2.25 +14.75 2.25 +15.25 2.25 +15.75 2.25 +16.25 2.25 +16.75 2.25 +17.25 2.25 +17.75 2.25 +18.25 2.25 +18.75 2.25 +19.25 2.25 +19.75 2.25 +20.25 2.25 +20.75 2.25 +21.25 2.25 +21.75 2.25 +22.25 2.25 +22.75 2.25 +23.25 2.25 +23.75 2.25 +24.25 2.25 +24.75 2.25 +0.25 2.766666667 +0.75 2.766666667 +1.25 2.766666667 +1.75 2.766666667 +2.25 2.766666667 +2.75 2.766666667 +3.25 2.766666667 +3.75 2.766666667 +4.25 2.766666667 +4.75 2.766666667 +5.25 2.766666667 +5.75 2.766666667 +6.25 2.766666667 +6.75 2.766666667 +7.25 2.766666667 +7.75 2.766666667 +8.25 2.766666667 +8.75 2.766666667 +9.25 2.766666667 +9.75 2.766666667 +10.25 2.766666667 +10.75 2.766666667 +11.25 2.766666667 +11.75 2.766666667 +12.25 2.766666667 +12.75 2.766666667 +13.25 2.766666667 +13.75 2.766666667 +14.25 2.766666667 +14.75 2.766666667 +15.25 2.766666667 +15.75 2.766666667 +16.25 2.766666667 +16.75 2.766666667 +17.25 2.766666667 +17.75 2.766666667 +18.25 2.766666667 +18.75 2.766666667 +19.25 2.766666667 +19.75 2.766666667 +20.25 2.766666667 +20.75 2.766666667 +21.25 2.766666667 +21.75 2.766666667 +22.25 2.766666667 +22.75 2.766666667 +23.25 2.766666667 +23.75 2.766666667 +24.25 2.766666667 +24.75 2.766666667 +0.25 3.3 +0.75 3.3 +1.25 3.3 +1.75 3.3 +2.25 3.3 +2.75 3.3 +3.25 3.3 +3.75 3.3 +4.25 3.3 +4.75 3.3 +5.25 3.3 +5.75 3.3 +6.25 3.3 +6.75 3.3 +7.25 3.3 +7.75 3.3 +8.25 3.3 +8.75 3.3 +9.25 3.3 +9.75 3.3 +10.25 3.3 +10.75 3.3 +11.25 3.3 +11.75 3.3 +12.25 3.3 +12.75 3.3 +13.25 3.3 +13.75 3.3 +14.25 3.3 +14.75 3.3 +15.25 3.3 +15.75 3.3 +16.25 3.3 +16.75 3.3 +17.25 3.3 +17.75 3.3 +18.25 3.3 +18.75 3.3 +19.25 3.3 +19.75 3.3 +20.25 3.3 +20.75 3.3 +21.25 3.3 +21.75 3.3 +22.25 3.3 +22.75 3.3 +23.25 3.3 +23.75 3.3 +24.25 3.3 +24.75 3.3 +0.25 3.833333333 +0.75 3.833333333 +1.25 3.833333333 +1.75 3.833333333 +2.25 3.833333333 +2.75 3.833333333 +3.25 3.833333333 +3.75 3.833333333 +4.25 3.833333333 +4.75 3.833333333 +5.25 3.833333333 +5.75 3.833333333 +6.25 3.833333333 +6.75 3.833333333 +7.25 3.833333333 +7.75 3.833333333 +8.25 3.833333333 +8.75 3.833333333 +9.25 3.833333333 +9.75 3.833333333 +10.25 3.833333333 +10.75 3.833333333 +11.25 3.833333333 +11.75 3.833333333 +12.25 3.833333333 +12.75 3.833333333 +13.25 3.833333333 +13.75 3.833333333 +14.25 3.833333333 +14.75 3.833333333 +15.25 3.833333333 +15.75 3.833333333 +16.25 3.833333333 +16.75 3.833333333 +17.25 3.833333333 +17.75 3.833333333 +18.25 3.833333333 +18.75 3.833333333 +19.25 3.833333333 +19.75 3.833333333 +20.25 3.833333333 +20.75 3.833333333 +21.25 3.833333333 +21.75 3.833333333 +22.25 3.833333333 +22.75 3.833333333 +23.25 3.833333333 +23.75 3.833333333 +24.25 3.833333333 +24.75 3.833333333 +1.75 2.25 +2.25 2.25 +1.75 1.75 +2.25 1.75 +2.5 2.05 +3 2.05 +3.5 2.05 +4 2.05 + +element vertices: +51 0 458 +52 51 458 +1 52 458 +0 1 458 +52 1 459 +53 52 459 +2 53 459 +1 2 459 +53 2 460 +54 53 460 +3 54 460 +2 3 460 +54 3 461 +55 54 461 +4 55 461 +3 4 461 +55 4 462 +56 55 462 +5 56 462 +4 5 462 +56 5 463 +57 56 463 +6 57 463 +5 6 463 +57 6 464 +58 57 464 +7 58 464 +6 7 464 +58 7 465 +59 58 465 +8 59 465 +7 8 465 +59 8 466 +60 59 466 +9 60 466 +8 9 466 +60 9 467 +61 60 467 +10 61 467 +9 10 467 +61 10 468 +62 61 468 +11 62 468 +10 11 468 +62 11 469 +63 62 469 +12 63 469 +11 12 469 +63 12 470 +64 63 470 +13 64 470 +12 13 470 +64 13 471 +65 64 471 +14 65 471 +13 14 471 +65 14 472 +66 65 472 +15 66 472 +14 15 472 +66 15 473 +67 66 473 +16 67 473 +15 16 473 +67 16 474 +68 67 474 +17 68 474 +16 17 474 +68 17 475 +69 68 475 +18 69 475 +17 18 475 +69 18 476 +70 69 476 +19 70 476 +18 19 476 +70 19 477 +71 70 477 +20 71 477 +19 20 477 +71 20 478 +72 71 478 +21 72 478 +20 21 478 +72 21 479 +73 72 479 +22 73 479 +21 22 479 +73 22 480 +74 73 480 +23 74 480 +22 23 480 +74 23 481 +75 74 481 +24 75 481 +23 24 481 +75 24 482 +76 75 482 +25 76 482 +24 25 482 +76 25 483 +77 76 483 +26 77 483 +25 26 483 +77 26 484 +78 77 484 +27 78 484 +26 27 484 +78 27 485 +79 78 485 +28 79 485 +27 28 485 +79 28 486 +80 79 486 +29 80 486 +28 29 486 +80 29 487 +81 80 487 +30 81 487 +29 30 487 +81 30 488 +82 81 488 +31 82 488 +30 31 488 +82 31 489 +83 82 489 +32 83 489 +31 32 489 +83 32 490 +84 83 490 +33 84 490 +32 33 490 +84 33 491 +85 84 491 +34 85 491 +33 34 491 +85 34 492 +86 85 492 +35 86 492 +34 35 492 +86 35 493 +87 86 493 +36 87 493 +35 36 493 +87 36 494 +88 87 494 +37 88 494 +36 37 494 +88 37 495 +89 88 495 +38 89 495 +37 38 495 +89 38 496 +90 89 496 +39 90 496 +38 39 496 +90 39 497 +91 90 497 +40 91 497 +39 40 497 +91 40 498 +92 91 498 +41 92 498 +40 41 498 +92 41 499 +93 92 499 +42 93 499 +41 42 499 +93 42 500 +94 93 500 +43 94 500 +42 43 500 +94 43 501 +95 94 501 +44 95 501 +43 44 501 +95 44 502 +96 95 502 +45 96 502 +44 45 502 +96 45 503 +97 96 503 +46 97 503 +45 46 503 +97 46 504 +98 97 504 +47 98 504 +46 47 504 +98 47 505 +99 98 505 +48 99 505 +47 48 505 +99 48 506 +100 99 506 +49 100 506 +48 49 506 +100 49 507 +101 100 507 +50 101 507 +49 50 507 +102 51 508 +103 102 508 +52 103 508 +51 52 508 +103 52 509 +104 103 509 +53 104 509 +52 53 509 +104 53 510 +105 104 510 +54 105 510 +53 54 510 +105 54 511 +106 105 511 +55 106 511 +54 55 511 +106 55 512 +107 106 512 +56 107 512 +55 56 512 +107 56 513 +108 107 513 +57 108 513 +56 57 513 +108 57 514 +109 108 514 +58 109 514 +57 58 514 +109 58 515 +110 109 515 +59 110 515 +58 59 515 +110 59 516 +111 110 516 +60 111 516 +59 60 516 +111 60 517 +112 111 517 +61 112 517 +60 61 517 +112 61 518 +113 112 518 +62 113 518 +61 62 518 +113 62 519 +114 113 519 +63 114 519 +62 63 519 +114 63 520 +115 114 520 +64 115 520 +63 64 520 +115 64 521 +116 115 521 +65 116 521 +64 65 521 +116 65 522 +117 116 522 +66 117 522 +65 66 522 +117 66 523 +118 117 523 +67 118 523 +66 67 523 +118 67 524 +119 118 524 +68 119 524 +67 68 524 +119 68 525 +120 119 525 +69 120 525 +68 69 525 +120 69 526 +121 120 526 +70 121 526 +69 70 526 +121 70 527 +122 121 527 +71 122 527 +70 71 527 +122 71 528 +123 122 528 +72 123 528 +71 72 528 +123 72 529 +124 123 529 +73 124 529 +72 73 529 +124 73 530 +125 124 530 +74 125 530 +73 74 530 +125 74 531 +126 125 531 +75 126 531 +74 75 531 +126 75 532 +127 126 532 +76 127 532 +75 76 532 +127 76 533 +128 127 533 +77 128 533 +76 77 533 +128 77 534 +129 128 534 +78 129 534 +77 78 534 +129 78 535 +130 129 535 +79 130 535 +78 79 535 +130 79 536 +131 130 536 +80 131 536 +79 80 536 +131 80 537 +132 131 537 +81 132 537 +80 81 537 +132 81 538 +133 132 538 +82 133 538 +81 82 538 +133 82 539 +134 133 539 +83 134 539 +82 83 539 +134 83 540 +135 134 540 +84 135 540 +83 84 540 +135 84 541 +136 135 541 +85 136 541 +84 85 541 +136 85 542 +137 136 542 +86 137 542 +85 86 542 +137 86 543 +138 137 543 +87 138 543 +86 87 543 +138 87 544 +139 138 544 +88 139 544 +87 88 544 +139 88 545 +140 139 545 +89 140 545 +88 89 545 +140 89 546 +141 140 546 +90 141 546 +89 90 546 +141 90 547 +142 141 547 +91 142 547 +90 91 547 +142 91 548 +143 142 548 +92 143 548 +91 92 548 +143 92 549 +144 143 549 +93 144 549 +92 93 549 +144 93 550 +145 144 550 +94 145 550 +93 94 550 +145 94 551 +146 145 551 +95 146 551 +94 95 551 +146 95 552 +147 146 552 +96 147 552 +95 96 552 +147 96 553 +148 147 553 +97 148 553 +96 97 553 +148 97 554 +149 148 554 +98 149 554 +97 98 554 +149 98 555 +150 149 555 +99 150 555 +98 99 555 +150 99 556 +151 150 556 +100 151 556 +99 100 556 +151 100 557 +152 151 557 +101 152 557 +100 101 557 +153 102 558 +154 153 558 +103 154 558 +102 103 558 +154 103 559 +155 154 559 +104 155 559 +103 104 559 +155 104 560 +156 155 560 +105 156 560 +104 105 560 +156 105 561 +157 156 561 +106 157 561 +105 106 561 +157 106 562 +158 157 562 +107 158 562 +106 107 562 +158 107 563 +159 158 563 +108 159 563 +107 108 563 +159 108 564 +160 159 564 +109 160 564 +108 109 564 +160 109 565 +161 160 565 +110 161 565 +109 110 565 +161 110 566 +162 161 566 +111 162 566 +110 111 566 +162 111 567 +163 162 567 +112 163 567 +111 112 567 +163 112 568 +164 163 568 +113 164 568 +112 113 568 +164 113 569 +165 164 569 +114 165 569 +113 114 569 +165 114 570 +166 165 570 +115 166 570 +114 115 570 +166 115 571 +167 166 571 +116 167 571 +115 116 571 +167 116 572 +168 167 572 +117 168 572 +116 117 572 +168 117 573 +169 168 573 +118 169 573 +117 118 573 +169 118 574 +170 169 574 +119 170 574 +118 119 574 +170 119 575 +171 170 575 +120 171 575 +119 120 575 +171 120 576 +172 171 576 +121 172 576 +120 121 576 +172 121 577 +173 172 577 +122 173 577 +121 122 577 +173 122 578 +174 173 578 +123 174 578 +122 123 578 +174 123 579 +175 174 579 +124 175 579 +123 124 579 +175 124 580 +176 175 580 +125 176 580 +124 125 580 +176 125 581 +177 176 581 +126 177 581 +125 126 581 +177 126 582 +178 177 582 +127 178 582 +126 127 582 +178 127 583 +179 178 583 +128 179 583 +127 128 583 +179 128 584 +180 179 584 +129 180 584 +128 129 584 +180 129 585 +181 180 585 +130 181 585 +129 130 585 +181 130 586 +182 181 586 +131 182 586 +130 131 586 +182 131 587 +183 182 587 +132 183 587 +131 132 587 +183 132 588 +184 183 588 +133 184 588 +132 133 588 +184 133 589 +185 184 589 +134 185 589 +133 134 589 +185 134 590 +186 185 590 +135 186 590 +134 135 590 +186 135 591 +187 186 591 +136 187 591 +135 136 591 +187 136 592 +188 187 592 +137 188 592 +136 137 592 +188 137 593 +189 188 593 +138 189 593 +137 138 593 +189 138 594 +190 189 594 +139 190 594 +138 139 594 +190 139 595 +191 190 595 +140 191 595 +139 140 595 +191 140 596 +192 191 596 +141 192 596 +140 141 596 +192 141 597 +193 192 597 +142 193 597 +141 142 597 +193 142 598 +194 193 598 +143 194 598 +142 143 598 +194 143 599 +195 194 599 +144 195 599 +143 144 599 +195 144 600 +196 195 600 +145 196 600 +144 145 600 +196 145 601 +197 196 601 +146 197 601 +145 146 601 +197 146 602 +198 197 602 +147 198 602 +146 147 602 +198 147 603 +199 198 603 +148 199 603 +147 148 603 +199 148 604 +200 199 604 +149 200 604 +148 149 604 +200 149 605 +201 200 605 +150 201 605 +149 150 605 +201 150 606 +202 201 606 +151 202 606 +150 151 606 +202 151 607 +203 202 607 +152 203 607 +151 152 607 +204 153 608 +205 204 608 +154 205 608 +153 154 608 +205 154 609 +206 205 609 +155 206 609 +154 155 609 +206 155 610 +207 206 610 +156 207 610 +155 156 610 +208 158 611 +209 208 611 +159 209 611 +158 159 611 +209 159 612 +210 209 612 +160 210 612 +159 160 612 +210 160 613 +211 210 613 +161 211 613 +160 161 613 +211 161 614 +212 211 614 +162 212 614 +161 162 614 +212 162 615 +213 212 615 +163 213 615 +162 163 615 +213 163 616 +214 213 616 +164 214 616 +163 164 616 +214 164 617 +215 214 617 +165 215 617 +164 165 617 +215 165 618 +216 215 618 +166 216 618 +165 166 618 +216 166 619 +217 216 619 +167 217 619 +166 167 619 +217 167 620 +218 217 620 +168 218 620 +167 168 620 +218 168 621 +219 218 621 +169 219 621 +168 169 621 +219 169 622 +220 219 622 +170 220 622 +169 170 622 +220 170 623 +221 220 623 +171 221 623 +170 171 623 +221 171 624 +222 221 624 +172 222 624 +171 172 624 +222 172 625 +223 222 625 +173 223 625 +172 173 625 +223 173 626 +224 223 626 +174 224 626 +173 174 626 +224 174 627 +225 224 627 +175 225 627 +174 175 627 +225 175 628 +226 225 628 +176 226 628 +175 176 628 +226 176 629 +227 226 629 +177 227 629 +176 177 629 +227 177 630 +228 227 630 +178 228 630 +177 178 630 +228 178 631 +229 228 631 +179 229 631 +178 179 631 +229 179 632 +230 229 632 +180 230 632 +179 180 632 +230 180 633 +231 230 633 +181 231 633 +180 181 633 +231 181 634 +232 231 634 +182 232 634 +181 182 634 +232 182 635 +233 232 635 +183 233 635 +182 183 635 +233 183 636 +234 233 636 +184 234 636 +183 184 636 +234 184 637 +235 234 637 +185 235 637 +184 185 637 +235 185 638 +236 235 638 +186 236 638 +185 186 638 +236 186 639 +237 236 639 +187 237 639 +186 187 639 +237 187 640 +238 237 640 +188 238 640 +187 188 640 +238 188 641 +239 238 641 +189 239 641 +188 189 641 +239 189 642 +240 239 642 +190 240 642 +189 190 642 +240 190 643 +241 240 643 +191 241 643 +190 191 643 +241 191 644 +242 241 644 +192 242 644 +191 192 644 +242 192 645 +243 242 645 +193 243 645 +192 193 645 +243 193 646 +244 243 646 +194 244 646 +193 194 646 +244 194 647 +245 244 647 +195 245 647 +194 195 647 +245 195 648 +246 245 648 +196 246 648 +195 196 648 +246 196 649 +247 246 649 +197 247 649 +196 197 649 +247 197 650 +248 247 650 +198 248 650 +197 198 650 +248 198 651 +249 248 651 +199 249 651 +198 199 651 +249 199 652 +250 249 652 +200 250 652 +199 200 652 +250 200 653 +251 250 653 +201 251 653 +200 201 653 +251 201 654 +252 251 654 +202 252 654 +201 202 654 +252 202 655 +253 252 655 +203 253 655 +202 203 655 +254 204 656 +255 254 656 +205 255 656 +204 205 656 +255 205 657 +256 255 657 +206 256 657 +205 206 657 +256 206 658 +257 256 658 +207 257 658 +206 207 658 +258 858 659 +259 258 659 +859 259 659 +858 859 659 +259 859 660 +260 259 660 +860 260 660 +859 860 660 +260 860 661 +261 260 661 +861 261 661 +860 861 661 +261 861 662 +262 261 662 +212 262 662 +861 212 662 +262 212 663 +263 262 663 +213 263 663 +212 213 663 +263 213 664 +264 263 664 +214 264 664 +213 214 664 +264 214 665 +265 264 665 +215 265 665 +214 215 665 +265 215 666 +266 265 666 +216 266 666 +215 216 666 +266 216 667 +267 266 667 +217 267 667 +216 217 667 +267 217 668 +268 267 668 +218 268 668 +217 218 668 +268 218 669 +269 268 669 +219 269 669 +218 219 669 +269 219 670 +270 269 670 +220 270 670 +219 220 670 +270 220 671 +271 270 671 +221 271 671 +220 221 671 +271 221 672 +272 271 672 +222 272 672 +221 222 672 +272 222 673 +273 272 673 +223 273 673 +222 223 673 +273 223 674 +274 273 674 +224 274 674 +223 224 674 +274 224 675 +275 274 675 +225 275 675 +224 225 675 +275 225 676 +276 275 676 +226 276 676 +225 226 676 +276 226 677 +277 276 677 +227 277 677 +226 227 677 +277 227 678 +278 277 678 +228 278 678 +227 228 678 +278 228 679 +279 278 679 +229 279 679 +228 229 679 +279 229 680 +280 279 680 +230 280 680 +229 230 680 +280 230 681 +281 280 681 +231 281 681 +230 231 681 +281 231 682 +282 281 682 +232 282 682 +231 232 682 +282 232 683 +283 282 683 +233 283 683 +232 233 683 +283 233 684 +284 283 684 +234 284 684 +233 234 684 +284 234 685 +285 284 685 +235 285 685 +234 235 685 +285 235 686 +286 285 686 +236 286 686 +235 236 686 +286 236 687 +287 286 687 +237 287 687 +236 237 687 +287 237 688 +288 287 688 +238 288 688 +237 238 688 +288 238 689 +289 288 689 +239 289 689 +238 239 689 +289 239 690 +290 289 690 +240 290 690 +239 240 690 +290 240 691 +291 290 691 +241 291 691 +240 241 691 +291 241 692 +292 291 692 +242 292 692 +241 242 692 +292 242 693 +293 292 693 +243 293 693 +242 243 693 +293 243 694 +294 293 694 +244 294 694 +243 244 694 +294 244 695 +295 294 695 +245 295 695 +244 245 695 +295 245 696 +296 295 696 +246 296 696 +245 246 696 +296 246 697 +297 296 697 +247 297 697 +246 247 697 +297 247 698 +298 297 698 +248 298 698 +247 248 698 +298 248 699 +299 298 699 +249 299 699 +248 249 699 +299 249 700 +300 299 700 +250 300 700 +249 250 700 +300 250 701 +301 300 701 +251 301 701 +250 251 701 +301 251 702 +302 301 702 +252 302 702 +251 252 702 +302 252 703 +303 302 703 +253 303 703 +252 253 703 +305 254 704 +306 305 704 +255 306 704 +254 255 704 +306 255 705 +307 306 705 +256 307 705 +255 256 705 +307 256 706 +308 307 706 +257 308 706 +256 257 706 +308 257 707 +309 308 707 +304 309 707 +257 304 707 +309 304 708 +310 309 708 +258 310 708 +304 258 708 +310 258 709 +311 310 709 +259 311 709 +258 259 709 +311 259 710 +312 311 710 +260 312 710 +259 260 710 +312 260 711 +313 312 711 +261 313 711 +260 261 711 +313 261 712 +314 313 712 +262 314 712 +261 262 712 +314 262 713 +315 314 713 +263 315 713 +262 263 713 +315 263 714 +316 315 714 +264 316 714 +263 264 714 +316 264 715 +317 316 715 +265 317 715 +264 265 715 +317 265 716 +318 317 716 +266 318 716 +265 266 716 +318 266 717 +319 318 717 +267 319 717 +266 267 717 +319 267 718 +320 319 718 +268 320 718 +267 268 718 +320 268 719 +321 320 719 +269 321 719 +268 269 719 +321 269 720 +322 321 720 +270 322 720 +269 270 720 +322 270 721 +323 322 721 +271 323 721 +270 271 721 +323 271 722 +324 323 722 +272 324 722 +271 272 722 +324 272 723 +325 324 723 +273 325 723 +272 273 723 +325 273 724 +326 325 724 +274 326 724 +273 274 724 +326 274 725 +327 326 725 +275 327 725 +274 275 725 +327 275 726 +328 327 726 +276 328 726 +275 276 726 +328 276 727 +329 328 727 +277 329 727 +276 277 727 +329 277 728 +330 329 728 +278 330 728 +277 278 728 +330 278 729 +331 330 729 +279 331 729 +278 279 729 +331 279 730 +332 331 730 +280 332 730 +279 280 730 +332 280 731 +333 332 731 +281 333 731 +280 281 731 +333 281 732 +334 333 732 +282 334 732 +281 282 732 +334 282 733 +335 334 733 +283 335 733 +282 283 733 +335 283 734 +336 335 734 +284 336 734 +283 284 734 +336 284 735 +337 336 735 +285 337 735 +284 285 735 +337 285 736 +338 337 736 +286 338 736 +285 286 736 +338 286 737 +339 338 737 +287 339 737 +286 287 737 +339 287 738 +340 339 738 +288 340 738 +287 288 738 +340 288 739 +341 340 739 +289 341 739 +288 289 739 +341 289 740 +342 341 740 +290 342 740 +289 290 740 +342 290 741 +343 342 741 +291 343 741 +290 291 741 +343 291 742 +344 343 742 +292 344 742 +291 292 742 +344 292 743 +345 344 743 +293 345 743 +292 293 743 +345 293 744 +346 345 744 +294 346 744 +293 294 744 +346 294 745 +347 346 745 +295 347 745 +294 295 745 +347 295 746 +348 347 746 +296 348 746 +295 296 746 +348 296 747 +349 348 747 +297 349 747 +296 297 747 +349 297 748 +350 349 748 +298 350 748 +297 298 748 +350 298 749 +351 350 749 +299 351 749 +298 299 749 +351 299 750 +352 351 750 +300 352 750 +299 300 750 +352 300 751 +353 352 751 +301 353 751 +300 301 751 +353 301 752 +354 353 752 +302 354 752 +301 302 752 +354 302 753 +355 354 753 +303 355 753 +302 303 753 +356 305 754 +357 356 754 +306 357 754 +305 306 754 +357 306 755 +358 357 755 +307 358 755 +306 307 755 +358 307 756 +359 358 756 +308 359 756 +307 308 756 +359 308 757 +360 359 757 +309 360 757 +308 309 757 +360 309 758 +361 360 758 +310 361 758 +309 310 758 +361 310 759 +362 361 759 +311 362 759 +310 311 759 +362 311 760 +363 362 760 +312 363 760 +311 312 760 +363 312 761 +364 363 761 +313 364 761 +312 313 761 +364 313 762 +365 364 762 +314 365 762 +313 314 762 +365 314 763 +366 365 763 +315 366 763 +314 315 763 +366 315 764 +367 366 764 +316 367 764 +315 316 764 +367 316 765 +368 367 765 +317 368 765 +316 317 765 +368 317 766 +369 368 766 +318 369 766 +317 318 766 +369 318 767 +370 369 767 +319 370 767 +318 319 767 +370 319 768 +371 370 768 +320 371 768 +319 320 768 +371 320 769 +372 371 769 +321 372 769 +320 321 769 +372 321 770 +373 372 770 +322 373 770 +321 322 770 +373 322 771 +374 373 771 +323 374 771 +322 323 771 +374 323 772 +375 374 772 +324 375 772 +323 324 772 +375 324 773 +376 375 773 +325 376 773 +324 325 773 +376 325 774 +377 376 774 +326 377 774 +325 326 774 +377 326 775 +378 377 775 +327 378 775 +326 327 775 +378 327 776 +379 378 776 +328 379 776 +327 328 776 +379 328 777 +380 379 777 +329 380 777 +328 329 777 +380 329 778 +381 380 778 +330 381 778 +329 330 778 +381 330 779 +382 381 779 +331 382 779 +330 331 779 +382 331 780 +383 382 780 +332 383 780 +331 332 780 +383 332 781 +384 383 781 +333 384 781 +332 333 781 +384 333 782 +385 384 782 +334 385 782 +333 334 782 +385 334 783 +386 385 783 +335 386 783 +334 335 783 +386 335 784 +387 386 784 +336 387 784 +335 336 784 +387 336 785 +388 387 785 +337 388 785 +336 337 785 +388 337 786 +389 388 786 +338 389 786 +337 338 786 +389 338 787 +390 389 787 +339 390 787 +338 339 787 +390 339 788 +391 390 788 +340 391 788 +339 340 788 +391 340 789 +392 391 789 +341 392 789 +340 341 789 +392 341 790 +393 392 790 +342 393 790 +341 342 790 +393 342 791 +394 393 791 +343 394 791 +342 343 791 +394 343 792 +395 394 792 +344 395 792 +343 344 792 +395 344 793 +396 395 793 +345 396 793 +344 345 793 +396 345 794 +397 396 794 +346 397 794 +345 346 794 +397 346 795 +398 397 795 +347 398 795 +346 347 795 +398 347 796 +399 398 796 +348 399 796 +347 348 796 +399 348 797 +400 399 797 +349 400 797 +348 349 797 +400 349 798 +401 400 798 +350 401 798 +349 350 798 +401 350 799 +402 401 799 +351 402 799 +350 351 799 +402 351 800 +403 402 800 +352 403 800 +351 352 800 +403 352 801 +404 403 801 +353 404 801 +352 353 801 +404 353 802 +405 404 802 +354 405 802 +353 354 802 +405 354 803 +406 405 803 +355 406 803 +354 355 803 +407 356 804 +408 407 804 +357 408 804 +356 357 804 +408 357 805 +409 408 805 +358 409 805 +357 358 805 +409 358 806 +410 409 806 +359 410 806 +358 359 806 +410 359 807 +411 410 807 +360 411 807 +359 360 807 +411 360 808 +412 411 808 +361 412 808 +360 361 808 +412 361 809 +413 412 809 +362 413 809 +361 362 809 +413 362 810 +414 413 810 +363 414 810 +362 363 810 +414 363 811 +415 414 811 +364 415 811 +363 364 811 +415 364 812 +416 415 812 +365 416 812 +364 365 812 +416 365 813 +417 416 813 +366 417 813 +365 366 813 +417 366 814 +418 417 814 +367 418 814 +366 367 814 +418 367 815 +419 418 815 +368 419 815 +367 368 815 +419 368 816 +420 419 816 +369 420 816 +368 369 816 +420 369 817 +421 420 817 +370 421 817 +369 370 817 +421 370 818 +422 421 818 +371 422 818 +370 371 818 +422 371 819 +423 422 819 +372 423 819 +371 372 819 +423 372 820 +424 423 820 +373 424 820 +372 373 820 +424 373 821 +425 424 821 +374 425 821 +373 374 821 +425 374 822 +426 425 822 +375 426 822 +374 375 822 +426 375 823 +427 426 823 +376 427 823 +375 376 823 +427 376 824 +428 427 824 +377 428 824 +376 377 824 +428 377 825 +429 428 825 +378 429 825 +377 378 825 +429 378 826 +430 429 826 +379 430 826 +378 379 826 +430 379 827 +431 430 827 +380 431 827 +379 380 827 +431 380 828 +432 431 828 +381 432 828 +380 381 828 +432 381 829 +433 432 829 +382 433 829 +381 382 829 +433 382 830 +434 433 830 +383 434 830 +382 383 830 +434 383 831 +435 434 831 +384 435 831 +383 384 831 +435 384 832 +436 435 832 +385 436 832 +384 385 832 +436 385 833 +437 436 833 +386 437 833 +385 386 833 +437 386 834 +438 437 834 +387 438 834 +386 387 834 +438 387 835 +439 438 835 +388 439 835 +387 388 835 +439 388 836 +440 439 836 +389 440 836 +388 389 836 +440 389 837 +441 440 837 +390 441 837 +389 390 837 +441 390 838 +442 441 838 +391 442 838 +390 391 838 +442 391 839 +443 442 839 +392 443 839 +391 392 839 +443 392 840 +444 443 840 +393 444 840 +392 393 840 +444 393 841 +445 444 841 +394 445 841 +393 394 841 +445 394 842 +446 445 842 +395 446 842 +394 395 842 +446 395 843 +447 446 843 +396 447 843 +395 396 843 +447 396 844 +448 447 844 +397 448 844 +396 397 844 +448 397 845 +449 448 845 +398 449 845 +397 398 845 +449 398 846 +450 449 846 +399 450 846 +398 399 846 +450 399 847 +451 450 847 +400 451 847 +399 400 847 +451 400 848 +452 451 848 +401 452 848 +400 401 848 +452 401 849 +453 452 849 +402 453 849 +401 402 849 +453 402 850 +454 453 850 +403 454 850 +402 403 850 +454 403 851 +455 454 851 +404 455 851 +403 404 851 +455 404 852 +456 455 852 +405 456 852 +404 405 852 +456 405 853 +457 456 853 +406 457 853 +405 406 853 +257 207 854 +304 257 854 +258 304 855 +858 258 855 +158 208 857 +157 158 857 +156 157 856 +207 156 856 + +element boundaries: +0 0 2 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 0 +0 0 4 +0 0 0 +0 0 0 +0 0 3 +0 0 4 +0 0 2 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 3 +0 0 0 +0 0 2 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 3 +0 0 0 +0 0 2 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 1 +0 0 0 +0 0 0 +0 0 0 +0 0 1 +0 0 0 +0 0 0 +0 0 0 +0 0 1 +0 0 0 +0 0 0 +0 0 0 +0 0 1 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 3 +0 0 0 +0 0 2 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 1 +0 0 0 +0 0 0 +0 0 0 +0 0 1 +0 0 0 +0 0 0 +0 0 0 +0 0 1 +0 0 0 +0 0 0 +0 0 0 +0 0 1 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 3 +0 0 0 +0 0 2 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 3 +0 0 0 +0 0 2 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 0 +0 0 3 +0 0 0 +0 0 2 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 0 +0 0 0 +0 0 0 +0 0 5 +0 0 3 +0 0 0 +1 0 0 +0 1 0 +1 0 0 +0 1 0 +1 0 0 +0 1 0 +1 0 0 +0 1 0 + +element neighbours: +3 1 -1 +0 2 203 +1 3 4 +2 0 -1 +7 5 2 +4 6 207 +5 7 8 +6 4 -1 +11 9 6 +8 10 211 +9 11 12 +10 8 -1 +15 13 10 +12 14 215 +13 15 16 +14 12 -1 +19 17 14 +16 18 219 +17 19 20 +18 16 -1 +23 21 18 +20 22 223 +21 23 24 +22 20 -1 +27 25 22 +24 26 227 +25 27 28 +26 24 -1 +31 29 26 +28 30 231 +29 31 32 +30 28 -1 +35 33 30 +32 34 235 +33 35 36 +34 32 -1 +39 37 34 +36 38 239 +37 39 40 +38 36 -1 +43 41 38 +40 42 243 +41 43 44 +42 40 -1 +47 45 42 +44 46 247 +45 47 48 +46 44 -1 +51 49 46 +48 50 251 +49 51 52 +50 48 -1 +55 53 50 +52 54 255 +53 55 56 +54 52 -1 +59 57 54 +56 58 259 +57 59 60 +58 56 -1 +63 61 58 +60 62 263 +61 63 64 +62 60 -1 +67 65 62 +64 66 267 +65 67 68 +66 64 -1 +71 69 66 +68 70 271 +69 71 72 +70 68 -1 +75 73 70 +72 74 275 +73 75 76 +74 72 -1 +79 77 74 +76 78 279 +77 79 80 +78 76 -1 +83 81 78 +80 82 283 +81 83 84 +82 80 -1 +87 85 82 +84 86 287 +85 87 88 +86 84 -1 +91 89 86 +88 90 291 +89 91 92 +90 88 -1 +95 93 90 +92 94 295 +93 95 96 +94 92 -1 +99 97 94 +96 98 299 +97 99 100 +98 96 -1 +103 101 98 +100 102 303 +101 103 104 +102 100 -1 +107 105 102 +104 106 307 +105 107 108 +106 104 -1 +111 109 106 +108 110 311 +109 111 112 +110 108 -1 +115 113 110 +112 114 315 +113 115 116 +114 112 -1 +119 117 114 +116 118 319 +117 119 120 +118 116 -1 +123 121 118 +120 122 323 +121 123 124 +122 120 -1 +127 125 122 +124 126 327 +125 127 128 +126 124 -1 +131 129 126 +128 130 331 +129 131 132 +130 128 -1 +135 133 130 +132 134 335 +133 135 136 +134 132 -1 +139 137 134 +136 138 339 +137 139 140 +138 136 -1 +143 141 138 +140 142 343 +141 143 144 +142 140 -1 +147 145 142 +144 146 347 +145 147 148 +146 144 -1 +151 149 146 +148 150 351 +149 151 152 +150 148 -1 +155 153 150 +152 154 355 +153 155 156 +154 152 -1 +159 157 154 +156 158 359 +157 159 160 +158 156 -1 +163 161 158 +160 162 363 +161 163 164 +162 160 -1 +167 165 162 +164 166 367 +165 167 168 +166 164 -1 +171 169 166 +168 170 371 +169 171 172 +170 168 -1 +175 173 170 +172 174 375 +173 175 176 +174 172 -1 +179 177 174 +176 178 379 +177 179 180 +178 176 -1 +183 181 178 +180 182 383 +181 183 184 +182 180 -1 +187 185 182 +184 186 387 +185 187 188 +186 184 -1 +191 189 186 +188 190 391 +189 191 192 +190 188 -1 +195 193 190 +192 194 395 +193 195 196 +194 192 -1 +199 197 194 +196 198 399 +197 199 -1 +198 196 -1 +203 201 -1 +200 202 403 +201 203 204 +202 200 1 +207 205 202 +204 206 407 +205 207 208 +206 204 5 +211 209 206 +208 210 411 +209 211 212 +210 208 9 +215 213 210 +212 214 415 +213 215 216 +214 212 13 +219 217 214 +216 218 419 +217 219 220 +218 216 17 +223 221 218 +220 222 423 +221 223 224 +222 220 21 +227 225 222 +224 226 427 +225 227 228 +226 224 25 +231 229 226 +228 230 431 +229 231 232 +230 228 29 +235 233 230 +232 234 435 +233 235 236 +234 232 33 +239 237 234 +236 238 439 +237 239 240 +238 236 37 +243 241 238 +240 242 443 +241 243 244 +242 240 41 +247 245 242 +244 246 447 +245 247 248 +246 244 45 +251 249 246 +248 250 451 +249 251 252 +250 248 49 +255 253 250 +252 254 455 +253 255 256 +254 252 53 +259 257 254 +256 258 459 +257 259 260 +258 256 57 +263 261 258 +260 262 463 +261 263 264 +262 260 61 +267 265 262 +264 266 467 +265 267 268 +266 264 65 +271 269 266 +268 270 471 +269 271 272 +270 268 69 +275 273 270 +272 274 475 +273 275 276 +274 272 73 +279 277 274 +276 278 479 +277 279 280 +278 276 77 +283 281 278 +280 282 483 +281 283 284 +282 280 81 +287 285 282 +284 286 487 +285 287 288 +286 284 85 +291 289 286 +288 290 491 +289 291 292 +290 288 89 +295 293 290 +292 294 495 +293 295 296 +294 292 93 +299 297 294 +296 298 499 +297 299 300 +298 296 97 +303 301 298 +300 302 503 +301 303 304 +302 300 101 +307 305 302 +304 306 507 +305 307 308 +306 304 105 +311 309 306 +308 310 511 +309 311 312 +310 308 109 +315 313 310 +312 314 515 +313 315 316 +314 312 113 +319 317 314 +316 318 519 +317 319 320 +318 316 117 +323 321 318 +320 322 523 +321 323 324 +322 320 121 +327 325 322 +324 326 527 +325 327 328 +326 324 125 +331 329 326 +328 330 531 +329 331 332 +330 328 129 +335 333 330 +332 334 535 +333 335 336 +334 332 133 +339 337 334 +336 338 539 +337 339 340 +338 336 137 +343 341 338 +340 342 543 +341 343 344 +342 340 141 +347 345 342 +344 346 547 +345 347 348 +346 344 145 +351 349 346 +348 350 551 +349 351 352 +350 348 149 +355 353 350 +352 354 555 +353 355 356 +354 352 153 +359 357 354 +356 358 559 +357 359 360 +358 356 157 +363 361 358 +360 362 563 +361 363 364 +362 360 161 +367 365 362 +364 366 567 +365 367 368 +366 364 165 +371 369 366 +368 370 571 +369 371 372 +370 368 169 +375 373 370 +372 374 575 +373 375 376 +374 372 173 +379 377 374 +376 378 579 +377 379 380 +378 376 177 +383 381 378 +380 382 583 +381 383 384 +382 380 181 +387 385 382 +384 386 587 +385 387 388 +386 384 185 +391 389 386 +388 390 591 +389 391 392 +390 388 189 +395 393 390 +392 394 595 +393 395 396 +394 392 193 +399 397 394 +396 398 599 +397 399 -1 +398 396 197 +403 401 -1 +400 402 603 +401 403 404 +402 400 201 +407 405 402 +404 406 607 +405 407 408 +406 404 205 +411 409 406 +408 410 611 +409 411 412 +410 408 209 +415 413 410 +412 414 1590 +413 415 416 +414 412 213 +419 417 414 +416 418 1589 +417 419 420 +418 416 217 +423 421 418 +420 422 615 +421 423 424 +422 420 221 +427 425 422 +424 426 619 +425 427 428 +426 424 225 +431 429 426 +428 430 623 +429 431 432 +430 428 229 +435 433 430 +432 434 627 +433 435 436 +434 432 233 +439 437 434 +436 438 631 +437 439 440 +438 436 237 +443 441 438 +440 442 635 +441 443 444 +442 440 241 +447 445 442 +444 446 639 +445 447 448 +446 444 245 +451 449 446 +448 450 643 +449 451 452 +450 448 249 +455 453 450 +452 454 647 +453 455 456 +454 452 253 +459 457 454 +456 458 651 +457 459 460 +458 456 257 +463 461 458 +460 462 655 +461 463 464 +462 460 261 +467 465 462 +464 466 659 +465 467 468 +466 464 265 +471 469 466 +468 470 663 +469 471 472 +470 468 269 +475 473 470 +472 474 667 +473 475 476 +474 472 273 +479 477 474 +476 478 671 +477 479 480 +478 476 277 +483 481 478 +480 482 675 +481 483 484 +482 480 281 +487 485 482 +484 486 679 +485 487 488 +486 484 285 +491 489 486 +488 490 683 +489 491 492 +490 488 289 +495 493 490 +492 494 687 +493 495 496 +494 492 293 +499 497 494 +496 498 691 +497 499 500 +498 496 297 +503 501 498 +500 502 695 +501 503 504 +502 500 301 +507 505 502 +504 506 699 +505 507 508 +506 504 305 +511 509 506 +508 510 703 +509 511 512 +510 508 309 +515 513 510 +512 514 707 +513 515 516 +514 512 313 +519 517 514 +516 518 711 +517 519 520 +518 516 317 +523 521 518 +520 522 715 +521 523 524 +522 520 321 +527 525 522 +524 526 719 +525 527 528 +526 524 325 +531 529 526 +528 530 723 +529 531 532 +530 528 329 +535 533 530 +532 534 727 +533 535 536 +534 532 333 +539 537 534 +536 538 731 +537 539 540 +538 536 337 +543 541 538 +540 542 735 +541 543 544 +542 540 341 +547 545 542 +544 546 739 +545 547 548 +546 544 345 +551 549 546 +548 550 743 +549 551 552 +550 548 349 +555 553 550 +552 554 747 +553 555 556 +554 552 353 +559 557 554 +556 558 751 +557 559 560 +558 556 357 +563 561 558 +560 562 755 +561 563 564 +562 560 361 +567 565 562 +564 566 759 +565 567 568 +566 564 365 +571 569 566 +568 570 763 +569 571 572 +570 568 369 +575 573 570 +572 574 767 +573 575 576 +574 572 373 +579 577 574 +576 578 771 +577 579 580 +578 576 377 +583 581 578 +580 582 775 +581 583 584 +582 580 381 +587 585 582 +584 586 779 +585 587 588 +586 584 385 +591 589 586 +588 590 783 +589 591 592 +590 588 389 +595 593 590 +592 594 787 +593 595 596 +594 592 393 +599 597 594 +596 598 791 +597 599 -1 +598 596 397 +603 601 -1 +600 602 795 +601 603 604 +602 600 401 +607 605 602 +604 606 799 +605 607 608 +606 604 405 +611 609 606 +608 610 803 +609 611 1591 +610 608 409 +615 613 1588 +612 614 -1 +613 615 616 +614 612 421 +619 617 614 +616 618 -1 +617 619 620 +618 616 425 +623 621 618 +620 622 -1 +621 623 624 +622 620 429 +627 625 622 +624 626 -1 +625 627 628 +626 624 433 +631 629 626 +628 630 823 +629 631 632 +630 628 437 +635 633 630 +632 634 827 +633 635 636 +634 632 441 +639 637 634 +636 638 831 +637 639 640 +638 636 445 +643 641 638 +640 642 835 +641 643 644 +642 640 449 +647 645 642 +644 646 839 +645 647 648 +646 644 453 +651 649 646 +648 650 843 +649 651 652 +650 648 457 +655 653 650 +652 654 847 +653 655 656 +654 652 461 +659 657 654 +656 658 851 +657 659 660 +658 656 465 +663 661 658 +660 662 855 +661 663 664 +662 660 469 +667 665 662 +664 666 859 +665 667 668 +666 664 473 +671 669 666 +668 670 863 +669 671 672 +670 668 477 +675 673 670 +672 674 867 +673 675 676 +674 672 481 +679 677 674 +676 678 871 +677 679 680 +678 676 485 +683 681 678 +680 682 875 +681 683 684 +682 680 489 +687 685 682 +684 686 879 +685 687 688 +686 684 493 +691 689 686 +688 690 883 +689 691 692 +690 688 497 +695 693 690 +692 694 887 +693 695 696 +694 692 501 +699 697 694 +696 698 891 +697 699 700 +698 696 505 +703 701 698 +700 702 895 +701 703 704 +702 700 509 +707 705 702 +704 706 899 +705 707 708 +706 704 513 +711 709 706 +708 710 903 +709 711 712 +710 708 517 +715 713 710 +712 714 907 +713 715 716 +714 712 521 +719 717 714 +716 718 911 +717 719 720 +718 716 525 +723 721 718 +720 722 915 +721 723 724 +722 720 529 +727 725 722 +724 726 919 +725 727 728 +726 724 533 +731 729 726 +728 730 923 +729 731 732 +730 728 537 +735 733 730 +732 734 927 +733 735 736 +734 732 541 +739 737 734 +736 738 931 +737 739 740 +738 736 545 +743 741 738 +740 742 935 +741 743 744 +742 740 549 +747 745 742 +744 746 939 +745 747 748 +746 744 553 +751 749 746 +748 750 943 +749 751 752 +750 748 557 +755 753 750 +752 754 947 +753 755 756 +754 752 561 +759 757 754 +756 758 951 +757 759 760 +758 756 565 +763 761 758 +760 762 955 +761 763 764 +762 760 569 +767 765 762 +764 766 959 +765 767 768 +766 764 573 +771 769 766 +768 770 963 +769 771 772 +770 768 577 +775 773 770 +772 774 967 +773 775 776 +774 772 581 +779 777 774 +776 778 971 +777 779 780 +778 776 585 +783 781 778 +780 782 975 +781 783 784 +782 780 589 +787 785 782 +784 786 979 +785 787 788 +786 784 593 +791 789 786 +788 790 983 +789 791 -1 +790 788 597 +795 793 -1 +792 794 987 +793 795 796 +794 792 601 +799 797 794 +796 798 991 +797 799 800 +798 796 605 +803 801 798 +800 802 995 +801 803 1584 +802 800 609 +807 805 1587 +804 806 1007 +805 807 808 +806 804 -1 +811 809 806 +808 810 1011 +809 811 812 +810 808 -1 +815 813 810 +812 814 1015 +813 815 816 +814 812 -1 +819 817 814 +816 818 1019 +817 819 820 +818 816 -1 +823 821 818 +820 822 1023 +821 823 824 +822 820 629 +827 825 822 +824 826 1027 +825 827 828 +826 824 633 +831 829 826 +828 830 1031 +829 831 832 +830 828 637 +835 833 830 +832 834 1035 +833 835 836 +834 832 641 +839 837 834 +836 838 1039 +837 839 840 +838 836 645 +843 841 838 +840 842 1043 +841 843 844 +842 840 649 +847 845 842 +844 846 1047 +845 847 848 +846 844 653 +851 849 846 +848 850 1051 +849 851 852 +850 848 657 +855 853 850 +852 854 1055 +853 855 856 +854 852 661 +859 857 854 +856 858 1059 +857 859 860 +858 856 665 +863 861 858 +860 862 1063 +861 863 864 +862 860 669 +867 865 862 +864 866 1067 +865 867 868 +866 864 673 +871 869 866 +868 870 1071 +869 871 872 +870 868 677 +875 873 870 +872 874 1075 +873 875 876 +874 872 681 +879 877 874 +876 878 1079 +877 879 880 +878 876 685 +883 881 878 +880 882 1083 +881 883 884 +882 880 689 +887 885 882 +884 886 1087 +885 887 888 +886 884 693 +891 889 886 +888 890 1091 +889 891 892 +890 888 697 +895 893 890 +892 894 1095 +893 895 896 +894 892 701 +899 897 894 +896 898 1099 +897 899 900 +898 896 705 +903 901 898 +900 902 1103 +901 903 904 +902 900 709 +907 905 902 +904 906 1107 +905 907 908 +906 904 713 +911 909 906 +908 910 1111 +909 911 912 +910 908 717 +915 913 910 +912 914 1115 +913 915 916 +914 912 721 +919 917 914 +916 918 1119 +917 919 920 +918 916 725 +923 921 918 +920 922 1123 +921 923 924 +922 920 729 +927 925 922 +924 926 1127 +925 927 928 +926 924 733 +931 929 926 +928 930 1131 +929 931 932 +930 928 737 +935 933 930 +932 934 1135 +933 935 936 +934 932 741 +939 937 934 +936 938 1139 +937 939 940 +938 936 745 +943 941 938 +940 942 1143 +941 943 944 +942 940 749 +947 945 942 +944 946 1147 +945 947 948 +946 944 753 +951 949 946 +948 950 1151 +949 951 952 +950 948 757 +955 953 950 +952 954 1155 +953 955 956 +954 952 761 +959 957 954 +956 958 1159 +957 959 960 +958 956 765 +963 961 958 +960 962 1163 +961 963 964 +962 960 769 +967 965 962 +964 966 1167 +965 967 968 +966 964 773 +971 969 966 +968 970 1171 +969 971 972 +970 968 777 +975 973 970 +972 974 1175 +973 975 976 +974 972 781 +979 977 974 +976 978 1179 +977 979 980 +978 976 785 +983 981 978 +980 982 1183 +981 983 -1 +982 980 789 +987 985 -1 +984 986 1187 +985 987 988 +986 984 793 +991 989 986 +988 990 1191 +989 991 992 +990 988 797 +995 993 990 +992 994 1195 +993 995 996 +994 992 801 +999 997 994 +996 998 1199 +997 999 1000 +998 996 1585 +1003 1001 998 +1000 1002 1203 +1001 1003 1004 +1002 1000 1586 +1007 1005 1002 +1004 1006 1207 +1005 1007 1008 +1006 1004 805 +1011 1009 1006 +1008 1010 1211 +1009 1011 1012 +1010 1008 809 +1015 1013 1010 +1012 1014 1215 +1013 1015 1016 +1014 1012 813 +1019 1017 1014 +1016 1018 1219 +1017 1019 1020 +1018 1016 817 +1023 1021 1018 +1020 1022 1223 +1021 1023 1024 +1022 1020 821 +1027 1025 1022 +1024 1026 1227 +1025 1027 1028 +1026 1024 825 +1031 1029 1026 +1028 1030 1231 +1029 1031 1032 +1030 1028 829 +1035 1033 1030 +1032 1034 1235 +1033 1035 1036 +1034 1032 833 +1039 1037 1034 +1036 1038 1239 +1037 1039 1040 +1038 1036 837 +1043 1041 1038 +1040 1042 1243 +1041 1043 1044 +1042 1040 841 +1047 1045 1042 +1044 1046 1247 +1045 1047 1048 +1046 1044 845 +1051 1049 1046 +1048 1050 1251 +1049 1051 1052 +1050 1048 849 +1055 1053 1050 +1052 1054 1255 +1053 1055 1056 +1054 1052 853 +1059 1057 1054 +1056 1058 1259 +1057 1059 1060 +1058 1056 857 +1063 1061 1058 +1060 1062 1263 +1061 1063 1064 +1062 1060 861 +1067 1065 1062 +1064 1066 1267 +1065 1067 1068 +1066 1064 865 +1071 1069 1066 +1068 1070 1271 +1069 1071 1072 +1070 1068 869 +1075 1073 1070 +1072 1074 1275 +1073 1075 1076 +1074 1072 873 +1079 1077 1074 +1076 1078 1279 +1077 1079 1080 +1078 1076 877 +1083 1081 1078 +1080 1082 1283 +1081 1083 1084 +1082 1080 881 +1087 1085 1082 +1084 1086 1287 +1085 1087 1088 +1086 1084 885 +1091 1089 1086 +1088 1090 1291 +1089 1091 1092 +1090 1088 889 +1095 1093 1090 +1092 1094 1295 +1093 1095 1096 +1094 1092 893 +1099 1097 1094 +1096 1098 1299 +1097 1099 1100 +1098 1096 897 +1103 1101 1098 +1100 1102 1303 +1101 1103 1104 +1102 1100 901 +1107 1105 1102 +1104 1106 1307 +1105 1107 1108 +1106 1104 905 +1111 1109 1106 +1108 1110 1311 +1109 1111 1112 +1110 1108 909 +1115 1113 1110 +1112 1114 1315 +1113 1115 1116 +1114 1112 913 +1119 1117 1114 +1116 1118 1319 +1117 1119 1120 +1118 1116 917 +1123 1121 1118 +1120 1122 1323 +1121 1123 1124 +1122 1120 921 +1127 1125 1122 +1124 1126 1327 +1125 1127 1128 +1126 1124 925 +1131 1129 1126 +1128 1130 1331 +1129 1131 1132 +1130 1128 929 +1135 1133 1130 +1132 1134 1335 +1133 1135 1136 +1134 1132 933 +1139 1137 1134 +1136 1138 1339 +1137 1139 1140 +1138 1136 937 +1143 1141 1138 +1140 1142 1343 +1141 1143 1144 +1142 1140 941 +1147 1145 1142 +1144 1146 1347 +1145 1147 1148 +1146 1144 945 +1151 1149 1146 +1148 1150 1351 +1149 1151 1152 +1150 1148 949 +1155 1153 1150 +1152 1154 1355 +1153 1155 1156 +1154 1152 953 +1159 1157 1154 +1156 1158 1359 +1157 1159 1160 +1158 1156 957 +1163 1161 1158 +1160 1162 1363 +1161 1163 1164 +1162 1160 961 +1167 1165 1162 +1164 1166 1367 +1165 1167 1168 +1166 1164 965 +1171 1169 1166 +1168 1170 1371 +1169 1171 1172 +1170 1168 969 +1175 1173 1170 +1172 1174 1375 +1173 1175 1176 +1174 1172 973 +1179 1177 1174 +1176 1178 1379 +1177 1179 1180 +1178 1176 977 +1183 1181 1178 +1180 1182 1383 +1181 1183 -1 +1182 1180 981 +1187 1185 -1 +1184 1186 1387 +1185 1187 1188 +1186 1184 985 +1191 1189 1186 +1188 1190 1391 +1189 1191 1192 +1190 1188 989 +1195 1193 1190 +1192 1194 1395 +1193 1195 1196 +1194 1192 993 +1199 1197 1194 +1196 1198 1399 +1197 1199 1200 +1198 1196 997 +1203 1201 1198 +1200 1202 1403 +1201 1203 1204 +1202 1200 1001 +1207 1205 1202 +1204 1206 1407 +1205 1207 1208 +1206 1204 1005 +1211 1209 1206 +1208 1210 1411 +1209 1211 1212 +1210 1208 1009 +1215 1213 1210 +1212 1214 1415 +1213 1215 1216 +1214 1212 1013 +1219 1217 1214 +1216 1218 1419 +1217 1219 1220 +1218 1216 1017 +1223 1221 1218 +1220 1222 1423 +1221 1223 1224 +1222 1220 1021 +1227 1225 1222 +1224 1226 1427 +1225 1227 1228 +1226 1224 1025 +1231 1229 1226 +1228 1230 1431 +1229 1231 1232 +1230 1228 1029 +1235 1233 1230 +1232 1234 1435 +1233 1235 1236 +1234 1232 1033 +1239 1237 1234 +1236 1238 1439 +1237 1239 1240 +1238 1236 1037 +1243 1241 1238 +1240 1242 1443 +1241 1243 1244 +1242 1240 1041 +1247 1245 1242 +1244 1246 1447 +1245 1247 1248 +1246 1244 1045 +1251 1249 1246 +1248 1250 1451 +1249 1251 1252 +1250 1248 1049 +1255 1253 1250 +1252 1254 1455 +1253 1255 1256 +1254 1252 1053 +1259 1257 1254 +1256 1258 1459 +1257 1259 1260 +1258 1256 1057 +1263 1261 1258 +1260 1262 1463 +1261 1263 1264 +1262 1260 1061 +1267 1265 1262 +1264 1266 1467 +1265 1267 1268 +1266 1264 1065 +1271 1269 1266 +1268 1270 1471 +1269 1271 1272 +1270 1268 1069 +1275 1273 1270 +1272 1274 1475 +1273 1275 1276 +1274 1272 1073 +1279 1277 1274 +1276 1278 1479 +1277 1279 1280 +1278 1276 1077 +1283 1281 1278 +1280 1282 1483 +1281 1283 1284 +1282 1280 1081 +1287 1285 1282 +1284 1286 1487 +1285 1287 1288 +1286 1284 1085 +1291 1289 1286 +1288 1290 1491 +1289 1291 1292 +1290 1288 1089 +1295 1293 1290 +1292 1294 1495 +1293 1295 1296 +1294 1292 1093 +1299 1297 1294 +1296 1298 1499 +1297 1299 1300 +1298 1296 1097 +1303 1301 1298 +1300 1302 1503 +1301 1303 1304 +1302 1300 1101 +1307 1305 1302 +1304 1306 1507 +1305 1307 1308 +1306 1304 1105 +1311 1309 1306 +1308 1310 1511 +1309 1311 1312 +1310 1308 1109 +1315 1313 1310 +1312 1314 1515 +1313 1315 1316 +1314 1312 1113 +1319 1317 1314 +1316 1318 1519 +1317 1319 1320 +1318 1316 1117 +1323 1321 1318 +1320 1322 1523 +1321 1323 1324 +1322 1320 1121 +1327 1325 1322 +1324 1326 1527 +1325 1327 1328 +1326 1324 1125 +1331 1329 1326 +1328 1330 1531 +1329 1331 1332 +1330 1328 1129 +1335 1333 1330 +1332 1334 1535 +1333 1335 1336 +1334 1332 1133 +1339 1337 1334 +1336 1338 1539 +1337 1339 1340 +1338 1336 1137 +1343 1341 1338 +1340 1342 1543 +1341 1343 1344 +1342 1340 1141 +1347 1345 1342 +1344 1346 1547 +1345 1347 1348 +1346 1344 1145 +1351 1349 1346 +1348 1350 1551 +1349 1351 1352 +1350 1348 1149 +1355 1353 1350 +1352 1354 1555 +1353 1355 1356 +1354 1352 1153 +1359 1357 1354 +1356 1358 1559 +1357 1359 1360 +1358 1356 1157 +1363 1361 1358 +1360 1362 1563 +1361 1363 1364 +1362 1360 1161 +1367 1365 1362 +1364 1366 1567 +1365 1367 1368 +1366 1364 1165 +1371 1369 1366 +1368 1370 1571 +1369 1371 1372 +1370 1368 1169 +1375 1373 1370 +1372 1374 1575 +1373 1375 1376 +1374 1372 1173 +1379 1377 1374 +1376 1378 1579 +1377 1379 1380 +1378 1376 1177 +1383 1381 1378 +1380 1382 1583 +1381 1383 -1 +1382 1380 1181 +1387 1385 -1 +1384 1386 -1 +1385 1387 1388 +1386 1384 1185 +1391 1389 1386 +1388 1390 -1 +1389 1391 1392 +1390 1388 1189 +1395 1393 1390 +1392 1394 -1 +1393 1395 1396 +1394 1392 1193 +1399 1397 1394 +1396 1398 -1 +1397 1399 1400 +1398 1396 1197 +1403 1401 1398 +1400 1402 -1 +1401 1403 1404 +1402 1400 1201 +1407 1405 1402 +1404 1406 -1 +1405 1407 1408 +1406 1404 1205 +1411 1409 1406 +1408 1410 -1 +1409 1411 1412 +1410 1408 1209 +1415 1413 1410 +1412 1414 -1 +1413 1415 1416 +1414 1412 1213 +1419 1417 1414 +1416 1418 -1 +1417 1419 1420 +1418 1416 1217 +1423 1421 1418 +1420 1422 -1 +1421 1423 1424 +1422 1420 1221 +1427 1425 1422 +1424 1426 -1 +1425 1427 1428 +1426 1424 1225 +1431 1429 1426 +1428 1430 -1 +1429 1431 1432 +1430 1428 1229 +1435 1433 1430 +1432 1434 -1 +1433 1435 1436 +1434 1432 1233 +1439 1437 1434 +1436 1438 -1 +1437 1439 1440 +1438 1436 1237 +1443 1441 1438 +1440 1442 -1 +1441 1443 1444 +1442 1440 1241 +1447 1445 1442 +1444 1446 -1 +1445 1447 1448 +1446 1444 1245 +1451 1449 1446 +1448 1450 -1 +1449 1451 1452 +1450 1448 1249 +1455 1453 1450 +1452 1454 -1 +1453 1455 1456 +1454 1452 1253 +1459 1457 1454 +1456 1458 -1 +1457 1459 1460 +1458 1456 1257 +1463 1461 1458 +1460 1462 -1 +1461 1463 1464 +1462 1460 1261 +1467 1465 1462 +1464 1466 -1 +1465 1467 1468 +1466 1464 1265 +1471 1469 1466 +1468 1470 -1 +1469 1471 1472 +1470 1468 1269 +1475 1473 1470 +1472 1474 -1 +1473 1475 1476 +1474 1472 1273 +1479 1477 1474 +1476 1478 -1 +1477 1479 1480 +1478 1476 1277 +1483 1481 1478 +1480 1482 -1 +1481 1483 1484 +1482 1480 1281 +1487 1485 1482 +1484 1486 -1 +1485 1487 1488 +1486 1484 1285 +1491 1489 1486 +1488 1490 -1 +1489 1491 1492 +1490 1488 1289 +1495 1493 1490 +1492 1494 -1 +1493 1495 1496 +1494 1492 1293 +1499 1497 1494 +1496 1498 -1 +1497 1499 1500 +1498 1496 1297 +1503 1501 1498 +1500 1502 -1 +1501 1503 1504 +1502 1500 1301 +1507 1505 1502 +1504 1506 -1 +1505 1507 1508 +1506 1504 1305 +1511 1509 1506 +1508 1510 -1 +1509 1511 1512 +1510 1508 1309 +1515 1513 1510 +1512 1514 -1 +1513 1515 1516 +1514 1512 1313 +1519 1517 1514 +1516 1518 -1 +1517 1519 1520 +1518 1516 1317 +1523 1521 1518 +1520 1522 -1 +1521 1523 1524 +1522 1520 1321 +1527 1525 1522 +1524 1526 -1 +1525 1527 1528 +1526 1524 1325 +1531 1529 1526 +1528 1530 -1 +1529 1531 1532 +1530 1528 1329 +1535 1533 1530 +1532 1534 -1 +1533 1535 1536 +1534 1532 1333 +1539 1537 1534 +1536 1538 -1 +1537 1539 1540 +1538 1536 1337 +1543 1541 1538 +1540 1542 -1 +1541 1543 1544 +1542 1540 1341 +1547 1545 1542 +1544 1546 -1 +1545 1547 1548 +1546 1544 1345 +1551 1549 1546 +1548 1550 -1 +1549 1551 1552 +1550 1548 1349 +1555 1553 1550 +1552 1554 -1 +1553 1555 1556 +1554 1552 1353 +1559 1557 1554 +1556 1558 -1 +1557 1559 1560 +1558 1556 1357 +1563 1561 1558 +1560 1562 -1 +1561 1563 1564 +1562 1560 1361 +1567 1565 1562 +1564 1566 -1 +1565 1567 1568 +1566 1564 1365 +1571 1569 1566 +1568 1570 -1 +1569 1571 1572 +1570 1568 1369 +1575 1573 1570 +1572 1574 -1 +1573 1575 1576 +1574 1572 1373 +1579 1577 1574 +1576 1578 -1 +1577 1579 1580 +1578 1576 1377 +1583 1581 1578 +1580 1582 -1 +1581 1583 -1 +1582 1580 1381 +-1 1585 802 +1584 -1 999 +-1 1587 1003 +1586 -1 804 +-1 1589 612 +1588 -1 417 +-1 1591 413 +1590 -1 610 + diff --git a/demo/src/navierStokes.cc b/demo/src/navierStokes.cc new file mode 100644 index 0000000000000000000000000000000000000000..867e9cd112d83372fcccf682a3fbc443cb17fb7a --- /dev/null +++ b/demo/src/navierStokes.cc @@ -0,0 +1,109 @@ +#include "AMDiS.h" +#include "NavierStokes_TaylorHood_RB.h" +#include "navierStokes.h" +#include "time/ExtendedRosenbrockAdaptInstationary.h" +#include "Refinement.h" +#include "MeshFunction_Level.h" + +#include "boost/date_time/posix_time/posix_time.hpp" + +using namespace AMDiS; +using namespace boost::posix_time; + + + +class NS_Channel : public NavierStokes_TaylorHood_RB +{ +public: + NS_Channel(std::string name_) : NavierStokes_TaylorHood_RB(name_) {} + + void solveInitialProblem(AdaptInfo *adaptInfo) + { + super::solveInitialProblem(adaptInfo); + + size_t nVertices = 0; + WorldVector<double> x; + Parameters::get("obstacle->num vertices",nVertices); + std::vector<WorldVector<double> > v(nVertices,x); + for (size_t i = 0; i < nVertices; i++) + Parameters::get("obstacle->vertex["+boost::lexical_cast<std::string>(i)+"]",v[i]); + v.push_back(v[0]); + + SignedDistRefinement refFunction(getMesh()); + RefinementLevelCoords2 refinement( + getFeSpace(), + &refFunction, + new Polygon(v)); + + // initial refinement + refinement.refine(10); + } +protected: + void fillBoundaryConditions() + { FUNCNAME("NS_DrivenCavity::fillBoundaryConditions()"); + + AbstractFunction<double, WorldVector<double> > *zero = new AMDiS::Const<double, WorldVector<double> >(0.0); + size_t dow = Global::getGeo(WORLD); + + // +------ 5 ------+ + // | | + // 2 # <--1 3 + // | | + // +------ 4 ------+ + + +// x[0] = 2.0; +// x[1] = 2.0; +#if DEBUG != 0 + DOFVector<double>* polygon = new DOFVector<double>(getFeSpace(0), "polygon"); + polygon->interpol(new Polygon(v)); + VtkWriter::writeFile(polygon, "polygon.vtu"); +#endif + + /// at rigid wall: no-slip boundary condition + for (size_t i = 0; i < dow; i++) { +// getProblem(0)->addImplicitDirichletBC(*(new Polygon(v)), i, i, *zero); +// getProblem(0)->addSingularDirichletBC(x, i, i, *zero); + getProblem(0)->addDirichletBC(1, i, i, zero); + getProblem(0)->addDirichletBC(4, i, i, zero); + getProblem(0)->addDirichletBC(5, i, i, zero); + } + + double H = 4.1; + double Um = 1.5; + Parameters::get("mesh->H",H); + Parameters::get("ns->Um",Um); + + /// at left wall: prescribed velocity + getProblem(0)->addDirichletBC(2, 0, 0, new InflowBC(H,Um)); + getProblem(0)->addDirichletBC(2, 1, 1, zero); + } +}; + + +int main(int argc, char** argv) +{ FUNCNAME("main"); + + AMDiS::init(argc, argv); + + NS_Channel nsProb("ns"); + nsProb.initialize(INIT_ALL | INIT_EXACT_SOLUTION); + + // Adapt-Infos + AdaptInfo adaptInfo("adapt", nsProb.getNumComponents()); + + // adaption loop - solve ch-prob and ns-prob +// AdaptInstationary adaptInstat("adapt", nsProb, adaptInfo, nsProb, adaptInfo); + ExtendedRosenbrockAdaptInstationary<NS_Channel> adaptInstat("adapt", nsProb, adaptInfo, nsProb, adaptInfo); + + ptime start_time = microsec_clock::local_time(); + nsProb.initTimeInterface(); + int error_code = adaptInstat.adapt(); + time_duration td = microsec_clock::local_time()-start_time; + + MSG("elapsed time= %d sec\n", td.total_seconds()); + + AMDiS::finalize(); + + return error_code; +}; diff --git a/demo/src/navierStokes.h b/demo/src/navierStokes.h new file mode 100644 index 0000000000000000000000000000000000000000..771f050117fdf6441c3da99b6cfe770659d5d7ec --- /dev/null +++ b/demo/src/navierStokes.h @@ -0,0 +1,48 @@ +/** \file navierStokes.h */ + +#ifndef NAVIER_STOKES_H +#define NAVIER_STOKES_H + +#include "AMDiS.h" +#include "GeometryTools.h" + +struct InflowBC : AbstractFunction<double, WorldVector<double> > +{ + InflowBC(double H_=4.1, double Um_=1.5) : H(H_), Um(Um_) {} + double operator()(const WorldVector<double> &x) const { + return 4.0 * Um * x[1] * (H - x[1]) / sqr(H); + } +protected: + double H; + double Um; +}; + + +class Polygon : public AbstractFunction<double, WorldVector<double> > +{ +public: + + Polygon(WorldVector<double> x0_, WorldVector<double> x1_, WorldVector<double> x2_, WorldVector<double> x3_) + { + vertices.push_back(x0_); + vertices.push_back(x1_); + vertices.push_back(x2_); + vertices.push_back(x3_); + vertices.push_back(x0_); + } + + Polygon(std::vector<WorldVector<double> > xi_) : vertices(xi_) { } + + double operator()(const WorldVector<double>& x) const + { + double result = 1.e15; + for (size_t i = 0; i < vertices.size()-1; i++) + result = std::min(result, meshconv::distance_point_line_2d(x.begin(), vertices[i].begin(), vertices[i+1].begin())); + return result * (meshconv::point_in_polygon(x.begin(), vertices) ? -1.0 : 1.0); + }; + +private: + std::vector<WorldVector<double> > vertices; +}; + +#endif diff --git a/demo/src/pfc.cc b/demo/src/pfc.cc index 76220ba8d6f9816fed352d928814f1b6575cdccd..387de288792577d9dc65dd0e1946eda33ce50be1 100644 --- a/demo/src/pfc.cc +++ b/demo/src/pfc.cc @@ -1,6 +1,7 @@ #include "AMDiS.h" #include "PhaseFieldCrystal_Base.h" #include "Helpers.h" +#include "getMaxima.h" #include "boost/date_time/posix_time/posix_time.hpp" @@ -87,6 +88,9 @@ int main(int argc, char** argv) MSG("elapsed time= %d sec\n", td.total_seconds()); + // calc maxima + getMaxima(&adaptInfo, pfcProb.getProblem(), "output"); + AMDiS::finalize(); return error_code;