diff --git a/bin/analyzegrid/CMakeLists.txt b/bin/analyzegrid/CMakeLists.txt index bfc2f7f0a275677bd0905103e9b7b5332c310841..d3322e2914a2272d0c70d7afd88b3d13d11d9e8d 100644 --- a/bin/analyzegrid/CMakeLists.txt +++ b/bin/analyzegrid/CMakeLists.txt @@ -1,4 +1,6 @@ if(consistent-edge-orientation_FOUND) add_executable(analyze_grid analyze_grid.cc) dune_symlink_to_source_files(FILES test_2d_structured.ini test_3d_structured.ini test_2d_unstructured.ini test_3d_unstructured.ini) + + add_executable(consistent_gmsh consistent_gmsh.cc) endif() diff --git a/bin/analyzegrid/analyze_grid.cc b/bin/analyzegrid/analyze_grid.cc index 5ceed213057f56917deafd37cc0f9ccf3722f703..38dc66c357dfe38c7c02b438f1798c99dbb1fdfa 100644 --- a/bin/analyzegrid/analyze_grid.cc +++ b/bin/analyzegrid/analyze_grid.cc @@ -42,6 +42,31 @@ int main(int argc, char** argv){ std::string filename = std::string(argv[2]); std::cout << "writing into: " << filename << std::endl; + bool gmsh = initree.hasKey("gmshFile"); + + if (gmsh){ + // TODO: For now just assume dim = 3 + const int dim = 3; + + // Setup grid (view)... + using Grid = Dune::UGGrid<3>; + using GV = Grid::LeafGridView; + IniGridFactory<Grid> factory(initree); + std::shared_ptr<Grid> grid = factory.getGrid(); + + // Note: As long as createConsistentGrid drops boundary faces we assume + // that the users already provides a consistent gmsh file. + // std::shared_ptr<Grid> grid_nonconsistent = factory.getGrid(); + // std::shared_ptr<Grid> grid = createConsistentGrid(grid_nonconsistent); + + GV gv = grid->leafGridView(); + + // Extract facemod/facedir intersection variation + using ES = Dune::PDELab::AllEntitySet<GV>; + ES es(gv); + std::cout << "Analyse 3d gmsh grid" << std::endl; + analyze_grid(es, filename); + } if (quadrilateral){ if (unstructured){ std::vector<int> tmp; @@ -50,8 +75,6 @@ int main(int argc, char** argv){ std::cout << "dim: " << dim << std::endl; if (dim == 2){ - std::cout << "Analyse 2d unstructured quadrilateral grid" << std::endl; - // Setup grid (view)... using Grid = Dune::UGGrid<2>; using GV = Grid::LeafGridView; @@ -63,11 +86,10 @@ int main(int argc, char** argv){ // Extract facemod/facedir intersection variation using ES = Dune::PDELab::AllEntitySet<GV>; ES es(gv); + std::cout << "Analyse 2d unstructured quadrilateral grid" << std::endl; analyze_grid(es, filename); } else if (dim == 3){ - std::cout << "Analyse 3d unstructured quadrilateral grid" << std::endl; - // Setup grid (view)... using Grid = Dune::UGGrid<3>; using GV = Grid::LeafGridView; @@ -79,6 +101,7 @@ int main(int argc, char** argv){ // Extract facemod/facedir intersection variation using ES = Dune::PDELab::AllEntitySet<GV>; ES es(gv); + std::cout << "Analyse 3d unstructured quadrilateral grid" << std::endl; analyze_grid(es, filename); } else{ @@ -92,8 +115,6 @@ int main(int argc, char** argv){ std::cout << "dim: " << dim << std::endl; if (dim == 2){ - std::cout << "Analyse 2d structured quadrilateral grid" << std::endl; - // For structured grids we already know the variation beforehand. This // is only implemented to cover all the cases. @@ -107,11 +128,10 @@ int main(int argc, char** argv){ // Extract facemod/facedir intersection variation using ES = Dune::PDELab::AllEntitySet<GV>; ES es(gv); + std::cout << "Analyse 2d structured quadrilateral grid" << std::endl; analyze_grid(es, filename); } if (dim == 3){ - std::cout << "Analyse 3d structured quadrilateral grid" << std::endl; - // For structured grids we already know the variation beforehand. This // is only implemented to cover all the cases. @@ -125,6 +145,7 @@ int main(int argc, char** argv){ // Extract facemod/facedir intersection variation using ES = Dune::PDELab::AllEntitySet<GV>; ES es(gv); + std::cout << "Analyse 3d structured quadrilateral grid" << std::endl; analyze_grid(es, filename); } else{ diff --git a/bin/analyzegrid/consistent_gmsh.cc b/bin/analyzegrid/consistent_gmsh.cc new file mode 100644 index 0000000000000000000000000000000000000000..a5e1d5bbd26390387fa7266b9f4f944051bf4f31 --- /dev/null +++ b/bin/analyzegrid/consistent_gmsh.cc @@ -0,0 +1,49 @@ +#include "config.h" + +#include <dune/common/parallel/mpihelper.hh> +#include <dune/common/parametertree.hh> +#include "dune/common/parametertreeparser.hh" +#include <dune/grid/uggrid.hh> +#include <dune/testtools/gridconstruction.hh> + +#include <dune/consistent-edge-orientation/createconsistentgrid.hh> + +int main(int argc, char** argv){ + try + { + if (argc != 3){ + std::cout << "Need ini file and output filename." << std::endl; + return 1; + } + + // Initialize basic stuff... + Dune::MPIHelper& mpihelper = Dune::MPIHelper::instance(argc, argv); + using RangeType = double; + Dune::ParameterTree initree; + Dune::ParameterTreeParser::readINITree(argv[1], initree); + + // TODO: For now just assume dim = 3 + const int dim = 3; + + // Setup grid (view)... + using Grid = Dune::UGGrid<3>; + using GV = Grid::LeafGridView; + IniGridFactory<Grid> factory(initree); + std::shared_ptr<Grid> grid = factory.getGrid(); + GV gv = grid->leafGridView(); + + std::string filename = std::string(argv[2]); + createConsistentGmshFile(gv, filename); + + return 0; + + } + catch (Dune::Exception& e) + { std::cerr << "Dune reported error: " << e << std::endl; + return 1; + } + catch (std::exception& e) + { std::cerr << "Unknown exception thrown!" << std::endl; + return 1; + } +} diff --git a/python/dune/perftool/pdelab/driver/gridfunctionspace.py b/python/dune/perftool/pdelab/driver/gridfunctionspace.py index a4ae3604c06af34a774ba0494751085fcffc1a9e..dce737323c4fd114149e4d3e1ebb4d52a83af174 100644 --- a/python/dune/perftool/pdelab/driver/gridfunctionspace.py +++ b/python/dune/perftool/pdelab/driver/gridfunctionspace.py @@ -83,10 +83,16 @@ def define_grid(name): # TODO: In principle this is only necessary if we use sum factorization in # one of the operators. So this could be turned off if that is not the case. if isQuadrilateral(get_trial_element().cell()) and get_option("grid_unstructured"): - include_file("dune/consistent-edge-orientation/createconsistentgrid.hh", filetag="driver") + # TODO: createConsistentGrid drops boundary faces. Instead of using + # this we assume that the user already provides a constistent gmsh + # file. This can be createt by createConsistentGmshFile, see + # dune-perftool/bin/analyzegrid/consistent_gmsh.cc + # include_file("dune/consistent-edge-orientation/createconsistentgrid.hh", filetag="driver") + # return ["IniGridFactory<{}> factory({});".format(_type, ini), + # "std::shared_ptr<{}> grid_nonconsistent = factory.getGrid();".format(_type), + # "std::shared_ptr<{}> grid = createConsistentGrid(grid_nonconsistent);".format(_type)] return ["IniGridFactory<{}> factory({});".format(_type, ini), - "std::shared_ptr<{}> grid_nonconsistent = factory.getGrid();".format(_type), - "std::shared_ptr<{}> grid = createConsistentGrid(grid_nonconsistent);".format(_type)] + "std::shared_ptr<{}> grid = factory.getGrid();".format(_type)] return ["IniGridFactory<{}> factory({});".format(_type, ini), "std::shared_ptr<{}> grid = factory.getGrid();".format(_type)] diff --git a/python/dune/perftool/sumfact/geometry.py b/python/dune/perftool/sumfact/geometry.py index c4471ee3ab329eef051bfd1a27f0eb0fce1bd4c3..ee0fbaae2866bae1017d8ef4ed25608332534139 100644 --- a/python/dune/perftool/sumfact/geometry.py +++ b/python/dune/perftool/sumfact/geometry.py @@ -152,7 +152,6 @@ def pymbolic_spatial_coordinate_multilinear(do_predicates, visitor): # Add a sum factorization kernel that implements the evaluation of # the basis functions at quadrature points (stage 1) - from dune.perftool.sumfact.realization import realize_sum_factorization_kernel var, _ = realize_sum_factorization_kernel(vsf) # The results of this function is already the right component of the @@ -260,12 +259,12 @@ def define_outer_normal(name, visitor): temporary_variable(name, shape=(world_dimension(),)) for i in range(world_dimension()): assignee = prim.Subscript(prim.Variable(name), i) - jit = pymbolic_jacobian_inverse(i, facedir_s, Restriction.POSITIVE, visitor) + ji = pymbolic_jacobian_inverse(facedir_s, i, Restriction.POSITIVE, visitor) # Note: 2*facemod_s-1 because of # -1 if facemod_s = 0 # +1 if facemod_s = 1 - expression = jit * (2 * facemod_s - 1) + expression = ji * (2 * facemod_s - 1) inames = default_quadrature_inames(visitor) + additional_inames(visitor) instruction(assignee=assignee, @@ -532,7 +531,7 @@ def _name_jacobian(i, j, restriction, visitor): def define_jacobian_inverse(name, restriction, visitor): - """Return jacobian inverse transposed of the geometry mapping (of a cell) + """Return jacobian inverse of the geometry mapping (of a cell) At the moment this only works for geometry mappings of cells and not for intersection. We only consider this case as it greatly simplifies code @@ -556,10 +555,12 @@ def define_jacobian_inverse(name, restriction, visitor): names_jacobian.append(name_jacobian) # Calculate the inverse of the jacobian of the geometry mapping and the - # determinant by calling a c++ function + # determinant by calling a c++ function. Note: The result will be column + # major -> fortran style. name_detjac = name_jacobian_determinant(visitor) temporary_variable(name_detjac, shape=()) - temporary_variable(name, shape=(dim, dim), managed=True) + ftags = ",".join(["f"] * 2) + temporary_variable(name, shape=(dim, dim), dim_tags=ftags, managed=True) include_file('dune/perftool/sumfact/invertgeometry.hh', filetag='operatorfile') code = "{} = invert_and_return_determinant({}, {});".format(name_detjac, ", ".join(names_jacobian), diff --git a/python/dune/perftool/sumfact/permutation.py b/python/dune/perftool/sumfact/permutation.py index b00be306a7f6c724b35107c30692c2f0d1b4f14a..1fed8e3cd9c722868574f486bc9aacceca1d2b0a 100644 --- a/python/dune/perftool/sumfact/permutation.py +++ b/python/dune/perftool/sumfact/permutation.py @@ -2,6 +2,8 @@ import itertools +from dune.perftool.options import get_option +from dune.perftool.sumfact.switch import get_facedir, get_facemod def sumfact_permutation_heuristic(permutations, stage): """Heuristic to choose a permutation @@ -74,6 +76,13 @@ def sumfact_permutation_strategy(sf): def permute_forward(t, perm): + """Forward permute t according to perm (not inplace) + + Example: + t = ('a', 'b', 'c') + perm = (1, 2, 0) + -> ('b', 'c', 'a') + """ tmp = [] for pos in perm: tmp.append(t[pos]) @@ -81,7 +90,55 @@ def permute_forward(t, perm): def permute_backward(t, perm): + """Backward permute t according to perm (not inplace) + + Inverse of permute_forward(t, perm). + + Example: + t = ('a', 'b', 'c') + perm = (1, 2, 0) + -> ('c', 'a', 'b') + """ tmp = [None] * len(t) for i, pos in enumerate(perm): tmp[pos] = t[i] return tuple(tmp) + + +def sumfact_quadrature_permutation_strategy(dim, restriction): + """Return order of direction for the quadrature points for stage 2 + + On intersection we need to make sure to go through the quadrature points of + self and neighbor in the same order. We do this by making the grid edge + consistent and complying to a convention of the directions depending on the + dimension, restriction, facedir and facemod. + + In order to derive those conventions you need to draw the cells with edge + orientation and see how to match the intersections. + """ + # Use a simpler convention for structured grids. In this case we can always + # go through the directions in the normal order. The same is true for 2D + # and sum factorization on volumes. + if (not get_option('grid_unstructured')) or dim == 2 or restriction == 0: + return tuple(range(dim)) + else: + def _order_on_self(restriction): + facedir = get_facedir(restriction) + facemod = get_facemod(restriction) + + # Here we specify the convention + if (facedir, facemod) in [(0, 0), (1, 1), (2, 0)]: + return tuple(range(dim)) + else: + l = list(range(dim)) + l.reverse() + return tuple(l) + + # On neighbor we need the reverse order + if restriction == 1: + return _order_on_self(restriction) + else: + assert restriction == 2 + l = list(_order_on_self(restriction)) + l.reverse() + return tuple(l) diff --git a/python/dune/perftool/sumfact/quadrature.py b/python/dune/perftool/sumfact/quadrature.py index 6fe70fa441979c4ae00a5f44acc8136def9f5fb0..398e793d49d5beb98d35e14a859936e15f03332c 100644 --- a/python/dune/perftool/sumfact/quadrature.py +++ b/python/dune/perftool/sumfact/quadrature.py @@ -177,6 +177,17 @@ def quadrature_weight(visitor): def define_quadrature_position(name, local_index): + # TODO: This whole function looks quite suspicious: + # - qps_per_dir not used + # - probably rename index to local_index + # + # This whole function is probably not used anymore since we use sum + # factorization for evaluation of global coordinate, jacobian inverse + # transposed and determinant of the jacobian of the geometry mapping. + # + # Needs to be fixed. Or removed. Added to my to do list ;) + assert False + qps_per_dir = quadrature_points_per_direction() local_qps_per_dir = local_quadrature_points_per_direction() qp_bound = local_qps_per_dir[index] diff --git a/python/dune/perftool/sumfact/realization.py b/python/dune/perftool/sumfact/realization.py index 794e67e1dc1a7947f82399554e9678114c017101..0acb21289b42a962b18d3989d11c90ca773252ff 100644 --- a/python/dune/perftool/sumfact/realization.py +++ b/python/dune/perftool/sumfact/realization.py @@ -190,7 +190,8 @@ def realize_sumfact_kernel_function(sf): # first matrix multiplication this can be taken from # * an input temporary (default) # * a global data structure (if FastDGGridOperator is in use) - # * a value from a global data structure, broadcasted to a vector type (vectorized + FastDGGridOperator) + # * a value from a global data structure, broadcasted to a vector type + # (vectorized + FastDGGridOperator) input_inames = (k_expr,) + tuple(prim.Variable(j) for j in out_inames[1:]) if l == 0 and sf.stage == 1 and sf.interface.direct_is_possible: # See comment below @@ -206,6 +207,9 @@ def realize_sumfact_kernel_function(sf): if l == 0: inp_shape = permute_backward(inp_shape, perm) input_inames = permute_backward(input_inames, perm) + if sf.stage == 1: + inp_shape = permute_backward(inp_shape, sf.quadrature_permutation) + input_inames = permute_backward(input_inames, sf.quadrature_permutation) # Get a temporary that interprets the base storage of the input # as a column-major matrix. In later iteration of the matrix loop @@ -215,7 +219,8 @@ def realize_sumfact_kernel_function(sf): dim_tags=ftags, ) - # The input temporary will only be read from, so we need to silence the loopy warning + # The input temporary will only be read from, so we need to silence + # the loopy warning silenced_warning('read_no_write({})'.format(inp)) input_summand = prim.Subscript(prim.Variable(inp), @@ -235,6 +240,9 @@ def realize_sumfact_kernel_function(sf): output_shape = tuple(out_shape[1:]) + (out_shape[0],) if l == len(matrix_sequence) - 1: output_shape = permute_backward(output_shape, perm) + if sf.stage == 3: + output_shape = permute_backward(output_shape, sf.quadrature_permutation) + out = buffer.get_temporary("buff_step{}_out".format(l), shape=output_shape + vec_shape, dim_tags=ftags, @@ -253,6 +261,8 @@ def realize_sumfact_kernel_function(sf): output_inames = tuple(prim.Variable(i) for i in out_inames[1:]) + (prim.Variable(out_inames[0]),) if l == len(matrix_sequence) - 1: output_inames = permute_backward(output_inames, perm) + if sf.stage == 3: + output_inames = permute_backward(output_inames, sf.quadrature_permutation) # Collect the key word arguments for the loopy instruction insn_args = {"depends_on": insn_dep} diff --git a/python/dune/perftool/sumfact/symbolic.py b/python/dune/perftool/sumfact/symbolic.py index fb283a0536318c1585a0963a33ca76105571cb84..c21176933d79139b4d7fa06aef15b77d5634370c 100644 --- a/python/dune/perftool/sumfact/symbolic.py +++ b/python/dune/perftool/sumfact/symbolic.py @@ -6,6 +6,7 @@ from dune.perftool.generation import (get_counted_variable, transform, ) from dune.perftool.pdelab.geometry import local_dimension, world_dimension +from dune.perftool.sumfact.permutation import permute_forward, sumfact_quadrature_permutation_strategy from dune.perftool.sumfact.quadrature import quadrature_inames from dune.perftool.sumfact.tabulation import BasisTabulationMatrixBase, BasisTabulationMatrixArray from dune.perftool.loopy.target import dtype_floatingpoint, type_floatingpoint @@ -224,6 +225,7 @@ class SumfactKernel(SumfactKernelBase, ImmutableRecord, prim.Variable): insn_dep=frozenset(), interface=SumfactKernelInterfaceBase(), predicates=frozenset(), + quadrature_permutation=None, ): """Create a sum factorization kernel @@ -296,6 +298,36 @@ class SumfactKernel(SumfactKernelBase, ImmutableRecord, prim.Variable): for a in SumfactKernel.init_arg_names: defaultdict[a] = eval(a) + # Not sure if this whole permuting would make sense if we would do sum + # factorized evaluation of intersections where len(matrix_sequence) + # would not be equal to world dim. + dim = len(matrix_sequence) + + # Get restriction for this sum factorization kernel. Note: For + # accumulation output we have a restriction for the test (index 0) and + # ansatz (index 1) space. We need the restriction corresponding to the + # test space since we are in stage 3 + restriction = interface.restriction + if isinstance(restriction, tuple): + assert interface.stage is 3 + assert len(restriction) is 2 + restriction = restriction[0] + + perm = sumfact_quadrature_permutation_strategy(dim, restriction) + permuted_matrix_sequence = [] + for i in perm: + for mat in matrix_sequence: + if mat.direction == i: + permuted_matrix_sequence.append(mat) + permuted_matrix_sequence = tuple(permuted_matrix_sequence) + + from dune.perftool.sumfact.switch import get_facedir, get_facemod + facedir = get_facedir(restriction) + facemod = get_facemod(restriction) + + defaultdict['matrix_sequence'] = permuted_matrix_sequence + defaultdict['quadrature_permutation'] = perm + # Call the base class constructors ImmutableRecord.__init__(self, **defaultdict) prim.Variable.__init__(self, "SUMFACT") @@ -325,8 +357,16 @@ class SumfactKernel(SumfactKernelBase, ImmutableRecord, prim.Variable): @property def function_name(self): """ The name of the function that implements this kernel """ - return "sfimpl_{}{}".format("_".join(str(m) for m in self.matrix_sequence), + name = "sfimpl_{}{}".format("_".join(str(m) for m in self.matrix_sequence), self.interface.function_name_suffix) + # On unstructured we need different permutation of the input to realize + # different permuation of quadrature points on self and neighbor. Mange + # the permutation of the quadrature points into the name to generate + # sperate functions. + if self.quadrature_permutation != tuple(range(len(self.matrix_sequence))): + name_quad_perm = "_qpperm_{}".format("_".join(str(a) for a in self.quadrature_permutation)) + name = name + name_quad_perm + return name @property def parallel_key(self): diff --git a/test/sumfact/poisson/CMakeLists.txt b/test/sumfact/poisson/CMakeLists.txt index a304dc7fcb6649ff6440688fcd0055d3645c65f4..ee95c2e25d3c1440cd88d7214ba4d80652415ab2 100644 --- a/test/sumfact/poisson/CMakeLists.txt +++ b/test/sumfact/poisson/CMakeLists.txt @@ -61,6 +61,13 @@ if(consistent-edge-orientation_FOUND) INIFILE poisson_dg_3d_unstructured.mini ANALYZE_GRID ) + dune_add_formcompiler_system_test(UFLFILE poisson_dg_3d.ufl + BASENAME sumfact_poisson_dg_3d_gmsh + INIFILE poisson_dg_3d_gmsh.mini + ANALYZE_GRID + ) + dune_symlink_to_source_files(FILES cube_hexa_2.msh cube_hexa_2_consistent.msh) + dune_symlink_to_source_files(FILES cube_hexa.msh cube_hexa_consistent.msh) endif() #============================================= diff --git a/test/sumfact/poisson/cube_hexa.msh b/test/sumfact/poisson/cube_hexa.msh new file mode 100644 index 0000000000000000000000000000000000000000..15101b483ece71544a5f223606e9d8f912ea5fa1 --- /dev/null +++ b/test/sumfact/poisson/cube_hexa.msh @@ -0,0 +1,1738 @@ +$MeshFormat +2.2 0 8 +$EndMeshFormat +$Nodes +729 +1 0 0 0 +2 1 0 0 +3 1 1 0 +4 0 1 0 +5 0 0 1 +6 1 0 1 +7 1 1 1 +8 0 1 1 +9 0.1249999999997731 0 0 +10 0.2499999999994109 0 0 +11 0.3749999999990461 0 0 +12 0.4999999999986921 0 0 +13 0.6249999999990109 0 0 +14 0.7499999999993406 0 0 +15 0.8749999999996703 0 0 +16 1 0.1249999999997731 0 +17 1 0.2499999999994109 0 +18 1 0.3749999999990461 0 +19 1 0.4999999999986921 0 +20 1 0.6249999999990109 0 +21 1 0.7499999999993406 0 +22 1 0.8749999999996703 0 +23 0.8749999999995012 1 0 +24 0.7500000000003471 1 0 +25 0.625000000001214 1 0 +26 0.5000000000020595 1 0 +27 0.3750000000015609 1 0 +28 0.2500000000010405 1 0 +29 0.1250000000005203 1 0 +30 0 0.8749999999995012 0 +31 0 0.7500000000003471 0 +32 0 0.625000000001214 0 +33 0 0.5000000000020595 0 +34 0 0.3750000000015609 0 +35 0 0.2500000000010405 0 +36 0 0.1250000000005203 0 +37 0.1249999999997731 0 1 +38 0.2499999999994109 0 1 +39 0.3749999999990461 0 1 +40 0.4999999999986921 0 1 +41 0.6249999999990109 0 1 +42 0.7499999999993406 0 1 +43 0.8749999999996703 0 1 +44 1 0.1249999999997731 1 +45 1 0.2499999999994109 1 +46 1 0.3749999999990461 1 +47 1 0.4999999999986921 1 +48 1 0.6249999999990109 1 +49 1 0.7499999999993406 1 +50 1 0.8749999999996703 1 +51 0.8749999999995012 1 1 +52 0.7500000000003471 1 1 +53 0.625000000001214 1 1 +54 0.5000000000020595 1 1 +55 0.3750000000015609 1 1 +56 0.2500000000010405 1 1 +57 0.1250000000005203 1 1 +58 0 0.8749999999995012 1 +59 0 0.7500000000003471 1 +60 0 0.625000000001214 1 +61 0 0.5000000000020595 1 +62 0 0.3750000000015609 1 +63 0 0.2500000000010405 1 +64 0 0.1250000000005203 1 +65 0 0 0.125 +66 0 0 0.25 +67 0 0 0.375 +68 0 0 0.5 +69 0 0 0.625 +70 0 0 0.75 +71 0 0 0.875 +72 1 0 0.125 +73 1 0 0.25 +74 1 0 0.375 +75 1 0 0.5 +76 1 0 0.625 +77 1 0 0.75 +78 1 0 0.875 +79 1 1 0.125 +80 1 1 0.25 +81 1 1 0.375 +82 1 1 0.5 +83 1 1 0.625 +84 1 1 0.75 +85 1 1 0.875 +86 0 1 0.125 +87 0 1 0.25 +88 0 1 0.375 +89 0 1 0.5 +90 0 1 0.625 +91 0 1 0.75 +92 0 1 0.875 +93 0.1249999999998665 0.1250000000004269 0 +94 0.1249999999999599 0.2500000000008368 0 +95 0.1250000000000533 0.3750000000012466 0 +96 0.1250000000001467 0.5000000000016386 0 +97 0.1250000000002401 0.6250000000009386 0 +98 0.1250000000003335 0.7500000000002213 0 +99 0.1250000000004269 0.8749999999995224 0 +100 0.2499999999996146 0.1250000000003335 0 +101 0.2499999999998183 0.2500000000006331 0 +102 0.250000000000022 0.3750000000009321 0 +103 0.2500000000002257 0.5000000000012177 0 +104 0.2500000000004294 0.6250000000006631 0 +105 0.2500000000006332 0.7500000000000954 0 +106 0.2500000000008369 0.8749999999995435 0 +107 0.3749999999993604 0.1250000000002401 0 +108 0.3749999999996748 0.2500000000004294 0 +109 0.3749999999999891 0.3750000000006178 0 +110 0.3750000000003035 0.5000000000007968 0 +111 0.3750000000006178 0.6250000000003878 0 +112 0.3750000000009321 0.7499999999999697 0 +113 0.3750000000012466 0.8749999999995646 0 +114 0.499999999999113 0.1250000000001467 0 +115 0.499999999999534 0.2500000000002257 0 +116 0.4999999999999548 0.3750000000003034 0 +117 0.5000000000003758 0.5000000000003758 0 +118 0.5000000000007967 0.6250000000001126 0 +119 0.5000000000012177 0.7499999999998439 0 +120 0.5000000000016385 0.8749999999995857 0 +121 0.6249999999992863 0.1250000000000533 0 +122 0.6249999999995618 0.250000000000022 0 +123 0.6249999999998369 0.3749999999999891 0 +124 0.6250000000001126 0.4999999999999548 0 +125 0.6250000000003878 0.6249999999998371 0 +126 0.6250000000006634 0.7499999999997179 0 +127 0.6250000000009387 0.8749999999996068 0 +128 0.7499999999994663 0.1249999999999599 0 +129 0.7499999999995922 0.2499999999998183 0 +130 0.7499999999997179 0.3749999999996748 0 +131 0.7499999999998439 0.499999999999534 0 +132 0.7499999999999697 0.6249999999995618 0 +133 0.7500000000000953 0.7499999999995921 0 +134 0.7500000000002213 0.8749999999996281 0 +135 0.8749999999996492 0.1249999999998666 0 +136 0.8749999999996279 0.2499999999996146 0 +137 0.8749999999996068 0.3749999999993605 0 +138 0.8749999999995857 0.499999999999113 0 +139 0.8749999999995646 0.6249999999992862 0 +140 0.8749999999995435 0.7499999999994665 0 +141 0.8749999999995224 0.8749999999996492 0 +142 0.1249999999997731 0 0.125 +143 0.1249999999997731 0 0.25 +144 0.1249999999997731 0 0.375 +145 0.1249999999997731 0 0.5 +146 0.1249999999997731 0 0.625 +147 0.1249999999997731 0 0.75 +148 0.1249999999997731 0 0.875 +149 0.2499999999994109 0 0.125 +150 0.2499999999994109 0 0.25 +151 0.2499999999994109 0 0.375 +152 0.2499999999994109 0 0.5 +153 0.2499999999994109 0 0.625 +154 0.2499999999994109 0 0.75 +155 0.2499999999994109 0 0.875 +156 0.3749999999990461 0 0.125 +157 0.3749999999990461 0 0.25 +158 0.3749999999990461 0 0.375 +159 0.3749999999990461 0 0.5 +160 0.3749999999990461 0 0.625 +161 0.3749999999990461 0 0.75 +162 0.3749999999990461 0 0.875 +163 0.4999999999986921 0 0.125 +164 0.4999999999986921 0 0.25 +165 0.4999999999986921 0 0.375 +166 0.4999999999986921 0 0.5 +167 0.4999999999986921 0 0.625 +168 0.4999999999986921 0 0.75 +169 0.4999999999986921 0 0.875 +170 0.6249999999990109 0 0.125 +171 0.6249999999990109 0 0.25 +172 0.6249999999990109 0 0.375 +173 0.6249999999990109 0 0.5 +174 0.6249999999990109 0 0.625 +175 0.6249999999990109 0 0.75 +176 0.6249999999990109 0 0.875 +177 0.7499999999993406 0 0.125 +178 0.7499999999993406 0 0.25 +179 0.7499999999993406 0 0.375 +180 0.7499999999993406 0 0.5 +181 0.7499999999993406 0 0.625 +182 0.7499999999993406 0 0.75 +183 0.7499999999993406 0 0.875 +184 0.8749999999996703 0 0.125 +185 0.8749999999996703 0 0.25 +186 0.8749999999996703 0 0.375 +187 0.8749999999996703 0 0.5 +188 0.8749999999996703 0 0.625 +189 0.8749999999996703 0 0.75 +190 0.8749999999996703 0 0.875 +191 1 0.1249999999997731 0.125 +192 1 0.1249999999997731 0.25 +193 1 0.1249999999997731 0.375 +194 1 0.1249999999997731 0.5 +195 1 0.1249999999997731 0.625 +196 1 0.1249999999997731 0.75 +197 1 0.1249999999997731 0.875 +198 1 0.2499999999994109 0.125 +199 1 0.2499999999994109 0.25 +200 1 0.2499999999994109 0.375 +201 1 0.2499999999994109 0.5 +202 1 0.2499999999994109 0.625 +203 1 0.2499999999994109 0.75 +204 1 0.2499999999994109 0.875 +205 1 0.3749999999990461 0.125 +206 1 0.3749999999990461 0.25 +207 1 0.3749999999990461 0.375 +208 1 0.3749999999990461 0.5 +209 1 0.3749999999990461 0.625 +210 1 0.3749999999990461 0.75 +211 1 0.3749999999990461 0.875 +212 1 0.4999999999986921 0.125 +213 1 0.4999999999986921 0.25 +214 1 0.4999999999986921 0.375 +215 1 0.4999999999986921 0.5 +216 1 0.4999999999986921 0.625 +217 1 0.4999999999986921 0.75 +218 1 0.4999999999986921 0.875 +219 1 0.6249999999990109 0.125 +220 1 0.6249999999990109 0.25 +221 1 0.6249999999990109 0.375 +222 1 0.6249999999990109 0.5 +223 1 0.6249999999990109 0.625 +224 1 0.6249999999990109 0.75 +225 1 0.6249999999990109 0.875 +226 1 0.7499999999993406 0.125 +227 1 0.7499999999993406 0.25 +228 1 0.7499999999993406 0.375 +229 1 0.7499999999993406 0.5 +230 1 0.7499999999993406 0.625 +231 1 0.7499999999993406 0.75 +232 1 0.7499999999993406 0.875 +233 1 0.8749999999996703 0.125 +234 1 0.8749999999996703 0.25 +235 1 0.8749999999996703 0.375 +236 1 0.8749999999996703 0.5 +237 1 0.8749999999996703 0.625 +238 1 0.8749999999996703 0.75 +239 1 0.8749999999996703 0.875 +240 0.8749999999995012 1 0.125 +241 0.8749999999995012 1 0.25 +242 0.8749999999995012 1 0.375 +243 0.8749999999995012 1 0.5 +244 0.8749999999995012 1 0.625 +245 0.8749999999995012 1 0.75 +246 0.8749999999995012 1 0.875 +247 0.7500000000003471 1 0.125 +248 0.7500000000003471 1 0.25 +249 0.7500000000003471 1 0.375 +250 0.7500000000003471 1 0.5 +251 0.7500000000003471 1 0.625 +252 0.7500000000003471 1 0.75 +253 0.7500000000003471 1 0.875 +254 0.625000000001214 1 0.125 +255 0.625000000001214 1 0.25 +256 0.625000000001214 1 0.375 +257 0.625000000001214 1 0.5 +258 0.625000000001214 1 0.625 +259 0.625000000001214 1 0.75 +260 0.625000000001214 1 0.875 +261 0.5000000000020595 1 0.125 +262 0.5000000000020595 1 0.25 +263 0.5000000000020595 1 0.375 +264 0.5000000000020595 1 0.5 +265 0.5000000000020595 1 0.625 +266 0.5000000000020595 1 0.75 +267 0.5000000000020595 1 0.875 +268 0.3750000000015609 1 0.125 +269 0.3750000000015609 1 0.25 +270 0.3750000000015609 1 0.375 +271 0.3750000000015609 1 0.5 +272 0.3750000000015609 1 0.625 +273 0.3750000000015609 1 0.75 +274 0.3750000000015609 1 0.875 +275 0.2500000000010405 1 0.125 +276 0.2500000000010405 1 0.25 +277 0.2500000000010405 1 0.375 +278 0.2500000000010405 1 0.5 +279 0.2500000000010405 1 0.625 +280 0.2500000000010405 1 0.75 +281 0.2500000000010405 1 0.875 +282 0.1250000000005203 1 0.125 +283 0.1250000000005203 1 0.25 +284 0.1250000000005203 1 0.375 +285 0.1250000000005203 1 0.5 +286 0.1250000000005203 1 0.625 +287 0.1250000000005203 1 0.75 +288 0.1250000000005203 1 0.875 +289 0 0.8749999999995012 0.125 +290 0 0.8749999999995012 0.25 +291 0 0.8749999999995012 0.375 +292 0 0.8749999999995012 0.5 +293 0 0.8749999999995012 0.625 +294 0 0.8749999999995012 0.75 +295 0 0.8749999999995012 0.875 +296 0 0.7500000000003471 0.125 +297 0 0.7500000000003471 0.25 +298 0 0.7500000000003471 0.375 +299 0 0.7500000000003471 0.5 +300 0 0.7500000000003471 0.625 +301 0 0.7500000000003471 0.75 +302 0 0.7500000000003471 0.875 +303 0 0.625000000001214 0.125 +304 0 0.625000000001214 0.25 +305 0 0.625000000001214 0.375 +306 0 0.625000000001214 0.5 +307 0 0.625000000001214 0.625 +308 0 0.625000000001214 0.75 +309 0 0.625000000001214 0.875 +310 0 0.5000000000020595 0.125 +311 0 0.5000000000020595 0.25 +312 0 0.5000000000020595 0.375 +313 0 0.5000000000020595 0.5 +314 0 0.5000000000020595 0.625 +315 0 0.5000000000020595 0.75 +316 0 0.5000000000020595 0.875 +317 0 0.3750000000015609 0.125 +318 0 0.3750000000015609 0.25 +319 0 0.3750000000015609 0.375 +320 0 0.3750000000015609 0.5 +321 0 0.3750000000015609 0.625 +322 0 0.3750000000015609 0.75 +323 0 0.3750000000015609 0.875 +324 0 0.2500000000010405 0.125 +325 0 0.2500000000010405 0.25 +326 0 0.2500000000010405 0.375 +327 0 0.2500000000010405 0.5 +328 0 0.2500000000010405 0.625 +329 0 0.2500000000010405 0.75 +330 0 0.2500000000010405 0.875 +331 0 0.1250000000005203 0.125 +332 0 0.1250000000005203 0.25 +333 0 0.1250000000005203 0.375 +334 0 0.1250000000005203 0.5 +335 0 0.1250000000005203 0.625 +336 0 0.1250000000005203 0.75 +337 0 0.1250000000005203 0.875 +338 0.1249999999998665 0.1250000000004269 1 +339 0.1249999999999599 0.2500000000008368 1 +340 0.1250000000000533 0.3750000000012466 1 +341 0.1250000000001467 0.5000000000016386 1 +342 0.1250000000002401 0.6250000000009386 1 +343 0.1250000000003335 0.7500000000002213 1 +344 0.1250000000004269 0.8749999999995224 1 +345 0.2499999999996146 0.1250000000003335 1 +346 0.2499999999998183 0.2500000000006331 1 +347 0.250000000000022 0.3750000000009321 1 +348 0.2500000000002257 0.5000000000012177 1 +349 0.2500000000004294 0.6250000000006631 1 +350 0.2500000000006332 0.7500000000000954 1 +351 0.2500000000008369 0.8749999999995435 1 +352 0.3749999999993604 0.1250000000002401 1 +353 0.3749999999996748 0.2500000000004294 1 +354 0.3749999999999891 0.3750000000006178 1 +355 0.3750000000003035 0.5000000000007968 1 +356 0.3750000000006178 0.6250000000003878 1 +357 0.3750000000009321 0.7499999999999697 1 +358 0.3750000000012466 0.8749999999995646 1 +359 0.499999999999113 0.1250000000001467 1 +360 0.499999999999534 0.2500000000002257 1 +361 0.4999999999999548 0.3750000000003034 1 +362 0.5000000000003758 0.5000000000003758 1 +363 0.5000000000007967 0.6250000000001126 1 +364 0.5000000000012177 0.7499999999998439 1 +365 0.5000000000016385 0.8749999999995857 1 +366 0.6249999999992863 0.1250000000000533 1 +367 0.6249999999995618 0.250000000000022 1 +368 0.6249999999998369 0.3749999999999891 1 +369 0.6250000000001126 0.4999999999999548 1 +370 0.6250000000003878 0.6249999999998371 1 +371 0.6250000000006634 0.7499999999997179 1 +372 0.6250000000009387 0.8749999999996068 1 +373 0.7499999999994663 0.1249999999999599 1 +374 0.7499999999995922 0.2499999999998183 1 +375 0.7499999999997179 0.3749999999996748 1 +376 0.7499999999998439 0.499999999999534 1 +377 0.7499999999999697 0.6249999999995618 1 +378 0.7500000000000953 0.7499999999995921 1 +379 0.7500000000002213 0.8749999999996281 1 +380 0.8749999999996492 0.1249999999998666 1 +381 0.8749999999996279 0.2499999999996146 1 +382 0.8749999999996068 0.3749999999993605 1 +383 0.8749999999995857 0.499999999999113 1 +384 0.8749999999995646 0.6249999999992862 1 +385 0.8749999999995435 0.7499999999994665 1 +386 0.8749999999995224 0.8749999999996492 1 +387 0.1249999999998665 0.1250000000004269 0.125 +388 0.1249999999998665 0.1250000000004269 0.25 +389 0.1249999999998665 0.1250000000004269 0.375 +390 0.1249999999998665 0.1250000000004269 0.5 +391 0.1249999999998665 0.1250000000004269 0.625 +392 0.1249999999998665 0.1250000000004269 0.75 +393 0.1249999999998665 0.1250000000004269 0.875 +394 0.1249999999999599 0.2500000000008368 0.125 +395 0.1249999999999599 0.2500000000008368 0.25 +396 0.1249999999999599 0.2500000000008368 0.375 +397 0.1249999999999599 0.2500000000008368 0.5 +398 0.1249999999999599 0.2500000000008368 0.625 +399 0.1249999999999599 0.2500000000008368 0.75 +400 0.1249999999999599 0.2500000000008368 0.875 +401 0.1250000000000533 0.3750000000012466 0.125 +402 0.1250000000000533 0.3750000000012466 0.25 +403 0.1250000000000533 0.3750000000012466 0.375 +404 0.1250000000000533 0.3750000000012466 0.5 +405 0.1250000000000533 0.3750000000012466 0.625 +406 0.1250000000000533 0.3750000000012466 0.75 +407 0.1250000000000533 0.3750000000012466 0.875 +408 0.1250000000001467 0.5000000000016386 0.125 +409 0.1250000000001467 0.5000000000016386 0.25 +410 0.1250000000001467 0.5000000000016386 0.375 +411 0.1250000000001467 0.5000000000016386 0.5 +412 0.1250000000001467 0.5000000000016386 0.625 +413 0.1250000000001467 0.5000000000016386 0.75 +414 0.1250000000001467 0.5000000000016386 0.875 +415 0.1250000000002401 0.6250000000009386 0.125 +416 0.1250000000002401 0.6250000000009386 0.25 +417 0.1250000000002401 0.6250000000009386 0.375 +418 0.1250000000002401 0.6250000000009386 0.5 +419 0.1250000000002401 0.6250000000009386 0.625 +420 0.1250000000002401 0.6250000000009386 0.75 +421 0.1250000000002401 0.6250000000009386 0.875 +422 0.1250000000003335 0.7500000000002213 0.125 +423 0.1250000000003335 0.7500000000002213 0.25 +424 0.1250000000003335 0.7500000000002213 0.375 +425 0.1250000000003335 0.7500000000002213 0.5 +426 0.1250000000003335 0.7500000000002213 0.625 +427 0.1250000000003335 0.7500000000002213 0.75 +428 0.1250000000003335 0.7500000000002213 0.875 +429 0.1250000000004269 0.8749999999995224 0.125 +430 0.1250000000004269 0.8749999999995224 0.25 +431 0.1250000000004269 0.8749999999995224 0.375 +432 0.1250000000004269 0.8749999999995224 0.5 +433 0.1250000000004269 0.8749999999995224 0.625 +434 0.1250000000004269 0.8749999999995224 0.75 +435 0.1250000000004269 0.8749999999995224 0.875 +436 0.2499999999996146 0.1250000000003335 0.125 +437 0.2499999999996146 0.1250000000003335 0.25 +438 0.2499999999996146 0.1250000000003335 0.375 +439 0.2499999999996146 0.1250000000003335 0.5 +440 0.2499999999996146 0.1250000000003335 0.625 +441 0.2499999999996146 0.1250000000003335 0.75 +442 0.2499999999996146 0.1250000000003335 0.875 +443 0.2499999999998183 0.2500000000006331 0.125 +444 0.2499999999998183 0.2500000000006331 0.25 +445 0.2499999999998183 0.2500000000006331 0.375 +446 0.2499999999998183 0.2500000000006331 0.5 +447 0.2499999999998183 0.2500000000006331 0.625 +448 0.2499999999998183 0.2500000000006331 0.75 +449 0.2499999999998183 0.2500000000006331 0.875 +450 0.250000000000022 0.3750000000009321 0.125 +451 0.250000000000022 0.3750000000009321 0.25 +452 0.250000000000022 0.3750000000009321 0.375 +453 0.250000000000022 0.3750000000009321 0.5 +454 0.250000000000022 0.3750000000009321 0.625 +455 0.250000000000022 0.3750000000009321 0.75 +456 0.250000000000022 0.3750000000009321 0.875 +457 0.2500000000002257 0.5000000000012177 0.125 +458 0.2500000000002257 0.5000000000012177 0.25 +459 0.2500000000002257 0.5000000000012177 0.375 +460 0.2500000000002257 0.5000000000012177 0.5 +461 0.2500000000002257 0.5000000000012177 0.625 +462 0.2500000000002257 0.5000000000012177 0.75 +463 0.2500000000002257 0.5000000000012177 0.875 +464 0.2500000000004294 0.6250000000006631 0.125 +465 0.2500000000004294 0.6250000000006631 0.25 +466 0.2500000000004294 0.6250000000006631 0.375 +467 0.2500000000004294 0.6250000000006631 0.5 +468 0.2500000000004294 0.6250000000006631 0.625 +469 0.2500000000004294 0.6250000000006631 0.75 +470 0.2500000000004294 0.6250000000006631 0.875 +471 0.2500000000006332 0.7500000000000954 0.125 +472 0.2500000000006332 0.7500000000000954 0.25 +473 0.2500000000006332 0.7500000000000954 0.375 +474 0.2500000000006332 0.7500000000000954 0.5 +475 0.2500000000006332 0.7500000000000954 0.625 +476 0.2500000000006332 0.7500000000000954 0.75 +477 0.2500000000006332 0.7500000000000954 0.875 +478 0.2500000000008369 0.8749999999995435 0.125 +479 0.2500000000008369 0.8749999999995435 0.25 +480 0.2500000000008369 0.8749999999995435 0.375 +481 0.2500000000008369 0.8749999999995435 0.5 +482 0.2500000000008369 0.8749999999995435 0.625 +483 0.2500000000008369 0.8749999999995435 0.75 +484 0.2500000000008369 0.8749999999995435 0.875 +485 0.3749999999993604 0.1250000000002401 0.125 +486 0.3749999999993604 0.1250000000002401 0.25 +487 0.3749999999993604 0.1250000000002401 0.375 +488 0.3749999999993604 0.1250000000002401 0.5 +489 0.3749999999993604 0.1250000000002401 0.625 +490 0.3749999999993604 0.1250000000002401 0.75 +491 0.3749999999993604 0.1250000000002401 0.875 +492 0.3749999999996748 0.2500000000004294 0.125 +493 0.3749999999996748 0.2500000000004294 0.25 +494 0.3749999999996748 0.2500000000004294 0.375 +495 0.3749999999996748 0.2500000000004294 0.5 +496 0.3749999999996748 0.2500000000004294 0.625 +497 0.3749999999996748 0.2500000000004294 0.75 +498 0.3749999999996748 0.2500000000004294 0.875 +499 0.3749999999999891 0.3750000000006178 0.125 +500 0.3749999999999891 0.3750000000006178 0.25 +501 0.3749999999999891 0.3750000000006178 0.375 +502 0.3749999999999891 0.3750000000006178 0.5 +503 0.3749999999999891 0.3750000000006178 0.625 +504 0.3749999999999891 0.3750000000006178 0.75 +505 0.3749999999999891 0.3750000000006178 0.875 +506 0.3750000000003035 0.5000000000007968 0.125 +507 0.3750000000003035 0.5000000000007968 0.25 +508 0.3750000000003035 0.5000000000007968 0.375 +509 0.3750000000003035 0.5000000000007968 0.5 +510 0.3750000000003035 0.5000000000007968 0.625 +511 0.3750000000003035 0.5000000000007968 0.75 +512 0.3750000000003035 0.5000000000007968 0.875 +513 0.3750000000006178 0.6250000000003878 0.125 +514 0.3750000000006178 0.6250000000003878 0.25 +515 0.3750000000006178 0.6250000000003878 0.375 +516 0.3750000000006178 0.6250000000003878 0.5 +517 0.3750000000006178 0.6250000000003878 0.625 +518 0.3750000000006178 0.6250000000003878 0.75 +519 0.3750000000006178 0.6250000000003878 0.875 +520 0.3750000000009321 0.7499999999999697 0.125 +521 0.3750000000009321 0.7499999999999697 0.25 +522 0.3750000000009321 0.7499999999999697 0.375 +523 0.3750000000009321 0.7499999999999697 0.5 +524 0.3750000000009321 0.7499999999999697 0.625 +525 0.3750000000009321 0.7499999999999697 0.75 +526 0.3750000000009321 0.7499999999999697 0.875 +527 0.3750000000012466 0.8749999999995646 0.125 +528 0.3750000000012466 0.8749999999995646 0.25 +529 0.3750000000012466 0.8749999999995646 0.375 +530 0.3750000000012466 0.8749999999995646 0.5 +531 0.3750000000012466 0.8749999999995646 0.625 +532 0.3750000000012466 0.8749999999995646 0.75 +533 0.3750000000012466 0.8749999999995646 0.875 +534 0.499999999999113 0.1250000000001467 0.125 +535 0.499999999999113 0.1250000000001467 0.25 +536 0.499999999999113 0.1250000000001467 0.375 +537 0.499999999999113 0.1250000000001467 0.5 +538 0.499999999999113 0.1250000000001467 0.625 +539 0.499999999999113 0.1250000000001467 0.75 +540 0.499999999999113 0.1250000000001467 0.875 +541 0.499999999999534 0.2500000000002257 0.125 +542 0.499999999999534 0.2500000000002257 0.25 +543 0.499999999999534 0.2500000000002257 0.375 +544 0.499999999999534 0.2500000000002257 0.5 +545 0.499999999999534 0.2500000000002257 0.625 +546 0.499999999999534 0.2500000000002257 0.75 +547 0.499999999999534 0.2500000000002257 0.875 +548 0.4999999999999548 0.3750000000003034 0.125 +549 0.4999999999999548 0.3750000000003034 0.25 +550 0.4999999999999548 0.3750000000003034 0.375 +551 0.4999999999999548 0.3750000000003034 0.5 +552 0.4999999999999548 0.3750000000003034 0.625 +553 0.4999999999999548 0.3750000000003034 0.75 +554 0.4999999999999548 0.3750000000003034 0.875 +555 0.5000000000003758 0.5000000000003758 0.125 +556 0.5000000000003758 0.5000000000003758 0.25 +557 0.5000000000003758 0.5000000000003758 0.375 +558 0.5000000000003758 0.5000000000003758 0.5 +559 0.5000000000003758 0.5000000000003758 0.625 +560 0.5000000000003758 0.5000000000003758 0.75 +561 0.5000000000003758 0.5000000000003758 0.875 +562 0.5000000000007967 0.6250000000001126 0.125 +563 0.5000000000007967 0.6250000000001126 0.25 +564 0.5000000000007967 0.6250000000001126 0.375 +565 0.5000000000007967 0.6250000000001126 0.5 +566 0.5000000000007967 0.6250000000001126 0.625 +567 0.5000000000007967 0.6250000000001126 0.75 +568 0.5000000000007967 0.6250000000001126 0.875 +569 0.5000000000012177 0.7499999999998439 0.125 +570 0.5000000000012177 0.7499999999998439 0.25 +571 0.5000000000012177 0.7499999999998439 0.375 +572 0.5000000000012177 0.7499999999998439 0.5 +573 0.5000000000012177 0.7499999999998439 0.625 +574 0.5000000000012177 0.7499999999998439 0.75 +575 0.5000000000012177 0.7499999999998439 0.875 +576 0.5000000000016385 0.8749999999995857 0.125 +577 0.5000000000016385 0.8749999999995857 0.25 +578 0.5000000000016385 0.8749999999995857 0.375 +579 0.5000000000016385 0.8749999999995857 0.5 +580 0.5000000000016385 0.8749999999995857 0.625 +581 0.5000000000016385 0.8749999999995857 0.75 +582 0.5000000000016385 0.8749999999995857 0.875 +583 0.6249999999992863 0.1250000000000533 0.125 +584 0.6249999999992863 0.1250000000000533 0.25 +585 0.6249999999992863 0.1250000000000533 0.375 +586 0.6249999999992863 0.1250000000000533 0.5 +587 0.6249999999992863 0.1250000000000533 0.625 +588 0.6249999999992863 0.1250000000000533 0.75 +589 0.6249999999992863 0.1250000000000533 0.875 +590 0.6249999999995618 0.250000000000022 0.125 +591 0.6249999999995618 0.250000000000022 0.25 +592 0.6249999999995618 0.250000000000022 0.375 +593 0.6249999999995618 0.250000000000022 0.5 +594 0.6249999999995618 0.250000000000022 0.625 +595 0.6249999999995618 0.250000000000022 0.75 +596 0.6249999999995618 0.250000000000022 0.875 +597 0.6249999999998369 0.3749999999999891 0.125 +598 0.6249999999998369 0.3749999999999891 0.25 +599 0.6249999999998369 0.3749999999999891 0.375 +600 0.6249999999998369 0.3749999999999891 0.5 +601 0.6249999999998369 0.3749999999999891 0.625 +602 0.6249999999998369 0.3749999999999891 0.75 +603 0.6249999999998369 0.3749999999999891 0.875 +604 0.6250000000001126 0.4999999999999548 0.125 +605 0.6250000000001126 0.4999999999999548 0.25 +606 0.6250000000001126 0.4999999999999548 0.375 +607 0.6250000000001126 0.4999999999999548 0.5 +608 0.6250000000001126 0.4999999999999548 0.625 +609 0.6250000000001126 0.4999999999999548 0.75 +610 0.6250000000001126 0.4999999999999548 0.875 +611 0.6250000000003878 0.6249999999998371 0.125 +612 0.6250000000003878 0.6249999999998371 0.25 +613 0.6250000000003878 0.6249999999998371 0.375 +614 0.6250000000003878 0.6249999999998371 0.5 +615 0.6250000000003878 0.6249999999998371 0.625 +616 0.6250000000003878 0.6249999999998371 0.75 +617 0.6250000000003878 0.6249999999998371 0.875 +618 0.6250000000006634 0.7499999999997179 0.125 +619 0.6250000000006634 0.7499999999997179 0.25 +620 0.6250000000006634 0.7499999999997179 0.375 +621 0.6250000000006634 0.7499999999997179 0.5 +622 0.6250000000006634 0.7499999999997179 0.625 +623 0.6250000000006634 0.7499999999997179 0.75 +624 0.6250000000006634 0.7499999999997179 0.875 +625 0.6250000000009387 0.8749999999996068 0.125 +626 0.6250000000009387 0.8749999999996068 0.25 +627 0.6250000000009387 0.8749999999996068 0.375 +628 0.6250000000009387 0.8749999999996068 0.5 +629 0.6250000000009387 0.8749999999996068 0.625 +630 0.6250000000009387 0.8749999999996068 0.75 +631 0.6250000000009387 0.8749999999996068 0.875 +632 0.7499999999994663 0.1249999999999599 0.125 +633 0.7499999999994663 0.1249999999999599 0.25 +634 0.7499999999994663 0.1249999999999599 0.375 +635 0.7499999999994663 0.1249999999999599 0.5 +636 0.7499999999994663 0.1249999999999599 0.625 +637 0.7499999999994663 0.1249999999999599 0.75 +638 0.7499999999994663 0.1249999999999599 0.875 +639 0.7499999999995922 0.2499999999998183 0.125 +640 0.7499999999995922 0.2499999999998183 0.25 +641 0.7499999999995922 0.2499999999998183 0.375 +642 0.7499999999995922 0.2499999999998183 0.5 +643 0.7499999999995922 0.2499999999998183 0.625 +644 0.7499999999995922 0.2499999999998183 0.75 +645 0.7499999999995922 0.2499999999998183 0.875 +646 0.7499999999997179 0.3749999999996748 0.125 +647 0.7499999999997179 0.3749999999996748 0.25 +648 0.7499999999997179 0.3749999999996748 0.375 +649 0.7499999999997179 0.3749999999996748 0.5 +650 0.7499999999997179 0.3749999999996748 0.625 +651 0.7499999999997179 0.3749999999996748 0.75 +652 0.7499999999997179 0.3749999999996748 0.875 +653 0.7499999999998439 0.499999999999534 0.125 +654 0.7499999999998439 0.499999999999534 0.25 +655 0.7499999999998439 0.499999999999534 0.375 +656 0.7499999999998439 0.499999999999534 0.5 +657 0.7499999999998439 0.499999999999534 0.625 +658 0.7499999999998439 0.499999999999534 0.75 +659 0.7499999999998439 0.499999999999534 0.875 +660 0.7499999999999697 0.6249999999995618 0.125 +661 0.7499999999999697 0.6249999999995618 0.25 +662 0.7499999999999697 0.6249999999995618 0.375 +663 0.7499999999999697 0.6249999999995618 0.5 +664 0.7499999999999697 0.6249999999995618 0.625 +665 0.7499999999999697 0.6249999999995618 0.75 +666 0.7499999999999697 0.6249999999995618 0.875 +667 0.7500000000000953 0.7499999999995921 0.125 +668 0.7500000000000953 0.7499999999995921 0.25 +669 0.7500000000000953 0.7499999999995921 0.375 +670 0.7500000000000953 0.7499999999995921 0.5 +671 0.7500000000000953 0.7499999999995921 0.625 +672 0.7500000000000953 0.7499999999995921 0.75 +673 0.7500000000000953 0.7499999999995921 0.875 +674 0.7500000000002213 0.8749999999996281 0.125 +675 0.7500000000002213 0.8749999999996281 0.25 +676 0.7500000000002213 0.8749999999996281 0.375 +677 0.7500000000002213 0.8749999999996281 0.5 +678 0.7500000000002213 0.8749999999996281 0.625 +679 0.7500000000002213 0.8749999999996281 0.75 +680 0.7500000000002213 0.8749999999996281 0.875 +681 0.8749999999996492 0.1249999999998666 0.125 +682 0.8749999999996492 0.1249999999998666 0.25 +683 0.8749999999996492 0.1249999999998666 0.375 +684 0.8749999999996492 0.1249999999998666 0.5 +685 0.8749999999996492 0.1249999999998666 0.625 +686 0.8749999999996492 0.1249999999998666 0.75 +687 0.8749999999996492 0.1249999999998666 0.875 +688 0.8749999999996279 0.2499999999996146 0.125 +689 0.8749999999996279 0.2499999999996146 0.25 +690 0.8749999999996279 0.2499999999996146 0.375 +691 0.8749999999996279 0.2499999999996146 0.5 +692 0.8749999999996279 0.2499999999996146 0.625 +693 0.8749999999996279 0.2499999999996146 0.75 +694 0.8749999999996279 0.2499999999996146 0.875 +695 0.8749999999996068 0.3749999999993605 0.125 +696 0.8749999999996068 0.3749999999993605 0.25 +697 0.8749999999996068 0.3749999999993605 0.375 +698 0.8749999999996068 0.3749999999993605 0.5 +699 0.8749999999996068 0.3749999999993605 0.625 +700 0.8749999999996068 0.3749999999993605 0.75 +701 0.8749999999996068 0.3749999999993605 0.875 +702 0.8749999999995857 0.499999999999113 0.125 +703 0.8749999999995857 0.499999999999113 0.25 +704 0.8749999999995857 0.499999999999113 0.375 +705 0.8749999999995857 0.499999999999113 0.5 +706 0.8749999999995857 0.499999999999113 0.625 +707 0.8749999999995857 0.499999999999113 0.75 +708 0.8749999999995857 0.499999999999113 0.875 +709 0.8749999999995646 0.6249999999992862 0.125 +710 0.8749999999995646 0.6249999999992862 0.25 +711 0.8749999999995646 0.6249999999992862 0.375 +712 0.8749999999995646 0.6249999999992862 0.5 +713 0.8749999999995646 0.6249999999992862 0.625 +714 0.8749999999995646 0.6249999999992862 0.75 +715 0.8749999999995646 0.6249999999992862 0.875 +716 0.8749999999995435 0.7499999999994665 0.125 +717 0.8749999999995435 0.7499999999994665 0.25 +718 0.8749999999995435 0.7499999999994665 0.375 +719 0.8749999999995435 0.7499999999994665 0.5 +720 0.8749999999995435 0.7499999999994665 0.625 +721 0.8749999999995435 0.7499999999994665 0.75 +722 0.8749999999995435 0.7499999999994665 0.875 +723 0.8749999999995224 0.8749999999996492 0.125 +724 0.8749999999995224 0.8749999999996492 0.25 +725 0.8749999999995224 0.8749999999996492 0.375 +726 0.8749999999995224 0.8749999999996492 0.5 +727 0.8749999999995224 0.8749999999996492 0.625 +728 0.8749999999995224 0.8749999999996492 0.75 +729 0.8749999999995224 0.8749999999996492 0.875 +$EndNodes +$Elements +1000 +1 15 2 0 1 1 +2 15 2 0 2 2 +3 15 2 0 3 3 +4 15 2 0 4 4 +5 15 2 0 5 5 +6 15 2 0 6 6 +7 15 2 0 10 7 +8 15 2 0 14 8 +9 1 2 0 5 1 9 +10 1 2 0 5 9 10 +11 1 2 0 5 10 11 +12 1 2 0 5 11 12 +13 1 2 0 5 12 13 +14 1 2 0 5 13 14 +15 1 2 0 5 14 15 +16 1 2 0 5 15 2 +17 1 2 0 6 2 16 +18 1 2 0 6 16 17 +19 1 2 0 6 17 18 +20 1 2 0 6 18 19 +21 1 2 0 6 19 20 +22 1 2 0 6 20 21 +23 1 2 0 6 21 22 +24 1 2 0 6 22 3 +25 1 2 0 7 3 23 +26 1 2 0 7 23 24 +27 1 2 0 7 24 25 +28 1 2 0 7 25 26 +29 1 2 0 7 26 27 +30 1 2 0 7 27 28 +31 1 2 0 7 28 29 +32 1 2 0 7 29 4 +33 1 2 0 8 4 30 +34 1 2 0 8 30 31 +35 1 2 0 8 31 32 +36 1 2 0 8 32 33 +37 1 2 0 8 33 34 +38 1 2 0 8 34 35 +39 1 2 0 8 35 36 +40 1 2 0 8 36 1 +41 1 2 0 12 5 37 +42 1 2 0 12 37 38 +43 1 2 0 12 38 39 +44 1 2 0 12 39 40 +45 1 2 0 12 40 41 +46 1 2 0 12 41 42 +47 1 2 0 12 42 43 +48 1 2 0 12 43 6 +49 1 2 0 13 6 44 +50 1 2 0 13 44 45 +51 1 2 0 13 45 46 +52 1 2 0 13 46 47 +53 1 2 0 13 47 48 +54 1 2 0 13 48 49 +55 1 2 0 13 49 50 +56 1 2 0 13 50 7 +57 1 2 0 14 7 51 +58 1 2 0 14 51 52 +59 1 2 0 14 52 53 +60 1 2 0 14 53 54 +61 1 2 0 14 54 55 +62 1 2 0 14 55 56 +63 1 2 0 14 56 57 +64 1 2 0 14 57 8 +65 1 2 0 15 8 58 +66 1 2 0 15 58 59 +67 1 2 0 15 59 60 +68 1 2 0 15 60 61 +69 1 2 0 15 61 62 +70 1 2 0 15 62 63 +71 1 2 0 15 63 64 +72 1 2 0 15 64 5 +73 1 2 0 17 1 65 +74 1 2 0 17 65 66 +75 1 2 0 17 66 67 +76 1 2 0 17 67 68 +77 1 2 0 17 68 69 +78 1 2 0 17 69 70 +79 1 2 0 17 70 71 +80 1 2 0 17 71 5 +81 1 2 0 18 2 72 +82 1 2 0 18 72 73 +83 1 2 0 18 73 74 +84 1 2 0 18 74 75 +85 1 2 0 18 75 76 +86 1 2 0 18 76 77 +87 1 2 0 18 77 78 +88 1 2 0 18 78 6 +89 1 2 0 22 3 79 +90 1 2 0 22 79 80 +91 1 2 0 22 80 81 +92 1 2 0 22 81 82 +93 1 2 0 22 82 83 +94 1 2 0 22 83 84 +95 1 2 0 22 84 85 +96 1 2 0 22 85 7 +97 1 2 0 26 4 86 +98 1 2 0 26 86 87 +99 1 2 0 26 87 88 +100 1 2 0 26 88 89 +101 1 2 0 26 89 90 +102 1 2 0 26 90 91 +103 1 2 0 26 91 92 +104 1 2 0 26 92 8 +105 3 2 0 10 1 9 93 36 +106 3 2 0 10 36 93 94 35 +107 3 2 0 10 35 94 95 34 +108 3 2 0 10 34 95 96 33 +109 3 2 0 10 33 96 97 32 +110 3 2 0 10 32 97 98 31 +111 3 2 0 10 31 98 99 30 +112 3 2 0 10 30 99 29 4 +113 3 2 0 10 9 10 100 93 +114 3 2 0 10 93 100 101 94 +115 3 2 0 10 94 101 102 95 +116 3 2 0 10 95 102 103 96 +117 3 2 0 10 96 103 104 97 +118 3 2 0 10 97 104 105 98 +119 3 2 0 10 98 105 106 99 +120 3 2 0 10 99 106 28 29 +121 3 2 0 10 10 11 107 100 +122 3 2 0 10 100 107 108 101 +123 3 2 0 10 101 108 109 102 +124 3 2 0 10 102 109 110 103 +125 3 2 0 10 103 110 111 104 +126 3 2 0 10 104 111 112 105 +127 3 2 0 10 105 112 113 106 +128 3 2 0 10 106 113 27 28 +129 3 2 0 10 11 12 114 107 +130 3 2 0 10 107 114 115 108 +131 3 2 0 10 108 115 116 109 +132 3 2 0 10 109 116 117 110 +133 3 2 0 10 110 117 118 111 +134 3 2 0 10 111 118 119 112 +135 3 2 0 10 112 119 120 113 +136 3 2 0 10 113 120 26 27 +137 3 2 0 10 12 13 121 114 +138 3 2 0 10 114 121 122 115 +139 3 2 0 10 115 122 123 116 +140 3 2 0 10 116 123 124 117 +141 3 2 0 10 117 124 125 118 +142 3 2 0 10 118 125 126 119 +143 3 2 0 10 119 126 127 120 +144 3 2 0 10 120 127 25 26 +145 3 2 0 10 13 14 128 121 +146 3 2 0 10 121 128 129 122 +147 3 2 0 10 122 129 130 123 +148 3 2 0 10 123 130 131 124 +149 3 2 0 10 124 131 132 125 +150 3 2 0 10 125 132 133 126 +151 3 2 0 10 126 133 134 127 +152 3 2 0 10 127 134 24 25 +153 3 2 0 10 14 15 135 128 +154 3 2 0 10 128 135 136 129 +155 3 2 0 10 129 136 137 130 +156 3 2 0 10 130 137 138 131 +157 3 2 0 10 131 138 139 132 +158 3 2 0 10 132 139 140 133 +159 3 2 0 10 133 140 141 134 +160 3 2 0 10 134 141 23 24 +161 3 2 0 10 15 2 16 135 +162 3 2 0 10 135 16 17 136 +163 3 2 0 10 136 17 18 137 +164 3 2 0 10 137 18 19 138 +165 3 2 0 10 138 19 20 139 +166 3 2 0 10 139 20 21 140 +167 3 2 0 10 140 21 22 141 +168 3 2 0 10 141 22 3 23 +169 3 2 0 19 1 9 142 65 +170 3 2 0 19 65 142 143 66 +171 3 2 0 19 66 143 144 67 +172 3 2 0 19 67 144 145 68 +173 3 2 0 19 68 145 146 69 +174 3 2 0 19 69 146 147 70 +175 3 2 0 19 70 147 148 71 +176 3 2 0 19 71 148 37 5 +177 3 2 0 19 9 10 149 142 +178 3 2 0 19 142 149 150 143 +179 3 2 0 19 143 150 151 144 +180 3 2 0 19 144 151 152 145 +181 3 2 0 19 145 152 153 146 +182 3 2 0 19 146 153 154 147 +183 3 2 0 19 147 154 155 148 +184 3 2 0 19 148 155 38 37 +185 3 2 0 19 10 11 156 149 +186 3 2 0 19 149 156 157 150 +187 3 2 0 19 150 157 158 151 +188 3 2 0 19 151 158 159 152 +189 3 2 0 19 152 159 160 153 +190 3 2 0 19 153 160 161 154 +191 3 2 0 19 154 161 162 155 +192 3 2 0 19 155 162 39 38 +193 3 2 0 19 11 12 163 156 +194 3 2 0 19 156 163 164 157 +195 3 2 0 19 157 164 165 158 +196 3 2 0 19 158 165 166 159 +197 3 2 0 19 159 166 167 160 +198 3 2 0 19 160 167 168 161 +199 3 2 0 19 161 168 169 162 +200 3 2 0 19 162 169 40 39 +201 3 2 0 19 12 13 170 163 +202 3 2 0 19 163 170 171 164 +203 3 2 0 19 164 171 172 165 +204 3 2 0 19 165 172 173 166 +205 3 2 0 19 166 173 174 167 +206 3 2 0 19 167 174 175 168 +207 3 2 0 19 168 175 176 169 +208 3 2 0 19 169 176 41 40 +209 3 2 0 19 13 14 177 170 +210 3 2 0 19 170 177 178 171 +211 3 2 0 19 171 178 179 172 +212 3 2 0 19 172 179 180 173 +213 3 2 0 19 173 180 181 174 +214 3 2 0 19 174 181 182 175 +215 3 2 0 19 175 182 183 176 +216 3 2 0 19 176 183 42 41 +217 3 2 0 19 14 15 184 177 +218 3 2 0 19 177 184 185 178 +219 3 2 0 19 178 185 186 179 +220 3 2 0 19 179 186 187 180 +221 3 2 0 19 180 187 188 181 +222 3 2 0 19 181 188 189 182 +223 3 2 0 19 182 189 190 183 +224 3 2 0 19 183 190 43 42 +225 3 2 0 19 15 2 72 184 +226 3 2 0 19 184 72 73 185 +227 3 2 0 19 185 73 74 186 +228 3 2 0 19 186 74 75 187 +229 3 2 0 19 187 75 76 188 +230 3 2 0 19 188 76 77 189 +231 3 2 0 19 189 77 78 190 +232 3 2 0 19 190 78 6 43 +233 3 2 0 23 2 16 191 72 +234 3 2 0 23 72 191 192 73 +235 3 2 0 23 73 192 193 74 +236 3 2 0 23 74 193 194 75 +237 3 2 0 23 75 194 195 76 +238 3 2 0 23 76 195 196 77 +239 3 2 0 23 77 196 197 78 +240 3 2 0 23 78 197 44 6 +241 3 2 0 23 16 17 198 191 +242 3 2 0 23 191 198 199 192 +243 3 2 0 23 192 199 200 193 +244 3 2 0 23 193 200 201 194 +245 3 2 0 23 194 201 202 195 +246 3 2 0 23 195 202 203 196 +247 3 2 0 23 196 203 204 197 +248 3 2 0 23 197 204 45 44 +249 3 2 0 23 17 18 205 198 +250 3 2 0 23 198 205 206 199 +251 3 2 0 23 199 206 207 200 +252 3 2 0 23 200 207 208 201 +253 3 2 0 23 201 208 209 202 +254 3 2 0 23 202 209 210 203 +255 3 2 0 23 203 210 211 204 +256 3 2 0 23 204 211 46 45 +257 3 2 0 23 18 19 212 205 +258 3 2 0 23 205 212 213 206 +259 3 2 0 23 206 213 214 207 +260 3 2 0 23 207 214 215 208 +261 3 2 0 23 208 215 216 209 +262 3 2 0 23 209 216 217 210 +263 3 2 0 23 210 217 218 211 +264 3 2 0 23 211 218 47 46 +265 3 2 0 23 19 20 219 212 +266 3 2 0 23 212 219 220 213 +267 3 2 0 23 213 220 221 214 +268 3 2 0 23 214 221 222 215 +269 3 2 0 23 215 222 223 216 +270 3 2 0 23 216 223 224 217 +271 3 2 0 23 217 224 225 218 +272 3 2 0 23 218 225 48 47 +273 3 2 0 23 20 21 226 219 +274 3 2 0 23 219 226 227 220 +275 3 2 0 23 220 227 228 221 +276 3 2 0 23 221 228 229 222 +277 3 2 0 23 222 229 230 223 +278 3 2 0 23 223 230 231 224 +279 3 2 0 23 224 231 232 225 +280 3 2 0 23 225 232 49 48 +281 3 2 0 23 21 22 233 226 +282 3 2 0 23 226 233 234 227 +283 3 2 0 23 227 234 235 228 +284 3 2 0 23 228 235 236 229 +285 3 2 0 23 229 236 237 230 +286 3 2 0 23 230 237 238 231 +287 3 2 0 23 231 238 239 232 +288 3 2 0 23 232 239 50 49 +289 3 2 0 23 22 3 79 233 +290 3 2 0 23 233 79 80 234 +291 3 2 0 23 234 80 81 235 +292 3 2 0 23 235 81 82 236 +293 3 2 0 23 236 82 83 237 +294 3 2 0 23 237 83 84 238 +295 3 2 0 23 238 84 85 239 +296 3 2 0 23 239 85 7 50 +297 3 2 0 27 3 23 240 79 +298 3 2 0 27 79 240 241 80 +299 3 2 0 27 80 241 242 81 +300 3 2 0 27 81 242 243 82 +301 3 2 0 27 82 243 244 83 +302 3 2 0 27 83 244 245 84 +303 3 2 0 27 84 245 246 85 +304 3 2 0 27 85 246 51 7 +305 3 2 0 27 23 24 247 240 +306 3 2 0 27 240 247 248 241 +307 3 2 0 27 241 248 249 242 +308 3 2 0 27 242 249 250 243 +309 3 2 0 27 243 250 251 244 +310 3 2 0 27 244 251 252 245 +311 3 2 0 27 245 252 253 246 +312 3 2 0 27 246 253 52 51 +313 3 2 0 27 24 25 254 247 +314 3 2 0 27 247 254 255 248 +315 3 2 0 27 248 255 256 249 +316 3 2 0 27 249 256 257 250 +317 3 2 0 27 250 257 258 251 +318 3 2 0 27 251 258 259 252 +319 3 2 0 27 252 259 260 253 +320 3 2 0 27 253 260 53 52 +321 3 2 0 27 25 26 261 254 +322 3 2 0 27 254 261 262 255 +323 3 2 0 27 255 262 263 256 +324 3 2 0 27 256 263 264 257 +325 3 2 0 27 257 264 265 258 +326 3 2 0 27 258 265 266 259 +327 3 2 0 27 259 266 267 260 +328 3 2 0 27 260 267 54 53 +329 3 2 0 27 26 27 268 261 +330 3 2 0 27 261 268 269 262 +331 3 2 0 27 262 269 270 263 +332 3 2 0 27 263 270 271 264 +333 3 2 0 27 264 271 272 265 +334 3 2 0 27 265 272 273 266 +335 3 2 0 27 266 273 274 267 +336 3 2 0 27 267 274 55 54 +337 3 2 0 27 27 28 275 268 +338 3 2 0 27 268 275 276 269 +339 3 2 0 27 269 276 277 270 +340 3 2 0 27 270 277 278 271 +341 3 2 0 27 271 278 279 272 +342 3 2 0 27 272 279 280 273 +343 3 2 0 27 273 280 281 274 +344 3 2 0 27 274 281 56 55 +345 3 2 0 27 28 29 282 275 +346 3 2 0 27 275 282 283 276 +347 3 2 0 27 276 283 284 277 +348 3 2 0 27 277 284 285 278 +349 3 2 0 27 278 285 286 279 +350 3 2 0 27 279 286 287 280 +351 3 2 0 27 280 287 288 281 +352 3 2 0 27 281 288 57 56 +353 3 2 0 27 29 4 86 282 +354 3 2 0 27 282 86 87 283 +355 3 2 0 27 283 87 88 284 +356 3 2 0 27 284 88 89 285 +357 3 2 0 27 285 89 90 286 +358 3 2 0 27 286 90 91 287 +359 3 2 0 27 287 91 92 288 +360 3 2 0 27 288 92 8 57 +361 3 2 0 31 4 30 289 86 +362 3 2 0 31 86 289 290 87 +363 3 2 0 31 87 290 291 88 +364 3 2 0 31 88 291 292 89 +365 3 2 0 31 89 292 293 90 +366 3 2 0 31 90 293 294 91 +367 3 2 0 31 91 294 295 92 +368 3 2 0 31 92 295 58 8 +369 3 2 0 31 30 31 296 289 +370 3 2 0 31 289 296 297 290 +371 3 2 0 31 290 297 298 291 +372 3 2 0 31 291 298 299 292 +373 3 2 0 31 292 299 300 293 +374 3 2 0 31 293 300 301 294 +375 3 2 0 31 294 301 302 295 +376 3 2 0 31 295 302 59 58 +377 3 2 0 31 31 32 303 296 +378 3 2 0 31 296 303 304 297 +379 3 2 0 31 297 304 305 298 +380 3 2 0 31 298 305 306 299 +381 3 2 0 31 299 306 307 300 +382 3 2 0 31 300 307 308 301 +383 3 2 0 31 301 308 309 302 +384 3 2 0 31 302 309 60 59 +385 3 2 0 31 32 33 310 303 +386 3 2 0 31 303 310 311 304 +387 3 2 0 31 304 311 312 305 +388 3 2 0 31 305 312 313 306 +389 3 2 0 31 306 313 314 307 +390 3 2 0 31 307 314 315 308 +391 3 2 0 31 308 315 316 309 +392 3 2 0 31 309 316 61 60 +393 3 2 0 31 33 34 317 310 +394 3 2 0 31 310 317 318 311 +395 3 2 0 31 311 318 319 312 +396 3 2 0 31 312 319 320 313 +397 3 2 0 31 313 320 321 314 +398 3 2 0 31 314 321 322 315 +399 3 2 0 31 315 322 323 316 +400 3 2 0 31 316 323 62 61 +401 3 2 0 31 34 35 324 317 +402 3 2 0 31 317 324 325 318 +403 3 2 0 31 318 325 326 319 +404 3 2 0 31 319 326 327 320 +405 3 2 0 31 320 327 328 321 +406 3 2 0 31 321 328 329 322 +407 3 2 0 31 322 329 330 323 +408 3 2 0 31 323 330 63 62 +409 3 2 0 31 35 36 331 324 +410 3 2 0 31 324 331 332 325 +411 3 2 0 31 325 332 333 326 +412 3 2 0 31 326 333 334 327 +413 3 2 0 31 327 334 335 328 +414 3 2 0 31 328 335 336 329 +415 3 2 0 31 329 336 337 330 +416 3 2 0 31 330 337 64 63 +417 3 2 0 31 36 1 65 331 +418 3 2 0 31 331 65 66 332 +419 3 2 0 31 332 66 67 333 +420 3 2 0 31 333 67 68 334 +421 3 2 0 31 334 68 69 335 +422 3 2 0 31 335 69 70 336 +423 3 2 0 31 336 70 71 337 +424 3 2 0 31 337 71 5 64 +425 3 2 0 32 5 37 338 64 +426 3 2 0 32 64 338 339 63 +427 3 2 0 32 63 339 340 62 +428 3 2 0 32 62 340 341 61 +429 3 2 0 32 61 341 342 60 +430 3 2 0 32 60 342 343 59 +431 3 2 0 32 59 343 344 58 +432 3 2 0 32 58 344 57 8 +433 3 2 0 32 37 38 345 338 +434 3 2 0 32 338 345 346 339 +435 3 2 0 32 339 346 347 340 +436 3 2 0 32 340 347 348 341 +437 3 2 0 32 341 348 349 342 +438 3 2 0 32 342 349 350 343 +439 3 2 0 32 343 350 351 344 +440 3 2 0 32 344 351 56 57 +441 3 2 0 32 38 39 352 345 +442 3 2 0 32 345 352 353 346 +443 3 2 0 32 346 353 354 347 +444 3 2 0 32 347 354 355 348 +445 3 2 0 32 348 355 356 349 +446 3 2 0 32 349 356 357 350 +447 3 2 0 32 350 357 358 351 +448 3 2 0 32 351 358 55 56 +449 3 2 0 32 39 40 359 352 +450 3 2 0 32 352 359 360 353 +451 3 2 0 32 353 360 361 354 +452 3 2 0 32 354 361 362 355 +453 3 2 0 32 355 362 363 356 +454 3 2 0 32 356 363 364 357 +455 3 2 0 32 357 364 365 358 +456 3 2 0 32 358 365 54 55 +457 3 2 0 32 40 41 366 359 +458 3 2 0 32 359 366 367 360 +459 3 2 0 32 360 367 368 361 +460 3 2 0 32 361 368 369 362 +461 3 2 0 32 362 369 370 363 +462 3 2 0 32 363 370 371 364 +463 3 2 0 32 364 371 372 365 +464 3 2 0 32 365 372 53 54 +465 3 2 0 32 41 42 373 366 +466 3 2 0 32 366 373 374 367 +467 3 2 0 32 367 374 375 368 +468 3 2 0 32 368 375 376 369 +469 3 2 0 32 369 376 377 370 +470 3 2 0 32 370 377 378 371 +471 3 2 0 32 371 378 379 372 +472 3 2 0 32 372 379 52 53 +473 3 2 0 32 42 43 380 373 +474 3 2 0 32 373 380 381 374 +475 3 2 0 32 374 381 382 375 +476 3 2 0 32 375 382 383 376 +477 3 2 0 32 376 383 384 377 +478 3 2 0 32 377 384 385 378 +479 3 2 0 32 378 385 386 379 +480 3 2 0 32 379 386 51 52 +481 3 2 0 32 43 6 44 380 +482 3 2 0 32 380 44 45 381 +483 3 2 0 32 381 45 46 382 +484 3 2 0 32 382 46 47 383 +485 3 2 0 32 383 47 48 384 +486 3 2 0 32 384 48 49 385 +487 3 2 0 32 385 49 50 386 +488 3 2 0 32 386 50 7 51 +489 5 2 0 1 1 9 93 36 65 142 387 331 +490 5 2 0 1 65 142 387 331 66 143 388 332 +491 5 2 0 1 66 143 388 332 67 144 389 333 +492 5 2 0 1 67 144 389 333 68 145 390 334 +493 5 2 0 1 68 145 390 334 69 146 391 335 +494 5 2 0 1 69 146 391 335 70 147 392 336 +495 5 2 0 1 70 147 392 336 71 148 393 337 +496 5 2 0 1 71 148 393 337 5 37 338 64 +497 5 2 0 1 36 93 94 35 331 387 394 324 +498 5 2 0 1 331 387 394 324 332 388 395 325 +499 5 2 0 1 332 388 395 325 333 389 396 326 +500 5 2 0 1 333 389 396 326 334 390 397 327 +501 5 2 0 1 334 390 397 327 335 391 398 328 +502 5 2 0 1 335 391 398 328 336 392 399 329 +503 5 2 0 1 336 392 399 329 337 393 400 330 +504 5 2 0 1 337 393 400 330 64 338 339 63 +505 5 2 0 1 35 94 95 34 324 394 401 317 +506 5 2 0 1 324 394 401 317 325 395 402 318 +507 5 2 0 1 325 395 402 318 326 396 403 319 +508 5 2 0 1 326 396 403 319 327 397 404 320 +509 5 2 0 1 327 397 404 320 328 398 405 321 +510 5 2 0 1 328 398 405 321 329 399 406 322 +511 5 2 0 1 329 399 406 322 330 400 407 323 +512 5 2 0 1 330 400 407 323 63 339 340 62 +513 5 2 0 1 34 95 96 33 317 401 408 310 +514 5 2 0 1 317 401 408 310 318 402 409 311 +515 5 2 0 1 318 402 409 311 319 403 410 312 +516 5 2 0 1 319 403 410 312 320 404 411 313 +517 5 2 0 1 320 404 411 313 321 405 412 314 +518 5 2 0 1 321 405 412 314 322 406 413 315 +519 5 2 0 1 322 406 413 315 323 407 414 316 +520 5 2 0 1 323 407 414 316 62 340 341 61 +521 5 2 0 1 33 96 97 32 310 408 415 303 +522 5 2 0 1 310 408 415 303 311 409 416 304 +523 5 2 0 1 311 409 416 304 312 410 417 305 +524 5 2 0 1 312 410 417 305 313 411 418 306 +525 5 2 0 1 313 411 418 306 314 412 419 307 +526 5 2 0 1 314 412 419 307 315 413 420 308 +527 5 2 0 1 315 413 420 308 316 414 421 309 +528 5 2 0 1 316 414 421 309 61 341 342 60 +529 5 2 0 1 32 97 98 31 303 415 422 296 +530 5 2 0 1 303 415 422 296 304 416 423 297 +531 5 2 0 1 304 416 423 297 305 417 424 298 +532 5 2 0 1 305 417 424 298 306 418 425 299 +533 5 2 0 1 306 418 425 299 307 419 426 300 +534 5 2 0 1 307 419 426 300 308 420 427 301 +535 5 2 0 1 308 420 427 301 309 421 428 302 +536 5 2 0 1 309 421 428 302 60 342 343 59 +537 5 2 0 1 31 98 99 30 296 422 429 289 +538 5 2 0 1 296 422 429 289 297 423 430 290 +539 5 2 0 1 297 423 430 290 298 424 431 291 +540 5 2 0 1 298 424 431 291 299 425 432 292 +541 5 2 0 1 299 425 432 292 300 426 433 293 +542 5 2 0 1 300 426 433 293 301 427 434 294 +543 5 2 0 1 301 427 434 294 302 428 435 295 +544 5 2 0 1 302 428 435 295 59 343 344 58 +545 5 2 0 1 30 99 29 4 289 429 282 86 +546 5 2 0 1 289 429 282 86 290 430 283 87 +547 5 2 0 1 290 430 283 87 291 431 284 88 +548 5 2 0 1 291 431 284 88 292 432 285 89 +549 5 2 0 1 292 432 285 89 293 433 286 90 +550 5 2 0 1 293 433 286 90 294 434 287 91 +551 5 2 0 1 294 434 287 91 295 435 288 92 +552 5 2 0 1 295 435 288 92 58 344 57 8 +553 5 2 0 1 9 10 100 93 142 149 436 387 +554 5 2 0 1 142 149 436 387 143 150 437 388 +555 5 2 0 1 143 150 437 388 144 151 438 389 +556 5 2 0 1 144 151 438 389 145 152 439 390 +557 5 2 0 1 145 152 439 390 146 153 440 391 +558 5 2 0 1 146 153 440 391 147 154 441 392 +559 5 2 0 1 147 154 441 392 148 155 442 393 +560 5 2 0 1 148 155 442 393 37 38 345 338 +561 5 2 0 1 93 100 101 94 387 436 443 394 +562 5 2 0 1 387 436 443 394 388 437 444 395 +563 5 2 0 1 388 437 444 395 389 438 445 396 +564 5 2 0 1 389 438 445 396 390 439 446 397 +565 5 2 0 1 390 439 446 397 391 440 447 398 +566 5 2 0 1 391 440 447 398 392 441 448 399 +567 5 2 0 1 392 441 448 399 393 442 449 400 +568 5 2 0 1 393 442 449 400 338 345 346 339 +569 5 2 0 1 94 101 102 95 394 443 450 401 +570 5 2 0 1 394 443 450 401 395 444 451 402 +571 5 2 0 1 395 444 451 402 396 445 452 403 +572 5 2 0 1 396 445 452 403 397 446 453 404 +573 5 2 0 1 397 446 453 404 398 447 454 405 +574 5 2 0 1 398 447 454 405 399 448 455 406 +575 5 2 0 1 399 448 455 406 400 449 456 407 +576 5 2 0 1 400 449 456 407 339 346 347 340 +577 5 2 0 1 95 102 103 96 401 450 457 408 +578 5 2 0 1 401 450 457 408 402 451 458 409 +579 5 2 0 1 402 451 458 409 403 452 459 410 +580 5 2 0 1 403 452 459 410 404 453 460 411 +581 5 2 0 1 404 453 460 411 405 454 461 412 +582 5 2 0 1 405 454 461 412 406 455 462 413 +583 5 2 0 1 406 455 462 413 407 456 463 414 +584 5 2 0 1 407 456 463 414 340 347 348 341 +585 5 2 0 1 96 103 104 97 408 457 464 415 +586 5 2 0 1 408 457 464 415 409 458 465 416 +587 5 2 0 1 409 458 465 416 410 459 466 417 +588 5 2 0 1 410 459 466 417 411 460 467 418 +589 5 2 0 1 411 460 467 418 412 461 468 419 +590 5 2 0 1 412 461 468 419 413 462 469 420 +591 5 2 0 1 413 462 469 420 414 463 470 421 +592 5 2 0 1 414 463 470 421 341 348 349 342 +593 5 2 0 1 97 104 105 98 415 464 471 422 +594 5 2 0 1 415 464 471 422 416 465 472 423 +595 5 2 0 1 416 465 472 423 417 466 473 424 +596 5 2 0 1 417 466 473 424 418 467 474 425 +597 5 2 0 1 418 467 474 425 419 468 475 426 +598 5 2 0 1 419 468 475 426 420 469 476 427 +599 5 2 0 1 420 469 476 427 421 470 477 428 +600 5 2 0 1 421 470 477 428 342 349 350 343 +601 5 2 0 1 98 105 106 99 422 471 478 429 +602 5 2 0 1 422 471 478 429 423 472 479 430 +603 5 2 0 1 423 472 479 430 424 473 480 431 +604 5 2 0 1 424 473 480 431 425 474 481 432 +605 5 2 0 1 425 474 481 432 426 475 482 433 +606 5 2 0 1 426 475 482 433 427 476 483 434 +607 5 2 0 1 427 476 483 434 428 477 484 435 +608 5 2 0 1 428 477 484 435 343 350 351 344 +609 5 2 0 1 99 106 28 29 429 478 275 282 +610 5 2 0 1 429 478 275 282 430 479 276 283 +611 5 2 0 1 430 479 276 283 431 480 277 284 +612 5 2 0 1 431 480 277 284 432 481 278 285 +613 5 2 0 1 432 481 278 285 433 482 279 286 +614 5 2 0 1 433 482 279 286 434 483 280 287 +615 5 2 0 1 434 483 280 287 435 484 281 288 +616 5 2 0 1 435 484 281 288 344 351 56 57 +617 5 2 0 1 10 11 107 100 149 156 485 436 +618 5 2 0 1 149 156 485 436 150 157 486 437 +619 5 2 0 1 150 157 486 437 151 158 487 438 +620 5 2 0 1 151 158 487 438 152 159 488 439 +621 5 2 0 1 152 159 488 439 153 160 489 440 +622 5 2 0 1 153 160 489 440 154 161 490 441 +623 5 2 0 1 154 161 490 441 155 162 491 442 +624 5 2 0 1 155 162 491 442 38 39 352 345 +625 5 2 0 1 100 107 108 101 436 485 492 443 +626 5 2 0 1 436 485 492 443 437 486 493 444 +627 5 2 0 1 437 486 493 444 438 487 494 445 +628 5 2 0 1 438 487 494 445 439 488 495 446 +629 5 2 0 1 439 488 495 446 440 489 496 447 +630 5 2 0 1 440 489 496 447 441 490 497 448 +631 5 2 0 1 441 490 497 448 442 491 498 449 +632 5 2 0 1 442 491 498 449 345 352 353 346 +633 5 2 0 1 101 108 109 102 443 492 499 450 +634 5 2 0 1 443 492 499 450 444 493 500 451 +635 5 2 0 1 444 493 500 451 445 494 501 452 +636 5 2 0 1 445 494 501 452 446 495 502 453 +637 5 2 0 1 446 495 502 453 447 496 503 454 +638 5 2 0 1 447 496 503 454 448 497 504 455 +639 5 2 0 1 448 497 504 455 449 498 505 456 +640 5 2 0 1 449 498 505 456 346 353 354 347 +641 5 2 0 1 102 109 110 103 450 499 506 457 +642 5 2 0 1 450 499 506 457 451 500 507 458 +643 5 2 0 1 451 500 507 458 452 501 508 459 +644 5 2 0 1 452 501 508 459 453 502 509 460 +645 5 2 0 1 453 502 509 460 454 503 510 461 +646 5 2 0 1 454 503 510 461 455 504 511 462 +647 5 2 0 1 455 504 511 462 456 505 512 463 +648 5 2 0 1 456 505 512 463 347 354 355 348 +649 5 2 0 1 103 110 111 104 457 506 513 464 +650 5 2 0 1 457 506 513 464 458 507 514 465 +651 5 2 0 1 458 507 514 465 459 508 515 466 +652 5 2 0 1 459 508 515 466 460 509 516 467 +653 5 2 0 1 460 509 516 467 461 510 517 468 +654 5 2 0 1 461 510 517 468 462 511 518 469 +655 5 2 0 1 462 511 518 469 463 512 519 470 +656 5 2 0 1 463 512 519 470 348 355 356 349 +657 5 2 0 1 104 111 112 105 464 513 520 471 +658 5 2 0 1 464 513 520 471 465 514 521 472 +659 5 2 0 1 465 514 521 472 466 515 522 473 +660 5 2 0 1 466 515 522 473 467 516 523 474 +661 5 2 0 1 467 516 523 474 468 517 524 475 +662 5 2 0 1 468 517 524 475 469 518 525 476 +663 5 2 0 1 469 518 525 476 470 519 526 477 +664 5 2 0 1 470 519 526 477 349 356 357 350 +665 5 2 0 1 105 112 113 106 471 520 527 478 +666 5 2 0 1 471 520 527 478 472 521 528 479 +667 5 2 0 1 472 521 528 479 473 522 529 480 +668 5 2 0 1 473 522 529 480 474 523 530 481 +669 5 2 0 1 474 523 530 481 475 524 531 482 +670 5 2 0 1 475 524 531 482 476 525 532 483 +671 5 2 0 1 476 525 532 483 477 526 533 484 +672 5 2 0 1 477 526 533 484 350 357 358 351 +673 5 2 0 1 106 113 27 28 478 527 268 275 +674 5 2 0 1 478 527 268 275 479 528 269 276 +675 5 2 0 1 479 528 269 276 480 529 270 277 +676 5 2 0 1 480 529 270 277 481 530 271 278 +677 5 2 0 1 481 530 271 278 482 531 272 279 +678 5 2 0 1 482 531 272 279 483 532 273 280 +679 5 2 0 1 483 532 273 280 484 533 274 281 +680 5 2 0 1 484 533 274 281 351 358 55 56 +681 5 2 0 1 11 12 114 107 156 163 534 485 +682 5 2 0 1 156 163 534 485 157 164 535 486 +683 5 2 0 1 157 164 535 486 158 165 536 487 +684 5 2 0 1 158 165 536 487 159 166 537 488 +685 5 2 0 1 159 166 537 488 160 167 538 489 +686 5 2 0 1 160 167 538 489 161 168 539 490 +687 5 2 0 1 161 168 539 490 162 169 540 491 +688 5 2 0 1 162 169 540 491 39 40 359 352 +689 5 2 0 1 107 114 115 108 485 534 541 492 +690 5 2 0 1 485 534 541 492 486 535 542 493 +691 5 2 0 1 486 535 542 493 487 536 543 494 +692 5 2 0 1 487 536 543 494 488 537 544 495 +693 5 2 0 1 488 537 544 495 489 538 545 496 +694 5 2 0 1 489 538 545 496 490 539 546 497 +695 5 2 0 1 490 539 546 497 491 540 547 498 +696 5 2 0 1 491 540 547 498 352 359 360 353 +697 5 2 0 1 108 115 116 109 492 541 548 499 +698 5 2 0 1 492 541 548 499 493 542 549 500 +699 5 2 0 1 493 542 549 500 494 543 550 501 +700 5 2 0 1 494 543 550 501 495 544 551 502 +701 5 2 0 1 495 544 551 502 496 545 552 503 +702 5 2 0 1 496 545 552 503 497 546 553 504 +703 5 2 0 1 497 546 553 504 498 547 554 505 +704 5 2 0 1 498 547 554 505 353 360 361 354 +705 5 2 0 1 109 116 117 110 499 548 555 506 +706 5 2 0 1 499 548 555 506 500 549 556 507 +707 5 2 0 1 500 549 556 507 501 550 557 508 +708 5 2 0 1 501 550 557 508 502 551 558 509 +709 5 2 0 1 502 551 558 509 503 552 559 510 +710 5 2 0 1 503 552 559 510 504 553 560 511 +711 5 2 0 1 504 553 560 511 505 554 561 512 +712 5 2 0 1 505 554 561 512 354 361 362 355 +713 5 2 0 1 110 117 118 111 506 555 562 513 +714 5 2 0 1 506 555 562 513 507 556 563 514 +715 5 2 0 1 507 556 563 514 508 557 564 515 +716 5 2 0 1 508 557 564 515 509 558 565 516 +717 5 2 0 1 509 558 565 516 510 559 566 517 +718 5 2 0 1 510 559 566 517 511 560 567 518 +719 5 2 0 1 511 560 567 518 512 561 568 519 +720 5 2 0 1 512 561 568 519 355 362 363 356 +721 5 2 0 1 111 118 119 112 513 562 569 520 +722 5 2 0 1 513 562 569 520 514 563 570 521 +723 5 2 0 1 514 563 570 521 515 564 571 522 +724 5 2 0 1 515 564 571 522 516 565 572 523 +725 5 2 0 1 516 565 572 523 517 566 573 524 +726 5 2 0 1 517 566 573 524 518 567 574 525 +727 5 2 0 1 518 567 574 525 519 568 575 526 +728 5 2 0 1 519 568 575 526 356 363 364 357 +729 5 2 0 1 112 119 120 113 520 569 576 527 +730 5 2 0 1 520 569 576 527 521 570 577 528 +731 5 2 0 1 521 570 577 528 522 571 578 529 +732 5 2 0 1 522 571 578 529 523 572 579 530 +733 5 2 0 1 523 572 579 530 524 573 580 531 +734 5 2 0 1 524 573 580 531 525 574 581 532 +735 5 2 0 1 525 574 581 532 526 575 582 533 +736 5 2 0 1 526 575 582 533 357 364 365 358 +737 5 2 0 1 113 120 26 27 527 576 261 268 +738 5 2 0 1 527 576 261 268 528 577 262 269 +739 5 2 0 1 528 577 262 269 529 578 263 270 +740 5 2 0 1 529 578 263 270 530 579 264 271 +741 5 2 0 1 530 579 264 271 531 580 265 272 +742 5 2 0 1 531 580 265 272 532 581 266 273 +743 5 2 0 1 532 581 266 273 533 582 267 274 +744 5 2 0 1 533 582 267 274 358 365 54 55 +745 5 2 0 1 12 13 121 114 163 170 583 534 +746 5 2 0 1 163 170 583 534 164 171 584 535 +747 5 2 0 1 164 171 584 535 165 172 585 536 +748 5 2 0 1 165 172 585 536 166 173 586 537 +749 5 2 0 1 166 173 586 537 167 174 587 538 +750 5 2 0 1 167 174 587 538 168 175 588 539 +751 5 2 0 1 168 175 588 539 169 176 589 540 +752 5 2 0 1 169 176 589 540 40 41 366 359 +753 5 2 0 1 114 121 122 115 534 583 590 541 +754 5 2 0 1 534 583 590 541 535 584 591 542 +755 5 2 0 1 535 584 591 542 536 585 592 543 +756 5 2 0 1 536 585 592 543 537 586 593 544 +757 5 2 0 1 537 586 593 544 538 587 594 545 +758 5 2 0 1 538 587 594 545 539 588 595 546 +759 5 2 0 1 539 588 595 546 540 589 596 547 +760 5 2 0 1 540 589 596 547 359 366 367 360 +761 5 2 0 1 115 122 123 116 541 590 597 548 +762 5 2 0 1 541 590 597 548 542 591 598 549 +763 5 2 0 1 542 591 598 549 543 592 599 550 +764 5 2 0 1 543 592 599 550 544 593 600 551 +765 5 2 0 1 544 593 600 551 545 594 601 552 +766 5 2 0 1 545 594 601 552 546 595 602 553 +767 5 2 0 1 546 595 602 553 547 596 603 554 +768 5 2 0 1 547 596 603 554 360 367 368 361 +769 5 2 0 1 116 123 124 117 548 597 604 555 +770 5 2 0 1 548 597 604 555 549 598 605 556 +771 5 2 0 1 549 598 605 556 550 599 606 557 +772 5 2 0 1 550 599 606 557 551 600 607 558 +773 5 2 0 1 551 600 607 558 552 601 608 559 +774 5 2 0 1 552 601 608 559 553 602 609 560 +775 5 2 0 1 553 602 609 560 554 603 610 561 +776 5 2 0 1 554 603 610 561 361 368 369 362 +777 5 2 0 1 117 124 125 118 555 604 611 562 +778 5 2 0 1 555 604 611 562 556 605 612 563 +779 5 2 0 1 556 605 612 563 557 606 613 564 +780 5 2 0 1 557 606 613 564 558 607 614 565 +781 5 2 0 1 558 607 614 565 559 608 615 566 +782 5 2 0 1 559 608 615 566 560 609 616 567 +783 5 2 0 1 560 609 616 567 561 610 617 568 +784 5 2 0 1 561 610 617 568 362 369 370 363 +785 5 2 0 1 118 125 126 119 562 611 618 569 +786 5 2 0 1 562 611 618 569 563 612 619 570 +787 5 2 0 1 563 612 619 570 564 613 620 571 +788 5 2 0 1 564 613 620 571 565 614 621 572 +789 5 2 0 1 565 614 621 572 566 615 622 573 +790 5 2 0 1 566 615 622 573 567 616 623 574 +791 5 2 0 1 567 616 623 574 568 617 624 575 +792 5 2 0 1 568 617 624 575 363 370 371 364 +793 5 2 0 1 119 126 127 120 569 618 625 576 +794 5 2 0 1 569 618 625 576 570 619 626 577 +795 5 2 0 1 570 619 626 577 571 620 627 578 +796 5 2 0 1 571 620 627 578 572 621 628 579 +797 5 2 0 1 572 621 628 579 573 622 629 580 +798 5 2 0 1 573 622 629 580 574 623 630 581 +799 5 2 0 1 574 623 630 581 575 624 631 582 +800 5 2 0 1 575 624 631 582 364 371 372 365 +801 5 2 0 1 120 127 25 26 576 625 254 261 +802 5 2 0 1 576 625 254 261 577 626 255 262 +803 5 2 0 1 577 626 255 262 578 627 256 263 +804 5 2 0 1 578 627 256 263 579 628 257 264 +805 5 2 0 1 579 628 257 264 580 629 258 265 +806 5 2 0 1 580 629 258 265 581 630 259 266 +807 5 2 0 1 581 630 259 266 582 631 260 267 +808 5 2 0 1 582 631 260 267 365 372 53 54 +809 5 2 0 1 13 14 128 121 170 177 632 583 +810 5 2 0 1 170 177 632 583 171 178 633 584 +811 5 2 0 1 171 178 633 584 172 179 634 585 +812 5 2 0 1 172 179 634 585 173 180 635 586 +813 5 2 0 1 173 180 635 586 174 181 636 587 +814 5 2 0 1 174 181 636 587 175 182 637 588 +815 5 2 0 1 175 182 637 588 176 183 638 589 +816 5 2 0 1 176 183 638 589 41 42 373 366 +817 5 2 0 1 121 128 129 122 583 632 639 590 +818 5 2 0 1 583 632 639 590 584 633 640 591 +819 5 2 0 1 584 633 640 591 585 634 641 592 +820 5 2 0 1 585 634 641 592 586 635 642 593 +821 5 2 0 1 586 635 642 593 587 636 643 594 +822 5 2 0 1 587 636 643 594 588 637 644 595 +823 5 2 0 1 588 637 644 595 589 638 645 596 +824 5 2 0 1 589 638 645 596 366 373 374 367 +825 5 2 0 1 122 129 130 123 590 639 646 597 +826 5 2 0 1 590 639 646 597 591 640 647 598 +827 5 2 0 1 591 640 647 598 592 641 648 599 +828 5 2 0 1 592 641 648 599 593 642 649 600 +829 5 2 0 1 593 642 649 600 594 643 650 601 +830 5 2 0 1 594 643 650 601 595 644 651 602 +831 5 2 0 1 595 644 651 602 596 645 652 603 +832 5 2 0 1 596 645 652 603 367 374 375 368 +833 5 2 0 1 123 130 131 124 597 646 653 604 +834 5 2 0 1 597 646 653 604 598 647 654 605 +835 5 2 0 1 598 647 654 605 599 648 655 606 +836 5 2 0 1 599 648 655 606 600 649 656 607 +837 5 2 0 1 600 649 656 607 601 650 657 608 +838 5 2 0 1 601 650 657 608 602 651 658 609 +839 5 2 0 1 602 651 658 609 603 652 659 610 +840 5 2 0 1 603 652 659 610 368 375 376 369 +841 5 2 0 1 124 131 132 125 604 653 660 611 +842 5 2 0 1 604 653 660 611 605 654 661 612 +843 5 2 0 1 605 654 661 612 606 655 662 613 +844 5 2 0 1 606 655 662 613 607 656 663 614 +845 5 2 0 1 607 656 663 614 608 657 664 615 +846 5 2 0 1 608 657 664 615 609 658 665 616 +847 5 2 0 1 609 658 665 616 610 659 666 617 +848 5 2 0 1 610 659 666 617 369 376 377 370 +849 5 2 0 1 125 132 133 126 611 660 667 618 +850 5 2 0 1 611 660 667 618 612 661 668 619 +851 5 2 0 1 612 661 668 619 613 662 669 620 +852 5 2 0 1 613 662 669 620 614 663 670 621 +853 5 2 0 1 614 663 670 621 615 664 671 622 +854 5 2 0 1 615 664 671 622 616 665 672 623 +855 5 2 0 1 616 665 672 623 617 666 673 624 +856 5 2 0 1 617 666 673 624 370 377 378 371 +857 5 2 0 1 126 133 134 127 618 667 674 625 +858 5 2 0 1 618 667 674 625 619 668 675 626 +859 5 2 0 1 619 668 675 626 620 669 676 627 +860 5 2 0 1 620 669 676 627 621 670 677 628 +861 5 2 0 1 621 670 677 628 622 671 678 629 +862 5 2 0 1 622 671 678 629 623 672 679 630 +863 5 2 0 1 623 672 679 630 624 673 680 631 +864 5 2 0 1 624 673 680 631 371 378 379 372 +865 5 2 0 1 127 134 24 25 625 674 247 254 +866 5 2 0 1 625 674 247 254 626 675 248 255 +867 5 2 0 1 626 675 248 255 627 676 249 256 +868 5 2 0 1 627 676 249 256 628 677 250 257 +869 5 2 0 1 628 677 250 257 629 678 251 258 +870 5 2 0 1 629 678 251 258 630 679 252 259 +871 5 2 0 1 630 679 252 259 631 680 253 260 +872 5 2 0 1 631 680 253 260 372 379 52 53 +873 5 2 0 1 14 15 135 128 177 184 681 632 +874 5 2 0 1 177 184 681 632 178 185 682 633 +875 5 2 0 1 178 185 682 633 179 186 683 634 +876 5 2 0 1 179 186 683 634 180 187 684 635 +877 5 2 0 1 180 187 684 635 181 188 685 636 +878 5 2 0 1 181 188 685 636 182 189 686 637 +879 5 2 0 1 182 189 686 637 183 190 687 638 +880 5 2 0 1 183 190 687 638 42 43 380 373 +881 5 2 0 1 128 135 136 129 632 681 688 639 +882 5 2 0 1 632 681 688 639 633 682 689 640 +883 5 2 0 1 633 682 689 640 634 683 690 641 +884 5 2 0 1 634 683 690 641 635 684 691 642 +885 5 2 0 1 635 684 691 642 636 685 692 643 +886 5 2 0 1 636 685 692 643 637 686 693 644 +887 5 2 0 1 637 686 693 644 638 687 694 645 +888 5 2 0 1 638 687 694 645 373 380 381 374 +889 5 2 0 1 129 136 137 130 639 688 695 646 +890 5 2 0 1 639 688 695 646 640 689 696 647 +891 5 2 0 1 640 689 696 647 641 690 697 648 +892 5 2 0 1 641 690 697 648 642 691 698 649 +893 5 2 0 1 642 691 698 649 643 692 699 650 +894 5 2 0 1 643 692 699 650 644 693 700 651 +895 5 2 0 1 644 693 700 651 645 694 701 652 +896 5 2 0 1 645 694 701 652 374 381 382 375 +897 5 2 0 1 130 137 138 131 646 695 702 653 +898 5 2 0 1 646 695 702 653 647 696 703 654 +899 5 2 0 1 647 696 703 654 648 697 704 655 +900 5 2 0 1 648 697 704 655 649 698 705 656 +901 5 2 0 1 649 698 705 656 650 699 706 657 +902 5 2 0 1 650 699 706 657 651 700 707 658 +903 5 2 0 1 651 700 707 658 652 701 708 659 +904 5 2 0 1 652 701 708 659 375 382 383 376 +905 5 2 0 1 131 138 139 132 653 702 709 660 +906 5 2 0 1 653 702 709 660 654 703 710 661 +907 5 2 0 1 654 703 710 661 655 704 711 662 +908 5 2 0 1 655 704 711 662 656 705 712 663 +909 5 2 0 1 656 705 712 663 657 706 713 664 +910 5 2 0 1 657 706 713 664 658 707 714 665 +911 5 2 0 1 658 707 714 665 659 708 715 666 +912 5 2 0 1 659 708 715 666 376 383 384 377 +913 5 2 0 1 132 139 140 133 660 709 716 667 +914 5 2 0 1 660 709 716 667 661 710 717 668 +915 5 2 0 1 661 710 717 668 662 711 718 669 +916 5 2 0 1 662 711 718 669 663 712 719 670 +917 5 2 0 1 663 712 719 670 664 713 720 671 +918 5 2 0 1 664 713 720 671 665 714 721 672 +919 5 2 0 1 665 714 721 672 666 715 722 673 +920 5 2 0 1 666 715 722 673 377 384 385 378 +921 5 2 0 1 133 140 141 134 667 716 723 674 +922 5 2 0 1 667 716 723 674 668 717 724 675 +923 5 2 0 1 668 717 724 675 669 718 725 676 +924 5 2 0 1 669 718 725 676 670 719 726 677 +925 5 2 0 1 670 719 726 677 671 720 727 678 +926 5 2 0 1 671 720 727 678 672 721 728 679 +927 5 2 0 1 672 721 728 679 673 722 729 680 +928 5 2 0 1 673 722 729 680 378 385 386 379 +929 5 2 0 1 134 141 23 24 674 723 240 247 +930 5 2 0 1 674 723 240 247 675 724 241 248 +931 5 2 0 1 675 724 241 248 676 725 242 249 +932 5 2 0 1 676 725 242 249 677 726 243 250 +933 5 2 0 1 677 726 243 250 678 727 244 251 +934 5 2 0 1 678 727 244 251 679 728 245 252 +935 5 2 0 1 679 728 245 252 680 729 246 253 +936 5 2 0 1 680 729 246 253 379 386 51 52 +937 5 2 0 1 15 2 16 135 184 72 191 681 +938 5 2 0 1 184 72 191 681 185 73 192 682 +939 5 2 0 1 185 73 192 682 186 74 193 683 +940 5 2 0 1 186 74 193 683 187 75 194 684 +941 5 2 0 1 187 75 194 684 188 76 195 685 +942 5 2 0 1 188 76 195 685 189 77 196 686 +943 5 2 0 1 189 77 196 686 190 78 197 687 +944 5 2 0 1 190 78 197 687 43 6 44 380 +945 5 2 0 1 135 16 17 136 681 191 198 688 +946 5 2 0 1 681 191 198 688 682 192 199 689 +947 5 2 0 1 682 192 199 689 683 193 200 690 +948 5 2 0 1 683 193 200 690 684 194 201 691 +949 5 2 0 1 684 194 201 691 685 195 202 692 +950 5 2 0 1 685 195 202 692 686 196 203 693 +951 5 2 0 1 686 196 203 693 687 197 204 694 +952 5 2 0 1 687 197 204 694 380 44 45 381 +953 5 2 0 1 136 17 18 137 688 198 205 695 +954 5 2 0 1 688 198 205 695 689 199 206 696 +955 5 2 0 1 689 199 206 696 690 200 207 697 +956 5 2 0 1 690 200 207 697 691 201 208 698 +957 5 2 0 1 691 201 208 698 692 202 209 699 +958 5 2 0 1 692 202 209 699 693 203 210 700 +959 5 2 0 1 693 203 210 700 694 204 211 701 +960 5 2 0 1 694 204 211 701 381 45 46 382 +961 5 2 0 1 137 18 19 138 695 205 212 702 +962 5 2 0 1 695 205 212 702 696 206 213 703 +963 5 2 0 1 696 206 213 703 697 207 214 704 +964 5 2 0 1 697 207 214 704 698 208 215 705 +965 5 2 0 1 698 208 215 705 699 209 216 706 +966 5 2 0 1 699 209 216 706 700 210 217 707 +967 5 2 0 1 700 210 217 707 701 211 218 708 +968 5 2 0 1 701 211 218 708 382 46 47 383 +969 5 2 0 1 138 19 20 139 702 212 219 709 +970 5 2 0 1 702 212 219 709 703 213 220 710 +971 5 2 0 1 703 213 220 710 704 214 221 711 +972 5 2 0 1 704 214 221 711 705 215 222 712 +973 5 2 0 1 705 215 222 712 706 216 223 713 +974 5 2 0 1 706 216 223 713 707 217 224 714 +975 5 2 0 1 707 217 224 714 708 218 225 715 +976 5 2 0 1 708 218 225 715 383 47 48 384 +977 5 2 0 1 139 20 21 140 709 219 226 716 +978 5 2 0 1 709 219 226 716 710 220 227 717 +979 5 2 0 1 710 220 227 717 711 221 228 718 +980 5 2 0 1 711 221 228 718 712 222 229 719 +981 5 2 0 1 712 222 229 719 713 223 230 720 +982 5 2 0 1 713 223 230 720 714 224 231 721 +983 5 2 0 1 714 224 231 721 715 225 232 722 +984 5 2 0 1 715 225 232 722 384 48 49 385 +985 5 2 0 1 140 21 22 141 716 226 233 723 +986 5 2 0 1 716 226 233 723 717 227 234 724 +987 5 2 0 1 717 227 234 724 718 228 235 725 +988 5 2 0 1 718 228 235 725 719 229 236 726 +989 5 2 0 1 719 229 236 726 720 230 237 727 +990 5 2 0 1 720 230 237 727 721 231 238 728 +991 5 2 0 1 721 231 238 728 722 232 239 729 +992 5 2 0 1 722 232 239 729 385 49 50 386 +993 5 2 0 1 141 22 3 23 723 233 79 240 +994 5 2 0 1 723 233 79 240 724 234 80 241 +995 5 2 0 1 724 234 80 241 725 235 81 242 +996 5 2 0 1 725 235 81 242 726 236 82 243 +997 5 2 0 1 726 236 82 243 727 237 83 244 +998 5 2 0 1 727 237 83 244 728 238 84 245 +999 5 2 0 1 728 238 84 245 729 239 85 246 +1000 5 2 0 1 729 239 85 246 386 50 7 51 +$EndElements diff --git a/test/sumfact/poisson/cube_hexa_2.msh b/test/sumfact/poisson/cube_hexa_2.msh new file mode 100644 index 0000000000000000000000000000000000000000..c2ecb0be0438733f3617ca4985870605b834aba8 --- /dev/null +++ b/test/sumfact/poisson/cube_hexa_2.msh @@ -0,0 +1,100 @@ +$MeshFormat +2.2 0 8 +$EndMeshFormat +$Nodes +27 +1 0 0 0 +2 1 0 0 +3 1 1 0 +4 0 1 0 +5 0 0 1 +6 1 0 1 +7 1 1 1 +8 0 1 1 +9 0.5 0 0 +10 1 0.5 0 +11 0.5 1 0 +12 0 0.5 0 +13 0.5 0 1 +14 1 0.5 1 +15 0.5 1 1 +16 0 0.5 1 +17 0 0 0.5 +18 1 0 0.5 +19 1 1 0.5 +20 0 1 0.5 +21 0.5 0.5 0 +22 0.5 0 0.5 +23 1 0.5 0.5 +24 0.5 1 0.5 +25 0 0.5 0.5 +26 0.5 0.5 1 +27 0.5 0.5 0.5 +$EndNodes +$Elements +64 +1 15 2 0 1 1 +2 15 2 0 2 2 +3 15 2 0 3 3 +4 15 2 0 4 4 +5 15 2 0 5 5 +6 15 2 0 6 6 +7 15 2 0 10 7 +8 15 2 0 14 8 +9 1 2 0 5 1 9 +10 1 2 0 5 9 2 +11 1 2 0 6 2 10 +12 1 2 0 6 10 3 +13 1 2 0 7 3 11 +14 1 2 0 7 11 4 +15 1 2 0 8 4 12 +16 1 2 0 8 12 1 +17 1 2 0 12 5 13 +18 1 2 0 12 13 6 +19 1 2 0 13 6 14 +20 1 2 0 13 14 7 +21 1 2 0 14 7 15 +22 1 2 0 14 15 8 +23 1 2 0 15 8 16 +24 1 2 0 15 16 5 +25 1 2 0 17 1 17 +26 1 2 0 17 17 5 +27 1 2 0 18 2 18 +28 1 2 0 18 18 6 +29 1 2 0 22 3 19 +30 1 2 0 22 19 7 +31 1 2 0 26 4 20 +32 1 2 0 26 20 8 +33 3 2 0 10 1 9 21 12 +34 3 2 0 10 12 21 11 4 +35 3 2 0 10 9 2 10 21 +36 3 2 0 10 21 10 3 11 +37 3 2 0 19 1 9 22 17 +38 3 2 0 19 17 22 13 5 +39 3 2 0 19 9 2 18 22 +40 3 2 0 19 22 18 6 13 +41 3 2 0 23 2 10 23 18 +42 3 2 0 23 18 23 14 6 +43 3 2 0 23 10 3 19 23 +44 3 2 0 23 23 19 7 14 +45 3 2 0 27 3 11 24 19 +46 3 2 0 27 19 24 15 7 +47 3 2 0 27 11 4 20 24 +48 3 2 0 27 24 20 8 15 +49 3 2 0 31 4 12 25 20 +50 3 2 0 31 20 25 16 8 +51 3 2 0 31 12 1 17 25 +52 3 2 0 31 25 17 5 16 +53 3 2 0 32 5 13 26 16 +54 3 2 0 32 16 26 15 8 +55 3 2 0 32 13 6 14 26 +56 3 2 0 32 26 14 7 15 +57 5 2 0 1 1 9 21 12 17 22 27 25 +58 5 2 0 1 17 22 27 25 5 13 26 16 +59 5 2 0 1 12 21 11 4 25 27 24 20 +60 5 2 0 1 25 27 24 20 16 26 15 8 +61 5 2 0 1 9 2 10 21 22 18 23 27 +62 5 2 0 1 22 18 23 27 13 6 14 26 +63 5 2 0 1 21 10 3 11 27 23 19 24 +64 5 2 0 1 27 23 19 24 26 14 7 15 +$EndElements diff --git a/test/sumfact/poisson/cube_hexa_2_consistent.msh b/test/sumfact/poisson/cube_hexa_2_consistent.msh new file mode 100644 index 0000000000000000000000000000000000000000..9516c8ae4871eb72b8a60b9d5db680ba8decf5b8 --- /dev/null +++ b/test/sumfact/poisson/cube_hexa_2_consistent.msh @@ -0,0 +1,44 @@ +$MeshFormat +2.0 0 8 +$EndMeshFormat +$Nodes +27 +1 0.5 0 0 +2 0.5 0.5 0 +3 0 0.5 0 +4 0 0 0 +5 0.5 0 0.5 +7 0 0.5 0.5 +8 0 0 0.5 +9 0.5 0 1 +10 0.5 0.5 1 +11 0 0.5 1 +12 0 0 1 +13 0.5 1 0 +14 0 1 0 +15 0.5 1 0.5 +16 0 1 0.5 +17 0.5 1 1 +18 0 1 1 +19 1 0 0 +20 1 0.5 0 +21 1 0 0.5 +22 1 0.5 0.5 +23 1 0 1 +24 1 0.5 1 +25 1 1 0 +26 1 1 0.5 +27 1 1 1 +6 0.5 0.5 0.5 +$EndNodes +$Elements +8 +1 5 0 2 3 4 1 6 7 8 5 +2 5 0 6 7 8 5 10 11 12 9 +3 5 0 2 13 14 3 6 15 16 7 +4 5 0 6 15 16 7 10 17 18 11 +5 5 0 2 1 19 20 6 5 21 22 +6 5 0 6 5 21 22 10 9 23 24 +7 5 0 2 20 25 13 6 22 26 15 +8 5 0 6 22 26 15 10 24 27 17 +$EndElements diff --git a/test/sumfact/poisson/cube_hexa_consistent.msh b/test/sumfact/poisson/cube_hexa_consistent.msh new file mode 100644 index 0000000000000000000000000000000000000000..b14f54f4a303156b44dc1e929c58c35c9a887637 --- /dev/null +++ b/test/sumfact/poisson/cube_hexa_consistent.msh @@ -0,0 +1,1250 @@ +$MeshFormat +2.0 0 8 +$EndMeshFormat +$Nodes +729 +1 0 0 0 +2 0.125 0 0 +3 0.125 0.125 0 +4 0 0.125 0 +5 0.125 0.25 0 +6 0 0.25 0 +7 0.125 0.375 0 +8 0 0.375 0 +9 0.125 0.5 0 +10 0 0.5 0 +11 0.125 0.625 0 +12 0 0.625 0 +13 0.125 0.75 0 +14 0 0.75 0 +15 0.125 0.875 0 +16 0 0.875 0 +17 0.125 1 0 +18 0 1 0 +19 0.25 0 0 +20 0.25 0.125 0 +21 0.25 0.25 0 +22 0.25 0.375 0 +23 0.25 0.5 0 +24 0.25 0.625 0 +25 0.25 0.75 0 +26 0.25 0.875 0 +27 0.25 1 0 +28 0.375 0 0 +29 0.375 0.125 0 +30 0.375 0.25 0 +31 0.375 0.375 0 +32 0.375 0.5 0 +33 0.375 0.625 0 +34 0.375 0.75 0 +35 0.375 0.875 0 +36 0.375 1 0 +37 0.5 0 0 +38 0.5 0.125 0 +39 0.5 0.25 0 +40 0.5 0.375 0 +41 0.5 0.5 0 +42 0.5 0.625 0 +43 0.5 0.75 0 +44 0.5 0.875 0 +45 0.5 1 0 +46 0.625 0 0 +47 0.625 0.125 0 +48 0.625 0.25 0 +49 0.625 0.375 0 +50 0.625 0.5 0 +51 0.625 0.625 0 +52 0.625 0.75 0 +53 0.625 0.875 0 +54 0.625 1 0 +55 0.75 0 0 +56 0.75 0.125 0 +57 0.75 0.25 0 +58 0.75 0.375 0 +59 0.75 0.5 0 +60 0.75 0.625 0 +61 0.75 0.75 0 +62 0.75 0.875 0 +63 0.75 1 0 +64 0.875 0 0 +65 0.875 0.125 0 +66 0.875 0.25 0 +67 0.875 0.375 0 +68 0.875 0.5 0 +69 0.875 0.625 0 +70 0.875 0.75 0 +71 0.875 0.875 0 +72 0.875 1 0 +73 1 0 0 +74 1 0.125 0 +75 1 0.25 0 +76 1 0.375 0 +77 1 0.5 0 +78 1 0.625 0 +79 1 0.75 0 +80 1 0.875 0 +81 1 1 0 +82 0.125 0 0.125 +83 0 0 0.125 +84 0.125 0 0.25 +85 0 0 0.25 +86 0.125 0 0.375 +87 0 0 0.375 +88 0.125 0 0.5 +89 0 0 0.5 +90 0.125 0 0.625 +91 0 0 0.625 +92 0.125 0 0.75 +93 0 0 0.75 +94 0.125 0 0.875 +95 0 0 0.875 +96 0.125 0 1 +97 0 0 1 +98 0.25 0 0.125 +99 0.25 0 0.25 +100 0.25 0 0.375 +101 0.25 0 0.5 +102 0.25 0 0.625 +103 0.25 0 0.75 +104 0.25 0 0.875 +105 0.25 0 1 +106 0.375 0 0.125 +107 0.375 0 0.25 +108 0.375 0 0.375 +109 0.375 0 0.5 +110 0.375 0 0.625 +111 0.375 0 0.75 +112 0.375 0 0.875 +113 0.375 0 1 +114 0.5 0 0.125 +115 0.5 0 0.25 +116 0.5 0 0.375 +117 0.5 0 0.5 +118 0.5 0 0.625 +119 0.5 0 0.75 +120 0.5 0 0.875 +121 0.5 0 1 +122 0.625 0 0.125 +123 0.625 0 0.25 +124 0.625 0 0.375 +125 0.625 0 0.5 +126 0.625 0 0.625 +127 0.625 0 0.75 +128 0.625 0 0.875 +129 0.625 0 1 +130 0.75 0 0.125 +131 0.75 0 0.25 +132 0.75 0 0.375 +133 0.75 0 0.5 +134 0.75 0 0.625 +135 0.75 0 0.75 +136 0.75 0 0.875 +137 0.75 0 1 +138 0.875 0 0.125 +139 0.875 0 0.25 +140 0.875 0 0.375 +141 0.875 0 0.5 +142 0.875 0 0.625 +143 0.875 0 0.75 +144 0.875 0 0.875 +145 0.875 0 1 +146 1 0 0.125 +147 1 0 0.25 +148 1 0 0.375 +149 1 0 0.5 +150 1 0 0.625 +151 1 0 0.75 +152 1 0 0.875 +153 1 0 1 +154 1 0.125 0.125 +155 1 0.125 0.25 +156 1 0.125 0.375 +157 1 0.125 0.5 +158 1 0.125 0.625 +159 1 0.125 0.75 +160 1 0.125 0.875 +161 1 0.125 1 +162 1 0.25 0.125 +163 1 0.25 0.25 +164 1 0.25 0.375 +165 1 0.25 0.5 +166 1 0.25 0.625 +167 1 0.25 0.75 +168 1 0.25 0.875 +169 1 0.25 1 +170 1 0.375 0.125 +171 1 0.375 0.25 +172 1 0.375 0.375 +173 1 0.375 0.5 +174 1 0.375 0.625 +175 1 0.375 0.75 +176 1 0.375 0.875 +177 1 0.375 1 +178 1 0.5 0.125 +179 1 0.5 0.25 +180 1 0.5 0.375 +181 1 0.5 0.5 +182 1 0.5 0.625 +183 1 0.5 0.75 +184 1 0.5 0.875 +185 1 0.5 1 +186 1 0.625 0.125 +187 1 0.625 0.25 +188 1 0.625 0.375 +189 1 0.625 0.5 +190 1 0.625 0.625 +191 1 0.625 0.75 +192 1 0.625 0.875 +193 1 0.625 1 +194 1 0.75 0.125 +195 1 0.75 0.25 +196 1 0.75 0.375 +197 1 0.75 0.5 +198 1 0.75 0.625 +199 1 0.75 0.75 +200 1 0.75 0.875 +201 1 0.75 1 +202 1 0.875 0.125 +203 1 0.875 0.25 +204 1 0.875 0.375 +205 1 0.875 0.5 +206 1 0.875 0.625 +207 1 0.875 0.75 +208 1 0.875 0.875 +209 1 0.875 1 +210 1 1 0.125 +211 1 1 0.25 +212 1 1 0.375 +213 1 1 0.5 +214 1 1 0.625 +215 1 1 0.75 +216 1 1 0.875 +217 1 1 1 +218 0.875 1 0.125 +219 0.875 1 0.25 +220 0.875 1 0.375 +221 0.875 1 0.5 +222 0.875 1 0.625 +223 0.875 1 0.75 +224 0.875 1 0.875 +225 0.875 1 1 +226 0.75 1 0.125 +227 0.75 1 0.25 +228 0.75 1 0.375 +229 0.75 1 0.5 +230 0.75 1 0.625 +231 0.75 1 0.75 +232 0.75 1 0.875 +233 0.75 1 1 +234 0.625 1 0.125 +235 0.625 1 0.25 +236 0.625 1 0.375 +237 0.625 1 0.5 +238 0.625 1 0.625 +239 0.625 1 0.75 +240 0.625 1 0.875 +241 0.625 1 1 +242 0.5 1 0.125 +243 0.5 1 0.25 +244 0.5 1 0.375 +245 0.5 1 0.5 +246 0.5 1 0.625 +247 0.5 1 0.75 +248 0.5 1 0.875 +249 0.5 1 1 +250 0.375 1 0.125 +251 0.375 1 0.25 +252 0.375 1 0.375 +253 0.375 1 0.5 +254 0.375 1 0.625 +255 0.375 1 0.75 +256 0.375 1 0.875 +257 0.375 1 1 +258 0.25 1 0.125 +259 0.25 1 0.25 +260 0.25 1 0.375 +261 0.25 1 0.5 +262 0.25 1 0.625 +263 0.25 1 0.75 +264 0.25 1 0.875 +265 0.25 1 1 +266 0.125 1 0.125 +267 0.125 1 0.25 +268 0.125 1 0.375 +269 0.125 1 0.5 +270 0.125 1 0.625 +271 0.125 1 0.75 +272 0.125 1 0.875 +273 0.125 1 1 +274 0 1 0.125 +275 0 1 0.25 +276 0 1 0.375 +277 0 1 0.5 +278 0 1 0.625 +279 0 1 0.75 +280 0 1 0.875 +281 0 1 1 +282 0 0.875 0.125 +283 0 0.875 0.25 +284 0 0.875 0.375 +285 0 0.875 0.5 +286 0 0.875 0.625 +287 0 0.875 0.75 +288 0 0.875 0.875 +289 0 0.875 1 +290 0 0.75 0.125 +291 0 0.75 0.25 +292 0 0.75 0.375 +293 0 0.75 0.5 +294 0 0.75 0.625 +295 0 0.75 0.75 +296 0 0.75 0.875 +297 0 0.75 1 +298 0 0.625 0.125 +299 0 0.625 0.25 +300 0 0.625 0.375 +301 0 0.625 0.5 +302 0 0.625 0.625 +303 0 0.625 0.75 +304 0 0.625 0.875 +305 0 0.625 1 +306 0 0.5 0.125 +307 0 0.5 0.25 +308 0 0.5 0.375 +309 0 0.5 0.5 +310 0 0.5 0.625 +311 0 0.5 0.75 +312 0 0.5 0.875 +313 0 0.5 1 +314 0 0.375 0.125 +315 0 0.375 0.25 +316 0 0.375 0.375 +317 0 0.375 0.5 +318 0 0.375 0.625 +319 0 0.375 0.75 +320 0 0.375 0.875 +321 0 0.375 1 +322 0 0.25 0.125 +323 0 0.25 0.25 +324 0 0.25 0.375 +325 0 0.25 0.5 +326 0 0.25 0.625 +327 0 0.25 0.75 +328 0 0.25 0.875 +329 0 0.25 1 +330 0 0.125 0.125 +331 0 0.125 0.25 +332 0 0.125 0.375 +333 0 0.125 0.5 +334 0 0.125 0.625 +335 0 0.125 0.75 +336 0 0.125 0.875 +337 0 0.125 1 +338 0.125 0.125 1 +339 0.125 0.25 1 +340 0.125 0.375 1 +341 0.125 0.5 1 +342 0.125 0.625 1 +343 0.125 0.75 1 +344 0.125 0.875 1 +345 0.25 0.125 1 +346 0.25 0.25 1 +347 0.25 0.375 1 +348 0.25 0.5 1 +349 0.25 0.625 1 +350 0.25 0.75 1 +351 0.25 0.875 1 +352 0.375 0.125 1 +353 0.375 0.25 1 +354 0.375 0.375 1 +355 0.375 0.5 1 +356 0.375 0.625 1 +357 0.375 0.75 1 +358 0.375 0.875 1 +359 0.5 0.125 1 +360 0.5 0.25 1 +361 0.5 0.375 1 +362 0.5 0.5 1 +363 0.5 0.625 1 +364 0.5 0.75 1 +365 0.5 0.875 1 +366 0.625 0.125 1 +367 0.625 0.25 1 +368 0.625 0.375 1 +369 0.625 0.5 1 +370 0.625 0.625 1 +371 0.625 0.75 1 +372 0.625 0.875 1 +373 0.75 0.125 1 +374 0.75 0.25 1 +375 0.75 0.375 1 +376 0.75 0.5 1 +377 0.75 0.625 1 +378 0.75 0.75 1 +379 0.75 0.875 1 +380 0.875 0.125 1 +381 0.875 0.25 1 +382 0.875 0.375 1 +383 0.875 0.5 1 +384 0.875 0.625 1 +385 0.875 0.75 1 +386 0.875 0.875 1 +387 0.125 0.125 0.125 +388 0.125 0.125 0.25 +389 0.125 0.125 0.375 +390 0.125 0.125 0.5 +391 0.125 0.125 0.625 +392 0.125 0.125 0.75 +393 0.125 0.125 0.875 +394 0.125 0.25 0.125 +395 0.125 0.25 0.25 +396 0.125 0.25 0.375 +397 0.125 0.25 0.5 +398 0.125 0.25 0.625 +399 0.125 0.25 0.75 +400 0.125 0.25 0.875 +401 0.125 0.375 0.125 +402 0.125 0.375 0.25 +403 0.125 0.375 0.375 +404 0.125 0.375 0.5 +405 0.125 0.375 0.625 +406 0.125 0.375 0.75 +407 0.125 0.375 0.875 +408 0.125 0.5 0.125 +409 0.125 0.5 0.25 +410 0.125 0.5 0.375 +411 0.125 0.5 0.5 +412 0.125 0.5 0.625 +413 0.125 0.5 0.75 +414 0.125 0.5 0.875 +415 0.125 0.625 0.125 +416 0.125 0.625 0.25 +417 0.125 0.625 0.375 +418 0.125 0.625 0.5 +419 0.125 0.625 0.625 +420 0.125 0.625 0.75 +421 0.125 0.625 0.875 +422 0.125 0.75 0.125 +423 0.125 0.75 0.25 +424 0.125 0.75 0.375 +425 0.125 0.75 0.5 +426 0.125 0.75 0.625 +427 0.125 0.75 0.75 +428 0.125 0.75 0.875 +429 0.125 0.875 0.125 +430 0.125 0.875 0.25 +431 0.125 0.875 0.375 +432 0.125 0.875 0.5 +433 0.125 0.875 0.625 +434 0.125 0.875 0.75 +435 0.125 0.875 0.875 +436 0.25 0.125 0.125 +437 0.25 0.125 0.25 +438 0.25 0.125 0.375 +439 0.25 0.125 0.5 +440 0.25 0.125 0.625 +441 0.25 0.125 0.75 +442 0.25 0.125 0.875 +443 0.25 0.25 0.125 +444 0.25 0.25 0.25 +445 0.25 0.25 0.375 +446 0.25 0.25 0.5 +447 0.25 0.25 0.625 +448 0.25 0.25 0.75 +449 0.25 0.25 0.875 +450 0.25 0.375 0.125 +451 0.25 0.375 0.25 +452 0.25 0.375 0.375 +453 0.25 0.375 0.5 +454 0.25 0.375 0.625 +455 0.25 0.375 0.75 +456 0.25 0.375 0.875 +457 0.25 0.5 0.125 +458 0.25 0.5 0.25 +459 0.25 0.5 0.375 +460 0.25 0.5 0.5 +461 0.25 0.5 0.625 +462 0.25 0.5 0.75 +463 0.25 0.5 0.875 +464 0.25 0.625 0.125 +465 0.25 0.625 0.25 +466 0.25 0.625 0.375 +467 0.25 0.625 0.5 +468 0.25 0.625 0.625 +469 0.25 0.625 0.75 +470 0.25 0.625 0.875 +471 0.25 0.75 0.125 +472 0.25 0.75 0.25 +473 0.25 0.75 0.375 +474 0.25 0.75 0.5 +475 0.25 0.75 0.625 +476 0.25 0.75 0.75 +477 0.25 0.75 0.875 +478 0.25 0.875 0.125 +479 0.25 0.875 0.25 +480 0.25 0.875 0.375 +481 0.25 0.875 0.5 +482 0.25 0.875 0.625 +483 0.25 0.875 0.75 +484 0.25 0.875 0.875 +485 0.375 0.125 0.125 +486 0.375 0.125 0.25 +487 0.375 0.125 0.375 +488 0.375 0.125 0.5 +489 0.375 0.125 0.625 +490 0.375 0.125 0.75 +491 0.375 0.125 0.875 +492 0.375 0.25 0.125 +493 0.375 0.25 0.25 +494 0.375 0.25 0.375 +495 0.375 0.25 0.5 +496 0.375 0.25 0.625 +497 0.375 0.25 0.75 +498 0.375 0.25 0.875 +499 0.375 0.375 0.125 +500 0.375 0.375 0.25 +501 0.375 0.375 0.375 +502 0.375 0.375 0.5 +503 0.375 0.375 0.625 +504 0.375 0.375 0.75 +505 0.375 0.375 0.875 +506 0.375 0.5 0.125 +507 0.375 0.5 0.25 +508 0.375 0.5 0.375 +509 0.375 0.5 0.5 +510 0.375 0.5 0.625 +511 0.375 0.5 0.75 +512 0.375 0.5 0.875 +513 0.375 0.625 0.125 +514 0.375 0.625 0.25 +515 0.375 0.625 0.375 +516 0.375 0.625 0.5 +517 0.375 0.625 0.625 +518 0.375 0.625 0.75 +519 0.375 0.625 0.875 +520 0.375 0.75 0.125 +521 0.375 0.75 0.25 +522 0.375 0.75 0.375 +523 0.375 0.75 0.5 +524 0.375 0.75 0.625 +525 0.375 0.75 0.75 +526 0.375 0.75 0.875 +527 0.375 0.875 0.125 +528 0.375 0.875 0.25 +529 0.375 0.875 0.375 +530 0.375 0.875 0.5 +531 0.375 0.875 0.625 +532 0.375 0.875 0.75 +533 0.375 0.875 0.875 +534 0.5 0.125 0.125 +535 0.5 0.125 0.25 +536 0.5 0.125 0.375 +537 0.5 0.125 0.5 +538 0.5 0.125 0.625 +539 0.5 0.125 0.75 +540 0.5 0.125 0.875 +541 0.5 0.25 0.125 +542 0.5 0.25 0.25 +543 0.5 0.25 0.375 +544 0.5 0.25 0.5 +545 0.5 0.25 0.625 +546 0.5 0.25 0.75 +547 0.5 0.25 0.875 +548 0.5 0.375 0.125 +549 0.5 0.375 0.25 +550 0.5 0.375 0.375 +551 0.5 0.375 0.5 +552 0.5 0.375 0.625 +553 0.5 0.375 0.75 +554 0.5 0.375 0.875 +555 0.5 0.5 0.125 +556 0.5 0.5 0.25 +557 0.5 0.5 0.375 +558 0.5 0.5 0.5 +559 0.5 0.5 0.625 +560 0.5 0.5 0.75 +561 0.5 0.5 0.875 +562 0.5 0.625 0.125 +563 0.5 0.625 0.25 +564 0.5 0.625 0.375 +565 0.5 0.625 0.5 +566 0.5 0.625 0.625 +567 0.5 0.625 0.75 +568 0.5 0.625 0.875 +569 0.5 0.75 0.125 +570 0.5 0.75 0.25 +571 0.5 0.75 0.375 +572 0.5 0.75 0.5 +573 0.5 0.75 0.625 +574 0.5 0.75 0.75 +575 0.5 0.75 0.875 +576 0.5 0.875 0.125 +577 0.5 0.875 0.25 +578 0.5 0.875 0.375 +579 0.5 0.875 0.5 +580 0.5 0.875 0.625 +581 0.5 0.875 0.75 +582 0.5 0.875 0.875 +583 0.625 0.125 0.125 +584 0.625 0.125 0.25 +585 0.625 0.125 0.375 +586 0.625 0.125 0.5 +587 0.625 0.125 0.625 +588 0.625 0.125 0.75 +589 0.625 0.125 0.875 +590 0.625 0.25 0.125 +591 0.625 0.25 0.25 +592 0.625 0.25 0.375 +593 0.625 0.25 0.5 +594 0.625 0.25 0.625 +595 0.625 0.25 0.75 +596 0.625 0.25 0.875 +597 0.625 0.375 0.125 +598 0.625 0.375 0.25 +599 0.625 0.375 0.375 +600 0.625 0.375 0.5 +601 0.625 0.375 0.625 +602 0.625 0.375 0.75 +603 0.625 0.375 0.875 +604 0.625 0.5 0.125 +605 0.625 0.5 0.25 +606 0.625 0.5 0.375 +607 0.625 0.5 0.5 +608 0.625 0.5 0.625 +609 0.625 0.5 0.75 +610 0.625 0.5 0.875 +611 0.625 0.625 0.125 +612 0.625 0.625 0.25 +613 0.625 0.625 0.375 +614 0.625 0.625 0.5 +615 0.625 0.625 0.625 +616 0.625 0.625 0.75 +617 0.625 0.625 0.875 +618 0.625 0.75 0.125 +619 0.625 0.75 0.25 +620 0.625 0.75 0.375 +621 0.625 0.75 0.5 +622 0.625 0.75 0.625 +623 0.625 0.75 0.75 +624 0.625 0.75 0.875 +625 0.625 0.875 0.125 +626 0.625 0.875 0.25 +627 0.625 0.875 0.375 +628 0.625 0.875 0.5 +629 0.625 0.875 0.625 +630 0.625 0.875 0.75 +631 0.625 0.875 0.875 +632 0.75 0.125 0.125 +633 0.75 0.125 0.25 +634 0.75 0.125 0.375 +635 0.75 0.125 0.5 +636 0.75 0.125 0.625 +637 0.75 0.125 0.75 +638 0.75 0.125 0.875 +639 0.75 0.25 0.125 +640 0.75 0.25 0.25 +641 0.75 0.25 0.375 +642 0.75 0.25 0.5 +643 0.75 0.25 0.625 +644 0.75 0.25 0.75 +645 0.75 0.25 0.875 +646 0.75 0.375 0.125 +647 0.75 0.375 0.25 +648 0.75 0.375 0.375 +649 0.75 0.375 0.5 +650 0.75 0.375 0.625 +651 0.75 0.375 0.75 +652 0.75 0.375 0.875 +653 0.75 0.5 0.125 +654 0.75 0.5 0.25 +655 0.75 0.5 0.375 +656 0.75 0.5 0.5 +657 0.75 0.5 0.625 +658 0.75 0.5 0.75 +659 0.75 0.5 0.875 +660 0.75 0.625 0.125 +661 0.75 0.625 0.25 +662 0.75 0.625 0.375 +663 0.75 0.625 0.5 +664 0.75 0.625 0.625 +665 0.75 0.625 0.75 +666 0.75 0.625 0.875 +667 0.75 0.75 0.125 +668 0.75 0.75 0.25 +669 0.75 0.75 0.375 +670 0.75 0.75 0.5 +671 0.75 0.75 0.625 +672 0.75 0.75 0.75 +673 0.75 0.75 0.875 +674 0.75 0.875 0.125 +675 0.75 0.875 0.25 +676 0.75 0.875 0.375 +677 0.75 0.875 0.5 +678 0.75 0.875 0.625 +679 0.75 0.875 0.75 +680 0.75 0.875 0.875 +681 0.875 0.125 0.125 +682 0.875 0.125 0.25 +683 0.875 0.125 0.375 +684 0.875 0.125 0.5 +685 0.875 0.125 0.625 +686 0.875 0.125 0.75 +687 0.875 0.125 0.875 +688 0.875 0.25 0.125 +689 0.875 0.25 0.25 +690 0.875 0.25 0.375 +691 0.875 0.25 0.5 +692 0.875 0.25 0.625 +693 0.875 0.25 0.75 +694 0.875 0.25 0.875 +695 0.875 0.375 0.125 +696 0.875 0.375 0.25 +697 0.875 0.375 0.375 +698 0.875 0.375 0.5 +699 0.875 0.375 0.625 +700 0.875 0.375 0.75 +701 0.875 0.375 0.875 +702 0.875 0.5 0.125 +703 0.875 0.5 0.25 +704 0.875 0.5 0.375 +705 0.875 0.5 0.5 +706 0.875 0.5 0.625 +707 0.875 0.5 0.75 +708 0.875 0.5 0.875 +709 0.875 0.625 0.125 +710 0.875 0.625 0.25 +711 0.875 0.625 0.375 +712 0.875 0.625 0.5 +713 0.875 0.625 0.625 +714 0.875 0.625 0.75 +715 0.875 0.625 0.875 +716 0.875 0.75 0.125 +717 0.875 0.75 0.25 +718 0.875 0.75 0.375 +719 0.875 0.75 0.5 +720 0.875 0.75 0.625 +721 0.875 0.75 0.75 +722 0.875 0.75 0.875 +723 0.875 0.875 0.125 +724 0.875 0.875 0.25 +725 0.875 0.875 0.375 +726 0.875 0.875 0.5 +727 0.875 0.875 0.625 +728 0.875 0.875 0.75 +729 0.875 0.875 0.875 +$EndNodes +$Elements +512 +1 5 0 2 3 4 1 82 387 330 83 +2 5 0 82 387 330 83 84 388 331 85 +3 5 0 84 388 331 85 86 389 332 87 +4 5 0 86 389 332 87 88 390 333 89 +5 5 0 88 390 333 89 90 391 334 91 +6 5 0 90 391 334 91 92 392 335 93 +7 5 0 92 392 335 93 94 393 336 95 +8 5 0 94 393 336 95 96 338 337 97 +9 5 0 3 5 6 4 387 394 322 330 +10 5 0 387 394 322 330 388 395 323 331 +11 5 0 388 395 323 331 389 396 324 332 +12 5 0 389 396 324 332 390 397 325 333 +13 5 0 390 397 325 333 391 398 326 334 +14 5 0 391 398 326 334 392 399 327 335 +15 5 0 392 399 327 335 393 400 328 336 +16 5 0 393 400 328 336 338 339 329 337 +17 5 0 5 7 8 6 394 401 314 322 +18 5 0 394 401 314 322 395 402 315 323 +19 5 0 395 402 315 323 396 403 316 324 +20 5 0 396 403 316 324 397 404 317 325 +21 5 0 397 404 317 325 398 405 318 326 +22 5 0 398 405 318 326 399 406 319 327 +23 5 0 399 406 319 327 400 407 320 328 +24 5 0 400 407 320 328 339 340 321 329 +25 5 0 7 9 10 8 401 408 306 314 +26 5 0 401 408 306 314 402 409 307 315 +27 5 0 402 409 307 315 403 410 308 316 +28 5 0 403 410 308 316 404 411 309 317 +29 5 0 404 411 309 317 405 412 310 318 +30 5 0 405 412 310 318 406 413 311 319 +31 5 0 406 413 311 319 407 414 312 320 +32 5 0 407 414 312 320 340 341 313 321 +33 5 0 9 11 12 10 408 415 298 306 +34 5 0 408 415 298 306 409 416 299 307 +35 5 0 409 416 299 307 410 417 300 308 +36 5 0 410 417 300 308 411 418 301 309 +37 5 0 411 418 301 309 412 419 302 310 +38 5 0 412 419 302 310 413 420 303 311 +39 5 0 413 420 303 311 414 421 304 312 +40 5 0 414 421 304 312 341 342 305 313 +41 5 0 11 13 14 12 415 422 290 298 +42 5 0 415 422 290 298 416 423 291 299 +43 5 0 416 423 291 299 417 424 292 300 +44 5 0 417 424 292 300 418 425 293 301 +45 5 0 418 425 293 301 419 426 294 302 +46 5 0 419 426 294 302 420 427 295 303 +47 5 0 420 427 295 303 421 428 296 304 +48 5 0 421 428 296 304 342 343 297 305 +49 5 0 13 15 16 14 422 429 282 290 +50 5 0 422 429 282 290 423 430 283 291 +51 5 0 423 430 283 291 424 431 284 292 +52 5 0 424 431 284 292 425 432 285 293 +53 5 0 425 432 285 293 426 433 286 294 +54 5 0 426 433 286 294 427 434 287 295 +55 5 0 427 434 287 295 428 435 288 296 +56 5 0 428 435 288 296 343 344 289 297 +57 5 0 15 17 18 16 429 266 274 282 +58 5 0 429 266 274 282 430 267 275 283 +59 5 0 430 267 275 283 431 268 276 284 +60 5 0 431 268 276 284 432 269 277 285 +61 5 0 432 269 277 285 433 270 278 286 +62 5 0 433 270 278 286 434 271 279 287 +63 5 0 434 271 279 287 435 272 280 288 +64 5 0 435 272 280 288 344 273 281 289 +65 5 0 2 19 20 3 82 98 436 387 +66 5 0 82 98 436 387 84 99 437 388 +67 5 0 84 99 437 388 86 100 438 389 +68 5 0 86 100 438 389 88 101 439 390 +69 5 0 88 101 439 390 90 102 440 391 +70 5 0 90 102 440 391 92 103 441 392 +71 5 0 92 103 441 392 94 104 442 393 +72 5 0 94 104 442 393 96 105 345 338 +73 5 0 3 20 21 5 387 436 443 394 +74 5 0 387 436 443 394 388 437 444 395 +75 5 0 388 437 444 395 389 438 445 396 +76 5 0 389 438 445 396 390 439 446 397 +77 5 0 390 439 446 397 391 440 447 398 +78 5 0 391 440 447 398 392 441 448 399 +79 5 0 392 441 448 399 393 442 449 400 +80 5 0 393 442 449 400 338 345 346 339 +81 5 0 5 21 22 7 394 443 450 401 +82 5 0 394 443 450 401 395 444 451 402 +83 5 0 395 444 451 402 396 445 452 403 +84 5 0 396 445 452 403 397 446 453 404 +85 5 0 397 446 453 404 398 447 454 405 +86 5 0 398 447 454 405 399 448 455 406 +87 5 0 399 448 455 406 400 449 456 407 +88 5 0 400 449 456 407 339 346 347 340 +89 5 0 7 22 23 9 401 450 457 408 +90 5 0 401 450 457 408 402 451 458 409 +91 5 0 402 451 458 409 403 452 459 410 +92 5 0 403 452 459 410 404 453 460 411 +93 5 0 404 453 460 411 405 454 461 412 +94 5 0 405 454 461 412 406 455 462 413 +95 5 0 406 455 462 413 407 456 463 414 +96 5 0 407 456 463 414 340 347 348 341 +97 5 0 9 23 24 11 408 457 464 415 +98 5 0 408 457 464 415 409 458 465 416 +99 5 0 409 458 465 416 410 459 466 417 +100 5 0 410 459 466 417 411 460 467 418 +101 5 0 411 460 467 418 412 461 468 419 +102 5 0 412 461 468 419 413 462 469 420 +103 5 0 413 462 469 420 414 463 470 421 +104 5 0 414 463 470 421 341 348 349 342 +105 5 0 11 24 25 13 415 464 471 422 +106 5 0 415 464 471 422 416 465 472 423 +107 5 0 416 465 472 423 417 466 473 424 +108 5 0 417 466 473 424 418 467 474 425 +109 5 0 418 467 474 425 419 468 475 426 +110 5 0 419 468 475 426 420 469 476 427 +111 5 0 420 469 476 427 421 470 477 428 +112 5 0 421 470 477 428 342 349 350 343 +113 5 0 13 25 26 15 422 471 478 429 +114 5 0 422 471 478 429 423 472 479 430 +115 5 0 423 472 479 430 424 473 480 431 +116 5 0 424 473 480 431 425 474 481 432 +117 5 0 425 474 481 432 426 475 482 433 +118 5 0 426 475 482 433 427 476 483 434 +119 5 0 427 476 483 434 428 477 484 435 +120 5 0 428 477 484 435 343 350 351 344 +121 5 0 15 26 27 17 429 478 258 266 +122 5 0 429 478 258 266 430 479 259 267 +123 5 0 430 479 259 267 431 480 260 268 +124 5 0 431 480 260 268 432 481 261 269 +125 5 0 432 481 261 269 433 482 262 270 +126 5 0 433 482 262 270 434 483 263 271 +127 5 0 434 483 263 271 435 484 264 272 +128 5 0 435 484 264 272 344 351 265 273 +129 5 0 19 28 29 20 98 106 485 436 +130 5 0 98 106 485 436 99 107 486 437 +131 5 0 99 107 486 437 100 108 487 438 +132 5 0 100 108 487 438 101 109 488 439 +133 5 0 101 109 488 439 102 110 489 440 +134 5 0 102 110 489 440 103 111 490 441 +135 5 0 103 111 490 441 104 112 491 442 +136 5 0 104 112 491 442 105 113 352 345 +137 5 0 20 29 30 21 436 485 492 443 +138 5 0 436 485 492 443 437 486 493 444 +139 5 0 437 486 493 444 438 487 494 445 +140 5 0 438 487 494 445 439 488 495 446 +141 5 0 439 488 495 446 440 489 496 447 +142 5 0 440 489 496 447 441 490 497 448 +143 5 0 441 490 497 448 442 491 498 449 +144 5 0 442 491 498 449 345 352 353 346 +145 5 0 21 30 31 22 443 492 499 450 +146 5 0 443 492 499 450 444 493 500 451 +147 5 0 444 493 500 451 445 494 501 452 +148 5 0 445 494 501 452 446 495 502 453 +149 5 0 446 495 502 453 447 496 503 454 +150 5 0 447 496 503 454 448 497 504 455 +151 5 0 448 497 504 455 449 498 505 456 +152 5 0 449 498 505 456 346 353 354 347 +153 5 0 22 31 32 23 450 499 506 457 +154 5 0 450 499 506 457 451 500 507 458 +155 5 0 451 500 507 458 452 501 508 459 +156 5 0 452 501 508 459 453 502 509 460 +157 5 0 453 502 509 460 454 503 510 461 +158 5 0 454 503 510 461 455 504 511 462 +159 5 0 455 504 511 462 456 505 512 463 +160 5 0 456 505 512 463 347 354 355 348 +161 5 0 23 32 33 24 457 506 513 464 +162 5 0 457 506 513 464 458 507 514 465 +163 5 0 458 507 514 465 459 508 515 466 +164 5 0 459 508 515 466 460 509 516 467 +165 5 0 460 509 516 467 461 510 517 468 +166 5 0 461 510 517 468 462 511 518 469 +167 5 0 462 511 518 469 463 512 519 470 +168 5 0 463 512 519 470 348 355 356 349 +169 5 0 24 33 34 25 464 513 520 471 +170 5 0 464 513 520 471 465 514 521 472 +171 5 0 465 514 521 472 466 515 522 473 +172 5 0 466 515 522 473 467 516 523 474 +173 5 0 467 516 523 474 468 517 524 475 +174 5 0 468 517 524 475 469 518 525 476 +175 5 0 469 518 525 476 470 519 526 477 +176 5 0 470 519 526 477 349 356 357 350 +177 5 0 25 34 35 26 471 520 527 478 +178 5 0 471 520 527 478 472 521 528 479 +179 5 0 472 521 528 479 473 522 529 480 +180 5 0 473 522 529 480 474 523 530 481 +181 5 0 474 523 530 481 475 524 531 482 +182 5 0 475 524 531 482 476 525 532 483 +183 5 0 476 525 532 483 477 526 533 484 +184 5 0 477 526 533 484 350 357 358 351 +185 5 0 26 35 36 27 478 527 250 258 +186 5 0 478 527 250 258 479 528 251 259 +187 5 0 479 528 251 259 480 529 252 260 +188 5 0 480 529 252 260 481 530 253 261 +189 5 0 481 530 253 261 482 531 254 262 +190 5 0 482 531 254 262 483 532 255 263 +191 5 0 483 532 255 263 484 533 256 264 +192 5 0 484 533 256 264 351 358 257 265 +193 5 0 28 37 38 29 106 114 534 485 +194 5 0 106 114 534 485 107 115 535 486 +195 5 0 107 115 535 486 108 116 536 487 +196 5 0 108 116 536 487 109 117 537 488 +197 5 0 109 117 537 488 110 118 538 489 +198 5 0 110 118 538 489 111 119 539 490 +199 5 0 111 119 539 490 112 120 540 491 +200 5 0 112 120 540 491 113 121 359 352 +201 5 0 29 38 39 30 485 534 541 492 +202 5 0 485 534 541 492 486 535 542 493 +203 5 0 486 535 542 493 487 536 543 494 +204 5 0 487 536 543 494 488 537 544 495 +205 5 0 488 537 544 495 489 538 545 496 +206 5 0 489 538 545 496 490 539 546 497 +207 5 0 490 539 546 497 491 540 547 498 +208 5 0 491 540 547 498 352 359 360 353 +209 5 0 30 39 40 31 492 541 548 499 +210 5 0 492 541 548 499 493 542 549 500 +211 5 0 493 542 549 500 494 543 550 501 +212 5 0 494 543 550 501 495 544 551 502 +213 5 0 495 544 551 502 496 545 552 503 +214 5 0 496 545 552 503 497 546 553 504 +215 5 0 497 546 553 504 498 547 554 505 +216 5 0 498 547 554 505 353 360 361 354 +217 5 0 31 40 41 32 499 548 555 506 +218 5 0 499 548 555 506 500 549 556 507 +219 5 0 500 549 556 507 501 550 557 508 +220 5 0 501 550 557 508 502 551 558 509 +221 5 0 502 551 558 509 503 552 559 510 +222 5 0 503 552 559 510 504 553 560 511 +223 5 0 504 553 560 511 505 554 561 512 +224 5 0 505 554 561 512 354 361 362 355 +225 5 0 32 41 42 33 506 555 562 513 +226 5 0 506 555 562 513 507 556 563 514 +227 5 0 507 556 563 514 508 557 564 515 +228 5 0 508 557 564 515 509 558 565 516 +229 5 0 509 558 565 516 510 559 566 517 +230 5 0 510 559 566 517 511 560 567 518 +231 5 0 511 560 567 518 512 561 568 519 +232 5 0 512 561 568 519 355 362 363 356 +233 5 0 33 42 43 34 513 562 569 520 +234 5 0 513 562 569 520 514 563 570 521 +235 5 0 514 563 570 521 515 564 571 522 +236 5 0 515 564 571 522 516 565 572 523 +237 5 0 516 565 572 523 517 566 573 524 +238 5 0 517 566 573 524 518 567 574 525 +239 5 0 518 567 574 525 519 568 575 526 +240 5 0 519 568 575 526 356 363 364 357 +241 5 0 34 43 44 35 520 569 576 527 +242 5 0 520 569 576 527 521 570 577 528 +243 5 0 521 570 577 528 522 571 578 529 +244 5 0 522 571 578 529 523 572 579 530 +245 5 0 523 572 579 530 524 573 580 531 +246 5 0 524 573 580 531 525 574 581 532 +247 5 0 525 574 581 532 526 575 582 533 +248 5 0 526 575 582 533 357 364 365 358 +249 5 0 35 44 45 36 527 576 242 250 +250 5 0 527 576 242 250 528 577 243 251 +251 5 0 528 577 243 251 529 578 244 252 +252 5 0 529 578 244 252 530 579 245 253 +253 5 0 530 579 245 253 531 580 246 254 +254 5 0 531 580 246 254 532 581 247 255 +255 5 0 532 581 247 255 533 582 248 256 +256 5 0 533 582 248 256 358 365 249 257 +257 5 0 37 46 47 38 114 122 583 534 +258 5 0 114 122 583 534 115 123 584 535 +259 5 0 115 123 584 535 116 124 585 536 +260 5 0 116 124 585 536 117 125 586 537 +261 5 0 117 125 586 537 118 126 587 538 +262 5 0 118 126 587 538 119 127 588 539 +263 5 0 119 127 588 539 120 128 589 540 +264 5 0 120 128 589 540 121 129 366 359 +265 5 0 38 47 48 39 534 583 590 541 +266 5 0 534 583 590 541 535 584 591 542 +267 5 0 535 584 591 542 536 585 592 543 +268 5 0 536 585 592 543 537 586 593 544 +269 5 0 537 586 593 544 538 587 594 545 +270 5 0 538 587 594 545 539 588 595 546 +271 5 0 539 588 595 546 540 589 596 547 +272 5 0 540 589 596 547 359 366 367 360 +273 5 0 39 48 49 40 541 590 597 548 +274 5 0 541 590 597 548 542 591 598 549 +275 5 0 542 591 598 549 543 592 599 550 +276 5 0 543 592 599 550 544 593 600 551 +277 5 0 544 593 600 551 545 594 601 552 +278 5 0 545 594 601 552 546 595 602 553 +279 5 0 546 595 602 553 547 596 603 554 +280 5 0 547 596 603 554 360 367 368 361 +281 5 0 40 49 50 41 548 597 604 555 +282 5 0 548 597 604 555 549 598 605 556 +283 5 0 549 598 605 556 550 599 606 557 +284 5 0 550 599 606 557 551 600 607 558 +285 5 0 551 600 607 558 552 601 608 559 +286 5 0 552 601 608 559 553 602 609 560 +287 5 0 553 602 609 560 554 603 610 561 +288 5 0 554 603 610 561 361 368 369 362 +289 5 0 41 50 51 42 555 604 611 562 +290 5 0 555 604 611 562 556 605 612 563 +291 5 0 556 605 612 563 557 606 613 564 +292 5 0 557 606 613 564 558 607 614 565 +293 5 0 558 607 614 565 559 608 615 566 +294 5 0 559 608 615 566 560 609 616 567 +295 5 0 560 609 616 567 561 610 617 568 +296 5 0 561 610 617 568 362 369 370 363 +297 5 0 42 51 52 43 562 611 618 569 +298 5 0 562 611 618 569 563 612 619 570 +299 5 0 563 612 619 570 564 613 620 571 +300 5 0 564 613 620 571 565 614 621 572 +301 5 0 565 614 621 572 566 615 622 573 +302 5 0 566 615 622 573 567 616 623 574 +303 5 0 567 616 623 574 568 617 624 575 +304 5 0 568 617 624 575 363 370 371 364 +305 5 0 43 52 53 44 569 618 625 576 +306 5 0 569 618 625 576 570 619 626 577 +307 5 0 570 619 626 577 571 620 627 578 +308 5 0 571 620 627 578 572 621 628 579 +309 5 0 572 621 628 579 573 622 629 580 +310 5 0 573 622 629 580 574 623 630 581 +311 5 0 574 623 630 581 575 624 631 582 +312 5 0 575 624 631 582 364 371 372 365 +313 5 0 44 53 54 45 576 625 234 242 +314 5 0 576 625 234 242 577 626 235 243 +315 5 0 577 626 235 243 578 627 236 244 +316 5 0 578 627 236 244 579 628 237 245 +317 5 0 579 628 237 245 580 629 238 246 +318 5 0 580 629 238 246 581 630 239 247 +319 5 0 581 630 239 247 582 631 240 248 +320 5 0 582 631 240 248 365 372 241 249 +321 5 0 46 55 56 47 122 130 632 583 +322 5 0 122 130 632 583 123 131 633 584 +323 5 0 123 131 633 584 124 132 634 585 +324 5 0 124 132 634 585 125 133 635 586 +325 5 0 125 133 635 586 126 134 636 587 +326 5 0 126 134 636 587 127 135 637 588 +327 5 0 127 135 637 588 128 136 638 589 +328 5 0 128 136 638 589 129 137 373 366 +329 5 0 47 56 57 48 583 632 639 590 +330 5 0 583 632 639 590 584 633 640 591 +331 5 0 584 633 640 591 585 634 641 592 +332 5 0 585 634 641 592 586 635 642 593 +333 5 0 586 635 642 593 587 636 643 594 +334 5 0 587 636 643 594 588 637 644 595 +335 5 0 588 637 644 595 589 638 645 596 +336 5 0 589 638 645 596 366 373 374 367 +337 5 0 48 57 58 49 590 639 646 597 +338 5 0 590 639 646 597 591 640 647 598 +339 5 0 591 640 647 598 592 641 648 599 +340 5 0 592 641 648 599 593 642 649 600 +341 5 0 593 642 649 600 594 643 650 601 +342 5 0 594 643 650 601 595 644 651 602 +343 5 0 595 644 651 602 596 645 652 603 +344 5 0 596 645 652 603 367 374 375 368 +345 5 0 49 58 59 50 597 646 653 604 +346 5 0 597 646 653 604 598 647 654 605 +347 5 0 598 647 654 605 599 648 655 606 +348 5 0 599 648 655 606 600 649 656 607 +349 5 0 600 649 656 607 601 650 657 608 +350 5 0 601 650 657 608 602 651 658 609 +351 5 0 602 651 658 609 603 652 659 610 +352 5 0 603 652 659 610 368 375 376 369 +353 5 0 50 59 60 51 604 653 660 611 +354 5 0 604 653 660 611 605 654 661 612 +355 5 0 605 654 661 612 606 655 662 613 +356 5 0 606 655 662 613 607 656 663 614 +357 5 0 607 656 663 614 608 657 664 615 +358 5 0 608 657 664 615 609 658 665 616 +359 5 0 609 658 665 616 610 659 666 617 +360 5 0 610 659 666 617 369 376 377 370 +361 5 0 51 60 61 52 611 660 667 618 +362 5 0 611 660 667 618 612 661 668 619 +363 5 0 612 661 668 619 613 662 669 620 +364 5 0 613 662 669 620 614 663 670 621 +365 5 0 614 663 670 621 615 664 671 622 +366 5 0 615 664 671 622 616 665 672 623 +367 5 0 616 665 672 623 617 666 673 624 +368 5 0 617 666 673 624 370 377 378 371 +369 5 0 52 61 62 53 618 667 674 625 +370 5 0 618 667 674 625 619 668 675 626 +371 5 0 619 668 675 626 620 669 676 627 +372 5 0 620 669 676 627 621 670 677 628 +373 5 0 621 670 677 628 622 671 678 629 +374 5 0 622 671 678 629 623 672 679 630 +375 5 0 623 672 679 630 624 673 680 631 +376 5 0 624 673 680 631 371 378 379 372 +377 5 0 53 62 63 54 625 674 226 234 +378 5 0 625 674 226 234 626 675 227 235 +379 5 0 626 675 227 235 627 676 228 236 +380 5 0 627 676 228 236 628 677 229 237 +381 5 0 628 677 229 237 629 678 230 238 +382 5 0 629 678 230 238 630 679 231 239 +383 5 0 630 679 231 239 631 680 232 240 +384 5 0 631 680 232 240 372 379 233 241 +385 5 0 55 64 65 56 130 138 681 632 +386 5 0 130 138 681 632 131 139 682 633 +387 5 0 131 139 682 633 132 140 683 634 +388 5 0 132 140 683 634 133 141 684 635 +389 5 0 133 141 684 635 134 142 685 636 +390 5 0 134 142 685 636 135 143 686 637 +391 5 0 135 143 686 637 136 144 687 638 +392 5 0 136 144 687 638 137 145 380 373 +393 5 0 56 65 66 57 632 681 688 639 +394 5 0 632 681 688 639 633 682 689 640 +395 5 0 633 682 689 640 634 683 690 641 +396 5 0 634 683 690 641 635 684 691 642 +397 5 0 635 684 691 642 636 685 692 643 +398 5 0 636 685 692 643 637 686 693 644 +399 5 0 637 686 693 644 638 687 694 645 +400 5 0 638 687 694 645 373 380 381 374 +401 5 0 57 66 67 58 639 688 695 646 +402 5 0 639 688 695 646 640 689 696 647 +403 5 0 640 689 696 647 641 690 697 648 +404 5 0 641 690 697 648 642 691 698 649 +405 5 0 642 691 698 649 643 692 699 650 +406 5 0 643 692 699 650 644 693 700 651 +407 5 0 644 693 700 651 645 694 701 652 +408 5 0 645 694 701 652 374 381 382 375 +409 5 0 58 67 68 59 646 695 702 653 +410 5 0 646 695 702 653 647 696 703 654 +411 5 0 647 696 703 654 648 697 704 655 +412 5 0 648 697 704 655 649 698 705 656 +413 5 0 649 698 705 656 650 699 706 657 +414 5 0 650 699 706 657 651 700 707 658 +415 5 0 651 700 707 658 652 701 708 659 +416 5 0 652 701 708 659 375 382 383 376 +417 5 0 59 68 69 60 653 702 709 660 +418 5 0 653 702 709 660 654 703 710 661 +419 5 0 654 703 710 661 655 704 711 662 +420 5 0 655 704 711 662 656 705 712 663 +421 5 0 656 705 712 663 657 706 713 664 +422 5 0 657 706 713 664 658 707 714 665 +423 5 0 658 707 714 665 659 708 715 666 +424 5 0 659 708 715 666 376 383 384 377 +425 5 0 60 69 70 61 660 709 716 667 +426 5 0 660 709 716 667 661 710 717 668 +427 5 0 661 710 717 668 662 711 718 669 +428 5 0 662 711 718 669 663 712 719 670 +429 5 0 663 712 719 670 664 713 720 671 +430 5 0 664 713 720 671 665 714 721 672 +431 5 0 665 714 721 672 666 715 722 673 +432 5 0 666 715 722 673 377 384 385 378 +433 5 0 61 70 71 62 667 716 723 674 +434 5 0 667 716 723 674 668 717 724 675 +435 5 0 668 717 724 675 669 718 725 676 +436 5 0 669 718 725 676 670 719 726 677 +437 5 0 670 719 726 677 671 720 727 678 +438 5 0 671 720 727 678 672 721 728 679 +439 5 0 672 721 728 679 673 722 729 680 +440 5 0 673 722 729 680 378 385 386 379 +441 5 0 62 71 72 63 674 723 218 226 +442 5 0 674 723 218 226 675 724 219 227 +443 5 0 675 724 219 227 676 725 220 228 +444 5 0 676 725 220 228 677 726 221 229 +445 5 0 677 726 221 229 678 727 222 230 +446 5 0 678 727 222 230 679 728 223 231 +447 5 0 679 728 223 231 680 729 224 232 +448 5 0 680 729 224 232 379 386 225 233 +449 5 0 64 73 74 65 138 146 154 681 +450 5 0 138 146 154 681 139 147 155 682 +451 5 0 139 147 155 682 140 148 156 683 +452 5 0 140 148 156 683 141 149 157 684 +453 5 0 141 149 157 684 142 150 158 685 +454 5 0 142 150 158 685 143 151 159 686 +455 5 0 143 151 159 686 144 152 160 687 +456 5 0 144 152 160 687 145 153 161 380 +457 5 0 65 74 75 66 681 154 162 688 +458 5 0 681 154 162 688 682 155 163 689 +459 5 0 682 155 163 689 683 156 164 690 +460 5 0 683 156 164 690 684 157 165 691 +461 5 0 684 157 165 691 685 158 166 692 +462 5 0 685 158 166 692 686 159 167 693 +463 5 0 686 159 167 693 687 160 168 694 +464 5 0 687 160 168 694 380 161 169 381 +465 5 0 66 75 76 67 688 162 170 695 +466 5 0 688 162 170 695 689 163 171 696 +467 5 0 689 163 171 696 690 164 172 697 +468 5 0 690 164 172 697 691 165 173 698 +469 5 0 691 165 173 698 692 166 174 699 +470 5 0 692 166 174 699 693 167 175 700 +471 5 0 693 167 175 700 694 168 176 701 +472 5 0 694 168 176 701 381 169 177 382 +473 5 0 67 76 77 68 695 170 178 702 +474 5 0 695 170 178 702 696 171 179 703 +475 5 0 696 171 179 703 697 172 180 704 +476 5 0 697 172 180 704 698 173 181 705 +477 5 0 698 173 181 705 699 174 182 706 +478 5 0 699 174 182 706 700 175 183 707 +479 5 0 700 175 183 707 701 176 184 708 +480 5 0 701 176 184 708 382 177 185 383 +481 5 0 68 77 78 69 702 178 186 709 +482 5 0 702 178 186 709 703 179 187 710 +483 5 0 703 179 187 710 704 180 188 711 +484 5 0 704 180 188 711 705 181 189 712 +485 5 0 705 181 189 712 706 182 190 713 +486 5 0 706 182 190 713 707 183 191 714 +487 5 0 707 183 191 714 708 184 192 715 +488 5 0 708 184 192 715 383 185 193 384 +489 5 0 69 78 79 70 709 186 194 716 +490 5 0 709 186 194 716 710 187 195 717 +491 5 0 710 187 195 717 711 188 196 718 +492 5 0 711 188 196 718 712 189 197 719 +493 5 0 712 189 197 719 713 190 198 720 +494 5 0 713 190 198 720 714 191 199 721 +495 5 0 714 191 199 721 715 192 200 722 +496 5 0 715 192 200 722 384 193 201 385 +497 5 0 70 79 80 71 716 194 202 723 +498 5 0 716 194 202 723 717 195 203 724 +499 5 0 717 195 203 724 718 196 204 725 +500 5 0 718 196 204 725 719 197 205 726 +501 5 0 719 197 205 726 720 198 206 727 +502 5 0 720 198 206 727 721 199 207 728 +503 5 0 721 199 207 728 722 200 208 729 +504 5 0 722 200 208 729 385 201 209 386 +505 5 0 71 80 81 72 723 202 210 218 +506 5 0 723 202 210 218 724 203 211 219 +507 5 0 724 203 211 219 725 204 212 220 +508 5 0 725 204 212 220 726 205 213 221 +509 5 0 726 205 213 221 727 206 214 222 +510 5 0 727 206 214 222 728 207 215 223 +511 5 0 728 207 215 223 729 208 216 224 +512 5 0 729 208 216 224 386 209 217 225 +$EndElements diff --git a/test/sumfact/poisson/poisson_dg_3d_gmsh.mini b/test/sumfact/poisson/poisson_dg_3d_gmsh.mini new file mode 100644 index 0000000000000000000000000000000000000000..23f7ebc020396f071fc661ee0fe6a44cf9ee28a1 --- /dev/null +++ b/test/sumfact/poisson/poisson_dg_3d_gmsh.mini @@ -0,0 +1,35 @@ +__name = sumfact_poisson_dg_3d_gmsh_{__exec_suffix} +__exec_suffix = {deg_suffix}_{diff_suffix}_{quadvec_suffix}_{gradvec_suffix} + +diff_suffix = numdiff, symdiff | expand num +quadvec_suffix = quadvec, nonquadvec | expand quad +gradvec_suffix = gradvec, nongradvec | expand grad +deg_suffix = deg{formcompiler.ufl_variants.degree} + +# {deg_suffix} == deg1 | exclude +{diff_suffix} == numdiff | exclude +{quadvec_suffix} == quadvec | exclude +{gradvec_suffix} == gradvec | exclude + +# gmshFile = cube_hexa_2.msh +# gmshFile = cube_hexa_2_consistent.msh +# gmshFile = cube_hexa.msh +gmshFile = cube_hexa_consistent.msh + +[wrapper.vtkcompare] +name = {__name} +extension = vtu + +[formcompiler] +compare_l2errorsquared = 1e-4, 5e-6 | expand deg +grid_unstructured = 1 + +[formcompiler.r] +numerical_jacobian = 1, 0 | expand num +sumfact = 1 +sumfact_regular_jacobians = 1 +vectorization_quadloop = 1, 0 | expand quad +vectorization_strategy = explicit, none | expand grad + +[formcompiler.ufl_variants] +degree = 1, 2 | expand deg diff --git a/test/sumfact/poisson/poisson_dg_3d_unstructured.mini b/test/sumfact/poisson/poisson_dg_3d_unstructured.mini index 012b3356b887eedc31a848b714790d2ac633445b..5c833f28f3dc61bb28fb163d43a4d09eb1243315 100644 --- a/test/sumfact/poisson/poisson_dg_3d_unstructured.mini +++ b/test/sumfact/poisson/poisson_dg_3d_unstructured.mini @@ -6,7 +6,7 @@ quadvec_suffix = quadvec, nonquadvec | expand quad gradvec_suffix = gradvec, nongradvec | expand grad deg_suffix = deg{formcompiler.ufl_variants.degree} -{deg_suffix} == deg1 | exclude +# {deg_suffix} == deg1 | exclude {diff_suffix} == numdiff | exclude {quadvec_suffix} == quadvec | exclude {gradvec_suffix} == gradvec | exclude