diff --git a/doc/source/tutorials/lbs/9pebs/9_pebs_dfem_9grp.ipynb b/doc/source/tutorials/lbs/9pebs/9_pebs_dfem_9grp.ipynb new file mode 100644 index 0000000000..a2ca519f18 --- /dev/null +++ b/doc/source/tutorials/lbs/9pebs/9_pebs_dfem_9grp.ipynb @@ -0,0 +1,306 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "id": "0292bb91-5e35-486f-b90e-39edc840ca87", + "metadata": {}, + "source": [ + "## OpenSn Tutoiral: 3D 9-Pebble BCC Lattice Cell\n", + "\n", + "This script performs a steady-state DFEM k-eigenvalue calculation for a nine-pebble geometry representing an infinite lattice. The corner pebbles are assigned homogenized multigroup cross sections that differ from those of the center pebble. The gap material consists of low-density helium. The unstructured tetrahedral mesh can be created using the jupyter notebook,\n", + "\n", + "opensn/doc/source/tutorials/lbs/9pebs/mesh/gmsh_9peb.ipynb\n", + "\n", + "![Visualization of 9-pebble mesh](images/9peb_mesh_visualization.png)\n", + "\n", + "This figure shows the nine-pebble mesh, with the helium gap material hidden for visual clarity. Each colored region corresponds to a distinct mesh subdomain or block, which can be assigned a specific cross-section representing the material properties for that region.\n", + "\n", + "The cell below begins by setting up MPI and importing all useful Python and OpenSn modules. An executable form of this script,\n", + "\n", + "opensn/doc/source/tutorials/lbs/9pebs/9_pebs_dfem_9grp.py\n", + "\n", + "is designed to be run in parallel using mpirun." + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "4003a12b-7314-4ad3-8254-2234c160ab2e", + "metadata": {}, + "outputs": [], + "source": [ + "import os\n", + "import sys\n", + "import re\n", + "\n", + "# ---------------------------------------------------------\n", + "# MPI + OpenSn Python bindings setup\n", + "# ---------------------------------------------------------\n", + "if \"opensn_console\" not in globals():\n", + " from mpi4py import MPI\n", + " size = MPI.COMM_WORLD.size\n", + " rank = MPI.COMM_WORLD.rank\n", + "\n", + " # Add OpenSn Python module path (adjust if necessary)\n", + " sys.path.append(os.path.abspath(os.path.join(os.path.dirname(\"__file__\"), \"../../../../../\")))\n", + "\n", + " # Import OpenSn components\n", + " from pyopensn.aquad import GLCProductQuadrature3DXYZ\n", + " from pyopensn.solver import DiscreteOrdinatesSolver\n", + " from pyopensn.fieldfunc import FieldFunctionGridBased\n", + " from pyopensn.solver import NonLinearKEigen\n" + ] + }, + { + "cell_type": "markdown", + "id": "99bd45b1-42db-447b-9823-490547c9c77d", + "metadata": {}, + "source": [ + "## Load the 9-Pebble Mesh\n", + "\n", + "This step is important not only to load the mesh file, but also to partition it for parallel execution. Each mpi process will be assigned a subset or partition of the mesh. The partitioner (parmetis) uses this information to decide how to split the mesh for parallel computation. The image below shows an example where the 9-pebble domain is partitioned between 16 mpi tasks.\n", + "\n", + "![Visualization of 9-pebble mesh partition](images/9peb_mesh_partitioned.png)\n", + "\n" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "51446690-1e17-40a3-8cf7-d75d3c6910c0", + "metadata": {}, + "outputs": [], + "source": [ + "mesh_file = \"mesh/msh_meshes/9_peb.msh\"\n", + "\n", + "# Load mesh and partition it for parallel execution\n", + "meshgen = FromFileMeshGenerator(\n", + " filename=mesh_file,\n", + " partitioner=PETScGraphPartitioner(type='parmetis')\n", + ")\n", + "\n", + "grid = meshgen.Execute()\n", + "\n", + "# Export mesh to Paraview format for inspection\n", + "grid.ExportToPVTU(\"9pebs\")\n" + ] + }, + { + "cell_type": "markdown", + "id": "a841ee56-dc91-4488-8ab5-823759b58520", + "metadata": {}, + "source": [ + "## Load Multigroup Cross Sections\n", + "\n", + "Here, there are three cross section files corresponding to three separate materials, the corner pebbles, the center pebble, and the helium gap." + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "6a3b0993-b3b9-421c-a624-cc4ed785fe94", + "metadata": {}, + "outputs": [], + "source": [ + "xss = []\n", + "for m in range(3):\n", + " xss.append(MultiGroupXS())\n", + "\n", + "# Store cross section data as a list\n", + "xss[0].LoadFromOpenSn(\"xs_files/xs_p_1.xs\") # corner pebbles\n", + "xss[1].LoadFromOpenSn(\"xs_files/xs_p_5.xs\") # center pebble\n", + "xss[2].LoadFromOpenSn(\"xs_files/xs_50.xs\") # helium gap\n", + "\n", + "# Store the number of energy groups\n", + "num_groups = xss[0].num_groups\n", + "print(\"Number of energy groups:\", num_groups)\n" + ] + }, + { + "cell_type": "markdown", + "id": "23a9b537-8283-4120-a2c7-c81c254e147a", + "metadata": {}, + "source": [ + "## Map Cross Sections to Mesh Block IDs\n", + "\n", + "Assign XS to each mesh region based on block IDs. Note that a single material may be replicated across several block IDs." + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "abb98609-7b2c-4ff7-9f8b-2380517b5e01", + "metadata": {}, + "outputs": [], + "source": [ + "\n", + "xs_map = []\n", + "xs_map.append({\"block_ids\": [0, 1, 2, 3, 5, 6, 7, 8], \"xs\": xss[0]}) # corner pebbles\n", + "xs_map.append({\"block_ids\": [4], \"xs\": xss[1]}) # center pebble\n", + "xs_map.append({\"block_ids\": [9], \"xs\": xss[2]}) # coolant gap\n" + ] + }, + { + "cell_type": "markdown", + "id": "28eaefad-0db2-416e-b767-71ec1ff4f9dc", + "metadata": {}, + "source": [ + "## Angular Quadrature:\n", + "\n", + "Declare Gauss-Legendre Chebyshev product quadrature" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "59988bef-708d-4784-b659-47baffa54d8c", + "metadata": {}, + "outputs": [], + "source": [ + "\n", + "# Choose S4 quadrature (can try S2 or S8)\n", + "pquad = GLCProductQuadrature3DXYZ(4, 8)" + ] + }, + { + "cell_type": "markdown", + "id": "fe4d9d6f-b0d4-46a8-9635-dd0c53153d0d", + "metadata": {}, + "source": [ + "## Discrete Ordiantes (SN) Solver Setup:\n", + "\n", + "Configure solver with mesh, groups, quadrature, and material map which were all set up in the previous sections." + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "f21b578a-3ef8-4865-beda-6a7a4d616c09", + "metadata": {}, + "outputs": [], + "source": [ + "\n", + "phys = DiscreteOrdinatesSolver(\n", + " mesh=grid,\n", + " num_groups=num_groups,\n", + " groupsets=[\n", + " {\n", + " \"groups_from_to\": (0, num_groups - 1),\n", + " \"angular_quadrature\": pquad,\n", + " \"inner_linear_method\": \"petsc_gmres\",\n", + " \"l_max_its\": 10,\n", + " \"l_abs_tol\": 1.0e-7,\n", + " \"angle_aggregation_type\": \"single\",\n", + " \"angle_aggregation_num_subsets\": 1,\n", + " },\n", + " ],\n", + " xs_map=xs_map,\n", + " options={\n", + " \"boundary_conditions\": [\n", + " {\"name\": \"xmin\", \"type\": \"reflecting\"},\n", + " {\"name\": \"xmax\", \"type\": \"reflecting\"},\n", + " {\"name\": \"ymin\", \"type\": \"reflecting\"},\n", + " {\"name\": \"ymax\", \"type\": \"reflecting\"},\n", + " {\"name\": \"zmin\", \"type\": \"reflecting\"},\n", + " {\"name\": \"zmax\", \"type\": \"reflecting\"}\n", + " ],\n", + " \"scattering_order\": 2,\n", + " \"verbose_outer_iterations\": True,\n", + " \"verbose_inner_iterations\": True,\n", + " \"power_field_function_on\": True,\n", + " \"power_default_kappa\": 1.0,\n", + " \"power_normalization\": 2000.0,\n", + " \"save_angular_flux\": True,\n", + " },\n", + " sweep_type=\"CBC\",\n", + ")\n" + ] + }, + { + "cell_type": "markdown", + "id": "789244d2-f913-4453-8eee-0c0c2734f291", + "metadata": {}, + "source": [ + "## Nonlinear k-Eigenvalue Solver\n", + "\n", + "Initialize and run the criticality calculation." + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "ff3d301d-746a-4871-a259-11a763858f59", + "metadata": {}, + "outputs": [], + "source": [ + "\n", + "k_solver = NonLinearKEigen(\n", + " lbs_solver=phys,\n", + " nl_max_its=50,\n", + " nl_abs_tol=1.0e-10,\n", + " nl_rel_tol=1.0e-10,\n", + " l_max_its=20,\n", + " num_initial_power_iterations=2,\n", + ")\n", + "\n", + "k_solver.Initialize()\n", + "k_solver.Execute()\n" + ] + }, + { + "cell_type": "markdown", + "id": "8bb2fc87-9482-4442-84fe-138bcde09296", + "metadata": {}, + "source": [ + "## Export Scalar Flux to VTK\n", + "\n", + "This code exports the scalar flux field functions computed by the solver to VTK format for visualization and post-processing. By calling GetScalarFieldFunctionList, only the scalar flux fields are retrieved, omitting angular fluxes to reduce file size and simplify analysis. ExportMultipleToVTK then writes these fields to files with the specified basename, creating one or more VTK files suitable for ParaView or other visualization tools. This allows the user to inspect the spatial distribution of neutron flux throughout the nine-pebble geometry, verify solver results, and perform further analyses such as integrating flux to compute reaction rates or visualizing cross-section assignments across different mesh blocks." + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "e9df8de5-e019-4823-a71e-b3fd8f4bf819", + "metadata": {}, + "outputs": [], + "source": [ + "\n", + "fflist = phys.GetScalarFieldFunctionList(only_scalar_flux=True)\n", + "vtk_basename = \"9pebs_dfem_9grp_out\"\n", + "\n", + "FieldFunctionGridBased.ExportMultipleToVTK(fflist, vtk_basename)\n" + ] + }, + { + "cell_type": "markdown", + "id": "ca7633e1-6f7d-479e-9ebc-78374b899886", + "metadata": {}, + "source": [ + "Exporting the scalar flux field to a VTK-compatible mesh enables detailed postprocessing and quantitative analysis in visualization tools such as ParaView. The image below shows the spatial distribution of thermal flux (group 7 of 9) in the nine-pebble problem.\n", + "\n", + "![Visualization g7 Flux Solution](images/9peb_flux_solution_g7.png)" + ] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 3 (ipykernel)", + "language": "python", + "name": "python3" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 3 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython3", + "version": "3.12.2" + } + }, + "nbformat": 4, + "nbformat_minor": 5 +} diff --git a/doc/source/tutorials/lbs/9pebs/images/9peb_flux_solution_g7.png b/doc/source/tutorials/lbs/9pebs/images/9peb_flux_solution_g7.png new file mode 100644 index 0000000000..8e482cc91e Binary files /dev/null and b/doc/source/tutorials/lbs/9pebs/images/9peb_flux_solution_g7.png differ diff --git a/doc/source/tutorials/lbs/9pebs/images/9peb_mesh_partitioned.png b/doc/source/tutorials/lbs/9pebs/images/9peb_mesh_partitioned.png new file mode 100644 index 0000000000..2f9e4ef3e6 Binary files /dev/null and b/doc/source/tutorials/lbs/9pebs/images/9peb_mesh_partitioned.png differ diff --git a/doc/source/tutorials/lbs/9pebs/images/9peb_mesh_visualization.png b/doc/source/tutorials/lbs/9pebs/images/9peb_mesh_visualization.png new file mode 100644 index 0000000000..4f4a01ca6b Binary files /dev/null and b/doc/source/tutorials/lbs/9pebs/images/9peb_mesh_visualization.png differ diff --git a/doc/source/tutorials/lbs/9pebs/index.rst b/doc/source/tutorials/lbs/9pebs/index.rst new file mode 100644 index 0000000000..508a2ae9e9 --- /dev/null +++ b/doc/source/tutorials/lbs/9pebs/index.rst @@ -0,0 +1,8 @@ + +KEigen +======================= + +.. toctree:: + :maxdepth: 1 + + pincell_example diff --git a/doc/source/tutorials/lbs/9pebs/mesh/gmsh_9peb.ipynb b/doc/source/tutorials/lbs/9pebs/mesh/gmsh_9peb.ipynb new file mode 100644 index 0000000000..f8f27f02b4 --- /dev/null +++ b/doc/source/tutorials/lbs/9pebs/mesh/gmsh_9peb.ipynb @@ -0,0 +1,587 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "id": "c8c1808d-7e7c-4a05-b64b-f3467f26985c", + "metadata": {}, + "source": [ + "### 9 pebbles\n", + "One pebble in center: 8 pebbles for each corner" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "b613da6a-ab38-4db2-bfe4-027279c0b4f5", + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Volumes: [(3, 1), (3, 2), (3, 3), (3, 4), (3, 5), (3, 6), (3, 7), (3, 8), (3, 9), (3, 10)]\n", + "Info : Meshing 1D...\n", + "Info : [ 10%] Meshing curve 14 (Circle)\n", + "Info : [ 10%] Meshing curve 40 (Circle)\n", + "Info : [ 10%] Meshing curve 42 (Circle)\n", + "Info : [ 10%] Meshing curve 43 (Circle)\n", + "Info : [ 10%] Meshing curve 44 (Line)\n", + "Info : [ 20%] Meshing curve 45 (Line)\n", + "Info : [ 20%] Meshing curve 46 (Line)\n", + "Info : [ 20%] Meshing curve 59 (Circle)\n", + "Info : [ 20%] Meshing curve 60 (Circle)\n", + "Info : [ 20%] Meshing curve 61 (Circle)\n", + "Info : [ 20%] Meshing curve 63 (Line)\n", + "Info : [ 30%] Meshing curve 64 (Line)\n", + "Info : [ 30%] Meshing curve 65 (Line)\n", + "Info : [ 30%] Meshing curve 78 (Circle)\n", + "Info : [ 30%] Meshing curve 79 (Circle)\n", + "Info : [ 30%] Meshing curve 80 (Circle)\n", + "Info : [ 30%] Meshing curve 82 (Line)\n", + "Info : [ 40%] Meshing curve 83 (Line)\n", + "Info : [ 40%] Meshing curve 84 (Line)\n", + "Info : [ 40%] Meshing curve 97 (Circle)\n", + "Info : [ 40%] Meshing curve 98 (Circle)\n", + "Info : [ 40%] Meshing curve 100 (Circle)\n", + "Info : [ 40%] Meshing curve 101 (Line)\n", + "Info : [ 50%] Meshing curve 102 (Line)\n", + "Info : [ 50%] Meshing curve 103 (Line)\n", + "Info : [ 50%] Meshing curve 116 (Circle)\n", + "Info : [ 50%] Meshing curve 117 (Circle)\n", + "Info : [ 50%] Meshing curve 118 (Circle)\n", + "Info : [ 50%] Meshing curve 120 (Line)\n", + "Info : [ 60%] Meshing curve 121 (Line)\n", + "Info : [ 60%] Meshing curve 122 (Line)\n", + "Info : [ 60%] Meshing curve 135 (Circle)\n", + "Info : [ 60%] Meshing curve 136 (Circle)\n", + "Info : [ 60%] Meshing curve 137 (Circle)\n", + "Info : [ 60%] Meshing curve 139 (Line)\n", + "Info : [ 70%] Meshing curve 140 (Line)\n", + "Info : [ 70%] Meshing curve 141 (Line)\n", + "Info : [ 70%] Meshing curve 154 (Circle)\n", + "Info : [ 70%] Meshing curve 155 (Circle)\n", + "Info : [ 70%] Meshing curve 156 (Circle)\n", + "Info : [ 70%] Meshing curve 158 (Line)\n", + "Info : [ 80%] Meshing curve 159 (Line)\n", + "Info : [ 80%] Meshing curve 160 (Line)\n", + "Info : [ 80%] Meshing curve 173 (Circle)\n", + "Info : [ 80%] Meshing curve 174 (Circle)\n", + "Info : [ 80%] Meshing curve 175 (Circle)\n", + "Info : [ 80%] Meshing curve 177 (Line)\n", + "Info : [ 90%] Meshing curve 178 (Line)\n", + "Info : [ 90%] Meshing curve 179 (Line)\n", + "Info : [ 90%] Meshing curve 180 (Line)\n", + "Info : [ 90%] Meshing curve 181 (Line)\n", + "Info : [ 90%] Meshing curve 182 (Line)\n", + "Info : [ 90%] Meshing curve 183 (Line)\n", + "Info : [ 90%] Meshing curve 184 (Line)\n", + "Info : [100%] Meshing curve 185 (Line)\n", + "Info : [100%] Meshing curve 186 (Line)\n", + "Info : [100%] Meshing curve 187 (Line)\n", + "Info : [100%] Meshing curve 188 (Line)\n", + "Info : [100%] Meshing curve 189 (Line)\n", + "Info : [100%] Meshing curve 190 (Line)\n", + "Info : [100%] Meshing curve 191 (Line)\n", + "Info : Done meshing 1D (Wall 0.00986587s, CPU 0.009547s)\n", + "Info : Meshing 2D...\n", + "Info : [ 0%] Meshing surface 5 (Sphere, Frontal-Delaunay)\n", + "Info : [ 10%] Meshing surface 16 (Sphere, Frontal-Delaunay)\n", + "Info : [ 10%] Meshing surface 17 (Plane, Frontal-Delaunay)\n", + "Info : [ 10%] Meshing surface 18 (Plane, Frontal-Delaunay)\n", + "Info : [ 20%] Meshing surface 19 (Plane, Frontal-Delaunay)\n", + "Info : [ 20%] Meshing surface 26 (Sphere, Frontal-Delaunay)\n", + "Info : [ 20%] Meshing surface 27 (Plane, Frontal-Delaunay)\n", + "Info : [ 20%] Meshing surface 28 (Plane, Frontal-Delaunay)\n", + "Info : [ 30%] Meshing surface 29 (Plane, Frontal-Delaunay)\n", + "Info : [ 30%] Meshing surface 36 (Sphere, Frontal-Delaunay)\n", + "Info : [ 30%] Meshing surface 37 (Plane, Frontal-Delaunay)\n", + "Info : [ 30%] Meshing surface 38 (Plane, Frontal-Delaunay)\n", + "Info : [ 40%] Meshing surface 39 (Plane, Frontal-Delaunay)\n", + "Info : [ 40%] Meshing surface 46 (Sphere, Frontal-Delaunay)\n", + "Info : [ 40%] Meshing surface 47 (Plane, Frontal-Delaunay)\n", + "Info : [ 40%] Meshing surface 48 (Plane, Frontal-Delaunay)\n", + "Info : [ 50%] Meshing surface 49 (Plane, Frontal-Delaunay)\n", + "Info : [ 50%] Meshing surface 56 (Sphere, Frontal-Delaunay)\n", + "Info : [ 50%] Meshing surface 57 (Plane, Frontal-Delaunay)\n", + "Info : [ 50%] Meshing surface 58 (Plane, Frontal-Delaunay)\n", + "Info : [ 60%] Meshing surface 59 (Plane, Frontal-Delaunay)\n", + "Info : [ 60%] Meshing surface 66 (Sphere, Frontal-Delaunay)\n", + "Info : [ 60%] Meshing surface 67 (Plane, Frontal-Delaunay)\n", + "Info : [ 60%] Meshing surface 68 (Plane, Frontal-Delaunay)\n", + "Info : [ 70%] Meshing surface 69 (Plane, Frontal-Delaunay)\n", + "Info : [ 70%] Meshing surface 76 (Sphere, Frontal-Delaunay)\n", + "Info : [ 70%] Meshing surface 77 (Plane, Frontal-Delaunay)\n", + "Info : [ 70%] Meshing surface 78 (Plane, Frontal-Delaunay)\n", + "Info : [ 80%] Meshing surface 79 (Plane, Frontal-Delaunay)\n", + "Info : [ 80%] Meshing surface 86 (Sphere, Frontal-Delaunay)\n", + "Info : [ 80%] Meshing surface 87 (Plane, Frontal-Delaunay)\n", + "Info : [ 80%] Meshing surface 88 (Plane, Frontal-Delaunay)\n", + "Info : [ 90%] Meshing surface 89 (Plane, Frontal-Delaunay)\n", + "Info : [ 90%] Meshing surface 90 (Plane, Frontal-Delaunay)\n", + "Info : [ 90%] Meshing surface 91 (Plane, Frontal-Delaunay)\n", + "Info : [ 90%] Meshing surface 92 (Plane, Frontal-Delaunay)\n", + "Info : [100%] Meshing surface 93 (Plane, Frontal-Delaunay)\n", + "Info : [100%] Meshing surface 94 (Plane, Frontal-Delaunay)\n", + "Info : [100%] Meshing surface 95 (Plane, Frontal-Delaunay)\n", + "Info : Done meshing 2D (Wall 0.0528853s, CPU 0.053003s)\n", + "Info : Meshing 3D...\n", + "Info : 3D Meshing 10 volumes with 1 connected component\n", + "Info : Tetrahedrizing 800 nodes...\n", + "Info : Done tetrahedrizing 808 nodes (Wall 0.00686339s, CPU 0.006449s)\n", + "Info : Reconstructing mesh...\n", + "Info : - Creating surface mesh\n", + "Info : - Identifying boundary edges\n", + "Info : - Recovering boundary\n", + "Info : Done reconstructing mesh (Wall 0.0164832s, CPU 0.014706s)\n", + "Info : Found volume 5\n", + "Info : Found volume 10\n", + "Info : Found volume 1\n", + "Info : Found volume 2\n", + "Info : Found volume 4\n", + "Info : Found volume 3\n", + "Info : Found volume 6\n", + "Info : Found volume 7\n", + "Info : Found volume 8\n", + "Info : Found volume 9\n", + "Info : It. 0 - 0 nodes created - worst tet radius 3.57505 (nodes removed 0 0)\n", + "Info : 3D refinement terminated (845 nodes total):\n", + "Info : - 1 Delaunay cavities modified for star shapeness\n", + "Info : - 1 nodes could not be inserted\n", + "Info : - 3516 tetrahedra created in 0.00345624 sec. (1017288 tets/s)\n", + "Info : Done meshing 3D (Wall 0.03856s, CPU 0.036096s)\n", + "Info : Optimizing mesh...\n", + "Info : Optimizing volume 1\n", + "Info : Optimization starts (volume = 13.6797) with worst = 0.0435143 / average = 0.755978:\n", + "Info : 0.00 < quality < 0.10 : 4 elements\n", + "Info : 0.10 < quality < 0.20 : 0 elements\n", + "Info : 0.20 < quality < 0.30 : 0 elements\n", + "Info : 0.30 < quality < 0.40 : 0 elements\n", + "Info : 0.40 < quality < 0.50 : 0 elements\n", + "Info : 0.50 < quality < 0.60 : 3 elements\n", + "Info : 0.60 < quality < 0.70 : 23 elements\n", + "Info : 0.70 < quality < 0.80 : 52 elements\n", + "Info : 0.80 < quality < 0.90 : 31 elements\n", + "Info : 0.90 < quality < 1.00 : 14 elements\n", + "Info : 4 edge swaps, 0 node relocations (volume = 13.6797): worst = 0.569653 / average = 0.77945 (Wall 8.1841e-05s, CPU 0.000105s)\n", + "Info : No ill-shaped tets in the mesh :-)\n", + "Info : 0.00 < quality < 0.10 : 0 elements\n", + "Info : 0.10 < quality < 0.20 : 0 elements\n", + "Info : 0.20 < quality < 0.30 : 0 elements\n", + "Info : 0.30 < quality < 0.40 : 0 elements\n", + "Info : 0.40 < quality < 0.50 : 0 elements\n", + "Info : 0.50 < quality < 0.60 : 3 elements\n", + "Info : 0.60 < quality < 0.70 : 22 elements\n", + "Info : 0.70 < quality < 0.80 : 52 elements\n", + "Info : 0.80 < quality < 0.90 : 32 elements\n", + "Info : 0.90 < quality < 1.00 : 14 elements\n", + "Info : Optimizing volume 2\n", + "Info : Optimization starts (volume = 13.6781) with worst = 0.0555913 / average = 0.762472:\n", + "Info : 0.00 < quality < 0.10 : 3 elements\n", + "Info : 0.10 < quality < 0.20 : 0 elements\n", + "Info : 0.20 < quality < 0.30 : 0 elements\n", + "Info : 0.30 < quality < 0.40 : 0 elements\n", + "Info : 0.40 < quality < 0.50 : 0 elements\n", + "Info : 0.50 < quality < 0.60 : 3 elements\n", + "Info : 0.60 < quality < 0.70 : 20 elements\n", + "Info : 0.70 < quality < 0.80 : 52 elements\n", + "Info : 0.80 < quality < 0.90 : 35 elements\n", + "Info : 0.90 < quality < 1.00 : 13 elements\n", + "Info : 3 edge swaps, 0 node relocations (volume = 13.6781): worst = 0.527942 / average = 0.780088 (Wall 9.7305e-05s, CPU 0.000111s)\n", + "Info : No ill-shaped tets in the mesh :-)\n", + "Info : 0.00 < quality < 0.10 : 0 elements\n", + "Info : 0.10 < quality < 0.20 : 0 elements\n", + "Info : 0.20 < quality < 0.30 : 0 elements\n", + "Info : 0.30 < quality < 0.40 : 0 elements\n", + "Info : 0.40 < quality < 0.50 : 0 elements\n", + "Info : 0.50 < quality < 0.60 : 3 elements\n", + "Info : 0.60 < quality < 0.70 : 18 elements\n", + "Info : 0.70 < quality < 0.80 : 54 elements\n", + "Info : 0.80 < quality < 0.90 : 35 elements\n", + "Info : 0.90 < quality < 1.00 : 13 elements\n", + "Info : Optimizing volume 3\n", + "Info : Optimization starts (volume = 13.6928) with worst = 0.196187 / average = 0.763788:\n", + "Info : 0.00 < quality < 0.10 : 0 elements\n", + "Info : 0.10 < quality < 0.20 : 1 elements\n", + "Info : 0.20 < quality < 0.30 : 0 elements\n", + "Info : 0.30 < quality < 0.40 : 1 elements\n", + "Info : 0.40 < quality < 0.50 : 0 elements\n", + "Info : 0.50 < quality < 0.60 : 2 elements\n", + "Info : 0.60 < quality < 0.70 : 27 elements\n", + "Info : 0.70 < quality < 0.80 : 51 elements\n", + "Info : 0.80 < quality < 0.90 : 33 elements\n", + "Info : 0.90 < quality < 1.00 : 14 elements\n", + "Info : 1 edge swaps, 0 node relocations (volume = 13.6928): worst = 0.315924 / average = 0.764203 (Wall 4.2622e-05s, CPU 5.7e-05s)\n", + "Info : No ill-shaped tets in the mesh :-)\n", + "Info : 0.00 < quality < 0.10 : 0 elements\n", + "Info : 0.10 < quality < 0.20 : 0 elements\n", + "Info : 0.20 < quality < 0.30 : 0 elements\n", + "Info : 0.30 < quality < 0.40 : 1 elements\n", + "Info : 0.40 < quality < 0.50 : 0 elements\n", + "Info : 0.50 < quality < 0.60 : 3 elements\n", + "Info : 0.60 < quality < 0.70 : 27 elements\n", + "Info : 0.70 < quality < 0.80 : 52 elements\n", + "Info : 0.80 < quality < 0.90 : 33 elements\n", + "Info : 0.90 < quality < 1.00 : 13 elements\n", + "Info : Optimizing volume 4\n", + "Info : Optimization starts (volume = 13.6928) with worst = 0.213054 / average = 0.763725:\n", + "Info : 0.00 < quality < 0.10 : 0 elements\n", + "Info : 0.10 < quality < 0.20 : 0 elements\n", + "Info : 0.20 < quality < 0.30 : 1 elements\n", + "Info : 0.30 < quality < 0.40 : 1 elements\n", + "Info : 0.40 < quality < 0.50 : 0 elements\n", + "Info : 0.50 < quality < 0.60 : 2 elements\n", + "Info : 0.60 < quality < 0.70 : 27 elements\n", + "Info : 0.70 < quality < 0.80 : 51 elements\n", + "Info : 0.80 < quality < 0.90 : 33 elements\n", + "Info : 0.90 < quality < 1.00 : 14 elements\n", + "Info : 1 edge swaps, 0 node relocations (volume = 13.6928): worst = 0.315924 / average = 0.764091 (Wall 4.1306e-05s, CPU 5.5e-05s)\n", + "Info : No ill-shaped tets in the mesh :-)\n", + "Info : 0.00 < quality < 0.10 : 0 elements\n", + "Info : 0.10 < quality < 0.20 : 0 elements\n", + "Info : 0.20 < quality < 0.30 : 0 elements\n", + "Info : 0.30 < quality < 0.40 : 1 elements\n", + "Info : 0.40 < quality < 0.50 : 0 elements\n", + "Info : 0.50 < quality < 0.60 : 3 elements\n", + "Info : 0.60 < quality < 0.70 : 27 elements\n", + "Info : 0.70 < quality < 0.80 : 52 elements\n", + "Info : 0.80 < quality < 0.90 : 33 elements\n", + "Info : 0.90 < quality < 1.00 : 13 elements\n", + "Info : Optimizing volume 5\n", + "Info : Optimization starts (volume = 109.175) with worst = 0.0826662 / average = 0.728682:\n", + "Info : 0.00 < quality < 0.10 : 1 elements\n", + "Info : 0.10 < quality < 0.20 : 3 elements\n", + "Info : 0.20 < quality < 0.30 : 12 elements\n", + "Info : 0.30 < quality < 0.40 : 10 elements\n", + "Info : 0.40 < quality < 0.50 : 20 elements\n", + "Info : 0.50 < quality < 0.60 : 40 elements\n", + "Info : 0.60 < quality < 0.70 : 133 elements\n", + "Info : 0.70 < quality < 0.80 : 218 elements\n", + "Info : 0.80 < quality < 0.90 : 177 elements\n", + "Info : 0.90 < quality < 1.00 : 38 elements\n", + "Info : 16 edge swaps, 0 node relocations (volume = 109.175): worst = 0.307981 / average = 0.744965 (Wall 0.000240412s, CPU 0.000253s)\n", + "Info : No ill-shaped tets in the mesh :-)\n", + "Info : 0.00 < quality < 0.10 : 0 elements\n", + "Info : 0.10 < quality < 0.20 : 0 elements\n", + "Info : 0.20 < quality < 0.30 : 0 elements\n", + "Info : 0.30 < quality < 0.40 : 9 elements\n", + "Info : 0.40 < quality < 0.50 : 15 elements\n", + "Info : 0.50 < quality < 0.60 : 37 elements\n", + "Info : 0.60 < quality < 0.70 : 140 elements\n", + "Info : 0.70 < quality < 0.80 : 219 elements\n", + "Info : 0.80 < quality < 0.90 : 178 elements\n", + "Info : 0.90 < quality < 1.00 : 38 elements\n", + "Info : Optimizing volume 6\n", + "Info : Optimization starts (volume = 13.6785) with worst = 0.0440853 / average = 0.743828:\n", + "Info : 0.00 < quality < 0.10 : 3 elements\n", + "Info : 0.10 < quality < 0.20 : 1 elements\n", + "Info : 0.20 < quality < 0.30 : 1 elements\n", + "Info : 0.30 < quality < 0.40 : 0 elements\n", + "Info : 0.40 < quality < 0.50 : 2 elements\n", + "Info : 0.50 < quality < 0.60 : 5 elements\n", + "Info : 0.60 < quality < 0.70 : 22 elements\n", + "Info : 0.70 < quality < 0.80 : 50 elements\n", + "Info : 0.80 < quality < 0.90 : 29 elements\n", + "Info : 0.90 < quality < 1.00 : 15 elements\n", + "Info : 4 edge swaps, 0 node relocations (volume = 13.6785): worst = 0.240644 / average = 0.766174 (Wall 0.000114599s, CPU 0.000129s)\n", + "Info : No ill-shaped tets in the mesh :-)\n", + "Info : 0.00 < quality < 0.10 : 0 elements\n", + "Info : 0.10 < quality < 0.20 : 0 elements\n", + "Info : 0.20 < quality < 0.30 : 1 elements\n", + "Info : 0.30 < quality < 0.40 : 0 elements\n", + "Info : 0.40 < quality < 0.50 : 2 elements\n", + "Info : 0.50 < quality < 0.60 : 4 elements\n", + "Info : 0.60 < quality < 0.70 : 21 elements\n", + "Info : 0.70 < quality < 0.80 : 52 elements\n", + "Info : 0.80 < quality < 0.90 : 29 elements\n", + "Info : 0.90 < quality < 1.00 : 15 elements\n", + "Info : Optimizing volume 7\n", + "Info : Optimization starts (volume = 13.6889) with worst = 0.137705 / average = 0.757659:\n", + "Info : 0.00 < quality < 0.10 : 0 elements\n", + "Info : 0.10 < quality < 0.20 : 1 elements\n", + "Info : 0.20 < quality < 0.30 : 2 elements\n", + "Info : 0.30 < quality < 0.40 : 0 elements\n", + "Info : 0.40 < quality < 0.50 : 0 elements\n", + "Info : 0.50 < quality < 0.60 : 1 elements\n", + "Info : 0.60 < quality < 0.70 : 30 elements\n", + "Info : 0.70 < quality < 0.80 : 50 elements\n", + "Info : 0.80 < quality < 0.90 : 32 elements\n", + "Info : 0.90 < quality < 1.00 : 14 elements\n", + "Info : 3 edge swaps, 0 node relocations (volume = 13.6889): worst = 0.577879 / average = 0.767897 (Wall 8.6809e-05s, CPU 8.6e-05s)\n", + "Info : No ill-shaped tets in the mesh :-)\n", + "Info : 0.00 < quality < 0.10 : 0 elements\n", + "Info : 0.10 < quality < 0.20 : 0 elements\n", + "Info : 0.20 < quality < 0.30 : 0 elements\n", + "Info : 0.30 < quality < 0.40 : 0 elements\n", + "Info : 0.40 < quality < 0.50 : 0 elements\n", + "Info : 0.50 < quality < 0.60 : 1 elements\n", + "Info : 0.60 < quality < 0.70 : 32 elements\n", + "Info : 0.70 < quality < 0.80 : 50 elements\n", + "Info : 0.80 < quality < 0.90 : 32 elements\n", + "Info : 0.90 < quality < 1.00 : 13 elements\n", + "Info : Optimizing volume 8\n", + "Info : Optimization starts (volume = 13.6889) with worst = 0.137705 / average = 0.718026:\n", + "Info : 0.00 < quality < 0.10 : 0 elements\n", + "Info : 0.10 < quality < 0.20 : 1 elements\n", + "Info : 0.20 < quality < 0.30 : 0 elements\n", + "Info : 0.30 < quality < 0.40 : 1 elements\n", + "Info : 0.40 < quality < 0.50 : 7 elements\n", + "Info : 0.50 < quality < 0.60 : 29 elements\n", + "Info : 0.60 < quality < 0.70 : 15 elements\n", + "Info : 0.70 < quality < 0.80 : 42 elements\n", + "Info : 0.80 < quality < 0.90 : 24 elements\n", + "Info : 0.90 < quality < 1.00 : 16 elements\n", + "Info : 1 edge swaps, 0 node relocations (volume = 13.6889): worst = 0.390684 / average = 0.723448 (Wall 6.3594e-05s, CPU 6.9e-05s)\n", + "Info : No ill-shaped tets in the mesh :-)\n", + "Info : 0.00 < quality < 0.10 : 0 elements\n", + "Info : 0.10 < quality < 0.20 : 0 elements\n", + "Info : 0.20 < quality < 0.30 : 0 elements\n", + "Info : 0.30 < quality < 0.40 : 1 elements\n", + "Info : 0.40 < quality < 0.50 : 6 elements\n", + "Info : 0.50 < quality < 0.60 : 30 elements\n", + "Info : 0.60 < quality < 0.70 : 15 elements\n", + "Info : 0.70 < quality < 0.80 : 42 elements\n", + "Info : 0.80 < quality < 0.90 : 24 elements\n", + "Info : 0.90 < quality < 1.00 : 16 elements\n", + "Info : Optimizing volume 9\n", + "Info : Optimization starts (volume = 13.6889) with worst = 0.137705 / average = 0.718026:\n", + "Info : 0.00 < quality < 0.10 : 0 elements\n", + "Info : 0.10 < quality < 0.20 : 1 elements\n", + "Info : 0.20 < quality < 0.30 : 0 elements\n", + "Info : 0.30 < quality < 0.40 : 1 elements\n", + "Info : 0.40 < quality < 0.50 : 7 elements\n", + "Info : 0.50 < quality < 0.60 : 29 elements\n", + "Info : 0.60 < quality < 0.70 : 15 elements\n", + "Info : 0.70 < quality < 0.80 : 42 elements\n", + "Info : 0.80 < quality < 0.90 : 24 elements\n", + "Info : 0.90 < quality < 1.00 : 16 elements\n", + "Info : 1 edge swaps, 0 node relocations (volume = 13.6889): worst = 0.390684 / average = 0.722083 (Wall 4.1958e-05s, CPU 5.4e-05s)\n", + "Info : No ill-shaped tets in the mesh :-)\n", + "Info : 0.00 < quality < 0.10 : 0 elements\n", + "Info : 0.10 < quality < 0.20 : 0 elements\n", + "Info : 0.20 < quality < 0.30 : 0 elements\n", + "Info : 0.30 < quality < 0.40 : 1 elements\n", + "Info : 0.40 < quality < 0.50 : 7 elements\n", + "Info : 0.50 < quality < 0.60 : 29 elements\n", + "Info : 0.60 < quality < 0.70 : 15 elements\n", + "Info : 0.70 < quality < 0.80 : 42 elements\n", + "Info : 0.80 < quality < 0.90 : 24 elements\n", + "Info : 0.90 < quality < 1.00 : 16 elements\n", + "Info : Optimizing volume 10\n", + "Info : Optimization starts (volume = 152.147) with worst = 0.0126257 / average = 0.749763:\n", + "Info : 0.00 < quality < 0.10 : 2 elements\n", + "Info : 0.10 < quality < 0.20 : 5 elements\n", + "Info : 0.20 < quality < 0.30 : 17 elements\n", + "Info : 0.30 < quality < 0.40 : 37 elements\n", + "Info : 0.40 < quality < 0.50 : 90 elements\n", + "Info : 0.50 < quality < 0.60 : 133 elements\n", + "Info : 0.60 < quality < 0.70 : 278 elements\n", + "Info : 0.70 < quality < 0.80 : 484 elements\n", + "Info : 0.80 < quality < 0.90 : 463 elements\n", + "Info : 0.90 < quality < 1.00 : 312 elements\n", + "Info : 23 edge swaps, 0 node relocations (volume = 152.147): worst = 0.309823 / average = 0.757786 (Wall 0.000409804s, CPU 0.000424s)\n", + "Info : No ill-shaped tets in the mesh :-)\n", + "Info : 0.00 < quality < 0.10 : 0 elements\n", + "Info : 0.10 < quality < 0.20 : 0 elements\n", + "Info : 0.20 < quality < 0.30 : 0 elements\n", + "Info : 0.30 < quality < 0.40 : 34 elements\n", + "Info : 0.40 < quality < 0.50 : 87 elements\n", + "Info : 0.50 < quality < 0.60 : 135 elements\n", + "Info : 0.60 < quality < 0.70 : 280 elements\n", + "Info : 0.70 < quality < 0.80 : 487 elements\n", + "Info : 0.80 < quality < 0.90 : 465 elements\n", + "Info : 0.90 < quality < 1.00 : 311 elements\n", + "Info : Done optimizing mesh (Wall 0.00358097s, CPU 0.003822s)\n", + "Info : 845 nodes 5419 elements\n", + "Info : Writing 'msh_meshes/9_peb.msh'...\n", + "Info : Done writing 'msh_meshes/9_peb.msh'\n", + "Info : Writing 'vtk_meshes/9_peb.vtk'...\n", + "Info : Done writing 'vtk_meshes/9_peb.vtk'\n", + "-------------------------------------------------------\n", + "Version : 4.11.1\n", + "License : GNU General Public License\n", + "Build OS : MacOSX-sdk\n", + "Build date : 20230730\n", + "Build host : Mac-1690686203827.local\n", + "Build options : 64Bit ALGLIB[contrib] ANN[contrib] Bamg Blossom Cairo DIntegration Dlopen DomHex Eigen[contrib] Fltk GMP Gmm[contrib] Hxt Jpeg Kbipack MathEx[contrib] Mesh Metis[contrib] Mpeg Netgen ONELAB ONELABMetamodel OpenCASCADE OpenCASCADE-CAF OpenGL OpenMP OptHom Parser Plugins Png Post QuadMeshingTools QuadTri Solver TetGen/BR Voro++[contrib] WinslowUntangler Zlib\n", + "FLTK version : 1.3.8\n", + "OCC version : 7.7.2\n", + "Packaged by : runner\n", + "Web site : https://gmsh.info\n", + "Issue tracker : https://gitlab.onelab.info/gmsh/gmsh/issues\n", + "-------------------------------------------------------\n" + ] + } + ], + "source": [ + "# Nine 3cm pebbles\n", + "# Side-length: 2*3.4641016151\n", + "\n", + "##########\n", + "# import #\n", + "##########\n", + "import gmsh\n", + "import sys\n", + "\n", + "######################\n", + "# Create Sphere Mesh #\n", + "######################\n", + "# Initialize the gmsh API\n", + "gmsh.initialize()\n", + "# Create a new model \"9pebs\"\n", + "gmsh.model.add(\"9pebs\")\n", + "\n", + "# Define radius of sphere and side length of bounding box\n", + "sphere_radius = 3.0 #cm\n", + "half_length = 3.5921476745623537 # Updated for PF = 0.61 #3.4641016151 # cm\n", + "side_length = half_length*2\n", + "\n", + "# Creates a sphere centered at (0,0,0) with given radius\n", + "p_1 = gmsh.model.occ.addSphere(-half_length, half_length, half_length, sphere_radius)\n", + "p_2 = gmsh.model.occ.addSphere( half_length, half_length, half_length, sphere_radius)\n", + "p_3 = gmsh.model.occ.addSphere( half_length, -half_length, half_length, sphere_radius)\n", + "p_4 = gmsh.model.occ.addSphere(-half_length, -half_length, half_length, sphere_radius)\n", + "p_5 = gmsh.model.occ.addSphere( 0, 0, 0, sphere_radius)\n", + "p_6 = gmsh.model.occ.addSphere(-half_length, half_length, -half_length, sphere_radius)\n", + "p_7 = gmsh.model.occ.addSphere( half_length, half_length, -half_length, sphere_radius)\n", + "p_8 = gmsh.model.occ.addSphere( half_length, -half_length, -half_length, sphere_radius)\n", + "p_9 = gmsh.model.occ.addSphere(-half_length, -half_length, -half_length, sphere_radius)\n", + "\n", + "# creates the bounding box.\n", + "bound_box = gmsh.model.occ.addBox(-half_length, -half_length, -half_length, side_length, side_length, side_length) # addBox(x, y, z, dx, dy, dz, tag=-1)\n", + "\n", + "\n", + "# gmsh.model.occ.cut(objectDimTags, toolDimTags, tag=-1, removeObject=True, removeTool=True)\n", + "# The function \"cut\" subtracts one volume from another volume.\n", + "# In this case, we are cutting the pebble sphere from the bounding box.\n", + "# objectDimTags: The entities to be cut. Each pair (dimension, tag) specifies an object\n", + "# toolDimTags: The cutting tools that remove material from objectDimTags.\n", + "# removeTool: boolean -> if set to True, removes the cutting tools: We don't want this because we want to mesh the pebble sphere also.\n", + "intersect_peb_1 = gmsh.model.occ.intersect([(3, p_1)], [(3, bound_box)], removeObject=True, removeTool=False)\n", + "cut_box_1 = gmsh.model.occ.cut([(3, bound_box)], [(3, p_1)], removeTool=False)\n", + "\n", + "intersect_peb_2 = gmsh.model.occ.intersect([(3, p_2)], [(3, bound_box)], removeObject=True, removeTool=False)\n", + "cut_box_2 = gmsh.model.occ.cut([(3, bound_box)], [(3, p_2)], removeTool=False)\n", + "\n", + "intersect_peb_3 = gmsh.model.occ.intersect([(3, p_3)], [(3, bound_box)], removeObject=True, removeTool=False)\n", + "cut_box_3 = gmsh.model.occ.cut([(3, bound_box)], [(3, p_3)], removeTool=False)\n", + "\n", + "intersect_peb_4 = gmsh.model.occ.intersect([(3, p_4)], [(3, bound_box)], removeObject=True, removeTool=False)\n", + "cut_box_4 = gmsh.model.occ.cut([(3, bound_box)], [(3, p_4)], removeTool=False)\n", + "\n", + "intersect_peb_5 = gmsh.model.occ.intersect([(3, p_5)], [(3, bound_box)], removeObject=True, removeTool=False)\n", + "cut_box_5 = gmsh.model.occ.cut([(3, bound_box)], [(3, p_5)], removeTool=False)\n", + "\n", + "intersect_peb_6 = gmsh.model.occ.intersect([(3, p_6)], [(3, bound_box)], removeObject=True, removeTool=False)\n", + "cut_box_6 = gmsh.model.occ.cut([(3, bound_box)], [(3, p_6)], removeTool=False)\n", + "\n", + "intersect_peb_7 = gmsh.model.occ.intersect([(3, p_7)], [(3, bound_box)], removeObject=True, removeTool=False)\n", + "cut_box_7 = gmsh.model.occ.cut([(3, bound_box)], [(3, p_7)], removeTool=False)\n", + "\n", + "intersect_peb_8 = gmsh.model.occ.intersect([(3, p_8)], [(3, bound_box)], removeObject=True, removeTool=False)\n", + "cut_box_8 = gmsh.model.occ.cut([(3, bound_box)], [(3, p_8)], removeTool=False)\n", + "\n", + "intersect_peb_9 = gmsh.model.occ.intersect([(3, p_9)], [(3, bound_box)], removeObject=True, removeTool=False)\n", + "cut_box_9 = gmsh.model.occ.cut([(3, bound_box)], [(3, p_9)], removeTool=False)\n", + "\n", + "\n", + "\n", + "# Synchronize the CAD kernel with the Gmsh model\n", + "# ChatGPT says: \n", + "# \"When working with OpenCASCADE(OCC)-based functions (gmsh.model.occ.*), changes happen only within the OCC kernel, \n", + "# and Gmsh doesn't immediately \"see\" them. Calling gmsh.model.occ.synchronize() commits these changes to the Gmsh model.\"\n", + "gmsh.model.occ.synchronize()\n", + "\n", + "# Define Surfaces!\n", + "# Get all 2D surfaces (faces)\n", + "surfaces = gmsh.model.getEntities(2)\n", + "# Determine the bounding surfaces of the cube based on location\n", + "tol = 1e-6\n", + "xmin = ymin = zmin = -half_length - tol\n", + "xmax = ymax = zmax = half_length + tol\n", + "\n", + "outer_surfaces = []\n", + "\n", + "for (dim, tag) in surfaces:\n", + " com = gmsh.model.occ.getCenterOfMass(dim, tag)\n", + " x, y, z = com\n", + " if (abs(x - xmin) < tol or abs(x - xmax) < tol or\n", + " abs(y - ymin) < tol or abs(y - ymax) < tol or\n", + " abs(z - zmin) < tol or abs(z - zmax) < tol):\n", + " outer_surfaces.append(tag)\n", + "\n", + "# Create a single physical group for all outer surfaces\n", + "gmsh.model.addPhysicalGroup(2, outer_surfaces, tag=100)\n", + "gmsh.model.setPhysicalName(2, 100, \"outer_boundary\")\n", + "\n", + "\n", + "# Identify the volumes: we expect two volumes\n", + "volEntities = gmsh.model.getEntities(3) # dim=3 specifies that we're looking for volumes (not polygons I suppose)\n", + "print(\"Volumes:\", volEntities)\n", + "\n", + "# Define physical groups\n", + "gmsh.model.addPhysicalGroup(3, [p_1], tag=1) # We specify \"physical groups\" (I think) to identify regions where we can assign properties to regions (i.e. Multi-group cross sections)\n", + "gmsh.model.setPhysicalName(3, 1, \"p_1\") # An integer ID is not enough, we want to assign a name to this region.\n", + "gmsh.model.addPhysicalGroup(3, [p_2], tag=2) \n", + "gmsh.model.setPhysicalName(3, 2, \"p_2\") \n", + "gmsh.model.addPhysicalGroup(3, [p_3], tag=3) \n", + "gmsh.model.setPhysicalName(3, 3, \"p_3\") \n", + "gmsh.model.addPhysicalGroup(3, [p_4], tag=4) \n", + "gmsh.model.setPhysicalName(3, 4, \"p_4\") \n", + "gmsh.model.addPhysicalGroup(3, [p_5], tag=5) \n", + "gmsh.model.setPhysicalName(3, 5, \"p_5\") \n", + "gmsh.model.addPhysicalGroup(3, [p_6], tag=6) \n", + "gmsh.model.setPhysicalName(3, 6, \"p_6\") \n", + "gmsh.model.addPhysicalGroup(3, [p_7], tag=7) \n", + "gmsh.model.setPhysicalName(3, 7, \"p_7\") \n", + "gmsh.model.addPhysicalGroup(3, [p_8], tag=8) \n", + "gmsh.model.setPhysicalName(3, 8, \"p_8\") \n", + "gmsh.model.addPhysicalGroup(3, [p_9], tag=9) \n", + "gmsh.model.setPhysicalName(3, 9, \"p_9\") \n", + "gmsh.model.addPhysicalGroup(3, [bound_box], tag=10) # Material ID: 6\n", + "gmsh.model.setPhysicalName(3, 10, \"gap\") # Region name: \"gap\"\n", + "\n", + "# Create a MathEval field to control the mesh size.\n", + "fieldTag = gmsh.model.mesh.field.add(\"MathEval\") # ChatGPT -> \"Fields in Gmsh are used to control the size, distribution, and behavior of the mesh\"\n", + "gmsh.model.mesh.field.setString(fieldTag, \"F\", \"1.0\")\n", + "\n", + "\n", + "# Use this field as the background mesh field.\n", + "gmsh.model.mesh.field.setAsBackgroundMesh(fieldTag) # Mesh field applied to the entire geometry\n", + "\n", + "# Generate the 3D mesh\n", + "gmsh.model.mesh.generate(3)\n", + "\n", + "# Save mesh in VTK format\n", + "gmsh.write(\"msh_meshes/9_peb.msh\")\n", + "gmsh.write(\"vtk_meshes/9_peb.vtk\")\n", + "\n", + "if \"-nopopup\" not in sys.argv:\n", + " gmsh.fltk.run()\n", + "\n", + "gmsh.finalize()\n", + "\n" + ] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 3 (ipykernel)", + "language": "python", + "name": "python3" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 3 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython3", + "version": "3.12.2" + } + }, + "nbformat": 4, + "nbformat_minor": 5 +} diff --git a/doc/source/tutorials/lbs/9pebs/mesh/msh_meshes/9_peb.msh b/doc/source/tutorials/lbs/9pebs/mesh/msh_meshes/9_peb.msh new file mode 100644 index 0000000000..3f27560ef8 --- /dev/null +++ b/doc/source/tutorials/lbs/9pebs/mesh/msh_meshes/9_peb.msh @@ -0,0 +1,5487 @@ +$MeshFormat +4.1 0 8 +$EndMeshFormat +$PhysicalNames +11 +2 100 "outer_boundary" +3 1 "p_1" +3 2 "p_2" +3 3 "p_3" +3 4 "p_4" +3 5 "p_5" +3 6 "p_6" +3 7 "p_7" +3 8 "p_8" +3 9 "p_9" +3 10 "gap" +$EndPhysicalNames +$Entities +34 71 39 10 +9 1.83697019872103e-16 -4.499279347985573e-32 3 0 +10 1.83697019872103e-16 -4.499279347985573e-32 -3 0 +21 -3.592147674562354 3.592147674562354 3.592147674562354 0 +27 -3.592147674562354 3.592147674562354 0.5921476745623537 0 +28 -0.5921476745623537 3.592147674562353 3.592147674562273 0 +29 -3.592147674562354 0.5921476745623528 3.592147674562354 0 +33 3.592147674562354 3.592147674562354 3.592147674562354 0 +37 3.592147674562354 3.592147674562354 0.5921476745623537 0 +38 0.5921476745623528 3.592147674562354 3.592147674562354 0 +39 3.592147674562354 0.5921476745623528 3.592147674562354 0 +43 3.592147674562354 -3.592147674562354 3.592147674562354 0 +46 3.592147674562354 -3.592147674562354 0.5921476745623537 0 +47 3.592147674562354 -0.5921476745623537 3.592147674562354 0 +48 0.5921476745623528 -3.592147674562354 3.592147674562354 0 +49 -3.592147674562354 -3.592147674562354 3.592147674562354 0 +54 -0.5921476745623537 -3.592147674562354 3.592147674562273 0 +55 -3.592147674562354 -0.5921476745623537 3.592147674562354 0 +56 -3.592147674562354 -3.592147674562354 0.5921476745623537 0 +57 -3.592147674562354 3.592147674562354 -3.592147674562354 0 +61 -0.5921476745623537 3.592147674562353 -3.592147674562434 0 +62 -3.592147674562354 3.592147674562354 -0.5921476745623537 0 +63 -3.592147674562354 0.5921476745623528 -3.592147674562354 0 +65 3.592147674562354 3.592147674562354 -3.592147674562354 0 +67 3.592147674562354 3.592147674562354 -0.5921476745623537 0 +68 3.592147674562354 0.5921476745623528 -3.592147674562354 0 +69 0.5921476745623528 3.592147674562354 -3.592147674562354 0 +71 3.592147674562354 -3.592147674562354 -3.592147674562354 0 +72 0.5921476745623528 -3.592147674562354 -3.592147674562354 0 +73 3.592147674562354 -3.592147674562354 -0.5921476745623537 0 +74 3.592147674562354 -0.5921476745623537 -3.592147674562354 0 +75 -3.592147674562354 -3.592147674562354 -3.592147674562354 0 +76 -3.592147674562354 -0.5921476745623537 -3.592147674562354 0 +77 -3.592147674562354 -3.592147674562354 -0.5921476745623537 0 +78 -0.5921476745623537 -3.592147674562354 -3.592147674562434 0 +13 -1.000000001849447e-07 -1.000000001847366e-07 2.9999999 1.000000001849447e-07 1.000000001847366e-07 3.0000001 0 2 9 -9 +14 -1.000000005024759e-07 -1.000000007347881e-07 -3.0000001 3.0000001 1e-07 3.0000001 0 2 10 -9 +15 -1.000000001849447e-07 -1.000000001847366e-07 -3.0000001 1.000000001849447e-07 1.000000001847366e-07 -2.9999999 0 2 10 -10 +40 -3.592147774562355 3.592147574562352 0.592147574562353 -0.592147574562353 3.592147774562354 3.592147774562274 0 2 27 -28 +41 -3.592147774562354 3.592147574562354 0.5921475745623538 -3.592147574562354 3.592147774562354 0.5921477745623537 0 2 27 -27 +42 -3.592147774562354 0.5921475745623539 0.5921475745623539 -3.592147574562354 3.592147774562353 3.592147774562354 0 2 29 -27 +43 -3.592147774562354 0.5921475745623539 3.592147574562354 -0.5921475745623539 3.592147774562353 3.592147774562354 0 2 29 -28 +44 -3.592147774562354 3.592147574562354 3.592147574562354 -0.5921475745623539 3.592147774562354 3.592147774562354 0 2 21 -28 +45 -3.592147774562354 3.592147574562354 0.5921475745623539 -3.592147574562354 3.592147774562354 3.592147774562354 0 2 27 -21 +46 -3.592147774562354 0.592147574562353 3.592147574562354 -3.592147574562354 3.592147774562354 3.592147774562354 0 2 29 -21 +59 0.5921475745623539 3.592147574562354 0.5921475745623539 3.592147774562354 3.592147774562354 3.592147774562353 0 2 37 -38 +60 0.5921475745623539 0.5921475745623539 3.592147574562354 3.592147774562353 3.592147774562354 3.592147774562354 0 2 38 -39 +61 3.592147574562354 0.5921475745623539 0.5921475745623539 3.592147774562354 3.592147774562353 3.592147774562354 0 2 39 -37 +62 3.592147574562354 3.592147574562354 0.5921475745623538 3.592147774562354 3.592147774562354 0.5921477745623537 0 2 37 -37 +63 3.592147574562354 3.592147574562354 0.5921475745623539 3.592147774562354 3.592147774562354 3.592147774562354 0 2 37 -33 +64 0.592147574562353 3.592147574562354 3.592147574562354 3.592147774562354 3.592147774562354 3.592147774562354 0 2 38 -33 +65 3.592147574562354 0.592147574562353 3.592147574562354 3.592147774562354 3.592147774562354 3.592147774562354 0 2 39 -33 +78 3.592147574562354 -3.592147774562354 0.5921475745623539 3.592147774562354 -0.5921475745623539 3.592147774562353 0 2 46 -47 +79 0.5921475745623539 -3.592147774562353 3.592147574562354 3.592147774562354 -0.5921475745623539 3.592147774562354 0 2 47 -48 +80 0.5921475745623539 -3.592147774562354 0.5921475745623539 3.592147774562354 -3.592147574562354 3.592147774562353 0 2 46 -48 +81 3.592147574562354 -3.592147774562354 0.5921475745623538 3.592147774562354 -3.592147574562354 0.5921477745623537 0 2 46 -46 +82 3.592147574562354 -3.592147774562354 3.592147574562354 3.592147774562354 -0.5921475745623539 3.592147774562354 0 2 43 -47 +83 3.592147574562354 -3.592147774562354 0.5921475745623539 3.592147774562354 -3.592147574562354 3.592147774562354 0 2 46 -43 +84 0.592147574562353 -3.592147774562354 3.592147574562354 3.592147774562354 -3.592147574562354 3.592147774562354 0 2 48 -43 +97 -3.592147774562354 -3.592147774562354 3.592147574562354 -0.5921475745623539 -0.5921475745623539 3.592147774562354 0 2 54 -55 +98 -3.592147774562354 -3.592147774562354 0.5921475745623539 -3.592147574562354 -0.5921475745623539 3.592147774562353 0 2 56 -55 +99 -3.592147774562354 -3.592147774562354 0.5921475745623538 -3.592147574562354 -3.592147574562354 0.5921477745623537 0 2 56 -56 +100 -3.592147774562355 -3.592147774562355 0.592147574562353 -0.592147574562353 -3.592147574562353 3.592147774562274 0 2 56 -54 +101 -3.592147774562354 -3.592147774562354 3.592147574562354 -3.592147574562354 -0.5921475745623539 3.592147774562354 0 2 49 -55 +102 -3.592147774562354 -3.592147774562354 3.592147574562354 -0.5921475745623539 -3.592147574562354 3.592147774562354 0 2 49 -54 +103 -3.592147774562354 -3.592147774562354 0.5921475745623539 -3.592147574562354 -3.592147574562354 3.592147774562354 0 2 56 -49 +116 -3.592147774562354 3.592147574562352 -3.592147774562435 -0.592147574562353 3.592147774562354 -0.592147574562353 0 2 61 -62 +117 -3.592147774562354 0.5921475745623539 -3.592147774562354 -0.5921475745623539 3.592147774562353 -3.592147574562354 0 2 63 -61 +118 -3.592147774562354 0.5921475745623539 -3.592147774562354 -3.592147574562354 3.592147774562354 -0.5921475745623539 0 2 62 -63 +119 -3.592147774562354 3.592147574562354 -0.5921477745623537 -3.592147574562354 3.592147774562354 -0.5921475745623538 0 2 62 -62 +120 -3.592147774562354 3.592147574562354 -3.592147774562354 -3.592147574562354 3.592147774562354 -0.5921475745623539 0 2 57 -62 +121 -3.592147774562354 3.592147574562354 -3.592147774562354 -0.5921475745623539 3.592147774562354 -3.592147574562354 0 2 57 -61 +122 -3.592147774562354 0.592147574562353 -3.592147774562354 -3.592147574562354 3.592147774562354 -3.592147574562354 0 2 63 -57 +135 3.592147574562354 0.5921475745623539 -3.592147774562354 3.592147774562354 3.592147774562354 -0.5921475745623539 0 2 67 -68 +136 0.5921475745623539 0.5921475745623539 -3.592147774562354 3.592147774562353 3.592147774562354 -3.592147574562354 0 2 69 -68 +137 0.5921475745623539 3.592147574562354 -3.592147774562434 3.592147774562353 3.592147774562354 -0.5921475745623539 0 2 69 -67 +138 3.592147574562354 3.592147574562354 -0.5921477745623537 3.592147774562354 3.592147774562354 -0.5921475745623538 0 2 67 -67 +139 3.592147574562354 3.592147574562354 -3.592147774562354 3.592147774562354 3.592147774562354 -0.5921475745623539 0 2 65 -67 +140 3.592147574562354 0.592147574562353 -3.592147774562354 3.592147774562354 3.592147774562354 -3.592147574562354 0 2 68 -65 +141 0.592147574562353 3.592147574562354 -3.592147774562354 3.592147774562354 3.592147774562354 -3.592147574562354 0 2 69 -65 +154 0.5921475745623539 -3.592147774562354 -3.592147774562354 3.592147774562353 -3.592147574562354 -0.5921475745623539 0 2 72 -73 +155 0.5921475745623539 -3.592147774562353 -3.592147774562354 3.592147774562354 -0.5921475745623539 -3.592147574562354 0 2 74 -72 +156 3.592147574562354 -3.592147774562353 -3.592147774562354 3.592147774562354 -0.5921475745623539 -0.5921475745623539 0 2 74 -73 +157 3.592147574562354 -3.592147774562354 -0.5921477745623537 3.592147774562354 -3.592147574562354 -0.5921475745623538 0 2 73 -73 +158 3.592147574562354 -3.592147774562354 -3.592147774562354 3.592147774562354 -3.592147574562354 -0.5921475745623539 0 2 71 -73 +159 0.592147574562353 -3.592147774562354 -3.592147774562354 3.592147774562354 -3.592147574562354 -3.592147574562354 0 2 72 -71 +160 3.592147574562354 -3.592147774562354 -3.592147774562354 3.592147774562354 -0.5921475745623539 -3.592147574562354 0 2 71 -74 +173 -3.592147774562354 -3.592147774562353 -3.592147774562354 -3.592147574562354 -0.5921475745623539 -0.5921475745623539 0 2 76 -77 +174 -3.592147774562354 -3.592147774562354 -3.592147774562354 -0.5921475745623539 -0.5921475745623539 -3.592147574562354 0 2 78 -76 +175 -3.592147774562354 -3.592147774562355 -3.592147774562435 -0.592147574562353 -3.592147574562353 -0.592147574562353 0 2 78 -77 +176 -3.592147774562354 -3.592147774562354 -0.5921477745623537 -3.592147574562354 -3.592147574562354 -0.5921475745623538 0 2 77 -77 +177 -3.592147774562354 -3.592147774562354 -3.592147774562354 -3.592147574562354 -0.5921475745623539 -3.592147574562354 0 2 75 -76 +178 -3.592147774562354 -3.592147774562354 -3.592147774562354 -3.592147574562354 -3.592147574562354 -0.5921475745623539 0 2 75 -77 +179 -3.592147774562354 -3.592147774562354 -3.592147774562354 -0.5921475745623539 -3.592147574562354 -3.592147574562354 0 2 75 -78 +180 -3.592147774562354 -0.5921477745623537 3.592147574562354 -3.592147574562354 0.5921477745623528 3.592147774562354 0 2 55 -29 +181 -3.592147774562354 3.592147574562354 -0.5921477745623537 -3.592147574562354 3.592147774562354 0.5921477745623537 0 2 62 -27 +182 -3.592147774562354 -0.5921477745623537 -3.592147774562354 -3.592147574562354 0.5921477745623528 -3.592147574562354 0 2 76 -63 +183 -3.592147774562354 -3.592147774562354 -0.5921477745623537 -3.592147574562354 -3.592147574562354 0.5921477745623537 0 2 77 -56 +184 -0.5921477745623537 3.592147574562354 3.592147574562354 0.5921477745623528 3.592147774562354 3.592147774562354 0 2 28 -38 +185 3.592147574562354 -0.5921477745623537 3.592147574562354 3.592147774562354 0.5921477745623528 3.592147774562354 0 2 47 -39 +186 -0.5921477745623537 -3.592147774562354 3.592147574562354 0.5921477745623528 -3.592147574562354 3.592147774562354 0 2 54 -48 +187 3.592147574562354 3.592147574562354 -0.5921477745623537 3.592147774562354 3.592147774562354 0.5921477745623537 0 2 67 -37 +188 -0.5921477745623537 3.592147574562354 -3.592147774562354 0.5921477745623528 3.592147774562354 -3.592147574562354 0 2 61 -69 +189 3.592147574562354 -0.5921477745623537 -3.592147774562354 3.592147774562354 0.5921477745623528 -3.592147574562354 0 2 74 -68 +190 -0.5921477745623537 -3.592147774562354 -3.592147774562354 0.5921477745623528 -3.592147574562354 -3.592147574562354 0 2 78 -72 +191 3.592147574562354 -3.592147774562354 -0.5921477745623537 3.592147774562354 -3.592147574562354 0.5921477745623537 0 2 73 -46 +5 -3.0000001 -3.0000001 -3.0000001 3.0000001 3.0000001 3.0000001 0 4 -13 -14 15 14 +16 -3.592147774562354 0.5921475745623539 0.5921475745623539 -0.5921475745623539 3.592147774562354 3.592147774562354 0 4 40 -43 42 41 +17 -3.592147774562355 3.592147574562352 0.592147574562353 -0.592147574562353 3.592147774562354 3.592147774562354 0 3 -40 45 44 +18 -3.592147774562354 0.592147574562353 0.5921475745623539 -3.592147574562354 3.592147774562354 3.592147774562354 0 3 -46 42 45 +19 -3.592147774562354 0.592147574562353 3.592147574562354 -0.5921475745623539 3.592147774562354 3.592147774562354 0 3 -46 43 -44 +26 0.5921475745623539 0.5921475745623539 0.5921475745623539 3.592147774562354 3.592147774562354 3.592147774562354 0 4 -59 62 -61 -60 +27 0.592147574562353 3.592147574562354 0.5921475745623539 3.592147774562354 3.592147774562354 3.592147774562354 0 3 -63 59 64 +28 0.592147574562353 0.592147574562353 3.592147574562354 3.592147774562354 3.592147774562354 3.592147774562354 0 3 60 65 -64 +29 3.592147574562354 0.592147574562353 0.5921475745623539 3.592147774562354 3.592147774562354 3.592147774562354 0 3 -65 61 63 +36 0.5921475745623539 -3.592147774562354 0.5921475745623539 3.592147774562354 -0.5921475745623539 3.592147774562354 0 4 -78 81 80 -79 +37 3.592147574562354 -3.592147774562354 0.5921475745623539 3.592147774562354 -0.5921475745623539 3.592147774562354 0 3 -82 -83 78 +38 0.592147574562353 -3.592147774562354 3.592147574562354 3.592147774562354 -0.5921475745623539 3.592147774562354 0 3 79 84 82 +39 0.592147574562353 -3.592147774562354 0.5921475745623539 3.592147774562354 -3.592147574562354 3.592147774562354 0 3 -83 80 84 +46 -3.592147774562354 -3.592147774562354 0.5921475745623539 -0.5921475745623539 -0.5921475745623539 3.592147774562354 0 4 -97 -100 99 98 +47 -3.592147774562354 -3.592147774562354 3.592147574562354 -0.5921475745623539 -0.5921475745623539 3.592147774562354 0 3 -101 102 97 +48 -3.592147774562354 -3.592147774562354 0.5921475745623539 -3.592147574562354 -0.5921475745623539 3.592147774562354 0 3 -101 -103 98 +49 -3.592147774562355 -3.592147774562355 0.592147574562353 -0.592147574562353 -3.592147574562353 3.592147774562354 0 3 -100 103 102 +56 -3.592147774562354 0.5921475745623539 -3.592147774562434 -0.5921475745623539 3.592147774562354 -0.5921475745623539 0 4 116 -119 118 117 +57 -3.592147774562354 3.592147574562352 -3.592147774562435 -0.592147574562353 3.592147774562354 -0.592147574562353 0 3 -116 -121 120 +58 -3.592147774562354 0.592147574562353 -3.592147774562354 -0.5921475745623539 3.592147774562354 -3.592147574562354 0 3 -122 117 -121 +59 -3.592147774562354 0.592147574562353 -3.592147774562354 -3.592147574562354 3.592147774562354 -0.5921475745623539 0 3 120 118 122 +66 0.5921475745623539 0.5921475745623539 -3.592147774562434 3.592147774562354 3.592147774562354 -0.5921475745623539 0 4 -135 -138 -137 136 +67 3.592147574562354 0.592147574562353 -3.592147774562354 3.592147774562354 3.592147774562354 -0.5921475745623539 0 3 139 135 140 +68 0.592147574562353 0.592147574562353 -3.592147774562354 3.592147774562354 3.592147774562354 -3.592147574562354 0 3 136 140 -141 +69 0.592147574562353 3.592147574562354 -3.592147774562434 3.592147774562354 3.592147774562354 -0.5921475745623539 0 3 -139 -141 137 +76 0.5921475745623539 -3.592147774562354 -3.592147774562354 3.592147774562354 -0.5921475745623539 -0.5921475745623539 0 4 154 -157 -156 155 +77 0.592147574562353 -3.592147774562354 -3.592147774562354 3.592147774562354 -3.592147574562354 -0.5921475745623539 0 3 -158 -159 154 +78 0.592147574562353 -3.592147774562354 -3.592147774562354 3.592147774562354 -0.5921475745623539 -3.592147574562354 0 3 155 159 160 +79 3.592147574562354 -3.592147774562354 -3.592147774562354 3.592147774562354 -0.5921475745623539 -0.5921475745623539 0 3 156 -158 160 +86 -3.592147774562354 -3.592147774562354 -3.592147774562434 -0.5921475745623539 -0.5921475745623539 -0.5921475745623539 0 4 173 -176 -175 174 +87 -3.592147774562354 -3.592147774562354 -3.592147774562354 -3.592147574562354 -0.5921475745623539 -0.5921475745623539 0 3 173 -178 177 +88 -3.592147774562354 -3.592147774562354 -3.592147774562354 -0.5921475745623539 -0.5921475745623539 -3.592147574562354 0 3 -177 179 174 +89 -3.592147774562354 -3.592147774562355 -3.592147774562435 -0.592147574562353 -3.592147574562353 -0.592147574562353 0 3 -175 -179 178 +90 -3.592147774562354 -3.592147774562354 -3.592147774562354 -3.592147574562354 3.592147774562354 3.592147774562354 0 8 -180 -98 -183 -173 182 -118 181 -42 +91 -3.592147774562354 -3.592147774562354 3.592147574562354 3.592147774562354 3.592147774562354 3.592147774562354 0 8 -180 -97 186 -79 185 -60 -184 -43 +92 -3.592147774562355 3.592147574562352 -3.592147774562435 3.592147774562355 3.592147774562354 3.592147774562354 0 8 -187 -137 -188 116 181 40 184 -59 +93 -3.592147774562354 -3.592147774562354 -3.592147774562354 3.592147774562354 3.592147774562354 -3.592147574562354 0 8 -182 -174 190 -155 189 -136 -188 -117 +94 -3.592147774562355 -3.592147774562355 -3.592147774562435 3.592147774562355 -3.592147574562353 3.592147774562354 0 8 -191 -154 -190 175 183 100 186 -80 +95 3.592147574562354 -3.592147774562354 -3.592147774562354 3.592147774562354 3.592147774562354 3.592147774562354 0 8 -185 -78 -191 -156 189 -135 187 -61 +1 -3.592147774562355 0.5921475745623521 0.592147574562353 -0.592147574562353 3.592147774562355 3.592147774562354 1 1 4 16 17 -18 19 +2 0.592147574562353 0.592147574562353 0.5921475745623539 3.592147774562354 3.592147774562354 3.592147774562354 1 2 4 26 27 28 29 +3 0.592147574562353 -3.592147774562354 0.5921475745623539 3.592147774562354 -0.5921475745623539 3.592147774562354 1 3 4 36 37 38 -39 +4 -3.592147774562355 -3.592147774562355 0.592147574562353 -0.592147574562353 -0.592147574562353 3.592147774562354 1 4 4 46 47 -48 -49 +5 -3.0000001 -3.0000001 -3.0000001 3.0000001 3.0000001 3.0000001 1 5 1 5 +6 -3.592147774562355 0.5921475745623521 -3.592147774562435 -0.592147574562353 3.592147774562354 -0.592147574562353 1 6 4 56 57 -58 -59 +7 0.592147574562353 0.592147574562353 -3.592147774562434 3.592147774562354 3.592147774562354 -0.5921475745623539 1 7 4 66 67 -68 69 +8 0.592147574562353 -3.592147774562354 -3.592147774562354 3.592147774562354 -0.5921475745623539 -0.5921475745623539 1 8 4 76 -77 -78 79 +9 -3.592147774562354 -3.592147774562355 -3.592147774562435 -0.592147574562353 -0.592147574562353 -0.592147574562353 1 9 4 86 -87 -88 -89 +10 -3.592147774562355 -3.592147774562355 -3.592147774562435 3.592147774562355 3.592147774562355 3.592147774562354 1 10 15 -90 91 -16 92 -56 -93 -86 -94 -46 -26 95 -36 -66 -76 -5 +$EndEntities +$Nodes +144 845 1 845 +0 9 0 1 +1 +1.83697019872103e-16 -4.499279347985573e-32 3 +0 10 0 1 +2 +1.83697019872103e-16 -4.499279347985573e-32 -3 +0 21 0 1 +3 +-3.592147674562354 3.592147674562354 3.592147674562354 +0 27 0 1 +4 +-3.592147674562354 3.592147674562354 0.5921476745623537 +0 28 0 1 +5 +-0.5921476745623537 3.592147674562353 3.592147674562273 +0 29 0 1 +6 +-3.592147674562354 0.5921476745623528 3.592147674562354 +0 33 0 1 +7 +3.592147674562354 3.592147674562354 3.592147674562354 +0 37 0 1 +8 +3.592147674562354 3.592147674562354 0.5921476745623537 +0 38 0 1 +9 +0.5921476745623528 3.592147674562354 3.592147674562354 +0 39 0 1 +10 +3.592147674562354 0.5921476745623528 3.592147674562354 +0 43 0 1 +11 +3.592147674562354 -3.592147674562354 3.592147674562354 +0 46 0 1 +12 +3.592147674562354 -3.592147674562354 0.5921476745623537 +0 47 0 1 +13 +3.592147674562354 -0.5921476745623537 3.592147674562354 +0 48 0 1 +14 +0.5921476745623528 -3.592147674562354 3.592147674562354 +0 49 0 1 +15 +-3.592147674562354 -3.592147674562354 3.592147674562354 +0 54 0 1 +16 +-0.5921476745623537 -3.592147674562354 3.592147674562273 +0 55 0 1 +17 +-3.592147674562354 -0.5921476745623537 3.592147674562354 +0 56 0 1 +18 +-3.592147674562354 -3.592147674562354 0.5921476745623537 +0 57 0 1 +19 +-3.592147674562354 3.592147674562354 -3.592147674562354 +0 61 0 1 +20 +-0.5921476745623537 3.592147674562353 -3.592147674562434 +0 62 0 1 +21 +-3.592147674562354 3.592147674562354 -0.5921476745623537 +0 63 0 1 +22 +-3.592147674562354 0.5921476745623528 -3.592147674562354 +0 65 0 1 +23 +3.592147674562354 3.592147674562354 -3.592147674562354 +0 67 0 1 +24 +3.592147674562354 3.592147674562354 -0.5921476745623537 +0 68 0 1 +25 +3.592147674562354 0.5921476745623528 -3.592147674562354 +0 69 0 1 +26 +0.5921476745623528 3.592147674562354 -3.592147674562354 +0 71 0 1 +27 +3.592147674562354 -3.592147674562354 -3.592147674562354 +0 72 0 1 +28 +0.5921476745623528 -3.592147674562354 -3.592147674562354 +0 73 0 1 +29 +3.592147674562354 -3.592147674562354 -0.5921476745623537 +0 74 0 1 +30 +3.592147674562354 -0.5921476745623537 -3.592147674562354 +0 75 0 1 +31 +-3.592147674562354 -3.592147674562354 -3.592147674562354 +0 76 0 1 +32 +-3.592147674562354 -0.5921476745623537 -3.592147674562354 +0 77 0 1 +33 +-3.592147674562354 -3.592147674562354 -0.5921476745623537 +0 78 0 1 +34 +-0.5921476745623537 -3.592147674562354 -3.592147674562434 +1 14 0 9 +35 +36 +37 +38 +39 +40 +41 +42 +43 +0.9270509831248417 -2.270620038260489e-16 -2.853169548885461 +1.763355756877419 -4.318975966835979e-16 -2.427050983124843 +2.427050983124843 -5.944560435702552e-16 -1.763355756877418 +2.853169548885461 -6.988249910934556e-16 -0.9270509831248404 +3 -7.347880794884119e-16 1.929747179611964e-15 +2.85316954888546 -6.988249910934553e-16 0.9270509831248441 +2.42705098312484 -5.944560435702543e-16 1.763355756877423 +1.763355756877416 -4.318975966835971e-16 2.427050983124845 +0.9270509831248432 -2.270620038260493e-16 2.85316954888546 +1 40 0 4 +44 +45 +46 +47 +-2.665096691437527 3.592147674562353 0.7389781256768879 +-1.828791917684961 3.592147674562353 1.165096691437492 +-1.16509669143754 3.592147674562353 1.828791917684895 +-0.7389781256769132 3.592147674562353 2.66509669143745 +1 42 0 4 +48 +49 +50 +51 +-3.592147674562354 0.7389781256768932 2.665096691437511 +-3.592147674562354 1.165096691437512 1.828791917684934 +-3.592147674562354 1.828791917684935 1.16509669143751 +-3.592147674562354 2.665096691437514 0.7389781256768919 +1 43 0 4 +52 +53 +54 +55 +-2.665096691437512 0.7389781256768928 3.592147674562354 +-1.828791917684935 1.165096691437511 3.592147674562354 +-1.165096691437512 1.828791917684934 3.592147674562354 +-0.7389781256768932 2.665096691437511 3.592147674562354 +1 44 0 2 +56 +57 +-2.592147674562354 3.592147674562354 3.592147674562354 +-1.592147674562354 3.592147674562354 3.592147674562354 +1 45 0 2 +58 +59 +-3.592147674562354 3.592147674562354 1.592147674562354 +-3.592147674562354 3.592147674562354 2.592147674562354 +1 46 0 2 +60 +61 +-3.592147674562354 1.592147674562353 3.592147674562354 +-3.592147674562354 2.592147674562353 3.592147674562354 +1 59 0 4 +62 +63 +64 +65 +2.665096691437512 3.592147674562354 0.7389781256768928 +1.828791917684935 3.592147674562354 1.165096691437511 +1.165096691437512 3.592147674562354 1.828791917684934 +0.7389781256768932 3.592147674562354 2.665096691437511 +1 60 0 4 +66 +67 +68 +69 +0.7389781256768928 2.665096691437512 3.592147674562354 +1.165096691437511 1.828791917684935 3.592147674562354 +1.828791917684934 1.165096691437512 3.592147674562354 +2.665096691437511 0.7389781256768932 3.592147674562354 +1 61 0 4 +70 +71 +72 +73 +3.592147674562354 0.7389781256768932 2.665096691437511 +3.592147674562354 1.165096691437512 1.828791917684934 +3.592147674562354 1.828791917684935 1.16509669143751 +3.592147674562354 2.665096691437514 0.7389781256768919 +1 63 0 2 +74 +75 +3.592147674562354 3.592147674562354 1.592147674562354 +3.592147674562354 3.592147674562354 2.592147674562354 +1 64 0 2 +76 +77 +1.592147674562353 3.592147674562354 3.592147674562354 +2.592147674562353 3.592147674562354 3.592147674562354 +1 65 0 2 +78 +79 +3.592147674562354 1.592147674562353 3.592147674562354 +3.592147674562354 2.592147674562353 3.592147674562354 +1 78 0 4 +80 +81 +82 +83 +3.592147674562354 -2.665096691437512 0.7389781256768928 +3.592147674562354 -1.828791917684935 1.165096691437511 +3.592147674562354 -1.165096691437512 1.828791917684934 +3.592147674562354 -0.7389781256768932 2.665096691437511 +1 79 0 4 +84 +85 +86 +87 +2.665096691437511 -0.7389781256768932 3.592147674562354 +1.828791917684934 -1.165096691437512 3.592147674562354 +1.16509669143751 -1.828791917684935 3.592147674562354 +0.7389781256768919 -2.665096691437514 3.592147674562354 +1 80 0 4 +88 +89 +90 +91 +2.665096691437512 -3.592147674562354 0.7389781256768928 +1.828791917684935 -3.592147674562354 1.165096691437511 +1.165096691437512 -3.592147674562354 1.828791917684934 +0.7389781256768932 -3.592147674562354 2.665096691437511 +1 82 0 2 +92 +93 +3.592147674562354 -2.592147674562354 3.592147674562354 +3.592147674562354 -1.592147674562354 3.592147674562354 +1 83 0 2 +94 +95 +3.592147674562354 -3.592147674562354 1.592147674562354 +3.592147674562354 -3.592147674562354 2.592147674562354 +1 84 0 2 +96 +97 +1.592147674562353 -3.592147674562354 3.592147674562354 +2.592147674562353 -3.592147674562354 3.592147674562354 +1 97 0 4 +98 +99 +100 +101 +-0.7389781256768932 -2.66509669143751 3.592147674562354 +-1.165096691437512 -1.828791917684933 3.592147674562354 +-1.828791917684936 -1.16509669143751 3.592147674562354 +-2.665096691437514 -0.7389781256768919 3.592147674562354 +1 98 0 4 +102 +103 +104 +105 +-3.592147674562354 -2.665096691437512 0.7389781256768928 +-3.592147674562354 -1.828791917684935 1.165096691437511 +-3.592147674562354 -1.165096691437512 1.828791917684934 +-3.592147674562354 -0.7389781256768932 2.665096691437511 +1 100 0 4 +106 +107 +108 +109 +-2.665096691437527 -3.592147674562354 0.7389781256768879 +-1.828791917684961 -3.592147674562354 1.165096691437492 +-1.16509669143754 -3.592147674562354 1.828791917684895 +-0.7389781256769132 -3.592147674562355 2.66509669143745 +1 101 0 2 +110 +111 +-3.592147674562354 -2.592147674562354 3.592147674562354 +-3.592147674562354 -1.592147674562354 3.592147674562354 +1 102 0 2 +112 +113 +-2.592147674562354 -3.592147674562354 3.592147674562354 +-1.592147674562354 -3.592147674562354 3.592147674562354 +1 103 0 2 +114 +115 +-3.592147674562354 -3.592147674562354 1.592147674562354 +-3.592147674562354 -3.592147674562354 2.592147674562354 +1 116 0 4 +116 +117 +118 +119 +-0.7389781256768733 3.592147674562353 -2.665096691437573 +-1.165096691437483 3.592147674562353 -1.828791917684974 +-1.828791917684908 3.592147674562353 -1.165096691437531 +-2.665096691437495 3.592147674562353 -0.7389781256768981 +1 117 0 4 +120 +121 +122 +123 +-2.665096691437512 0.7389781256768928 -3.592147674562354 +-1.828791917684935 1.165096691437511 -3.592147674562354 +-1.165096691437512 1.828791917684934 -3.592147674562354 +-0.7389781256768932 2.665096691437511 -3.592147674562354 +1 118 0 4 +124 +125 +126 +127 +-3.592147674562354 2.66509669143751 -0.7389781256768932 +-3.592147674562354 1.828791917684933 -1.165096691437512 +-3.592147674562354 1.16509669143751 -1.828791917684936 +-3.592147674562354 0.7389781256768919 -2.665096691437514 +1 120 0 2 +128 +129 +-3.592147674562354 3.592147674562354 -2.592147674562354 +-3.592147674562354 3.592147674562354 -1.592147674562354 +1 121 0 2 +130 +131 +-2.592147674562354 3.592147674562354 -3.592147674562354 +-1.592147674562354 3.592147674562354 -3.592147674562354 +1 122 0 2 +132 +133 +-3.592147674562354 1.592147674562353 -3.592147674562354 +-3.592147674562354 2.592147674562353 -3.592147674562354 +1 135 0 4 +134 +135 +136 +137 +3.592147674562354 2.66509669143751 -0.7389781256768932 +3.592147674562354 1.828791917684933 -1.165096691437512 +3.592147674562354 1.16509669143751 -1.828791917684936 +3.592147674562354 0.7389781256768919 -2.665096691437514 +1 136 0 4 +138 +139 +140 +141 +0.7389781256768928 2.665096691437512 -3.592147674562354 +1.165096691437511 1.828791917684935 -3.592147674562354 +1.828791917684934 1.165096691437512 -3.592147674562354 +2.665096691437511 0.7389781256768932 -3.592147674562354 +1 137 0 4 +142 +143 +144 +145 +0.7389781256768733 3.592147674562354 -2.665096691437573 +1.165096691437483 3.592147674562354 -1.828791917684974 +1.828791917684908 3.592147674562354 -1.165096691437531 +2.665096691437496 3.592147674562354 -0.7389781256768981 +1 139 0 2 +146 +147 +3.592147674562354 3.592147674562354 -2.592147674562354 +3.592147674562354 3.592147674562354 -1.592147674562354 +1 140 0 2 +148 +149 +3.592147674562354 1.592147674562353 -3.592147674562354 +3.592147674562354 2.592147674562353 -3.592147674562354 +1 141 0 2 +150 +151 +1.592147674562353 3.592147674562354 -3.592147674562354 +2.592147674562353 3.592147674562354 -3.592147674562354 +1 154 0 4 +152 +153 +154 +155 +0.7389781256768928 -3.592147674562354 -2.665096691437512 +1.165096691437511 -3.592147674562354 -1.828791917684935 +1.828791917684934 -3.592147674562354 -1.165096691437512 +2.665096691437511 -3.592147674562354 -0.7389781256768932 +1 155 0 4 +156 +157 +158 +159 +2.665096691437511 -0.7389781256768932 -3.592147674562354 +1.828791917684934 -1.165096691437512 -3.592147674562354 +1.16509669143751 -1.828791917684935 -3.592147674562354 +0.7389781256768919 -2.665096691437514 -3.592147674562354 +1 156 0 4 +160 +161 +162 +163 +3.592147674562354 -0.7389781256768928 -2.665096691437512 +3.592147674562354 -1.165096691437511 -1.828791917684935 +3.592147674562354 -1.828791917684934 -1.165096691437512 +3.592147674562354 -2.665096691437511 -0.7389781256768932 +1 158 0 2 +164 +165 +3.592147674562354 -3.592147674562354 -2.592147674562354 +3.592147674562354 -3.592147674562354 -1.592147674562354 +1 159 0 2 +166 +167 +1.592147674562353 -3.592147674562354 -3.592147674562354 +2.592147674562353 -3.592147674562354 -3.592147674562354 +1 160 0 2 +168 +169 +3.592147674562354 -2.592147674562354 -3.592147674562354 +3.592147674562354 -1.592147674562354 -3.592147674562354 +1 173 0 4 +170 +171 +172 +173 +-3.592147674562354 -0.7389781256768928 -2.665096691437512 +-3.592147674562354 -1.165096691437511 -1.828791917684935 +-3.592147674562354 -1.828791917684934 -1.165096691437512 +-3.592147674562354 -2.665096691437511 -0.7389781256768932 +1 174 0 4 +174 +175 +176 +177 +-0.7389781256768932 -2.66509669143751 -3.592147674562354 +-1.165096691437512 -1.828791917684933 -3.592147674562354 +-1.828791917684936 -1.16509669143751 -3.592147674562354 +-2.665096691437514 -0.7389781256768919 -3.592147674562354 +1 175 0 4 +178 +179 +180 +181 +-0.7389781256768733 -3.592147674562355 -2.665096691437573 +-1.165096691437483 -3.592147674562354 -1.828791917684974 +-1.828791917684908 -3.592147674562354 -1.165096691437531 +-2.665096691437495 -3.592147674562354 -0.7389781256768981 +1 177 0 2 +182 +183 +-3.592147674562354 -2.592147674562354 -3.592147674562354 +-3.592147674562354 -1.592147674562354 -3.592147674562354 +1 178 0 2 +184 +185 +-3.592147674562354 -3.592147674562354 -2.592147674562354 +-3.592147674562354 -3.592147674562354 -1.592147674562354 +1 179 0 2 +186 +187 +-2.592147674562354 -3.592147674562354 -3.592147674562354 +-1.592147674562354 -3.592147674562354 -3.592147674562354 +1 180 0 1 +188 +-3.592147674562354 -4.440892098500626e-16 3.592147674562354 +1 181 0 1 +189 +-3.592147674562354 3.592147674562354 0 +1 182 0 1 +190 +-3.592147674562354 -4.440892098500626e-16 -3.592147674562354 +1 183 0 1 +191 +-3.592147674562354 -3.592147674562354 0 +1 184 0 1 +192 +-4.440892098500626e-16 3.592147674562354 3.592147674562354 +1 185 0 1 +193 +3.592147674562354 -4.440892098500626e-16 3.592147674562354 +1 186 0 1 +194 +-4.440892098500626e-16 -3.592147674562354 3.592147674562354 +1 187 0 1 +195 +3.592147674562354 3.592147674562354 0 +1 188 0 1 +196 +-4.440892098500626e-16 3.592147674562354 -3.592147674562354 +1 189 0 1 +197 +3.592147674562354 -4.440892098500626e-16 -3.592147674562354 +1 190 0 1 +198 +-4.440892098500626e-16 -3.592147674562354 -3.592147674562354 +1 191 0 1 +199 +3.592147674562354 -3.592147674562354 0 +2 5 0 153 +200 +201 +202 +203 +204 +205 +206 +207 +208 +209 +210 +211 +212 +213 +214 +215 +216 +217 +218 +219 +220 +221 +222 +223 +224 +225 +226 +227 +228 +229 +230 +231 +232 +233 +234 +235 +236 +237 +238 +239 +240 +241 +242 +243 +244 +245 +246 +247 +248 +249 +250 +251 +252 +253 +254 +255 +256 +257 +258 +259 +260 +261 +262 +263 +264 +265 +266 +267 +268 +269 +270 +271 +272 +273 +274 +275 +276 +277 +278 +279 +280 +281 +282 +283 +284 +285 +286 +287 +288 +289 +290 +291 +292 +293 +294 +295 +296 +297 +298 +299 +300 +301 +302 +303 +304 +305 +306 +307 +308 +309 +310 +311 +312 +313 +314 +315 +316 +317 +318 +319 +320 +321 +322 +323 +324 +325 +326 +327 +328 +329 +330 +331 +332 +333 +334 +335 +336 +337 +338 +339 +340 +341 +342 +343 +344 +345 +346 +347 +348 +349 +350 +351 +352 +1.166838787500599 0.7024625318605942 2.673019572565105 +1.070781057746596 0.8416616250340471 -2.673019572565104 +2.838482652576577 0.8503501713195637 -0.4689571591932767 +2.573881434239385 0.8094516837062918 1.311458094726192 +2.031949497300402 0.9323702430707214 -2.000466688114773 +1.337852541506205 1.67701779582885 -2.09708413986152 +2.855405524791977 -0.8035762062148382 -0.448246104046453 +2.042643871524517 -0.8094516837091589 2.042643871521352 +1.311458094877398 -0.8094516838063265 -2.573881434130882 +0.4518976621510945 -0.8094516837059823 2.853169548885577 +1.998741555422382 1.622535643507644 -1.540263055511931 +1.30714931836422 2.26767573097023 -1.465949330184569 +1.928055290669739 2.151596842377602 -0.8082288190798201 +1.16617888203556 2.681852458703512 -0.6690995477806139 +1.702654943476126 2.469552603946462 0.04770827808210448 +0.8858845481026278 2.858583170488223 0.2090713391034341 +1.374497495457103 2.500625183692466 0.9260833254438661 +0.4835082840691534 2.771660798303231 1.041208988813294 +0.9255161843027009 2.2774110362192 1.71954027713827 +-0.000815087109725736 2.417754339898333 1.77608087865834 +-0.4549274445605066 2.7777911913596 1.037746171947004 +-0.9278225410985775 2.29768813030445 1.691086806786517 +-1.347247742730645 2.530403380659054 0.8842976031045271 +-1.779629276202755 1.965675349231977 1.40322487887858 +-2.105873894092762 2.054048982867079 0.5883688623315353 +-2.442184194124107 1.379094321989446 1.064816985696065 +-2.650118469938542 1.390657207991366 0.2072308499047835 +-2.844476671081695 0.6300501800183462 0.7155342328154854 +-2.93237370901657 0.5987269935927815 -0.206664989322696 +-2.980352534860807 -0.151993435770751 0.3072405628027776 +-2.920950651881632 -0.2243614256160597 -0.6463042936324949 +-2.834989936358431 -0.9768920945113065 -0.09227077775621093 +-2.608939566360919 1.3133123931727 -0.6845764362111011 +-2.617274556353147 -1.007394492567605 -1.065377976593703 +-2.379512550790599 -1.746152341400185 -0.537468158361655 +-2.481002636870314 -1.618812282889862 0.4733631889019928 +-1.90199262103857 -2.319822705720963 0.02909782700305693 +-1.718295498581493 -2.308990988998208 -0.8461803544633136 +-1.131444889304175 -2.757806058502443 -0.3381393294999612 +0.5019140337743767 1.528185496876763 -2.532337139844506 +-2.482687906307913 0.5301089400931862 1.598513456779676 +0.6707903642281456 2.648907060286367 -1.238277704403867 +-1.246167100064491 -2.669938192585199 0.5643559218186138 +-0.4245461267476231 -2.960951293872024 0.2291899246938634 +-0.2351304651205915 -2.911606601712137 -0.6835646723159934 +0.4676628787487339 -2.961326939366222 -0.1087850634247585 +0.3145916622615459 -2.869635584294085 0.8162252749018868 +1.134228255430645 -2.743077405624355 0.4348018069602655 +1.323702211785154 -2.639156535522901 -0.5316627084175189 +-1.37752644044332 1.666140148769328 2.079999497725534 +-2.597961299535923 -0.1987976838857904 -1.486968919310442 +-2.145673895543021 -0.9401038392988776 -1.874110003527207 +-1.989579730402101 -0.09146814759944336 -2.243480794289942 +-1.451798304274744 -0.818660594076047 -2.494409051341025 +-1.583777724926063 2.547692367396823 0.03019468721476255 +0.9771065481078853 -2.525449013461533 1.291266848506095 +0.09556398753088594 -2.524439985829067 1.618045204014566 +-1.209374999841602 0.05280625134280748 -2.744926886016683 +-1.666908628222298 0.7809685558803129 -2.368861274934494 +1.791267213710201 -2.219391864053113 0.9304093308118426 +0.8109007242736155 -2.064817218147347 2.019645976654022 +-0.1545530516457111 -1.939501258339237 2.283516635176436 +0.4179917879759447 1.809860327204153 2.355777676521523 +-1.517647609505193 -1.623176570294661 -2.015451203827479 +-0.7492815938880876 -1.44047829896044 -2.522617561836229 +-0.7757590789356855 -2.174480531013614 -1.915680628834464 +0.05474544755307197 -1.931721224506581 -2.29465815466335 +-0.8058149192615778 -2.288840203392654 1.764050123785866 +-1.043264953886426 -1.57763457142371 2.328662147036594 +-1.642622517583171 -1.8742343228275 1.670041007839049 +-1.81269537037492 -1.076436381740264 2.134343039505894 +-1.141694240022222 -0.7132575859152187 2.681006877729797 +1.812999145104664 1.976275904852242 1.344383742742455 +-0.8145505080982617 0.9043271877240585 -2.742024764166064 +-1.209835168400195 1.590212785217959 -2.237749352146351 +-1.984514813237198 1.421186173112576 -1.744113189387592 +-1.452684300129341 2.141419507793544 -1.517903427691954 +-0.6057001324341369 2.230650369152494 -1.912413731431856 +-0.8973089063588836 2.650919629345366 -1.080491483224396 +2.042643871523493 0.8094516837062629 2.042643871523523 +-1.837011929089651 -0.1671637938719402 2.365891679346428 +-1.091288075111598 0.2181009277555645 2.785950882989753 +2.438048869826848 1.747689155558955 -0.03873401447898107 +0.09617423890534907 -1.325979589050378 -2.68933609747699 +-2.465416357520354 -0.4731085144180506 1.642525652054242 +0.1135513175060471 -2.597566571555039 -1.49658077116888 +0.8011599015639801 -2.145367651004464 -1.937921684214711 +1.414787517226541 1.451146291570769 2.211911101640472 +2.559093393915486 0.7773340893512963 -1.358967517915903 +-1.779751042884875 0.8427424401873584 2.263243558448197 +2.84012264168358 0.8690925264860683 0.4223524128051719 +1.283472743252833 -0.7891460775092202 2.594213982245654 +2.851315291806211 -0.8023191292949974 0.4756943572078404 +2.572630738197097 -1.543213365464245 -0.007974555535052724 +2.428893108072551 -1.507333470545419 -0.9101780474891534 +2.44787515309311 -1.545407139474611 0.7871619961166093 +2.042198352126001 -2.158680449909419 -0.4117336587561384 +1.848087730725162 -2.048023005922409 -1.179056193213734 +2.084019161006539 -1.307283493388405 -1.716937390958625 +2.558131361944778 -0.7378051823969826 1.382608928027991 +2.436843593969496 1.544768462797274 0.8218781508695675 +1.519306051373513 -1.446797565824784 2.144408059531419 +2.085330399345064 -1.489089214779644 1.560131544452037 +-2.296576528792613 -1.354348827867009 1.375309237900384 +-2.756232640134091 -0.7658263423522825 0.9037099351112511 +-0.842417460376378 2.865786821220311 0.2785654640708347 +0.0203184990210467 2.98114877327212 0.3351703301541861 +0.3717083437582854 2.929065119641586 -0.5314230255450911 +-0.3173071521925619 1.583594883646623 -2.528150196419446 +-0.394986631276182 -1.252966576045579 2.697046592186667 +-2.079285262131575 -1.641577198580751 -1.407763154718222 +0.585307165582949 -2.808274073684216 -0.8780730316951506 +-0.5435214715442019 -2.731181689488794 1.115899184053632 +-2.317092279106719 0.5903579881870692 -1.811783876704875 +-0.9084413578236628 -2.61809309665394 -1.149053017335791 +-2.752097011040806 0.5210564489578312 -1.074459035431188 +-1.902513387185268 -2.138983504424202 0.897325124680014 +-2.186097816481397 2.033111927973896 -0.295689406481757 +0.8127807774151296 -1.640497510721404 -2.376584802859207 +-0.8328565263881603 -0.7941696881976671 -2.770495355130679 +-0.4639661218016813 1.923543276316281 2.254931640196379 +-0.06323151615798894 1.292265198347688 2.706668142292988 +2.45830302377172 1.496284700601122 -0.8472770137693034 +0.4302082286077329 2.276010303659844 -1.906488389073592 +2.587057577661003 -0.7439104756806251 -1.324284823608356 +-2.1244352203289 1.824661440549974 -1.075818582288016 +0.1826531976391183 0.6646013926015553 -2.91975046850656 +-0.6380455974880214 -0.3535051464757061 2.909971121324289 +0.7471537977627589 -1.5042744996002 2.485743235400884 +-1.577601951351406 2.438541597154005 -0.7514565603159153 +-0.7304328144979321 0.02191449973562313 -2.90963703203784 +1.477673024771865 -2.022690910322751 1.650819164281567 +1.444994921245449 -1.671219949529129 -2.029535305894114 +-0.9615063239320277 1.12424773931543 2.609898965416643 +2.075365465700437 -0.6210606794477924 -2.075365465700404 +1.999385410462859 2.128482877744701 0.6870361122793334 +0.3037621832680939 0.6437063599107553 2.914338802923683 +-2.79377814864436 -0.02702805406859542 1.092736537528858 +-0.4508163383972806 2.925620616174739 -0.4873485808406661 +-0.3323610264073347 0.7000451399926086 2.898184423065274 +1.129509042406343 -2.42864175198173 -1.351261841262979 +2.167667788634294 -2.056218409233815 0.2705221093464837 +-2.06031728688486 1.317105355503595 1.737908559121324 +-1.429478435665269 -2.152096735269548 -1.524818364268391 +-1.290827183951406 -2.40807980968034 1.238917596687521 +0.563320509497508 1.182728983135507 2.698855712340087 +2.122628582519578 1.34041139216879 1.642481415546324 +-0.1479085616745359 2.684951453268857 -1.330097271244785 +0.4597785581196268 -0.6013445050558898 -2.902927567772321 +-2.443923270424074 1.105337545500063 -1.343677029195482 +-1.099229766789748 2.761706273758291 -0.4057984441658821 +0.6582836210390843 -1.123714152215257 -2.702578246486643 +1.70374188937915 -2.467340365860318 0.09744276972170568 +2 16 0 15 +353 +354 +355 +356 +357 +358 +359 +360 +361 +362 +363 +364 +365 +366 +367 +-2.782695990856043 1.549503803038844 1.549503803038865 +-1.018266240262577 2.782695990759611 2.280689580014218 +-2.322656114688805 0.9950455901951707 2.789993111659756 +-2.280689579839125 2.78269599085887 1.018266240320575 +-1.585357560073541 1.515561259037636 2.77938553515859 +-1.537963656109497 2.783807264394211 1.560666983139943 +-2.049429550634824 2.134647789068088 1.471837585474033 +-1.019304267397318 2.232611976494843 2.86266875924849 +-2.803001597054942 2.308674931310569 0.9979336923156308 +-2.936172820467987 3.029654600584152 0.7192922078668169 +-1.490336155558607 2.074316402190283 2.082652206800844 +-2.143006770221075 1.555915503200771 1.932705242997754 +-0.7455227117199676 2.922548420932785 2.922548420932761 +-2.887607326556211 1.019130243693919 2.219848800964265 +-3.019881370300154 0.707931318835457 2.997330651971756 +2 17 0 6 +368 +369 +370 +371 +372 +373 +-2.071721230861963 3.592147674562354 2.07172123086192 +-1.589211509210295 3.592147674562354 2.594777048051896 +-2.613664536873653 3.592147674562354 1.573583070303509 +-2.338641899662283 3.592147674562354 2.841867416216015 +-2.885997533532256 3.592147674562354 2.290487956384799 +-3.048699857576836 3.592147674562354 3.048699857576835 +2 18 0 6 +374 +375 +376 +377 +378 +379 +-3.592147674562354 2.071721230861952 2.071721230861951 +-3.592147674562354 2.613664536873638 1.573583070303525 +-3.592147674562354 1.589211509210274 2.594777048051929 +-3.592147674562354 2.338641899662254 2.841867416216033 +-3.592147674562354 2.88599753353224 2.290487956384813 +-3.592147674562354 3.048699857576814 3.048699857576813 +2 19 0 6 +380 +381 +382 +383 +384 +385 +-2.071721230861952 2.071721230861951 3.592147674562354 +-2.613664536873638 1.573583070303525 3.592147674562354 +-1.589211509210275 2.594777048051929 3.592147674562354 +-2.338641899662255 2.841867416216036 3.592147674562354 +-2.885997533532239 2.290487956384813 3.592147674562354 +-3.048699857576814 3.048699857576813 3.592147674562354 +2 26 0 15 +386 +387 +388 +389 +390 +391 +392 +393 +394 +395 +396 +397 +398 +399 +400 +1.549503803038849 2.782695990856044 1.54950380303886 +2.782695990783318 1.018266240278273 2.28068957996878 +0.9950455901953661 2.322656114688331 2.789993111659872 +2.782695990858345 2.280689579838473 1.018266240321072 +1.515561259039879 1.585357560071675 2.779385535157467 +2.783807264417812 1.53796365611982 1.560666983120115 +2.134647789073942 2.049429550620491 1.471837585480436 +2.232611976494292 1.01930426740043 2.862668759238537 +2.308674931310491 2.803001597054868 0.9979336923156912 +2.886605357099309 2.89405694127026 0.7610914465159055 +2.074316402188821 1.49033615555972 2.082652206800763 +1.555915503202508 2.143006770220183 1.9327052429964 +2.92254842090954 0.7455227117308958 2.922548420909545 +1.019130243693954 2.887607326555836 2.219848800964391 +0.7079313188353855 3.019881370300331 2.997330651971932 +2 27 0 6 +401 +402 +403 +404 +405 +406 +2.071721230861951 3.592147674562354 2.071721230861951 +1.573583070303525 3.592147674562354 2.613664536873637 +2.594777048051929 3.592147674562354 1.589211509210274 +2.841867416216034 3.592147674562354 2.338641899662254 +2.290487956384813 3.592147674562354 2.88599753353224 +3.048699857576814 3.592147674562354 3.048699857576814 +2 28 0 6 +407 +408 +409 +410 +411 +412 +2.071721230861951 2.071721230861951 3.592147674562354 +2.613664536873638 1.573583070303525 3.592147674562354 +1.589211509210274 2.594777048051928 3.592147674562354 +2.338641899662254 2.841867416216033 3.592147674562354 +2.88599753353224 2.290487956384813 3.592147674562354 +3.048699857576814 3.048699857576814 3.592147674562354 +2 29 0 6 +413 +414 +415 +416 +417 +418 +3.592147674562354 2.071721230861952 2.071721230861951 +3.592147674562354 2.613664536873638 1.573583070303525 +3.592147674562354 1.589211509210274 2.594777048051929 +3.592147674562354 2.338641899662254 2.841867416216033 +3.592147674562354 2.88599753353224 2.290487956384813 +3.592147674562354 3.048699857576814 3.048699857576813 +2 36 0 17 +419 +420 +421 +422 +423 +424 +425 +426 +427 +428 +429 +430 +431 +432 +433 +434 +435 +2.80929305630978 -1.533810110967204 1.554898273783499 +1.549503803034235 -2.78269599085252 1.549503803044871 +1.549503803038942 -1.549503803038943 2.782695990855598 +2.79724563505727 -2.797245635057271 0.8107340959693898 +0.9592335210475769 -2.575348023642348 2.575348023642607 +2.575348023642487 -0.9592335210476826 2.575348023642194 +2.035870581770412 -2.876526444436664 1.129245348962645 +2.153484612516968 -2.123120928372589 1.407609047555428 +1.446214356648071 -2.116831652043867 2.102713947011635 +2.132187662858356 -1.43294522989695 2.106765041307134 +3.003670491146332 -2.107173664067061 1.052748770736436 +2.030312997470084 -1.11326981151634 2.947293441783724 +1.117996484912465 -2.033197199227184 2.922633443043494 +2.937776624869967 -1.074099091221559 2.098417491008477 +1.128995168951977 -2.930938949369638 2.012374666869683 +3.019881370300962 -0.7079313188351306 2.997330651972561 +0.7041580015129476 -3.018620529062065 3.01718483947275 +2 37 0 6 +436 +437 +438 +439 +440 +441 +3.592147674562354 -2.071721230861952 2.071721230861951 +3.592147674562354 -2.594777048051928 1.589211509210274 +3.592147674562354 -1.573583070303525 2.613664536873638 +3.592147674562354 -2.841867416216035 2.338641899662254 +3.592147674562354 -2.290487956384814 2.88599753353224 +3.592147674562354 -3.048699857576816 3.048699857576816 +2 38 0 6 +442 +443 +444 +445 +446 +447 +2.071721230861951 -2.071721230861953 3.592147674562354 +2.594777048051928 -1.589211509210276 3.592147674562354 +1.573583070303525 -2.613664536873639 3.592147674562354 +2.841867416216032 -2.338641899662258 3.592147674562354 +2.290487956384812 -2.885997533532241 3.592147674562354 +3.048699857576814 -3.048699857576814 3.592147674562354 +2 39 0 6 +448 +449 +450 +451 +452 +453 +2.071721230861951 -3.592147674562354 2.071721230861951 +1.573583070303525 -3.592147674562354 2.613664536873637 +2.594777048051929 -3.592147674562354 1.589211509210274 +2.841867416216034 -3.592147674562354 2.338641899662254 +2.290487956384813 -3.592147674562354 2.88599753353224 +3.048699857576814 -3.592147674562354 3.048699857576814 +2 46 0 17 +454 +455 +456 +457 +458 +459 +460 +461 +462 +463 +464 +465 +466 +467 +468 +469 +470 +-2.782695990856038 -1.549503803038855 1.549503803038855 +-1.533810110967838 -2.809293056310386 1.554898273782627 +-1.549503803038944 -1.549503803038942 2.782695990855598 +-0.959233521047655 -2.575348023642396 2.575348023642356 +-2.575348023642364 -0.9592335210476111 2.575348023642501 +-2.797245635057679 -2.797245635057678 0.8107340959691567 +-2.876526444439617 -2.035870581772754 1.129245348960306 +-2.123120928373712 -2.153484612518196 1.407609047553863 +-2.116831652044697 -1.446214356649564 2.102713947008661 +-1.432945229897022 -2.132187662858414 2.106765041306974 +-2.10717366406757 -3.003670491147014 1.052748770735981 +-1.074099091221852 -2.937776624869612 2.098417491008139 +-2.930938949366774 -1.128995168954801 2.012374666866478 +-1.113269811516391 -2.030312997469983 2.947293441783774 +-2.033197199227212 -1.117996484912452 2.922633443043479 +-0.707931318835116 -3.019881370301011 2.997330651972587 +-3.018620529062022 -0.704158001512964 3.017184839472712 +2 47 0 6 +471 +472 +473 +474 +475 +476 +-2.071721230861952 -2.071721230861951 3.592147674562354 +-2.613664536873638 -1.573583070303525 3.592147674562354 +-1.589211509210275 -2.594777048051927 3.592147674562354 +-2.338641899662256 -2.841867416216035 3.592147674562354 +-2.88599753353224 -2.290487956384814 3.592147674562354 +-3.048699857576815 -3.048699857576815 3.592147674562354 +2 48 0 6 +477 +478 +479 +480 +481 +482 +-3.592147674562354 -2.071721230861951 2.071721230861952 +-3.592147674562354 -2.594777048051928 1.589211509210274 +-3.592147674562354 -1.573583070303525 2.613664536873638 +-3.592147674562354 -2.841867416216035 2.338641899662256 +-3.592147674562354 -2.290487956384813 2.88599753353224 +-3.592147674562354 -3.048699857576817 3.048699857576816 +2 49 0 6 +483 +484 +485 +486 +487 +488 +-2.071721230861963 -3.592147674562354 2.07172123086192 +-1.589211509210295 -3.592147674562354 2.594777048051896 +-2.613664536873653 -3.592147674562354 1.573583070303509 +-2.338641899662283 -3.592147674562354 2.841867416216015 +-2.885997533532256 -3.592147674562354 2.290487956384799 +-3.048699857576836 -3.592147674562354 3.048699857576835 +2 56 0 15 +489 +490 +491 +492 +493 +494 +495 +496 +497 +498 +499 +500 +501 +502 +503 +-1.549503803038822 2.782695990855999 -1.549503803038904 +-2.847998580961539 0.9962945922759885 -2.28532369954603 +-1.014420513002601 2.277820977954445 -2.799766443625602 +-2.782695990853391 2.280689579829846 -1.018266240327026 +-1.524292919074852 1.565281079289619 -2.807425739622825 +-2.828108588096881 1.512087629444103 -1.569873093043893 +-2.054029202786717 1.482479198928469 -2.114483812346329 +-2.308674931310323 2.803001597054871 -0.9979336923157729 +-2.031688382512645 2.202691851021787 -1.439386816435069 +-1.007292629651348 2.857145521396991 -2.258619935203175 +-1.437467011119725 2.160224409849294 -2.073283953697336 +-2.246715235267021 1.007012325335437 -2.880187860057885 +-2.886605357102475 2.894056941258141 -0.7610914465181047 +-2.922548420761018 0.7455227118007892 -2.922548420760938 +-0.7479571321807903 2.917395809567792 -2.9173958095678 +2 57 0 6 +504 +505 +506 +507 +508 +509 +-2.071721230861941 3.592147674562354 -2.071721230861983 +-2.594777048051917 3.592147674562354 -1.589211509210288 +-1.57358307030351 3.592147674562354 -2.613664536873665 +-2.84186741621603 3.592147674562354 -2.338641899662267 +-2.290487956384806 3.592147674562354 -2.88599753353225 +-3.048699857576823 3.592147674562354 -3.048699857576823 +2 58 0 6 +510 +511 +512 +513 +514 +515 +-2.071721230861952 2.071721230861951 -3.592147674562354 +-2.613664536873638 1.573583070303525 -3.592147674562354 +-1.589211509210275 2.594777048051929 -3.592147674562354 +-2.338641899662255 2.841867416216036 -3.592147674562354 +-2.885997533532239 2.290487956384813 -3.592147674562354 +-3.048699857576814 3.048699857576813 -3.592147674562354 +2 59 0 6 +516 +517 +518 +519 +520 +521 +-3.592147674562354 2.071721230861951 -2.071721230861952 +-3.592147674562354 2.594777048051925 -1.589211509210276 +-3.592147674562354 1.573583070303525 -2.613664536873639 +-3.592147674562354 2.841867416216032 -2.338641899662258 +-3.592147674562354 2.290487956384812 -2.885997533532241 +-3.592147674562354 3.048699857576814 -3.048699857576814 +2 66 0 17 +522 +523 +524 +525 +526 +527 +528 +529 +530 +531 +532 +533 +534 +535 +536 +537 +538 +1.536676294943638 1.557195634400713 -2.795901731520691 +2.782695990848861 1.549503803029373 -1.549503803051182 +1.536127828521712 2.847961974574723 -1.538134520607167 +0.959233521047675 2.575348023642658 -2.575348023642042 +2.575348023642404 0.959233521047508 -2.575348023642728 +2.797245635059904 2.797245635059905 -0.8107340959678844 +3.018397535679533 2.139799903007257 -1.030605809085452 +2.184078656931329 2.221667672594235 -1.325183580053342 +2.113699630804211 1.406023559900192 -2.165594928534259 +1.432512681359686 2.151129959498295 -2.08900275772763 +2.125441013539916 3.002565331525734 -1.042409603447624 +1.128995168958233 2.930938949363586 -2.012374666862462 +1.117996484912315 2.033197199227469 -2.922633443043384 +2.030312997469952 1.113269811516381 -2.947293441783888 +2.946609792128102 1.11108618282038 -2.034067215620199 +0.7041580015129991 3.018620529061914 -3.017184839472641 +3.019881370301156 0.7079313188350516 -2.99733065197276 +2 67 0 6 +539 +540 +541 +542 +543 +544 +3.592147674562354 2.071721230861951 -2.071721230861952 +3.592147674562354 2.594777048051925 -1.589211509210276 +3.592147674562354 1.573583070303525 -2.613664536873639 +3.592147674562354 2.841867416216032 -2.338641899662258 +3.592147674562354 2.290487956384812 -2.885997533532241 +3.592147674562354 3.048699857576814 -3.048699857576814 +2 68 0 6 +545 +546 +547 +548 +549 +550 +2.071721230861951 2.071721230861951 -3.592147674562354 +2.594777048051929 1.589211509210274 -3.592147674562354 +1.573583070303526 2.613664536873637 -3.592147674562354 +2.841867416216034 2.338641899662254 -3.592147674562354 +2.290487956384813 2.885997533532238 -3.592147674562354 +3.048699857576811 3.048699857576811 -3.592147674562354 +2 69 0 6 +551 +552 +553 +554 +555 +556 +2.071721230861942 3.592147674562354 -2.071721230861983 +1.573583070303511 3.592147674562354 -2.613664536873664 +2.594777048051918 3.592147674562354 -1.589211509210288 +2.841867416216032 3.592147674562354 -2.338641899662264 +2.290487956384806 3.592147674562354 -2.885997533532248 +3.048699857576818 3.592147674562354 -3.048699857576818 +2 76 0 17 +557 +558 +559 +560 +561 +562 +563 +564 +565 +566 +567 +568 +569 +570 +571 +572 +573 +1.557195634400837 -1.536676294943542 -2.795901731520621 +1.54950380303276 -2.782695990851531 -1.549503803046737 +2.847961974576985 -1.536127828520454 -1.538134520607606 +0.9592335210475449 -2.575348023642388 -2.57534802364265 +2.575348023642517 -0.9592335210476604 -2.575348023642221 +2.797245635058387 -2.797245635058387 -0.8107340959687517 +2.139799903009629 -3.018397535683967 -1.030605809083114 +2.221667672598512 -2.184078656935648 -1.325183580048074 +1.406023559901287 -2.113699630804237 -2.165594928532553 +2.151129959498506 -1.432512681359007 -2.089002757728401 +3.002565331536621 -2.125441013537643 -1.042409603446414 +2.930938949365019 -1.128995168956697 -2.012374666864257 +2.033197199227373 -1.117996484912358 -2.92263344304345 +1.113269811516379 -2.030312997469958 -2.947293441783879 +1.111086182822476 -2.946609792125824 -2.034067215617806 +3.018620529061957 -0.704158001512988 -3.017184839472655 +0.707931318835068 -3.019881370301114 -2.99733065197272 +2 77 0 6 +574 +575 +576 +577 +578 +579 +2.071721230861951 -3.592147674562354 -2.071721230861952 +1.573583070303525 -3.592147674562354 -2.613664536873638 +2.594777048051929 -3.592147674562354 -1.589211509210275 +2.841867416216036 -3.592147674562354 -2.338641899662255 +2.290487956384813 -3.592147674562354 -2.885997533532239 +3.048699857576813 -3.592147674562354 -3.048699857576814 +2 78 0 6 +580 +581 +582 +583 +584 +585 +2.071721230861951 -2.071721230861953 -3.592147674562354 +2.594777048051928 -1.589211509210276 -3.592147674562354 +1.573583070303525 -2.613664536873639 -3.592147674562354 +2.841867416216032 -2.338641899662258 -3.592147674562354 +2.290487956384812 -2.885997533532241 -3.592147674562354 +3.048699857576814 -3.048699857576814 -3.592147674562354 +2 79 0 6 +586 +587 +588 +589 +590 +591 +3.592147674562354 -2.071721230861951 -2.071721230861952 +3.592147674562354 -2.594777048051928 -1.589211509210275 +3.592147674562354 -1.573583070303525 -2.613664536873638 +3.592147674562354 -2.841867416216035 -2.338641899662256 +3.592147674562354 -2.290487956384814 -2.88599753353224 +3.592147674562354 -3.048699857576815 -3.048699857576815 +2 86 0 17 +592 +593 +594 +595 +596 +597 +598 +599 +600 +601 +602 +603 +604 +605 +606 +607 +608 +-1.536676294943453 -1.557195634400953 -2.795901731520554 +-2.782695990854226 -1.549503803036182 -1.549503803042247 +-1.536127828519159 -2.847961974579226 -1.53813452060809 +-0.9592335210476133 -2.57534802364246 -2.5753480236424 +-2.575348023642403 -0.9592335210475822 -2.575348023642538 +-2.797245635056845 -2.797245635056846 -0.8107340959696328 +-3.018397535688377 -2.139799903011912 -1.030605809080831 +-2.184078656939949 -2.221667672602747 -1.325183580042841 +-2.113699630804277 -1.406023559902374 -2.165594928530846 +-1.432512681358324 -2.151129959498684 -2.089002757729213 +-2.125441013535329 -3.002565331547482 -1.042409603445233 +-1.12899516895515 -2.930938949366486 -2.012374666866055 +-1.117996484912387 -2.033197199227313 -2.922633443043483 +-2.030312997469991 -1.113269811516364 -2.947293441783856 +-2.946609792123535 -1.111086182824605 -2.034067215615363 +-3.019881370301072 -0.7079313188350849 -2.997330651972677 +-0.7041580015129685 -3.018620529061995 -3.017184839472714 +2 87 0 6 +609 +610 +611 +612 +613 +614 +-3.592147674562354 -2.071721230861951 -2.071721230861952 +-3.592147674562354 -2.594777048051928 -1.589211509210275 +-3.592147674562354 -1.573583070303525 -2.613664536873638 +-3.592147674562354 -2.841867416216035 -2.338641899662256 +-3.592147674562354 -2.290487956384814 -2.88599753353224 +-3.592147674562354 -3.048699857576815 -3.048699857576815 +2 88 0 6 +615 +616 +617 +618 +619 +620 +-2.071721230861953 -2.071721230861951 -3.592147674562354 +-1.589211509210275 -2.594777048051927 -3.592147674562354 +-2.613664536873638 -1.573583070303525 -3.592147674562354 +-2.338641899662256 -2.841867416216034 -3.592147674562354 +-2.88599753353224 -2.290487956384814 -3.592147674562354 +-3.048699857576815 -3.048699857576815 -3.592147674562354 +2 89 0 6 +621 +622 +623 +624 +625 +626 +-2.071721230861941 -3.592147674562354 -2.071721230861983 +-2.594777048051917 -3.592147674562354 -1.589211509210288 +-1.57358307030351 -3.592147674562354 -2.613664536873665 +-2.84186741621603 -3.592147674562354 -2.338641899662267 +-2.290487956384806 -3.592147674562354 -2.88599753353225 +-3.048699857576823 -3.592147674562354 -3.048699857576823 +2 90 0 29 +627 +628 +629 +630 +631 +632 +633 +634 +635 +636 +637 +638 +639 +640 +641 +642 +643 +644 +645 +646 +647 +648 +649 +650 +651 +652 +653 +654 +655 +-3.592147674562354 -0.9169864319818095 0.9221673782604944 +-3.592147674562354 0.9120792263157926 -0.9546359275938987 +-3.592147674562354 0.9598168738725819 0.9120792263157953 +-3.592147674562354 -0.9221673782604936 -0.9221673782604949 +-3.592147674562354 1.667633937798202 0.3912212927634418 +-3.592147674562354 -1.63917183758541 0.3544939373043721 +-3.592147674562354 -0.3868582140957599 1.661359021862855 +-3.592147674562354 0.39122129276344 -1.666770446751757 +-3.592147674562354 -0.6905826718832992 0.003117023538593777 +-3.592147674562354 0.6792380650834748 0.01554283883605434 +-3.592147674562354 0.007236014911530919 -0.6902298068981589 +-3.592147674562354 0.01629389772626366 0.6898338260041208 +-3.592147674562354 2.347526385870982 1.77635683940025e-15 +-3.592147674562354 -2.347526385870983 4.440892098500626e-16 +-3.592147674562354 4.440892098500626e-16 2.347526385870984 +-3.592147674562354 -8.881784197001252e-16 -2.347526385870985 +-3.592147674562354 0.2506333042762985 3.062524012448834 +-3.592147674562354 -0.2752556005858979 3.072188780110338 +-3.592147674562354 3.067045750092873 -0.2506333042762976 +-3.592147674562354 3.072424189192181 0.275255600585897 +-3.592147674562354 -3.062515683026558 -0.2506333042762976 +-3.592147674562354 -3.07212796665486 0.2752556005858975 +-3.592147674562354 -0.2506333042762985 -3.067045750092874 +-3.592147674562354 0.2752556005858979 -3.072424189192182 +-3.592147674562354 -1.637397350659283 -0.4322758495449248 +-3.592147674562354 0.4155545624362484 1.684114511529366 +-3.592147674562354 1.638234177024552 -0.4302895850897914 +-3.592147674562354 -0.4046308146166546 -1.686763771150641 +-3.592147674562354 -4.440892098500626e-16 -4.884981308350689e-15 +2 91 0 29 +656 +657 +658 +659 +660 +661 +662 +663 +664 +665 +666 +667 +668 +669 +670 +671 +672 +673 +674 +675 +676 +677 +678 +679 +680 +681 +682 +683 +684 +-0.922167378260494 0.9221673782604953 3.592147674562354 +0.8119854520994063 -0.8528604330698228 3.592147674562354 +0.9503532198006117 0.9503532198006117 3.592147674562354 +-0.9880728535300074 -0.940144759980639 3.592147674562354 +1.663276818187732 0.3779821385878033 3.592147674562354 +0.3578493946305765 1.643144074230505 3.592147674562354 +-1.64763450681734 -0.3566341017948642 3.592147674562354 +0.3411020711246024 -1.677873397995921 3.592147674562354 +0.0155428388360539 0.6792380650834775 3.592147674562354 +-0.6887262448832185 -0.02271116903861348 3.592147674562354 +0.1604152887191668 -0.1718613832041478 3.592147674562354 +4.440892098500626e-16 2.347526385870985 3.592147674562354 +2.347526385870985 4.440892098500626e-16 3.592147674562354 +-2.347526385870986 4.440892098500626e-16 3.592147674562354 +0 -2.347526385870984 3.592147674562354 +3.067045750092873 0.2506333042762976 3.592147674562354 +3.07242418919218 -0.2752556005858979 3.592147674562354 +-0.2506333042762985 -3.062524012448834 3.592147674562354 +0.275255600585897 -3.072188780110339 3.592147674562354 +-0.2506333042762958 3.062515683026557 3.592147674562354 +0.2752556005858979 3.072127966654862 3.592147674562354 +-3.067045750092874 0.2506333042762985 3.592147674562354 +-3.072424189192182 -0.2752556005858975 3.592147674562354 +-0.4303971266135367 1.636519687545626 3.592147674562354 +1.694171729934125 -0.3657158705866106 3.592147674562354 +-0.4560207193834014 -1.717844070007088 3.592147674562354 +-1.638183808067837 0.4296884316889122 3.592147674562354 +0.8728324788041824 0.1262279366661354 3.592147674562354 +-0.1504173721579618 -0.8830488397940868 3.592147674562354 +2 92 0 29 +685 +686 +687 +688 +689 +690 +691 +692 +693 +694 +695 +696 +697 +698 +699 +700 +701 +702 +703 +704 +705 +706 +707 +708 +709 +710 +711 +712 +713 +0.9221673782604944 3.592147674562354 0.9221673782604936 +-0.9120792263157718 3.592147674562354 -0.954635927593928 +-0.954635927593924 3.592147674562354 0.9120792263157749 +0.9221673782604491 3.592147674562354 -0.922167378260522 +0.3544939373043721 3.592147674562354 1.639788616904295 +1.639788616904295 3.592147674562354 0.3544939373043712 +-1.666770446751775 3.592147674562354 0.391221292763428 +-0.3912212927634187 3.592147674562354 -1.666770446751796 +0.01554283883605656 3.592147674562354 0.6792380650834793 +-0.6792380650835304 3.592147674562354 0.01554283883603125 +-0.01554283883599394 3.592147674562354 -0.6792380650835148 +0.6792380650834775 3.592147674562354 0.01554283883605256 +3.996802888650564e-15 3.592147674562354 2.34752638587098 +0 3.592147674562354 -2.347526385871032 +2.347526385870981 3.592147674562354 4.440892098500626e-16 +-2.347526385871017 3.592147674562354 2.220446049250313e-15 +-0.2506333042762674 3.592147674562354 3.067045750092846 +0.2752556005858975 3.592147674562354 3.072424189192234 +3.06251568302655 3.592147674562354 -0.2506333042762976 +3.072127966654857 3.592147674562354 0.2752556005858975 +-3.067045750092875 3.592147674562354 -0.2506333042762976 +-3.072424189192184 3.592147674562354 0.2752556005858984 +0.2506333042763238 3.592147674562354 -3.062524012448903 +-0.2752556005859463 3.592147674562354 -3.072188780110382 +-0.4346618893111196 3.592147674562354 1.634919257870493 +1.636120228340102 3.592147674562354 -0.4307965858190461 +-1.683250455524083 3.592147674562354 -0.4086577688514832 +0.4264242815977206 3.592147674562354 -1.639332350941051 +2.664535259100376e-15 3.592147674562354 0.007771419418011849 +2 93 0 29 +714 +715 +716 +717 +718 +719 +720 +721 +722 +723 +724 +725 +726 +727 +728 +729 +730 +731 +732 +733 +734 +735 +736 +737 +738 +739 +740 +741 +742 +0.9221673782604944 0.9221673782604944 -3.592147674562354 +-0.9546359275938996 -0.9120792263157953 -3.592147674562354 +0.9120792263157944 -0.9546359275938987 -3.592147674562354 +-0.9221673782604944 0.9221673782604944 -3.592147674562354 +1.639788616904301 0.354493937304373 -3.592147674562354 +0.3544939373043721 1.639788616904301 -3.592147674562354 +0.3912212927634418 -1.666770446751757 -3.592147674562354 +-1.666770446751757 -0.3912212927634418 -3.592147674562354 +0.6792380650834775 0.01554283883605434 -3.592147674562354 +0.01554283883605123 -0.6792380650834788 -3.592147674562354 +-0.6792380650834802 -0.01554283883605123 -3.592147674562354 +0.01554283883605567 0.6792380650834784 -3.592147674562354 +4.440892098500626e-16 2.347526385870985 -3.592147674562354 +2.347526385870982 4.440892098500626e-16 -3.592147674562354 +-2.347526385870986 4.440892098500626e-16 -3.592147674562354 +0 -2.347526385870984 -3.592147674562354 +3.062515683026557 0.2506333042762972 -3.592147674562354 +3.072127966654857 -0.2752556005858975 -3.592147674562354 +-0.2506333042762967 3.067045750092873 -3.592147674562354 +0.2752556005858979 3.07242418919218 -3.592147674562354 +-3.062524012448835 0.250633304276298 -3.592147674562354 +-3.072188780110339 -0.2752556005858966 -3.592147674562354 +-0.250633304276298 -3.067045750092873 -3.592147674562354 +0.275255600585897 -3.072424189192181 -3.592147674562354 +1.634919257870517 -0.4346618893111089 -3.592147674562354 +-0.4307965858190372 1.636120228340124 -3.592147674562354 +-0.408657768851469 -1.683250455524077 -3.592147674562354 +-1.639332350941013 0.4264242815977197 -3.592147674562354 +0.007771419418026504 1.332267629550188e-15 -3.592147674562354 +2 94 0 29 +743 +744 +745 +746 +747 +748 +749 +750 +751 +752 +753 +754 +755 +756 +757 +758 +759 +760 +761 +762 +763 +764 +765 +766 +767 +768 +769 +770 +771 +0.9221673782604953 -3.592147674562354 0.9221673782604944 +-0.9120792263157722 -3.592147674562354 -0.954635927593928 +-0.9546359275939236 -3.592147674562354 0.9120792263157749 +0.9221673782604953 -3.592147674562354 -0.922167378260494 +0.3544939373043721 -3.592147674562354 1.639788616904295 +1.639788616904301 -3.592147674562354 0.354493937304373 +-1.666770446751774 -3.592147674562354 0.3912212927634289 +-0.3912212927634173 -3.592147674562354 -1.66677044675179 +0.01554283883605923 -3.592147674562354 0.6792380650834793 +-0.6792380650835304 -3.592147674562354 0.01554283883603125 +-0.01554283883599661 -3.592147674562354 -0.6792380650835157 +0.6792380650834784 -3.592147674562354 0.01554283883605656 +3.996802888650564e-15 -3.592147674562354 2.34752638587098 +2.347526385870981 -3.592147674562354 4.440892098500626e-16 +-2.347526385871017 -3.592147674562354 2.220446049250313e-15 +-3.108624468950438e-15 -3.592147674562354 -2.34752638587103 +-0.2506333042763003 -3.592147674562354 3.067045750092866 +0.2752556005858962 -3.592147674562354 3.07242418919218 +-3.06252401244883 -3.592147674562354 -0.2506333042762976 +-3.07218878011034 -3.592147674562354 0.2752556005858975 +3.067045750092873 -3.592147674562354 -0.2506333042762967 +3.072424189192181 -3.592147674562354 0.2752556005858975 +0.2506333042762905 -3.592147674562354 -3.062524012448843 +-0.2752556005859486 -3.592147674562354 -3.072188780110379 +-0.4346618893111196 -3.592147674562354 1.634919257870493 +-1.683250455524083 -3.592147674562354 -0.4086577688514832 +1.636120228340127 -3.592147674562354 -0.4307965858190359 +0.4264242815977348 -3.592147674562354 -1.639332350941025 +2.220446049250313e-15 -3.592147674562354 0.007771419418012737 +2 95 0 29 +772 +773 +774 +775 +776 +777 +778 +779 +780 +781 +782 +783 +784 +785 +786 +787 +788 +789 +790 +791 +792 +793 +794 +795 +796 +797 +798 +799 +800 +3.592147674562354 -0.9120792263157962 0.9546359275938978 +3.592147674562354 0.9120792263157926 -0.9546359275938987 +3.592147674562354 0.9546359275938969 0.9120792263157953 +3.592147674562354 -0.9546359275938983 -0.9120792263157971 +3.592147674562354 1.666770446751755 0.3912212927634418 +3.592147674562354 -1.666770446751755 -0.3912212927634435 +3.592147674562354 -0.3912212927634422 1.666770446751756 +3.592147674562354 0.39122129276344 -1.666770446751757 +3.592147674562354 -0.6792380650834757 -0.01554283883605523 +3.592147674562354 0.6792380650834748 0.01554283883605434 +3.592147674562354 -0.01554283883605345 0.6792380650834793 +3.592147674562354 0.01554283883605079 -0.6792380650834806 +3.592147674562354 2.347526385870982 1.77635683940025e-15 +3.592147674562354 -2.347526385870982 4.440892098500626e-16 +3.592147674562354 8.881784197001252e-16 2.347526385870983 +3.592147674562354 -8.881784197001252e-16 -2.347526385870985 +3.592147674562354 0.2506333042762985 3.062524012448834 +3.592147674562354 -0.2752556005858962 3.072188780110337 +3.592147674562354 3.062524012448834 -0.2506333042762985 +3.592147674562354 3.072188780110338 0.2752556005858975 +3.592147674562354 -3.062515683026558 -0.2506333042762972 +3.592147674562354 -3.072127966654859 0.2752556005858979 +3.592147674562354 -0.2506333042762985 -3.067045750092874 +3.592147674562354 0.2752556005858979 -3.072424189192182 +3.592147674562354 -1.683250455524076 0.4086577688514672 +3.592147674562354 0.4086577688514672 1.683250455524076 +3.592147674562354 1.683250455524074 -0.4086577688514681 +3.592147674562354 -0.408657768851469 -1.683250455524078 +3.592147674562354 -8.881784197001252e-16 2.220446049250313e-15 +3 1 0 1 +801 +-2.513672588902969 2.513672588902981 2.45683769321783 +3 2 0 1 +802 +2.515774446909849 2.47862094033712 2.466594179341314 +3 3 0 1 +803 +2.462224826920431 -2.507602658973128 2.507602658973129 +3 4 0 1 +804 +-2.507602658973333 -2.462224826920637 2.507602658973333 +3 5 0 37 +805 +806 +807 +808 +809 +810 +811 +812 +813 +814 +815 +816 +817 +818 +819 +820 +821 +822 +823 +824 +825 +826 +827 +828 +829 +830 +831 +832 +833 +834 +835 +836 +837 +838 +839 +840 +841 +1.554312234475219e-15 -3.108624468950438e-15 1.419003803349028e-15 +-0.6019864961787025 1.391909826950624 0.01421291156576526 +1.006603691314102 -1.12936116893761 -0.01550935646413376 +-0.5019926463935688 -1.347583383343913 0.4898840005364453 +-1.322855755624372 0.3534258958944756 -0.7472358503778858 +0.3070237825724366 -0.7220138794524131 -1.290084408382666 +1.223808392357601 0.647425471067514 0.629835345510554 +-1.338967536002216 -0.1059284552942445 0.6921693291835738 +0.09015925392041726 0.1570666845917202 1.5 +0.8345516095812316 0.9867025438300996 -0.7992056723532944 +-0.9656140774036193 -1.045959196935288 -0.7985976444216489 +1.53713786669156 -0.1486229853453004 -0.7514178340406058 +-0.2201157037324852 0.5238242838664088 -1.557408483694158 +1.052653810669786 -0.6719682189262483 1.204916737430563 +0.2137237325399406 1.368248494660413 1.031804805193109 +0.0769499743953895 -1.716263371964097 -0.5871890476875993 +-1.066264988558218 1.021728776429337 1.085726161965663 +-0.7416201940781084 -0.7139075926203002 1.544192658943194 +0.4950665273694786 1.78452219653519 -0.02005376479665621 +1.790140040098087 -0.335705381112391 0.3693499864585799 +-0.1456429702252552 1.560644040148022 -1.04005326434094 +0.4044645900152823 -1.587147559583325 1.09601575518038 +-1.046564589837998 -0.1929243186628883 -1.683867314704843 +0.8994981966452579 0.2103773350309524 -1.669524846886719 +-1.695243273794634 0.820688454727679 0.1908015577106557 +-1.806194907278228 -0.5347051170120237 -0.2429842671349494 +1.774507202517597 0.8854259404968068 -0.2714747857769649 +-1.205886999817622 1.387516409935493 -0.8381845461490378 +1.13373426451237 0.2976517619960219 1.6270859736673 +-1.538158171249166 -1.153157965191507 0.6526384426726586 +0.9949077596289774 -1.453249593559713 -1.006973565740701 +-0.3467400380985516 1.014242149126385 1.785120287447713 +0.1115262663363888 -0.6008253535851218 0.8171425590526059 +0.2399370357619798 -0.8608061216926041 1.800263457991186 +-1.030721693791127 0.2479472598218396 1.760731648721755 +-0.1446987474713413 -0.3181002486424134 -2.073379611371378 +-0.4423602639115819 -0.1923519722634701 -0.8812423693782934 +3 6 0 1 +842 +-2.512220859968156 2.466809325468398 -2.512220859968155 +3 7 0 1 +843 +2.460025324332078 2.505386986978524 -2.505386986978518 +3 8 0 1 +844 +2.750429768354949 -2.750429768354947 -2.750429768354948 +3 9 0 1 +845 +-2.750429768354948 -2.750429768354957 -2.750429768354948 +3 10 0 0 +$EndNodes +$Elements +10 3463 1 3463 +3 1 4 123 +1 360 382 354 801 +2 354 382 369 801 +3 361 356 375 801 +4 360 354 363 801 +5 370 375 356 801 +6 355 376 381 801 +7 385 379 373 801 +8 376 355 366 801 +9 360 380 382 801 +10 380 360 357 801 +11 381 376 384 801 +12 357 381 380 801 +13 374 366 353 801 +14 353 361 374 801 +15 363 357 360 801 +16 356 358 368 801 +17 375 374 361 801 +18 374 376 366 801 +19 369 382 371 801 +20 356 368 370 801 +21 358 363 354 801 +22 355 381 357 801 +23 359 356 361 801 +24 370 372 375 801 +25 378 375 372 801 +26 384 376 377 801 +27 364 366 355 801 +28 371 382 383 801 +29 373 379 372 801 +30 356 359 358 801 +31 353 359 361 801 +32 360 365 354 382 +33 364 353 366 801 +34 365 369 354 382 +35 379 385 377 801 +36 385 373 383 801 +37 355 381 376 367 +38 355 357 364 801 +39 355 376 366 367 +40 384 380 381 801 +41 372 379 378 801 +42 385 384 377 801 +43 373 371 383 801 +44 383 382 380 801 +45 359 363 358 801 +46 385 383 384 801 +47 368 372 370 801 +48 378 374 375 801 +49 377 376 374 801 +50 371 368 369 801 +51 357 363 364 801 +52 376 381 60 367 +53 377 378 379 801 +54 373 372 371 801 +55 353 364 359 801 +56 375 362 58 370 +57 57 371 369 382 +58 379 372 59 373 +59 381 60 384 376 +60 58 375 370 372 +61 354 46 358 369 +62 375 58 378 372 +63 368 358 46 369 +64 57 369 365 382 +65 374 49 353 366 +66 371 57 383 382 +67 54 380 360 357 +68 360 380 54 382 +69 373 56 385 383 +70 385 377 61 379 +71 357 53 380 381 +72 380 384 383 801 +73 49 374 376 366 +74 60 377 384 376 +75 53 357 355 381 +76 59 378 372 379 +77 363 359 364 801 +78 377 374 378 801 +79 368 371 372 801 +80 384 377 61 385 +81 50 49 353 374 +82 358 46 45 368 +83 383 373 56 371 +84 380 54 53 357 +85 362 58 51 375 +86 52 60 381 367 +87 376 60 48 367 +88 55 57 365 382 +89 369 47 57 365 +90 58 51 4 362 +91 362 44 58 370 +92 60 52 6 367 +93 5 57 47 365 +94 57 5 55 365 +95 6 48 60 367 +96 58 4 44 362 +97 361 362 51 375 +98 366 376 48 367 +99 385 61 3 379 +100 3 56 385 373 +101 379 59 3 373 +102 379 3 385 373 +103 381 355 52 367 +104 383 56 57 371 +105 370 362 44 356 +106 370 44 45 356 +107 369 354 47 365 +108 354 47 46 369 +109 360 55 365 382 +110 377 61 60 384 +111 361 51 50 375 +112 378 58 59 372 +113 52 53 355 381 +114 360 54 55 382 +115 48 49 376 366 +116 375 356 362 370 +117 362 356 375 361 +118 369 358 801 368 +119 369 801 358 354 +120 374 361 50 375 +121 374 50 361 353 +122 368 356 45 358 +123 368 45 356 370 +3 2 4 123 +124 408 387 393 802 +125 403 394 389 802 +126 387 408 415 802 +127 389 414 403 802 +128 406 418 412 802 +129 388 402 409 802 +130 402 388 399 802 +131 396 393 387 802 +132 401 399 386 802 +133 409 402 405 802 +134 393 407 408 802 +135 399 401 402 802 +136 403 386 394 802 +137 401 386 403 802 +138 390 407 393 802 +139 415 408 411 802 +140 394 392 389 802 +141 409 407 390 802 +142 414 417 403 802 +143 390 393 396 802 +144 394 386 392 802 +145 409 405 410 802 +146 397 399 388 802 +147 397 386 399 802 +148 390 388 409 802 +149 412 410 406 802 +150 392 391 389 802 +151 387 391 396 802 +152 415 387 398 408 +153 408 387 398 393 +154 415 411 416 802 +155 418 416 412 802 +156 406 404 418 802 +157 403 417 404 802 +158 400 409 388 402 +159 390 397 388 802 +160 388 400 402 399 +161 406 410 405 802 +162 402 401 405 802 +163 411 412 416 802 +164 418 404 417 802 +165 411 408 407 802 +166 409 410 407 802 +167 414 413 417 802 +168 411 410 412 802 +169 406 405 404 802 +170 404 401 403 802 +171 72 414 389 391 +172 392 396 391 802 +173 414 72 413 391 +174 392 386 397 802 +175 418 417 416 802 +176 397 390 396 802 +177 400 76 409 402 +178 415 416 413 802 +179 394 403 63 386 +180 403 401 63 386 +181 408 78 415 411 +182 402 76 409 405 +183 74 417 403 414 +184 415 71 387 391 +185 71 415 413 391 +186 401 399 64 386 +187 68 393 390 407 +188 415 398 78 408 +189 393 68 408 407 +190 414 403 74 395 +191 412 418 79 416 +192 418 406 75 404 +193 77 410 406 412 +194 399 64 402 401 +195 411 407 410 802 +196 417 404 74 403 +197 405 76 409 410 +198 411 78 415 416 +199 409 67 390 407 +200 404 405 401 802 +201 388 67 390 409 +202 417 413 416 802 +203 396 392 397 802 +204 410 405 77 406 +205 411 412 79 416 +206 404 75 418 417 +207 401 64 63 386 +208 72 71 413 391 +209 67 68 390 407 +210 66 76 409 400 +211 414 74 73 395 +212 398 78 70 415 +213 400 65 76 402 +214 398 69 78 408 +215 403 62 74 395 +216 66 9 76 400 +217 73 74 8 395 +218 62 8 74 395 +219 10 78 70 398 +220 78 10 69 398 +221 76 9 65 400 +222 402 400 65 399 +223 409 388 66 400 +224 418 79 7 412 +225 406 75 7 418 +226 406 7 77 412 +227 7 406 418 412 +228 73 72 414 389 +229 387 398 70 415 +230 387 70 71 415 +231 77 76 405 410 +232 79 78 411 416 +233 404 74 75 417 +234 408 398 69 393 +235 403 63 62 394 +236 389 414 73 395 +237 394 62 403 395 +238 66 67 388 409 +239 69 68 408 393 +240 65 64 402 399 +241 389 403 395 394 +242 389 395 403 414 +243 415 391 802 413 +244 802 391 415 387 +245 414 391 802 389 +246 802 391 414 413 +3 3 4 130 +247 453 447 441 803 +248 429 436 419 803 +249 438 443 424 803 +250 437 436 429 803 +251 432 419 436 803 +252 440 443 438 803 +253 423 444 449 803 +254 432 436 438 803 +255 426 425 450 803 +256 424 432 438 803 +257 429 419 426 803 +258 444 452 449 803 +259 437 450 439 803 +260 446 452 444 803 +261 424 428 432 803 +262 432 428 419 803 +263 430 424 443 803 +264 453 452 447 803 +265 445 443 440 803 +266 421 430 442 803 +267 425 450 422 426 +268 430 428 424 803 +269 439 450 451 803 +270 442 430 443 803 +271 439 453 441 803 +272 447 445 441 803 +273 423 431 444 803 +274 423 449 433 803 +275 431 442 444 803 +276 421 442 431 803 +277 449 448 433 803 +278 425 448 450 803 +279 420 448 425 803 +280 420 433 448 803 +281 428 430 421 803 +282 423 433 427 803 +283 423 427 431 803 +284 449 435 444 423 +285 452 446 447 803 +286 427 421 431 803 +287 420 427 433 803 +288 436 440 438 803 +289 445 440 441 803 +290 439 451 453 803 +291 438 424 443 434 +292 420 425 426 803 +293 437 439 436 803 +294 446 444 442 803 +295 448 449 452 803 +296 441 440 439 803 +297 448 451 450 803 +298 442 443 445 803 +299 419 428 426 803 +300 427 428 421 803 +301 451 452 453 803 +302 446 445 447 803 +303 435 444 96 449 +304 437 450 94 439 +305 447 453 97 452 +306 438 440 93 443 +307 449 444 96 452 +308 446 96 444 452 +309 427 420 426 803 +310 439 450 94 451 +311 443 93 438 434 +312 92 441 447 445 +313 95 439 453 441 +314 436 439 440 803 +315 440 445 93 443 +316 450 94 422 437 +317 426 428 427 803 +318 446 447 97 452 +319 81 419 429 436 +320 451 448 452 803 +321 442 445 446 803 +322 82 419 436 432 +323 92 441 445 440 +324 90 89 420 448 +325 81 82 419 436 +326 451 95 439 453 +327 425 420 89 448 +328 438 82 436 432 +329 429 437 81 436 +330 420 433 90 448 +331 89 450 425 448 +332 431 421 86 442 +333 438 432 424 434 +334 442 421 85 430 +335 449 90 433 448 +336 444 96 87 435 +337 435 96 91 449 +338 422 88 94 450 +339 422 94 80 437 +340 421 85 86 442 +341 86 444 431 442 +342 87 96 14 435 +343 435 87 444 423 +344 91 14 96 435 +345 449 91 435 423 +346 85 442 430 443 +347 443 84 93 434 +348 12 88 94 422 +349 94 80 12 422 +350 93 83 438 434 +351 13 93 84 434 +352 443 424 84 434 +353 13 83 93 434 +354 430 424 84 443 +355 431 444 87 423 +356 449 433 91 423 +357 97 11 453 447 +358 441 11 92 447 +359 453 11 95 441 +360 432 438 83 434 +361 441 453 11 447 +362 439 94 95 451 +363 93 92 445 440 +364 446 97 96 452 +365 450 88 89 425 +366 450 422 88 425 +367 91 90 433 449 +368 437 80 422 429 +369 438 83 82 432 +370 87 86 444 431 +371 81 80 437 429 +372 430 84 85 443 +373 450 803 422 426 +374 450 422 803 437 +375 429 422 803 426 +376 429 803 422 437 +3 4 4 130 +377 482 476 488 804 +378 464 483 455 804 +379 483 464 485 804 +380 457 484 473 804 +381 483 465 455 804 +382 487 485 478 804 +383 472 479 458 804 +384 478 461 460 804 +385 465 483 484 804 +386 457 465 484 804 +387 464 455 461 804 +388 479 472 481 804 +389 484 486 473 804 +390 481 472 475 804 +391 457 463 465 804 +392 463 455 465 804 +393 467 457 473 804 +394 476 482 481 804 +395 478 480 487 804 +396 467 471 456 804 +397 478 460 461 459 +398 467 463 457 804 +399 474 473 486 804 +400 471 467 473 804 +401 488 480 482 804 +402 486 488 476 804 +403 468 472 458 804 +404 466 458 479 804 +405 468 471 472 804 +406 471 468 456 804 +407 479 477 466 804 +408 477 478 460 804 +409 460 454 477 804 +410 477 454 466 804 +411 467 456 463 804 +412 466 462 458 804 +413 468 458 462 804 +414 458 472 470 479 +415 475 476 481 804 +416 462 456 468 804 +417 466 454 462 804 +418 483 485 487 804 +419 487 480 488 804 +420 474 486 476 804 +421 457 484 469 473 +422 461 454 460 804 +423 486 484 483 804 +424 475 472 471 804 +425 481 477 479 804 +426 487 488 486 804 +427 478 477 480 804 +428 471 473 474 804 +429 455 463 461 804 +430 462 463 456 804 +431 482 480 481 804 +432 474 476 475 804 +433 479 472 470 111 +434 113 484 486 473 +435 110 481 476 482 +436 485 487 478 114 +437 111 479 472 481 +438 481 472 111 475 +439 462 454 461 804 +440 113 469 484 473 +441 113 486 474 473 +442 115 482 488 480 +443 476 486 112 488 +444 483 487 486 804 +445 478 487 480 114 +446 462 461 463 804 +447 478 485 114 459 +448 110 481 475 476 +449 107 464 483 455 +450 471 474 475 804 +451 480 477 481 804 +452 483 455 465 108 +453 488 487 115 480 +454 471 100 99 456 +455 104 477 103 454 +456 455 483 107 108 +457 486 476 112 474 +458 467 471 99 456 +459 471 468 100 456 +460 483 464 107 485 +461 460 103 477 454 +462 484 483 465 108 +463 477 104 466 454 +464 484 465 457 469 +465 472 468 100 471 +466 477 478 103 460 +467 479 104 466 477 +468 99 471 467 473 +469 485 114 459 106 +470 472 101 470 111 +471 470 105 479 111 +472 459 114 478 102 +473 101 17 470 111 +474 472 470 101 458 +475 105 470 17 111 +476 470 105 458 479 +477 98 469 113 473 +478 18 114 459 102 +479 459 114 18 106 +480 484 109 469 113 +481 16 98 469 113 +482 98 457 469 473 +483 469 109 16 113 +484 467 457 98 473 +485 465 484 109 469 +486 101 468 472 458 +487 466 458 105 479 +488 476 15 110 482 +489 476 112 15 488 +490 15 482 488 115 +491 482 476 15 488 +492 112 486 474 113 +493 487 115 480 114 +494 111 110 481 475 +495 478 103 460 102 +496 459 478 460 102 +497 100 101 468 472 +498 104 466 105 479 +499 467 98 99 473 +500 464 485 459 106 +501 465 484 108 109 +502 464 107 485 106 +503 478 804 459 461 +504 478 459 804 485 +505 464 459 804 461 +506 464 804 459 485 +3 5 4 636 +507 335 811 216 823 +508 810 334 828 835 +509 293 206 294 824 +510 222 220 221 821 +511 816 828 334 835 +512 238 236 808 815 +513 222 806 220 821 +514 293 294 807 824 +515 243 808 246 820 +516 221 220 819 821 +517 811 819 216 823 +518 265 810 264 815 +519 247 807 820 826 +520 246 820 808 826 +521 806 819 220 821 +522 245 820 246 826 +523 247 820 245 826 +524 206 816 294 824 +525 264 815 810 840 +526 335 811 823 831 +527 238 815 808 820 +528 274 258 817 832 +529 805 808 812 815 +530 805 813 819 821 +531 294 816 807 824 +532 247 245 246 826 +533 810 805 816 828 +534 238 237 236 815 +535 814 816 805 828 +536 805 806 825 841 +537 809 817 258 832 +538 811 805 819 823 +539 808 815 236 834 +540 810 815 265 820 +541 805 819 806 821 +542 805 815 812 841 +543 805 816 811 824 +544 812 813 805 821 +545 243 246 245 820 +546 810 828 816 835 +547 220 305 306 806 +548 811 814 805 823 +549 272 216 335 811 +550 258 817 809 827 +551 344 269 267 808 +552 805 828 810 841 +553 805 814 828 841 +554 264 827 815 840 +555 265 266 264 810 +556 805 814 811 816 +557 815 830 812 841 +558 216 811 272 819 +559 807 816 805 824 +560 806 805 812 829 +561 811 805 824 837 +562 806 819 805 823 +563 806 809 805 829 +564 335 823 214 831 +565 814 817 239 828 +566 239 817 814 825 +567 806 812 805 821 +568 809 812 805 829 +569 807 824 805 837 +570 220 806 306 823 +571 251 827 250 830 +572 805 806 814 825 +573 213 241 211 823 +574 266 810 265 820 +575 811 824 818 837 +576 220 819 806 823 +577 227 337 229 812 +578 808 815 805 820 +579 238 236 242 808 +580 805 820 815 841 +581 350 278 338 825 +582 243 238 808 820 +583 295 807 259 818 +584 812 815 808 830 +585 338 306 305 806 +586 305 220 222 806 +587 806 805 814 823 +588 41 811 203 824 +589 807 808 805 820 +590 350 338 806 825 +591 254 224 317 829 +592 251 815 827 830 +593 247 807 248 820 +594 807 810 805 816 +595 241 814 211 823 +596 243 312 246 808 +597 259 295 341 807 +598 208 334 810 835 +599 224 806 222 821 +600 807 816 294 835 +601 295 293 341 807 +602 331 818 259 826 +603 208 810 334 828 +604 222 224 254 806 +605 263 815 264 827 +606 274 275 258 832 +607 807 805 810 820 +608 275 809 258 832 +609 254 317 806 829 +610 230 315 228 830 +611 254 806 224 829 +612 269 267 808 822 +613 234 236 815 834 +614 269 808 344 834 +615 294 296 293 807 +616 220 306 217 823 +617 805 808 807 837 +618 807 818 824 837 +619 350 338 305 806 +620 806 825 809 832 +621 809 250 827 830 +622 203 811 41 833 +623 259 818 807 826 +624 241 814 823 825 +625 238 242 243 808 +626 805 813 811 819 +627 275 313 258 809 +628 811 287 819 833 +629 807 818 295 824 +630 808 822 267 826 +631 805 810 820 841 +632 41 811 824 833 +633 254 329 350 832 +634 812 821 806 829 +635 239 814 205 825 +636 819 833 345 836 +637 272 811 287 819 +638 245 247 248 820 +639 805 809 806 841 +640 234 815 830 834 +641 220 217 819 823 +642 241 211 814 825 +643 812 830 809 841 +644 293 296 341 807 +645 295 259 302 818 +646 239 205 814 828 +647 352 248 247 807 +648 331 302 259 818 +649 805 813 812 837 +650 817 828 814 841 +651 312 246 808 826 +652 254 806 317 832 +653 809 825 817 832 +654 314 815 238 820 +655 306 806 338 823 +656 805 812 809 841 +657 806 821 224 829 +658 236 237 234 815 +659 233 251 250 830 +660 247 259 807 826 +661 247 259 352 807 +662 259 341 352 807 +663 265 264 263 815 +664 261 267 822 826 +665 41 203 40 824 +666 287 345 819 833 +667 812 830 808 834 +668 315 232 228 829 +669 212 823 814 831 +670 350 305 254 806 +671 227 812 229 829 +672 269 822 808 834 +673 304 229 337 812 +674 300 811 335 831 +675 287 811 346 833 +676 294 296 807 835 +677 272 346 287 811 +678 214 823 212 831 +679 314 237 238 815 +680 284 812 822 834 +681 242 808 236 834 +682 317 806 829 832 +683 810 820 266 835 +684 817 825 274 832 +685 315 809 232 829 +686 811 816 814 831 +687 204 816 814 828 +688 284 822 812 839 +689 229 829 812 830 +690 315 349 232 809 +691 284 822 270 834 +692 254 350 806 832 +693 258 809 313 827 +694 230 250 315 830 +695 298 816 334 835 +696 203 811 290 824 +697 41 824 818 833 +698 293 807 295 824 +699 312 344 267 808 +700 264 810 283 840 +701 254 305 222 806 +702 228 829 229 830 +703 349 313 275 809 +704 250 809 315 830 +705 335 216 214 823 +706 267 312 808 826 +707 338 823 806 825 +708 341 296 352 807 +709 273 274 258 817 +710 811 813 805 837 +711 264 266 283 810 +712 808 830 815 834 +713 312 243 242 808 +714 262 287 345 819 +715 206 324 294 816 +716 313 809 250 827 +717 813 821 812 839 +718 242 344 808 834 +719 805 812 808 837 +720 40 203 290 824 +721 1 813 43 838 +722 239 205 323 825 +723 253 263 264 827 +724 299 818 41 824 +725 216 819 217 823 +726 39 38 816 831 +727 221 220 219 819 +728 232 809 349 832 +729 1 336 43 813 +730 344 312 242 808 +731 262 819 345 836 +732 261 826 822 838 +733 813 345 833 836 +734 243 244 238 820 +735 813 833 819 836 +736 816 824 39 831 +737 808 820 807 826 +738 251 815 263 827 +739 210 814 204 831 +740 248 352 296 807 +741 336 43 813 833 +742 223 224 222 821 +743 229 812 304 830 +744 248 820 807 835 +745 252 313 250 827 +746 43 833 291 838 +747 336 813 345 833 +748 218 272 287 819 +749 240 337 227 812 +750 815 827 810 840 +751 208 351 810 840 +752 221 821 819 836 +753 38 206 39 816 +754 346 811 203 833 +755 811 818 813 837 +756 272 335 300 811 +757 302 295 818 824 +758 216 272 218 819 +759 39 816 206 824 +760 267 256 312 826 +761 43 813 833 838 +762 806 823 814 825 +763 290 824 811 831 +764 313 315 250 809 +765 313 349 315 809 +766 210 204 288 831 +767 269 268 267 822 +768 242 316 344 834 +769 265 815 314 820 +770 300 203 346 811 +771 314 265 343 815 +772 300 290 203 811 +773 299 207 41 818 +774 307 823 338 825 +775 285 266 265 820 +776 247 255 259 826 +777 813 822 812 837 +778 349 809 275 832 +779 314 343 237 815 +780 268 261 267 822 +781 241 823 307 825 +782 343 265 263 815 +783 213 212 214 823 +784 263 251 310 815 +785 302 818 299 824 +786 214 282 335 831 +787 234 815 310 830 +788 241 323 211 825 +789 813 818 811 833 +790 207 41 818 833 +791 273 817 258 827 +792 326 239 817 828 +793 232 829 809 832 +794 300 346 272 811 +795 227 229 228 829 +796 339 813 1 839 +797 351 318 208 810 +798 345 813 336 836 +799 248 807 296 835 +800 302 299 295 824 +801 284 280 822 839 +802 350 329 278 832 +803 337 284 304 812 +804 810 827 815 841 +805 809 827 815 830 +806 221 819 219 836 +807 298 294 816 835 +808 301 818 331 826 +809 274 817 308 825 +810 208 298 334 835 +811 284 270 303 834 +812 240 821 812 829 +813 351 283 318 810 +814 290 811 300 831 +815 279 346 203 833 +816 320 221 219 836 +817 311 245 248 820 +818 291 833 818 838 +819 819 821 813 836 +820 230 228 229 830 +821 284 280 270 822 +822 813 812 822 839 +823 810 828 817 841 +824 814 823 811 831 +825 298 294 324 816 +826 318 332 810 835 +827 810 816 807 835 +828 240 812 227 829 +829 233 234 310 830 +830 310 234 237 815 +831 334 816 37 828 +832 329 254 317 832 +833 300 335 282 831 +834 290 39 824 831 +835 39 202 38 831 +836 219 220 217 819 +837 208 810 332 835 +838 317 829 232 832 +839 304 830 812 834 +840 1 327 813 838 +841 1 813 327 839 +842 212 814 210 831 +843 218 217 216 819 +844 253 251 263 827 +845 277 274 825 832 +846 815 820 810 841 +847 37 816 204 828 +848 211 210 212 814 +849 241 307 347 825 +850 336 345 200 833 +851 239 308 817 825 +852 222 221 223 821 +853 343 310 237 815 +854 260 301 331 826 +855 336 200 43 833 +856 812 821 240 839 +857 255 331 259 826 +858 204 814 205 828 +859 284 812 240 839 +860 327 822 813 838 +861 810 840 827 841 +862 283 810 351 840 +863 316 242 236 834 +864 347 307 338 825 +865 252 250 251 827 +866 284 304 812 834 +867 813 822 327 839 +868 808 826 807 837 +869 240 225 821 829 +870 325 232 349 832 +871 811 824 816 831 +872 208 332 298 835 +873 337 240 284 812 +874 288 204 37 816 +875 812 822 808 837 +876 215 216 217 823 +877 38 288 816 831 +878 224 821 225 829 +879 294 297 296 835 +880 818 826 301 838 +881 248 311 820 835 +882 209 1 43 838 +883 810 828 208 840 +884 270 822 269 834 +885 343 263 310 815 +886 336 1 339 813 +887 243 245 244 820 +888 277 276 274 832 +889 314 238 244 820 +890 306 338 307 823 +891 261 256 267 826 +892 239 201 205 828 +893 811 819 813 833 +894 223 225 224 821 +895 279 203 41 833 +896 260 301 826 838 +897 292 206 293 824 +898 322 210 288 831 +899 318 332 208 810 +900 322 212 210 831 +901 317 232 325 832 +902 213 307 241 823 +903 273 258 257 827 +904 258 313 252 827 +905 326 201 239 828 +906 210 205 204 814 +907 207 42 41 833 +908 312 256 246 826 +909 240 225 342 821 +910 287 262 218 819 +911 277 274 308 825 +912 38 324 206 816 +913 339 836 813 839 +914 255 247 246 826 +915 214 212 282 831 +916 302 331 301 818 +917 299 41 40 824 +918 342 240 821 839 +919 283 266 318 810 +920 210 211 205 814 +921 234 235 236 834 +922 290 39 40 824 +923 298 324 334 816 +924 326 828 817 840 +925 339 336 813 836 +926 807 820 810 835 +927 36 334 37 828 +928 292 293 295 824 +929 339 1 327 839 +930 225 240 227 829 +931 291 818 301 838 +932 37 334 324 816 +933 208 348 351 840 +934 215 214 216 823 +935 285 265 314 820 +936 808 822 812 834 +937 818 824 811 833 +938 308 274 273 817 +939 231 830 304 834 +940 342 289 240 839 +941 302 207 299 818 +942 809 830 815 841 +943 35 828 326 840 +944 249 821 221 836 +945 38 202 288 831 +946 42 291 43 833 +947 303 270 269 834 +948 338 278 347 825 +949 232 317 226 829 +950 326 817 273 840 +951 230 233 250 830 +952 215 217 306 823 +953 808 822 826 837 +954 812 829 809 830 +955 43 291 209 838 +956 821 836 249 839 +957 231 304 235 834 +958 291 207 301 818 +959 317 224 226 829 +960 284 303 304 834 +961 333 249 836 839 +962 200 345 287 833 +963 269 344 316 834 +964 229 304 231 830 +965 269 270 268 822 +966 288 37 38 816 +967 2 35 326 840 +968 260 328 301 838 +969 219 819 262 836 +970 826 837 822 838 +971 821 813 836 839 +972 809 829 806 832 +973 290 202 39 831 +974 299 292 295 824 +975 292 39 206 824 +976 228 232 226 829 +977 248 340 311 835 +978 218 219 217 819 +979 807 826 818 837 +980 325 349 275 832 +981 285 340 286 835 +982 274 276 275 832 +983 219 262 320 836 +984 326 239 308 817 +985 321 345 336 836 +986 280 284 240 839 +987 809 815 827 841 +988 327 281 339 839 +989 239 323 308 825 +990 225 223 342 821 +991 287 346 279 833 +992 333 289 249 839 +993 221 249 223 821 +994 264 319 827 840 +995 37 324 38 816 +996 268 309 261 822 +997 244 245 311 820 +998 320 249 221 836 +999 348 208 828 840 +1000 273 257 330 827 +1001 818 833 813 838 +1002 261 822 309 838 +1003 253 264 319 827 +1004 213 214 215 823 +1005 36 37 204 828 +1006 248 296 297 835 +1007 208 334 36 828 +1008 302 301 207 818 +1009 348 2 330 840 +1010 291 301 328 838 +1011 226 227 228 829 +1012 261 260 826 838 +1013 271 327 822 839 +1014 817 828 810 840 +1015 226 224 225 829 +1016 219 218 262 819 +1017 314 244 285 820 +1018 348 828 35 840 +1019 306 307 215 823 +1020 201 204 205 828 +1021 1 209 327 838 +1022 308 273 326 817 +1023 309 822 327 838 +1024 290 300 282 831 +1025 270 280 271 822 +1026 271 822 280 839 +1027 330 326 273 840 +1028 231 230 229 830 +1029 329 317 325 832 +1030 215 307 213 823 +1031 809 827 817 841 +1032 282 212 322 831 +1033 257 258 252 827 +1034 253 252 251 827 +1035 201 326 35 828 +1036 298 297 294 835 +1037 226 225 227 829 +1038 339 321 336 836 +1039 292 40 39 824 +1040 299 40 292 824 +1041 318 286 332 835 +1042 231 234 233 830 +1043 308 323 277 825 +1044 276 278 329 832 +1045 348 208 35 828 +1046 270 271 268 822 +1047 256 255 246 826 +1048 244 311 285 820 +1049 241 347 323 825 +1050 333 836 339 839 +1051 261 260 256 826 +1052 330 2 326 840 +1053 249 342 223 821 +1054 810 817 840 841 +1055 260 331 255 826 +1056 321 262 345 836 +1057 277 347 278 825 +1058 271 327 309 822 +1059 316 236 235 834 +1060 231 233 230 830 +1061 279 41 42 833 +1062 333 339 281 839 +1063 322 288 202 831 +1064 818 813 837 838 +1065 261 328 260 838 +1066 268 271 309 822 +1067 280 240 289 839 +1068 271 281 327 839 +1069 822 837 813 838 +1070 309 327 209 838 +1071 276 325 275 832 +1072 204 201 36 828 +1073 303 235 304 834 +1074 35 208 36 828 +1075 43 200 42 833 +1076 252 253 257 827 +1077 269 316 303 834 +1078 329 325 276 832 +1079 340 248 297 835 +1080 333 249 320 836 +1081 827 840 817 841 +1082 260 255 256 826 +1083 277 323 347 825 +1084 818 837 826 838 +1085 202 290 282 831 +1086 200 287 279 833 +1087 2 348 35 840 +1088 281 271 280 839 +1089 320 262 321 836 +1090 348 283 351 840 +1091 201 35 36 828 +1092 264 283 319 840 +1093 200 279 42 833 +1094 282 322 202 831 +1095 289 333 281 839 +1096 286 340 332 835 +1097 303 316 235 834 +1098 319 257 253 827 +1099 320 321 333 836 +1100 209 291 328 838 +1101 261 309 328 838 +1102 340 297 332 835 +1103 298 332 297 835 +1104 333 321 339 836 +1105 319 330 257 827 +1106 330 827 319 840 +1107 289 281 280 839 +1108 309 209 328 838 +1109 348 319 283 840 +1110 348 330 319 840 +1111 825 841 809 817 +1112 809 841 825 806 +1113 825 350 832 806 +1114 825 832 350 278 +1115 211 823 212 213 +1116 212 823 211 814 +1117 835 266 318 286 +1118 835 318 266 810 +1119 829 315 830 228 +1120 829 830 315 809 +1121 831 204 816 814 +1122 831 816 204 288 +1123 841 814 825 805 +1124 841 825 814 817 +1125 825 205 211 814 +1126 825 211 205 323 +1127 310 830 251 233 +1128 251 830 310 815 +1129 839 342 249 821 +1130 839 249 342 289 +1131 833 207 291 818 +1132 833 291 207 42 +1133 835 285 266 286 +1134 835 266 285 820 +1135 832 277 278 825 +1136 832 278 277 276 +1137 835 311 285 340 +1138 835 285 311 820 +1139 834 231 234 235 +1140 834 234 231 830 +1141 273 840 827 330 +1142 827 840 273 817 +3 6 4 125 +1143 491 498 506 842 +1144 491 506 512 842 +1145 491 499 498 842 +1146 500 511 490 842 +1147 490 511 518 842 +1148 504 498 489 842 +1149 499 489 498 842 +1150 509 521 515 842 +1151 506 498 504 842 +1152 491 510 493 842 +1153 500 490 495 842 +1154 500 493 510 842 +1155 512 506 508 842 +1156 512 510 491 842 +1157 492 496 497 842 +1158 500 510 511 842 +1159 504 505 497 489 +1160 504 489 497 842 +1161 497 505 496 489 +1162 518 511 520 842 +1163 517 507 505 842 +1164 504 497 505 842 +1165 499 491 493 842 +1166 514 520 511 842 +1167 498 503 491 506 +1168 505 497 496 842 +1169 500 495 493 842 +1170 500 502 490 511 +1171 512 508 513 842 +1172 517 519 507 842 +1173 520 515 521 842 +1174 518 490 502 511 +1175 503 512 491 506 +1176 517 492 496 501 +1177 515 513 509 842 +1178 509 519 521 842 +1179 517 496 505 501 +1180 490 494 495 842 +1181 497 494 492 842 +1182 520 514 515 842 +1183 508 506 504 842 +1184 509 513 508 842 +1185 519 509 507 842 +1186 507 504 505 842 +1187 499 493 495 842 +1188 520 516 518 842 +1189 514 511 510 842 +1190 512 513 510 842 +1191 517 516 519 842 +1192 509 508 507 842 +1193 498 117 504 489 +1194 499 497 489 842 +1195 498 117 506 504 +1196 513 515 514 842 +1197 520 521 519 842 +1198 493 121 500 510 +1199 118 496 505 489 +1200 497 495 494 842 +1201 505 504 118 489 +1202 507 505 129 517 +1203 520 521 515 133 +1204 512 506 131 508 +1205 132 518 511 520 +1206 121 511 500 510 +1207 514 132 511 520 +1208 129 519 507 517 +1209 506 503 512 131 +1210 518 502 132 511 +1211 505 129 517 501 +1212 125 517 492 494 +1213 125 516 517 494 +1214 509 130 513 515 +1215 509 519 128 521 +1216 508 504 507 842 +1217 508 131 512 513 +1218 495 497 499 842 +1219 516 126 518 494 +1220 490 518 126 494 +1221 515 514 520 133 +1222 514 510 513 842 +1223 520 519 516 842 +1224 130 513 508 509 +1225 118 504 117 489 +1226 519 509 128 507 +1227 121 493 122 510 +1228 125 126 516 494 +1229 124 517 129 501 +1230 503 123 512 131 +1231 502 132 511 120 +1232 518 132 502 127 +1233 506 116 503 131 +1234 505 129 501 119 +1235 124 129 21 501 +1236 20 123 503 131 +1237 116 20 503 131 +1238 501 129 21 119 +1239 132 502 22 120 +1240 132 22 502 127 +1241 515 521 19 133 +1242 521 509 19 128 +1243 509 19 130 515 +1244 19 509 521 515 +1245 507 128 519 129 +1246 512 491 123 503 +1247 492 517 125 124 +1248 122 512 491 123 +1249 513 508 131 130 +1250 498 116 503 506 +1251 505 496 118 119 +1252 511 500 502 120 +1253 492 517 124 501 +1254 520 514 132 133 +1255 518 502 490 127 +1256 116 498 117 506 +1257 496 505 501 119 +1258 126 518 490 127 +1259 500 511 121 120 +1260 496 517 842 492 +1261 842 517 496 505 +1262 510 491 122 512 +1263 510 122 491 493 +1264 494 518 842 490 +1265 494 842 518 516 +1266 494 517 842 516 +1267 494 842 517 492 +3 7 4 129 +1268 544 556 550 843 +1269 532 524 551 843 +1270 526 541 546 843 +1271 532 551 553 843 +1272 546 541 543 843 +1273 525 547 552 843 +1274 552 547 549 843 +1275 553 554 540 843 +1276 523 529 528 843 +1277 549 555 552 843 +1278 546 535 526 843 +1279 526 536 541 843 +1280 556 555 550 843 +1281 536 539 541 843 +1282 539 536 523 843 +1283 543 548 546 843 +1284 532 529 524 843 +1285 535 545 522 843 +1286 554 542 540 843 +1287 535 546 545 843 +1288 550 548 544 843 +1289 544 554 556 843 +1290 525 534 547 843 +1291 552 533 525 843 +1292 534 545 547 843 +1293 545 534 522 843 +1294 533 552 551 843 +1295 551 524 533 843 +1296 535 530 526 843 +1297 526 530 536 843 +1298 525 537 552 547 +1299 549 550 555 843 +1300 533 531 525 843 +1301 525 531 534 843 +1302 543 541 539 843 +1303 543 544 548 843 +1304 535 522 530 843 +1305 530 523 536 843 +1306 524 531 533 843 +1307 531 522 534 843 +1308 554 544 542 843 +1309 542 539 540 843 +1310 541 526 546 538 +1311 552 555 551 843 +1312 549 547 545 843 +1313 542 544 543 843 +1314 548 545 546 843 +1315 551 554 553 843 +1316 549 548 550 843 +1317 556 554 555 843 +1318 150 552 537 547 +1319 553 147 540 554 +1320 550 556 151 555 +1321 541 546 543 148 +1322 552 150 549 547 +1323 549 555 150 552 +1324 554 147 540 542 +1325 529 531 524 843 +1326 550 548 149 544 +1327 544 554 146 556 +1328 523 530 529 843 +1329 543 539 542 843 +1330 148 541 546 538 +1331 546 548 543 148 +1332 522 531 530 843 +1333 523 528 135 539 +1334 532 524 144 551 +1335 147 540 527 553 +1336 549 550 151 555 +1337 545 548 549 843 +1338 555 554 551 843 +1339 531 529 530 843 +1340 548 543 149 544 +1341 532 144 553 551 +1342 540 135 528 539 +1343 539 135 523 136 +1344 522 139 545 140 +1345 554 544 146 542 +1346 140 545 522 535 +1347 545 534 139 522 +1348 524 144 551 143 +1349 136 523 539 536 +1350 551 524 143 533 +1351 545 534 547 139 +1352 143 552 551 533 +1353 545 140 546 535 +1354 150 537 138 547 +1355 142 150 552 537 +1356 527 134 540 147 +1357 147 527 145 553 +1358 539 136 536 541 +1359 537 26 150 138 +1360 525 138 537 547 +1361 142 26 150 537 +1362 537 552 142 525 +1363 541 536 526 538 +1364 134 527 24 147 +1365 527 145 24 147 +1366 546 141 148 538 +1367 141 546 535 526 +1368 137 541 148 538 +1369 138 525 534 547 +1370 552 533 142 525 +1371 550 149 23 544 +1372 556 23 550 151 +1373 544 146 23 556 +1374 148 141 25 538 +1375 23 550 544 556 +1376 141 546 526 538 +1377 25 137 148 538 +1378 137 536 541 538 +1379 147 554 146 542 +1380 149 543 548 148 +1381 151 549 555 150 +1382 553 527 145 532 +1383 547 534 138 139 +1384 533 143 552 142 +1385 527 134 528 540 +1386 546 140 141 535 +1387 144 553 145 532 +1388 536 136 137 541 +1389 135 528 134 540 +1390 539 528 843 523 +1391 539 843 528 540 +1392 553 527 843 540 +1393 528 527 843 529 +1394 528 843 527 540 +1395 843 532 527 553 +1396 843 527 532 529 +3 8 4 134 +1397 582 560 575 844 +1398 588 561 581 844 +1399 557 580 566 844 +1400 566 580 581 844 +1401 586 566 588 844 +1402 566 561 588 844 +1403 569 581 566 580 +1404 566 586 559 844 +1405 566 581 561 844 +1406 557 569 566 580 +1407 576 564 587 844 +1408 574 575 565 844 +1409 559 586 564 844 +1410 582 575 584 844 +1411 558 564 574 844 +1412 558 574 565 844 +1413 588 568 566 586 +1414 584 575 578 844 +1415 582 565 560 844 +1416 565 575 560 844 +1417 565 582 580 844 +1418 581 590 588 844 +1419 576 574 564 844 +1420 565 580 557 844 +1421 564 586 587 844 +1422 565 571 575 574 +1423 565 582 570 580 +1424 559 566 568 586 +1425 566 559 564 844 +1426 571 565 558 574 +1427 576 587 589 844 +1428 567 564 559 586 +1429 581 569 566 561 +1430 565 557 566 844 +1431 581 583 590 844 +1432 564 567 587 586 +1433 570 557 565 580 +1434 576 589 577 844 +1435 558 565 564 844 +1436 566 568 588 561 +1437 582 584 580 844 +1438 574 578 575 844 +1439 586 588 590 844 +1440 571 560 565 575 +1441 570 582 565 560 +1442 576 564 562 587 +1443 564 565 566 844 +1444 576 577 574 844 +1445 587 586 589 844 +1446 581 580 583 844 +1447 588 572 581 561 +1448 567 562 564 587 +1449 564 563 576 562 +1450 575 560 582 573 +1451 169 572 581 588 +1452 576 587 165 589 +1453 581 169 588 590 +1454 575 584 166 582 +1455 584 575 166 578 +1456 589 165 576 577 +1457 167 579 585 844 +1458 579 164 591 844 +1459 585 591 168 844 +1460 166 575 582 573 +1461 581 169 590 583 +1462 154 558 563 574 +1463 162 567 559 586 +1464 587 562 576 165 +1465 579 591 585 844 +1466 574 577 578 844 +1467 580 584 583 844 +1468 586 590 589 844 +1469 578 579 167 844 +1470 584 167 585 844 +1471 590 168 591 844 +1472 567 162 587 586 +1473 563 576 154 574 +1474 154 153 558 574 +1475 158 157 557 580 +1476 579 577 164 844 +1477 589 591 164 844 +1478 585 168 583 844 +1479 557 570 158 580 +1480 157 569 557 580 +1481 162 559 161 586 +1482 153 571 558 574 +1483 568 161 559 586 +1484 582 158 570 580 +1485 584 578 167 844 +1486 569 157 581 580 +1487 161 568 588 586 +1488 581 156 572 169 +1489 588 169 572 160 +1490 562 155 576 165 +1491 587 163 562 165 +1492 571 153 575 574 +1493 581 572 156 561 +1494 590 583 168 844 +1495 156 30 572 169 +1496 169 30 572 160 +1497 561 588 572 160 +1498 575 571 560 573 +1499 589 164 577 844 +1500 163 29 562 165 +1501 29 155 562 165 +1502 166 582 159 573 +1503 582 570 159 560 +1504 156 569 581 561 +1505 573 166 575 152 +1506 561 568 588 160 +1507 585 167 27 579 +1508 579 164 27 591 +1509 27 168 585 591 +1510 159 28 166 573 +1511 27 585 579 591 +1512 28 166 573 152 +1513 560 159 582 573 +1514 573 575 571 152 +1515 577 579 578 844 +1516 590 591 589 844 +1517 584 585 583 844 +1518 165 589 164 577 +1519 169 590 583 168 +1520 578 166 584 167 +1521 567 163 562 587 +1522 588 568 161 160 +1523 156 569 157 581 +1524 576 563 155 562 +1525 158 570 159 582 +1526 163 567 162 587 +1527 575 153 571 152 +1528 154 576 563 155 +1529 574 563 564 558 +1530 574 564 563 576 +3 9 4 134 +1531 596 611 617 845 +1532 595 616 623 845 +1533 615 601 592 845 +1534 615 616 601 845 +1535 623 621 601 845 +1536 601 595 623 845 +1537 615 601 616 604 +1538 621 594 601 845 +1539 601 616 595 845 +1540 615 601 604 592 +1541 599 622 610 845 +1542 609 611 600 845 +1543 621 599 594 845 +1544 619 617 611 845 +1545 599 609 593 845 +1546 609 600 593 845 +1547 603 623 621 601 +1548 611 613 619 845 +1549 596 617 600 845 +1550 596 600 611 845 +1551 615 600 617 845 +1552 623 616 625 845 +1553 610 609 599 845 +1554 615 592 600 845 +1555 599 621 622 845 +1556 606 600 609 611 +1557 617 600 615 605 +1558 621 594 603 601 +1559 599 601 594 845 +1560 600 609 593 606 +1561 610 622 612 845 +1562 621 599 602 594 +1563 616 601 595 604 +1564 601 600 592 845 +1565 625 616 618 845 +1566 610 598 599 593 +1567 602 599 621 622 +1568 600 592 615 605 +1569 599 609 610 593 +1570 612 622 624 845 +1571 599 593 600 845 +1572 623 603 595 601 +1573 623 625 621 845 +1574 613 611 609 845 +1575 615 617 619 845 +1576 596 600 606 611 +1577 596 600 617 605 +1578 599 597 610 622 +1579 599 600 601 845 +1580 610 612 609 845 +1581 618 616 615 845 +1582 624 622 621 845 +1583 595 608 623 616 +1584 597 599 602 622 +1585 610 598 597 599 +1586 611 596 617 607 +1587 608 187 623 616 +1588 610 185 612 622 +1589 623 187 625 616 +1590 183 611 619 617 +1591 619 613 611 183 +1592 624 612 185 622 +1593 182 614 620 845 +1594 614 184 626 845 +1595 186 620 626 845 +1596 617 183 611 607 +1597 618 625 187 616 +1598 180 621 602 594 +1599 185 610 597 622 +1600 614 626 620 845 +1601 613 609 612 845 +1602 618 615 619 845 +1603 624 621 625 845 +1604 182 620 619 845 +1605 182 613 614 845 +1606 625 186 626 845 +1607 622 621 602 180 +1608 609 172 593 171 +1609 592 175 615 176 +1610 614 612 184 845 +1611 186 618 620 845 +1612 624 626 184 845 +1613 615 592 176 605 +1614 604 175 615 592 +1615 594 180 621 179 +1616 171 593 609 606 +1617 621 594 179 603 +1618 176 617 615 605 +1619 619 613 182 845 +1620 175 615 616 604 +1621 179 623 621 603 +1622 608 187 616 174 +1623 178 608 187 623 +1624 173 185 610 597 +1625 185 597 181 622 +1626 609 171 606 611 +1627 616 595 608 174 +1628 625 618 186 845 +1629 187 608 34 174 +1630 34 178 608 187 +1631 178 608 623 595 +1632 611 606 596 607 +1633 624 184 612 845 +1634 33 185 597 181 +1635 33 185 173 597 +1636 183 617 177 607 +1637 605 617 596 177 +1638 595 616 604 174 +1639 183 170 611 607 +1640 603 178 623 595 +1641 620 31 614 182 +1642 620 626 31 186 +1643 626 614 31 184 +1644 32 183 177 607 +1645 31 614 626 620 +1646 183 32 170 607 +1647 177 617 596 607 +1648 611 170 606 607 +1649 618 619 620 845 +1650 613 612 614 845 +1651 626 624 625 845 +1652 612 184 624 185 +1653 618 625 186 187 +1654 613 619 182 183 +1655 602 181 597 622 +1656 179 178 623 603 +1657 616 175 604 174 +1658 610 598 173 597 +1659 605 176 617 177 +1660 622 602 181 180 +1661 606 171 170 611 +1662 598 172 610 173 +1663 610 593 172 598 +1664 172 593 610 609 +3 10 4 1799 +1665 756 352 341 296 +1666 367 470 641 280 +1667 254 700 317 362 +1668 435 469 261 755 +1669 537 503 726 323 +1670 400 219 365 697 +1671 219 365 320 400 +1672 700 501 317 362 +1673 707 708 732 537 +1674 538 37 572 787 +1675 317 362 639 361 +1676 670 469 261 435 +1677 434 41 398 786 +1678 362 501 317 639 +1679 470 643 367 641 +1680 398 42 434 668 +1681 673 670 674 435 +1682 362 646 645 639 +1683 459 647 761 648 +1684 643 470 644 641 +1685 573 266 608 758 +1686 793 763 422 792 +1687 677 367 669 470 +1688 254 700 362 356 +1689 470 367 669 280 +1690 763 793 422 764 +1691 761 459 648 762 +1692 678 677 669 470 +1693 726 732 537 733 +1694 765 608 766 758 +1695 759 435 760 755 +1696 537 698 503 323 +1697 731 572 730 727 +1698 262 219 320 400 +1699 676 667 675 400 +1700 785 793 422 792 +1701 761 757 459 762 +1702 434 41 42 398 +1703 362 317 254 224 +1704 323 726 537 525 +1705 266 595 608 758 +1706 572 36 538 727 +1707 642 252 502 313 +1708 764 756 422 763 +1709 640 647 459 648 +1710 296 341 756 562 +1711 702 400 701 697 +1712 280 470 641 458 +1713 490 642 502 313 +1714 729 266 608 573 +1715 642 649 650 502 +1716 790 791 395 784 +1717 703 791 395 790 +1718 705 700 706 362 +1719 670 673 469 435 +1720 726 732 503 537 +1721 791 703 395 704 +1722 254 362 224 356 +1723 734 650 649 502 +1724 365 675 667 400 +1725 642 252 607 502 +1726 699 704 395 703 +1727 367 677 643 470 +1728 645 705 706 362 +1729 537 732 707 733 +1730 434 672 671 668 +1731 678 643 677 470 +1732 572 730 794 731 +1733 766 736 608 765 +1734 639 226 317 361 +1735 400 365 701 697 +1736 362 317 224 361 +1737 645 706 646 362 +1738 759 674 435 673 +1739 252 642 250 313 +1740 756 341 352 422 +1741 362 645 501 639 +1742 700 317 496 254 +1743 759 469 435 755 +1744 674 759 435 760 +1745 470 678 643 644 +1746 676 701 400 675 +1747 785 341 429 293 +1748 785 341 293 562 +1749 538 37 36 572 +1750 700 705 501 362 +1751 701 676 400 702 +1752 789 434 788 786 +1753 729 608 736 573 +1754 708 503 732 537 +1755 675 365 701 400 +1756 572 794 538 787 +1757 538 730 572 727 +1758 705 645 501 362 +1759 698 503 277 498 +1760 573 736 729 737 +1761 765 573 608 758 +1762 252 607 502 728 +1763 352 756 563 296 +1764 697 400 219 399 +1765 457 261 755 256 +1766 435 261 670 423 +1767 434 671 398 668 +1768 234 640 236 235 +1769 422 756 341 562 +1770 398 788 434 786 +1771 308 726 503 323 +1772 341 785 422 562 +1773 501 317 496 700 +1774 734 607 735 728 +1775 469 457 261 755 +1776 794 795 538 787 +1777 252 250 642 596 +1778 641 280 458 284 +1779 759 469 673 435 +1780 424 41 42 434 +1781 296 293 341 562 +1782 790 395 703 527 +1783 527 790 395 784 +1784 525 323 726 239 +1785 561 572 37 787 +1786 758 266 595 285 +1787 224 317 226 361 +1788 214 699 395 527 +1789 761 597 647 459 +1790 434 672 788 671 +1791 422 763 562 792 +1792 573 560 266 758 +1793 703 395 699 527 +1794 642 252 596 607 +1795 789 434 672 788 +1796 502 607 734 728 +1797 698 277 503 323 +1798 649 642 607 502 +1799 671 788 434 398 +1800 234 640 597 236 +1801 562 785 422 792 +1802 597 761 757 459 +1803 266 608 595 729 +1804 640 459 236 235 +1805 758 560 266 285 +1806 667 388 262 400 +1807 395 282 527 784 +1808 573 608 736 765 +1809 538 730 794 572 +1810 562 293 785 567 +1811 572 37 36 561 +1812 755 465 457 256 +1813 252 502 500 728 +1814 434 42 424 668 +1815 387 398 41 786 +1816 393 42 398 668 +1817 640 647 597 459 +1818 763 756 422 562 +1819 280 641 240 284 +1820 434 424 41 786 +1821 737 573 736 765 +1822 538 730 795 794 +1823 649 607 735 734 +1824 459 597 640 236 +1825 727 572 36 561 +1826 367 366 355 280 +1827 785 341 422 429 +1828 498 277 698 347 +1829 431 670 261 423 +1830 595 603 758 285 +1831 219 354 365 697 +1832 649 607 734 502 +1833 308 726 323 239 +1834 367 355 669 280 +1835 360 320 365 667 +1836 663 670 261 431 +1837 491 308 726 503 +1838 728 500 252 741 +1839 692 498 698 347 +1840 424 432 41 786 +1841 282 214 395 527 +1842 430 424 42 668 +1843 250 642 596 606 +1844 757 597 459 236 +1845 226 317 653 639 +1846 756 563 296 562 +1847 501 492 317 639 +1848 313 642 250 634 +1849 490 642 313 634 +1850 534 525 726 239 +1851 653 317 492 639 +1852 423 433 755 256 +1853 266 595 604 729 +1854 758 603 750 285 +1855 293 777 785 567 +1856 605 596 252 728 +1857 323 537 698 525 +1858 352 756 769 563 +1859 282 395 389 784 +1860 562 296 293 567 +1861 699 395 394 214 +1862 538 526 37 787 +1863 503 277 308 323 +1864 639 631 226 361 +1865 252 721 605 728 +1866 262 320 321 667 +1867 430 42 680 668 +1868 277 698 347 323 +1869 261 457 469 670 +1870 700 254 691 356 +1871 570 266 729 560 +1872 37 568 561 787 +1873 526 37 36 538 +1874 429 796 785 293 +1875 257 728 252 741 +1876 758 571 560 285 +1877 602 597 757 236 +1878 768 602 757 236 +1879 711 700 496 254 +1880 527 214 699 532 +1881 607 252 596 728 +1882 261 467 670 681 +1883 604 740 266 729 +1884 467 457 261 670 +1885 214 395 335 282 +1886 767 465 755 256 +1887 679 321 320 667 +1888 641 458 466 284 +1889 728 252 721 257 +1890 538 36 526 727 +1891 41 42 398 279 +1892 598 640 597 234 +1893 536 37 526 787 +1894 266 729 560 573 +1895 278 350 118 329 +1896 669 470 280 458 +1897 741 500 252 258 +1898 729 720 266 570 +1899 651 640 234 235 +1900 320 360 679 667 +1901 387 41 398 279 +1902 432 778 41 786 +1903 352 756 422 425 +1904 236 464 757 749 +1905 786 41 387 797 +1906 219 709 354 697 +1907 308 491 277 503 +1908 756 352 748 425 +1909 526 36 535 727 +1910 668 393 42 660 +1911 432 424 434 786 +1912 569 727 36 561 +1913 561 37 36 334 +1914 459 464 757 236 +1915 252 257 741 258 +1916 710 699 214 532 +1917 642 250 654 606 +1918 502 252 500 313 +1919 266 265 595 285 +1920 37 568 787 799 +1921 256 433 755 747 +1922 469 465 457 755 +1923 603 595 265 285 +1924 389 776 282 784 +1925 282 527 784 528 +1926 598 640 234 651 +1927 279 42 398 393 +1928 633 652 641 284 +1929 365 354 219 320 +1930 653 232 317 226 +1931 232 317 492 653 +1932 640 651 632 235 +1933 602 180 594 238 +1934 239 726 534 719 +1935 465 267 457 256 +1936 395 335 394 214 +1937 536 779 37 787 +1938 686 350 118 278 +1939 681 670 261 663 +1940 641 466 633 284 +1941 219 697 399 689 +1942 500 490 502 313 +1943 220 221 46 358 +1944 389 335 395 282 +1945 257 721 728 741 +1946 459 640 460 235 +1947 313 349 490 315 +1948 569 36 727 738 +1949 237 602 238 236 +1950 571 758 770 285 +1951 262 321 661 667 +1952 323 698 533 525 +1953 366 240 355 280 +1954 535 36 718 727 +1955 690 699 394 214 +1956 784 798 282 528 +1957 361 224 362 356 +1958 350 711 118 329 +1959 267 261 457 256 +1960 652 240 641 284 +1961 632 460 640 235 +1962 682 669 355 280 +1963 118 489 278 329 +1964 252 721 253 605 +1965 468 669 280 458 +1966 633 337 652 284 +1967 358 221 46 354 +1968 121 258 273 493 +1969 347 698 712 323 +1970 538 536 526 787 +1971 606 642 596 607 +1972 261 328 663 431 +1973 744 594 180 238 +1974 470 105 641 458 +1975 787 561 572 160 +1976 386 399 64 217 +1977 174 595 608 729 +1978 670 87 435 423 +1979 252 253 721 257 +1980 573 571 560 758 +1981 248 352 769 563 +1982 533 712 698 323 +1983 221 219 709 354 +1984 709 221 46 220 +1985 178 608 595 758 +1986 537 142 698 525 +1987 329 711 496 254 +1988 470 669 101 458 +1989 238 602 768 236 +1990 563 248 352 296 +1991 138 537 726 525 +1992 219 221 709 220 +1993 755 435 91 423 +1994 572 727 156 561 +1995 286 266 560 285 +1996 429 341 295 293 +1997 740 264 604 266 +1998 398 387 279 393 +1999 46 221 709 354 +2000 465 312 267 256 +2001 685 386 64 217 +2002 757 768 236 749 +2003 315 126 490 634 +2004 296 294 293 567 +2005 424 42 41 207 +2006 237 602 594 238 +2007 687 220 46 358 +2008 89 743 420 247 +2009 329 496 317 254 +2010 652 240 366 641 +2011 350 254 305 691 +2012 266 740 283 729 +2013 103 460 627 454 +2014 180 602 768 238 +2015 567 294 777 162 +2016 773 523 322 135 +2017 258 741 121 273 +2018 567 559 294 162 +2019 118 711 350 686 +2020 266 283 720 729 +2021 464 236 242 749 +2022 315 490 313 634 +2023 388 262 661 667 +2024 203 387 41 797 +2025 534 139 714 239 +2026 391 300 72 774 +2027 627 304 460 235 +2028 217 219 399 689 +2029 714 139 534 522 +2030 254 222 691 356 +2031 258 500 252 313 +2032 559 775 294 162 +2033 468 280 669 662 +2034 289 333 357 682 +2035 240 652 337 284 +2036 295 81 772 419 +2037 225 631 50 361 +2038 321 679 661 667 +2039 631 225 226 361 +2040 491 308 739 726 +2041 308 726 239 719 +2042 260 423 261 256 +2043 627 460 304 454 +2044 653 631 226 639 +2045 682 53 357 355 +2046 84 424 430 668 +2047 467 98 457 670 +2048 729 159 570 560 +2049 422 341 352 425 +2050 767 312 465 256 +2051 212 214 527 532 +2052 351 716 557 158 +2053 105 641 458 466 +2054 568 561 787 160 +2055 629 353 49 240 +2056 174 604 595 729 +2057 87 670 431 423 +2058 117 498 692 278 +2059 103 460 632 627 +2060 234 231 651 235 +2061 258 121 500 493 +2062 347 698 692 712 +2063 527 282 214 212 +2064 233 172 598 651 +2065 570 351 557 158 +2066 118 489 686 278 +2067 776 300 389 282 +2068 333 53 357 682 +2069 596 605 177 728 +2070 333 656 53 682 +2071 434 424 84 668 +2072 469 457 98 670 +2073 122 493 308 491 +2074 347 692 498 278 +2075 595 603 178 758 +2076 142 533 698 525 +2077 233 598 172 593 +2078 631 50 629 225 +2079 526 535 141 727 +2080 273 121 493 717 +2081 433 755 91 423 +2082 468 101 669 458 +2083 525 726 138 534 +2084 569 156 727 561 +2085 333 289 281 682 +2086 247 420 89 425 +2087 777 293 294 567 +2088 627 460 632 235 +2089 294 775 777 162 +2090 774 300 72 776 +2091 217 399 219 218 +2092 242 236 238 749 +2093 225 629 631 226 +2094 217 399 64 689 +2095 710 214 213 532 +2096 322 798 773 135 +2097 419 295 81 429 +2098 391 203 300 774 +2099 680 291 42 430 +2100 741 121 500 258 +2101 273 741 121 717 +2102 498 117 489 278 +2103 216 690 394 214 +2104 315 313 250 634 +2105 234 233 598 651 +2106 170 642 606 607 +2107 137 536 538 787 +2108 139 714 239 719 +2109 357 656 53 333 +2110 227 353 629 240 +2111 455 745 312 108 +2112 340 311 571 746 +2113 748 352 756 769 +2114 594 744 314 238 +2115 311 770 571 746 +2116 83 432 434 786 +2117 786 778 41 797 +2118 312 745 455 344 +2119 42 680 668 660 +2120 250 233 230 630 +2121 747 312 256 246 +2122 206 775 777 294 +2123 217 220 219 689 +2124 768 744 180 238 +2125 688 524 241 143 +2126 329 350 711 254 +2127 571 286 560 285 +2128 729 159 560 573 +2129 746 153 571 558 +2130 357 289 682 355 +2131 629 50 353 225 +2132 758 573 571 152 +2133 394 63 690 216 +2134 420 743 255 247 +2135 456 99 467 659 +2136 455 745 242 344 +2137 775 559 294 324 +2138 266 720 283 570 +2139 745 312 242 344 +2140 206 780 292 293 +2141 685 63 386 216 +2142 465 455 312 108 +2143 419 772 295 292 +2144 528 322 523 135 +2145 558 248 746 154 +2146 219 709 697 689 +2147 227 629 225 226 +2148 798 322 528 135 +2149 596 252 250 251 +2150 307 241 347 712 +2151 469 109 465 755 +2152 596 177 607 728 +2153 660 393 42 200 +2154 639 362 51 361 +2155 315 250 230 637 +2156 217 685 386 216 +2157 293 206 777 294 +2158 236 768 238 749 +2159 538 526 141 727 +2160 715 592 253 176 +2161 355 289 682 280 +2162 217 386 399 218 +2163 294 206 775 324 +2164 266 283 318 570 +2165 597 602 237 236 +2166 592 253 176 605 +2167 264 740 283 266 +2168 81 295 772 796 +2169 690 63 685 216 +2170 299 772 419 292 +2171 715 175 264 592 +2172 745 242 107 455 +2173 651 233 172 630 +2174 642 654 250 634 +2175 464 107 242 455 +2176 740 720 283 729 +2177 90 255 743 747 +2178 632 651 231 235 +2179 777 780 206 293 +2180 257 741 258 273 +2181 259 352 341 425 +2182 534 714 522 201 +2183 773 523 288 322 +2184 41 299 778 432 +2185 48 366 367 641 +2186 46 709 220 687 +2187 309 681 261 663 +2188 222 220 687 358 +2189 232 653 228 226 +2190 630 230 250 654 +2191 362 700 44 356 +2192 743 255 90 420 +2193 42 43 291 680 +2194 652 629 49 240 +2195 780 796 292 293 +2196 250 315 634 637 +2197 300 389 72 776 +2198 534 139 239 719 +2199 627 304 104 454 +2200 691 222 45 356 +2201 215 690 685 216 +2202 212 213 214 532 +2203 232 349 315 628 +2204 81 295 796 429 +2205 462 461 269 463 +2206 347 498 277 278 +2207 154 558 248 563 +2208 571 340 746 558 +2209 122 308 739 491 +2210 65 400 697 399 +2211 682 669 280 662 +2212 282 322 212 528 +2213 289 281 682 280 +2214 787 779 37 799 +2215 720 351 570 158 +2216 333 360 54 679 +2217 770 153 571 746 +2218 37 568 799 324 +2219 254 711 700 691 +2220 213 524 144 532 +2221 254 350 711 691 +2222 118 496 489 329 +2223 496 118 711 329 +2224 714 534 239 201 +2225 292 772 295 796 +2226 629 227 636 226 +2227 217 64 685 689 +2228 419 295 299 292 +2229 630 250 233 654 +2230 798 322 282 528 +2231 50 353 225 361 +2232 654 233 630 171 +2233 254 222 305 691 +2234 769 247 352 248 +2235 255 743 747 246 +2236 355 240 289 280 +2237 220 709 219 689 +2238 351 720 716 158 +2239 172 233 593 630 +2240 559 775 161 324 +2241 304 466 104 454 +2242 721 715 253 176 +2243 561 568 37 334 +2244 670 467 98 681 +2245 656 333 281 682 +2246 256 755 767 747 +2247 774 203 300 290 +2248 767 312 256 747 +2249 747 256 255 246 +2250 785 796 777 293 +2251 629 353 227 225 +2252 712 307 241 688 +2253 767 312 747 246 +2254 603 314 750 285 +2255 213 524 688 144 +2256 143 524 241 533 +2257 233 234 231 651 +2258 176 253 721 605 +2259 70 398 387 786 +2260 633 304 104 627 +2261 776 774 300 290 +2262 467 309 456 659 +2263 798 290 776 282 +2264 202 798 773 322 +2265 42 424 430 207 +2266 433 256 255 747 +2267 717 493 273 122 +2268 255 90 433 747 +2269 776 798 282 784 +2270 429 295 796 293 +2271 503 491 277 498 +2272 212 527 282 528 +2273 340 248 311 746 +2274 355 669 52 367 +2275 203 41 40 797 +2276 422 341 259 429 +2277 66 388 667 400 +2278 175 604 264 592 +2279 634 315 126 628 +2280 394 389 335 395 +2281 391 71 203 774 +2282 85 657 291 421 +2283 285 244 753 311 +2284 798 202 290 282 +2285 524 241 213 688 +2286 780 777 796 293 +2287 743 89 748 247 +2288 203 71 387 797 +2289 397 287 218 272 +2290 390 345 67 388 +2291 712 241 347 323 +2292 261 467 681 309 +2293 45 222 687 358 +2294 215 214 213 710 +2295 304 104 466 633 +2296 365 354 47 697 +2297 250 596 251 606 +2298 345 661 67 388 +2299 352 247 748 425 +2300 285 770 571 311 +2301 222 254 224 356 +2302 430 291 42 207 +2303 217 685 215 306 +2304 659 456 100 271 +2305 594 314 237 238 +2306 754 247 769 248 +2307 636 228 653 226 +2308 492 125 232 653 +2309 690 699 214 710 +2310 687 305 222 691 +2311 264 595 604 266 +2312 459 460 640 102 +2313 422 756 88 425 +2314 424 41 432 207 +2315 309 456 659 271 +2316 316 236 459 235 +2317 234 597 237 236 +2318 651 231 233 630 +2319 772 82 299 419 +2320 309 467 99 659 +2321 393 398 69 668 +2322 421 657 291 328 +2323 42 43 680 660 +2324 312 767 751 246 +2325 748 247 352 769 +2326 460 304 454 235 +2327 43 657 291 680 +2328 175 264 740 715 +2329 216 215 690 214 +2330 38 799 775 324 +2331 726 123 491 503 +2332 653 628 232 228 +2333 636 631 629 226 +2334 395 73 389 784 +2335 105 641 466 633 +2336 160 787 568 799 +2337 687 222 45 691 +2338 171 654 233 606 +2339 300 72 389 391 +2340 36 718 727 738 +2341 233 593 630 171 +2342 421 86 657 328 +2343 603 750 178 758 +2344 561 36 569 334 +2345 797 40 203 774 +2346 688 241 712 143 +2347 717 273 739 122 +2348 685 217 215 216 +2349 399 262 219 218 +2350 652 227 629 240 +2351 240 353 49 366 +2352 241 712 533 323 +2353 365 55 360 667 +2354 99 681 309 659 +2355 430 680 84 668 +2356 63 394 386 216 +2357 743 90 89 420 +2358 288 202 773 322 +2359 206 38 775 324 +2360 232 125 494 628 +2361 285 753 770 311 +2362 464 236 316 242 +2363 386 64 63 685 +2364 638 227 652 337 +2365 349 232 494 628 +2366 726 308 739 719 +2367 312 767 465 108 +2368 654 230 250 637 +2369 663 87 670 431 +2370 698 116 503 498 +2371 387 71 203 391 +2372 689 217 220 693 +2373 523 210 529 530 +2374 314 603 265 285 +2375 658 345 67 390 +2376 501 124 492 639 +2377 245 311 248 746 +2378 523 773 288 136 +2379 756 748 88 425 +2380 102 460 640 632 +2381 288 779 37 536 +2382 604 174 740 729 +2383 533 698 712 142 +2384 778 40 41 797 +2385 746 340 248 558 +2386 630 593 172 171 +2387 749 238 242 752 +2388 104 103 627 454 +2389 696 215 213 710 +2390 681 467 99 309 +2391 767 312 745 108 +2392 522 714 139 140 +2393 776 389 73 784 +2394 691 44 700 356 +2395 770 758 750 285 +2396 300 776 290 282 +2397 433 423 260 256 +2398 634 654 250 637 +2399 144 710 213 532 +2400 456 468 100 271 +2401 720 729 159 570 +2402 531 529 210 530 +2403 307 347 695 712 +2404 558 746 153 154 +2405 328 309 261 663 +2406 571 285 311 340 +2407 244 753 750 285 +2408 654 250 233 606 +2409 799 38 37 324 +2410 692 686 117 278 +2411 89 748 247 425 +2412 320 365 354 360 +2413 161 568 559 324 +2414 236 459 464 316 +2415 502 500 728 120 +2416 255 433 90 420 +2417 357 360 54 333 +2418 45 222 358 356 +2419 338 686 695 347 +2420 295 292 796 293 +2421 299 40 41 778 +2422 775 799 161 324 +2423 247 754 245 248 +2424 177 605 721 728 +2425 354 709 47 697 +2426 248 754 245 746 +2427 125 492 232 494 +2428 797 203 71 774 +2429 490 502 642 127 +2430 433 91 755 747 +2431 686 489 117 278 +2432 326 273 739 725 +2433 242 107 464 749 +2434 653 492 124 639 +2435 723 720 351 283 +2436 362 51 646 639 +2437 727 156 569 738 +2438 718 201 36 535 +2439 685 696 215 306 +2440 717 739 273 725 +2441 334 568 37 324 +2442 248 769 746 154 +2443 202 798 322 282 +2444 224 226 225 361 +2445 714 239 719 326 +2446 669 101 468 662 +2447 119 501 496 700 +2448 628 232 125 653 +2449 41 299 432 207 +2450 687 220 305 693 +2451 786 387 70 797 +2452 696 213 688 710 +2453 353 49 50 629 +2454 779 38 773 288 +2455 328 421 86 431 +2456 181 597 757 602 +2457 163 562 785 567 +2458 660 68 393 200 +2459 123 491 739 726 +2460 779 38 37 799 +2461 726 138 534 719 +2462 749 242 107 745 +2463 633 337 304 627 +2464 248 769 754 746 +2465 657 85 86 421 +2466 633 466 304 284 +2467 648 761 191 647 +2468 703 195 791 790 +2469 216 394 335 214 +2470 693 220 305 306 +2471 54 656 333 679 +2472 776 781 290 798 +2473 291 657 85 680 +2474 638 337 633 627 +2475 257 741 273 717 +2476 788 672 193 671 +2477 678 643 188 677 +2478 214 215 690 710 +2479 135 773 523 136 +2480 117 118 489 686 +2481 394 62 395 699 +2482 202 290 781 798 +2483 314 244 750 285 +2484 121 493 717 122 +2485 660 42 43 200 +2486 779 37 38 288 +2487 631 636 653 226 +2488 772 299 82 778 +2489 701 192 675 676 +2490 763 793 199 792 +2491 759 673 194 674 +2492 527 699 145 532 +2493 347 338 686 278 +2494 536 137 779 787 +2495 345 262 661 388 +2496 40 299 772 778 +2497 227 353 240 225 +2498 247 748 754 769 +2499 344 455 312 267 +2500 170 642 654 606 +2501 706 705 645 189 +2502 685 217 693 306 +2503 662 659 100 271 +2504 80 785 422 429 +2505 459 757 464 106 +2506 761 648 191 762 +2507 704 195 791 703 +2508 672 788 193 789 +2509 643 678 188 644 +2510 335 389 300 282 +2511 72 71 391 774 +2512 687 46 45 358 +2513 628 315 232 228 +2514 718 141 535 727 +2515 355 289 240 366 +2516 241 712 143 533 +2517 307 696 213 688 +2518 233 593 171 606 +2519 68 393 200 390 +2520 710 213 688 144 +2521 770 750 753 285 +2522 592 715 175 176 +2523 557 716 157 158 +2524 639 51 631 361 +2525 240 652 227 337 +2526 154 248 769 563 +2527 699 394 62 690 +2528 54 656 357 333 +2529 604 175 264 740 +2530 570 318 266 560 +2531 468 280 662 271 +2532 636 227 228 226 +2533 793 199 764 763 +2534 701 192 676 702 +2535 304 337 633 284 +2536 674 194 759 760 +2537 328 86 657 663 +2538 220 305 222 687 +2539 189 645 706 646 +2540 686 692 695 347 +2541 264 265 595 266 +2542 662 100 468 271 +2543 288 773 779 136 +2544 40 203 774 290 +2545 711 350 686 694 +2546 500 495 490 313 +2547 593 310 599 600 +2548 82 299 419 432 +2549 393 279 42 200 +2550 222 221 220 358 +2551 81 82 772 419 +2552 107 745 455 108 +2553 493 274 308 491 +2554 495 500 258 313 +2555 772 299 40 292 +2556 751 767 747 246 +2557 524 213 212 532 +2558 496 119 700 711 +2559 234 598 233 593 +2560 527 134 784 528 +2561 456 100 99 659 +2562 53 54 656 357 +2563 596 252 253 605 +2564 279 387 41 203 +2565 712 695 307 688 +2566 571 758 152 770 +2567 744 314 244 750 +2568 197 795 730 794 +2569 736 198 737 765 +2570 649 190 650 734 +2571 196 707 732 733 +2572 730 197 794 731 +2573 190 649 735 734 +2574 198 736 766 765 +2575 732 708 707 196 +2576 273 739 308 326 +2577 220 709 689 693 +2578 744 314 750 179 +2579 213 215 696 307 +2580 242 749 752 745 +2581 608 736 174 729 +2582 638 227 629 652 +2583 68 200 658 390 +2584 489 277 498 278 +2585 136 288 523 536 +2586 337 638 304 627 +2587 714 522 201 140 +2588 524 688 144 143 +2589 572 794 787 160 +2590 635 632 231 627 +2591 253 252 596 251 +2592 598 640 173 597 +2593 563 756 155 562 +2594 367 669 52 677 +2595 312 242 243 745 +2596 238 242 752 243 +2597 304 635 231 627 +2598 668 69 393 660 +2599 221 219 354 320 +2600 179 314 750 603 +2601 498 698 116 692 +2602 674 670 87 435 +2603 496 492 317 501 +2604 67 661 345 658 +2605 792 562 785 163 +2606 757 761 597 181 +2607 790 527 134 784 +2608 242 752 243 745 +2609 696 685 693 306 +2610 83 778 432 786 +2611 719 308 739 326 +2612 255 743 246 247 +2613 105 644 470 641 +2614 707 142 698 537 +2615 686 692 347 278 +2616 669 52 682 355 +2617 754 246 743 247 +2618 231 230 233 630 +2619 314 744 594 179 +2620 557 208 716 351 +2621 422 259 341 425 +2622 676 66 667 400 +2623 793 80 785 422 +2624 757 459 762 106 +2625 248 563 297 296 +2626 388 661 66 667 +2627 49 652 240 366 +2628 568 161 799 324 +2629 228 655 230 229 +2630 281 682 280 662 +2631 719 239 308 326 +2632 421 291 85 430 +2633 679 360 55 667 +2634 127 642 490 634 +2635 255 259 420 247 +2636 778 299 82 432 +2637 455 242 464 344 +2638 228 653 628 636 +2639 267 455 312 465 +2640 680 85 291 430 +2641 669 678 470 101 +2642 326 739 719 725 +2643 733 537 726 138 +2644 647 173 640 597 +2645 178 766 608 758 +2646 703 699 145 527 +2647 763 155 756 562 +2648 727 731 572 156 +2649 745 243 312 751 +2650 333 321 679 339 +2651 661 336 321 345 +2652 333 360 320 249 +2653 760 435 91 755 +2654 654 630 230 637 +2655 321 664 336 661 +2656 648 459 640 102 +2657 764 756 88 422 +2658 526 36 37 204 +2659 718 714 201 140 +2660 563 558 248 297 +2661 718 36 35 738 +2662 352 259 247 425 +2663 714 239 326 201 +2664 109 767 465 755 +2665 680 43 657 683 +2666 161 775 559 162 +2667 180 594 744 179 +2668 325 496 317 329 +2669 658 67 68 390 +2670 304 229 635 627 +2671 311 753 770 746 +2672 799 783 38 775 +2673 38 202 773 288 +2674 456 309 268 271 +2675 663 328 86 431 +2676 201 36 35 718 +2677 48 367 643 641 +2678 316 459 460 235 +2679 244 238 744 752 +2680 120 500 728 741 +2681 740 715 319 723 +2682 256 260 433 255 +2683 43 666 657 683 +2684 345 661 336 658 +2685 738 569 36 208 +2686 338 686 350 694 +2687 201 718 140 535 +2688 70 788 398 786 +2689 179 594 314 603 +2690 490 495 349 313 +2691 73 395 791 784 +2692 467 309 268 456 +2693 716 157 208 557 +2694 398 671 69 668 +2695 305 694 687 693 +2696 723 319 740 283 +2697 745 312 767 751 +2698 136 779 288 536 +2699 679 656 333 339 +2700 228 628 315 637 +2701 259 420 247 425 +2702 311 245 753 746 +2703 310 263 601 600 +2704 775 783 38 206 +2705 124 501 645 639 +2706 119 705 501 700 +2707 692 695 347 712 +2708 502 650 642 127 +2709 782 40 797 774 +2710 163 785 777 567 +2711 673 469 98 670 +2712 664 336 661 658 +2713 662 280 281 271 +2714 728 734 502 120 +2715 200 660 68 658 +2716 657 209 291 328 +2717 522 201 140 535 +2718 310 598 234 593 +2719 768 238 749 752 +2720 303 454 304 235 +2721 675 55 365 667 +2722 366 48 652 641 +2723 274 308 273 493 +2724 229 635 231 304 +2725 317 492 325 232 +2726 123 726 732 503 +2727 738 35 722 716 +2728 253 592 715 319 +2729 242 316 464 344 +2730 262 661 321 345 +2731 658 200 345 390 +2732 47 701 365 697 +2733 243 771 244 245 +2734 349 495 275 313 +2735 399 697 65 689 +2736 264 715 592 319 +2737 230 655 635 229 +2738 426 422 259 429 +2739 503 708 116 698 +2740 725 2 724 742 +2741 65 400 702 697 +2742 535 201 36 204 +2743 319 715 721 257 +2744 672 434 84 668 +2745 630 231 230 635 +2746 106 757 464 749 +2747 757 181 602 768 +2748 699 710 145 532 +2749 363 320 354 360 +2750 704 395 62 699 +2751 326 714 35 722 +2752 526 535 36 204 +2753 722 714 35 718 +2754 218 217 386 216 +2755 333 357 360 249 +2756 338 713 306 307 +2757 204 37 526 536 +2758 319 721 253 257 +2759 337 652 638 633 +2760 305 694 693 306 +2761 683 200 43 336 +2762 494 125 126 628 +2763 257 715 721 330 +2764 499 277 498 489 +2765 336 683 200 658 +2766 208 35 738 716 +2767 628 634 315 637 +2768 773 783 38 779 +2769 785 80 796 429 +2770 597 647 33 761 +2771 790 703 24 527 +2772 264 740 319 283 +2773 35 714 201 718 +2774 310 601 599 600 +2775 247 754 246 245 +2776 496 492 325 317 +2777 562 763 29 792 +2778 251 233 250 606 +2779 262 345 390 388 +2780 351 348 723 283 +2781 296 564 294 567 +2782 774 782 40 290 +2783 388 262 399 218 +2784 420 331 427 426 +2785 343 594 599 601 +2786 753 744 244 750 +2787 205 525 534 239 +2788 304 229 638 337 +2789 287 200 393 390 +2790 278 489 277 276 +2791 328 209 309 663 +2792 636 638 655 229 +2793 238 768 744 752 +2794 655 228 636 229 +2795 457 261 267 268 +2796 213 241 307 688 +2797 215 696 690 710 +2798 657 209 328 663 +2799 321 679 664 661 +2800 222 224 359 356 +2801 778 772 40 782 +2802 287 388 397 218 +2803 721 253 715 319 +2804 796 292 772 780 +2805 320 354 221 363 +2806 267 269 463 268 +2807 729 159 573 737 +2808 244 752 771 243 +2809 213 211 524 212 +2810 261 457 467 268 +2811 288 37 204 536 +2812 751 243 312 246 +2813 351 557 318 570 +2814 730 538 141 727 +2815 200 393 279 287 +2816 657 666 43 209 +2817 454 304 466 284 +2818 499 277 489 276 +2819 208 738 157 716 +2820 273 330 717 725 +2821 709 220 687 693 +2822 388 287 262 218 +2823 394 63 62 690 +2824 648 459 18 762 +2825 793 764 12 422 +2826 177 735 607 728 +2827 237 599 597 234 +2828 434 789 83 786 +2829 2 724 742 723 +2830 275 497 495 494 +2831 173 640 598 651 +2832 738 36 35 208 +2833 756 769 563 155 +2834 230 655 228 637 +2835 743 751 246 754 +2836 304 454 303 284 +2837 748 754 743 247 +2838 638 229 304 627 +2839 720 723 740 283 +2840 784 134 798 528 +2841 685 217 689 693 +2842 683 1 666 336 +2843 43 657 209 291 +2844 232 492 325 494 +2845 426 331 427 428 +2846 636 228 227 229 +2847 782 772 40 292 +2848 389 73 72 776 +2849 691 45 44 356 +2850 469 759 109 755 +2851 709 47 46 354 +2852 228 315 230 637 +2853 683 43 1 336 +2854 752 244 238 243 +2855 264 740 715 319 +2856 125 492 124 653 +2857 687 694 305 691 +2858 271 662 659 327 +2859 754 771 751 245 +2860 630 651 231 635 +2861 208 157 569 557 +2862 326 35 714 201 +2863 427 301 421 428 +2864 263 601 343 310 +2865 706 700 44 362 +2866 251 596 253 605 +2867 285 571 286 340 +2868 681 309 684 663 +2869 202 783 38 773 +2870 744 314 238 244 +2871 739 122 491 123 +2872 310 599 598 593 +2873 495 490 349 494 +2874 321 679 339 664 +2875 263 592 601 600 +2876 741 257 330 717 +2877 209 684 309 663 +2878 430 301 291 207 +2879 775 777 780 206 +2880 782 39 40 290 +2881 229 638 635 627 +2882 713 694 338 306 +2883 359 224 225 361 +2884 779 783 38 799 +2885 696 690 685 215 +2886 271 281 662 327 +2887 309 261 467 268 +2888 631 51 50 361 +2889 39 38 783 206 +2890 326 719 714 725 +2891 713 696 306 307 +2892 303 454 466 284 +2893 717 330 724 725 +2894 800 39 781 783 +2895 800 782 39 780 +2896 781 782 39 800 +2897 783 800 39 780 +2898 607 170 642 649 +2899 360 363 320 249 +2900 338 350 686 278 +2901 666 1 43 209 +2902 564 296 562 567 +2903 535 522 201 204 +2904 693 217 220 306 +2905 754 753 245 746 +2906 751 747 743 246 +2907 137 538 795 787 +2908 66 67 661 388 +2909 758 765 573 152 +2910 223 222 224 359 +2911 89 88 748 425 +2912 660 200 683 658 +2913 364 342 289 366 +2914 323 525 205 239 +2915 43 660 200 683 +2916 103 460 102 632 +2917 724 721 715 330 +2918 660 43 680 683 +2919 458 462 466 284 +2920 695 338 347 307 +2921 318 286 266 560 +2922 40 778 782 797 +2923 36 569 334 208 +2924 664 683 666 336 +2925 277 499 274 276 +2926 359 358 222 356 +2927 359 225 353 361 +2928 458 462 284 270 +2929 387 71 70 797 +2930 330 257 273 717 +2931 38 39 783 202 +2932 203 346 387 391 +2933 602 599 597 237 +2934 637 630 230 635 +2935 682 281 656 665 +2936 798 781 202 773 +2937 636 227 638 229 +2938 645 501 21 705 +2939 336 664 683 658 +2940 783 775 780 206 +2941 681 309 659 684 +2942 604 264 592 263 +2943 432 299 419 428 +2944 463 269 462 268 +2945 752 753 244 771 +2946 157 569 738 208 +2947 268 462 463 456 +2948 455 464 316 344 +2949 782 774 781 290 +2950 234 233 310 593 +2951 310 598 599 234 +2952 426 419 302 428 +2953 39 782 781 290 +2954 230 655 637 635 +2955 522 205 201 204 +2956 390 287 262 388 +2957 637 228 628 636 +2958 753 244 744 752 +2959 284 280 458 270 +2960 303 466 462 284 +2961 281 339 656 665 +2962 117 498 116 692 +2963 264 592 263 253 +2964 262 345 287 390 +2965 664 666 1 336 +2966 495 258 274 275 +2967 339 679 656 664 +2968 52 53 682 355 +2969 780 39 783 206 +2970 119 496 118 711 +2971 349 275 495 494 +2972 782 40 39 292 +2973 635 231 230 229 +2974 327 662 659 665 +2975 292 782 772 780 +2976 90 91 433 747 +2977 268 269 462 270 +2978 139 534 138 719 +2979 638 227 636 629 +2980 393 396 279 287 +2981 280 468 270 271 +2982 599 343 601 310 +2983 428 432 299 207 +2984 287 393 396 390 +2985 299 302 419 428 +2986 534 531 525 205 +2987 468 101 100 662 +2988 319 715 257 330 +2989 289 342 240 366 +2990 274 273 258 493 +2991 292 39 782 780 +2992 491 274 308 277 +2993 599 598 597 234 +2994 303 316 460 235 +2995 394 386 216 272 +2996 490 126 127 634 +2997 563 564 296 562 +2998 791 395 8 704 +2999 330 273 326 725 +3000 718 35 722 738 +3001 276 499 274 275 +3002 751 246 754 245 +3003 330 724 741 717 +3004 281 682 662 665 +3005 529 212 322 528 +3006 656 333 339 281 +3007 718 141 140 535 +3008 666 43 1 683 +3009 783 39 781 202 +3010 462 268 270 456 +3011 69 68 393 660 +3012 327 281 662 665 +3013 426 302 331 428 +3014 305 338 350 694 +3015 65 64 399 689 +3016 684 327 1 209 +3017 694 305 338 306 +3018 268 270 456 271 +3019 713 338 695 307 +3020 771 753 244 245 +3021 323 525 531 205 +3022 566 332 565 564 +3023 781 783 202 773 +3024 558 340 248 297 +3025 365 701 5 675 +3026 10 788 671 398 +3027 497 325 492 494 +3028 563 564 297 296 +3029 462 303 284 270 +3030 500 495 258 493 +3031 360 55 54 679 +3032 274 491 499 277 +3033 231 651 632 635 +3034 745 752 243 751 +3035 104 105 466 633 +3036 223 221 222 358 +3037 303 460 454 235 +3038 330 724 2 723 +3039 213 211 241 524 +3040 774 776 781 290 +3041 121 500 120 741 +3042 289 355 364 366 +3043 427 260 331 255 +3044 496 325 497 329 +3045 179 750 178 603 +3046 712 143 533 142 +3047 430 428 301 207 +3048 503 732 20 708 +3049 291 421 301 430 +3050 655 228 637 636 +3051 665 339 656 664 +3052 666 684 1 209 +3053 157 569 156 738 +3054 695 696 307 688 +3055 22 650 734 502 +3056 269 303 462 270 +3057 216 335 394 272 +3058 463 457 267 268 +3059 326 2 330 725 +3060 714 326 725 722 +3061 696 215 306 307 +3062 780 292 39 206 +3063 592 601 604 263 +3064 386 394 392 272 +3065 353 240 342 366 +3066 363 223 221 249 +3067 245 753 244 311 +3068 208 298 566 334 +3069 264 604 601 263 +3070 456 270 468 271 +3071 460 461 316 303 +3072 605 721 176 177 +3073 570 159 720 158 +3074 249 357 363 364 +3075 387 279 346 203 +3076 345 336 200 658 +3077 299 302 428 207 +3078 713 693 694 306 +3079 600 251 253 605 +3080 319 723 348 283 +3081 86 87 663 431 +3082 495 258 275 313 +3083 269 461 455 463 +3084 389 392 300 391 +3085 363 360 357 249 +3086 210 288 523 322 +3087 564 565 558 297 +3088 388 390 287 397 +3089 425 422 259 426 +3090 529 524 212 532 +3091 491 499 277 498 +3092 327 659 309 684 +3093 349 232 325 494 +3094 271 659 309 327 +3095 495 274 493 499 +3096 338 695 686 694 +3097 209 657 684 663 +3098 205 522 530 204 +3099 348 35 208 716 +3100 251 593 233 606 +3101 527 212 529 528 +3102 522 535 530 204 +3103 160 568 161 799 +3104 778 82 83 432 +3105 569 566 334 208 +3106 175 604 174 740 +3107 361 359 224 356 +3108 396 287 272 346 +3109 696 713 695 307 +3110 767 465 108 109 +3111 348 722 35 716 +3112 566 332 557 565 +3113 461 464 459 316 +3114 343 314 594 603 +3115 325 497 275 494 +3116 211 241 533 323 +3117 593 233 310 251 +3118 357 289 355 364 +3119 136 137 779 536 +3120 2 724 330 725 +3121 280 468 458 270 +3122 227 638 229 337 +3123 152 153 571 770 +3124 295 299 302 419 +3125 777 163 567 162 +3126 753 771 754 245 +3127 145 710 144 532 +3128 211 524 529 531 +3129 491 493 274 499 +3130 225 223 224 359 +3131 346 396 387 391 +3132 529 322 523 528 +3133 274 258 495 493 +3134 652 49 48 366 +3135 684 309 327 209 +3136 331 259 302 426 +3137 363 221 320 249 +3138 684 1 327 666 +3139 693 713 696 306 +3140 569 566 208 557 +3141 396 287 390 397 +3142 399 397 388 218 +3143 339 321 664 336 +3144 489 496 497 329 +3145 665 1 339 664 +3146 260 427 433 255 +3147 259 331 255 420 +3148 1 665 666 664 +3149 251 253 263 600 +3150 338 695 694 713 +3151 298 208 566 332 +3152 326 35 2 722 +3153 301 427 331 428 +3154 655 638 635 229 +3155 84 85 680 430 +3156 467 99 98 681 +3157 460 461 459 316 +3158 497 495 499 275 +3159 300 346 203 391 +3160 354 358 221 363 +3161 461 316 303 269 +3162 211 241 524 533 +3163 564 332 565 297 +3164 358 223 221 363 +3165 389 392 335 300 +3166 332 208 566 557 +3167 464 107 106 749 +3168 222 223 358 359 +3169 601 265 264 263 +3170 353 240 225 342 +3171 327 1 339 665 +3172 563 769 154 155 +3173 598 172 173 651 +3174 596 251 600 605 +3175 170 654 171 606 +3176 665 684 327 666 +3177 344 455 267 463 +3178 684 327 659 665 +3179 364 353 342 366 +3180 290 781 39 202 +3181 259 295 302 426 +3182 346 300 392 391 +3183 343 314 603 265 +3184 346 392 396 391 +3185 565 340 558 297 +3186 33 173 647 597 +3187 792 29 562 163 +3188 703 145 24 527 +3189 761 33 597 181 +3190 24 790 527 134 +3191 180 602 181 768 +3192 763 29 155 562 +3193 394 389 392 335 +3194 316 461 464 455 +3195 426 295 419 429 +3196 793 12 80 422 +3197 648 18 459 102 +3198 762 459 18 106 +3199 764 88 12 422 +3200 276 497 499 275 +3201 684 666 657 209 +3202 531 211 323 205 +3203 349 325 275 494 +3204 396 287 397 272 +3205 386 218 216 272 +3206 594 314 343 237 +3207 428 301 421 430 +3208 568 334 298 324 +3209 454 462 461 303 +3210 646 706 4 362 +3211 335 392 394 272 +3212 523 288 210 204 +3213 208 557 318 351 +3214 752 771 243 751 +3215 6 367 677 643 +3216 205 522 531 530 +3217 301 427 260 331 +3218 497 325 276 329 +3219 495 274 499 275 +3220 208 348 716 351 +3221 455 316 461 344 +3222 276 278 489 329 +3223 211 524 212 529 +3224 225 342 359 364 +3225 455 344 269 463 +3226 796 80 81 429 +3227 342 225 353 364 +3228 665 327 1 666 +3229 302 301 331 428 +3230 558 563 564 297 +3231 463 455 267 465 +3232 528 134 798 135 +3233 339 281 327 665 +3234 564 298 332 297 +3235 289 357 249 364 +3236 396 279 387 393 +3237 225 342 223 359 +3238 599 237 310 234 +3239 340 565 332 297 +3240 328 301 421 427 +3241 325 492 496 497 +3242 325 497 276 275 +3243 323 533 211 531 +3244 529 212 527 532 +3245 332 566 298 564 +3246 287 345 200 390 +3247 359 353 225 364 +3248 8 73 395 791 +3249 537 707 26 733 +3250 21 124 501 645 +3251 395 62 8 704 +3252 249 333 357 289 +3253 564 559 294 567 +3254 267 269 344 463 +3255 523 288 204 536 +3256 34 766 736 608 +3257 572 794 30 731 +3258 705 501 21 119 +3259 742 722 2 348 +3260 599 602 594 237 +3261 431 328 421 427 +3262 457 463 467 268 +3263 223 358 359 363 +3264 392 272 300 346 +3265 592 263 253 600 +3266 20 503 708 116 +3267 344 461 455 269 +3268 400 676 9 702 +3269 569 566 561 334 +3270 22 650 502 127 +3271 20 123 732 503 +3272 14 674 435 760 +3273 470 17 678 644 +3274 365 5 55 675 +3275 502 734 22 120 +3276 788 10 70 398 +3277 218 386 397 272 +3278 671 69 10 398 +3279 701 5 47 365 +3280 4 51 646 362 +3281 251 263 310 600 +3282 559 298 294 324 +3283 566 559 298 564 +3284 722 742 723 348 +3285 557 332 208 318 +3286 469 16 759 673 +3287 561 566 568 334 +3288 529 210 211 531 +3289 533 323 525 531 +3290 428 302 301 207 +3291 600 253 592 605 +3292 348 723 722 716 +3293 343 599 594 237 +3294 212 210 211 529 +3295 287 396 279 346 +3296 237 343 599 310 +3297 387 279 396 346 +3298 297 564 294 296 +3299 489 497 499 276 +3300 301 421 291 328 +3301 454 466 462 303 +3302 332 286 565 340 +3303 723 742 2 348 +3304 523 210 322 529 +3305 676 9 66 400 +3306 731 30 572 156 +3307 674 87 14 435 +3308 26 537 733 138 +3309 470 678 17 101 +3310 766 34 178 608 +3311 644 470 17 105 +3312 91 14 435 760 +3313 463 467 268 456 +3314 594 343 603 601 +3315 276 489 497 329 +3316 289 249 342 364 +3317 424 428 430 207 +3318 5 701 192 675 +3319 788 193 10 671 +3320 643 677 6 188 +3321 16 759 673 194 +3322 601 343 603 265 +3323 259 331 420 426 +3324 26 142 707 537 +3325 397 399 386 218 +3326 419 302 295 426 +3327 367 52 6 677 +3328 608 736 34 174 +3329 253 264 592 319 +3330 30 794 572 160 +3331 400 9 65 702 +3332 303 462 461 269 +3333 530 523 210 204 +3334 706 44 4 362 +3335 265 601 343 263 +3336 705 21 645 189 +3337 367 6 48 643 +3338 431 328 427 260 +3339 13 789 434 672 +3340 795 730 25 197 +3341 198 28 737 765 +3342 650 190 22 734 +3343 196 26 707 733 +3344 294 298 559 564 +3345 204 530 523 536 +3346 310 593 251 600 +3347 298 559 568 324 +3348 794 197 30 731 +3349 649 735 32 190 +3350 736 766 34 198 +3351 708 732 20 196 +3352 25 538 730 795 +3353 737 28 573 765 +3354 607 735 32 649 +3355 428 424 432 207 +3356 210 211 531 205 +3357 676 192 9 702 +3358 13 789 672 193 +3359 14 194 674 760 +3360 17 188 678 644 +3361 189 706 4 646 +3362 334 566 568 298 +3363 195 791 8 704 +3364 249 223 342 364 +3365 427 328 301 260 +3366 427 423 431 260 +3367 434 84 13 672 +3368 98 16 469 673 +3369 270 468 462 456 +3370 267 457 463 465 +3371 212 322 210 529 +3372 571 565 286 340 +3373 600 593 251 606 +3374 530 210 205 204 +3375 272 335 392 300 +3376 340 571 565 558 +3377 83 13 789 434 +3378 461 460 454 303 +3379 759 109 16 469 +3380 363 359 223 364 +3381 24 195 703 790 +3382 191 33 647 761 +3383 557 318 570 565 +3384 433 423 427 260 +3385 722 35 2 348 +3386 199 29 763 792 +3387 531 210 205 530 +3388 735 32 177 607 +3389 573 159 28 737 +3390 1 339 664 336 +3391 386 392 397 272 +3392 538 141 25 730 +3393 392 396 397 272 +3394 316 461 344 269 +3395 286 318 565 560 +3396 538 25 137 795 +3397 249 363 223 364 +3398 294 298 564 297 +3399 272 392 396 346 +3400 28 573 765 152 +3401 565 318 570 560 +3402 531 534 522 205 +3403 32 170 607 649 +3404 533 524 211 531 +3405 648 18 191 762 +3406 199 764 12 793 +3407 462 458 468 270 +3408 425 259 420 426 +3409 318 557 332 565 +3410 359 342 223 364 +3411 560 565 286 571 +3412 603 595 601 265 +3413 600 251 596 606 +3414 535 526 530 204 +3415 559 568 566 298 +3416 526 530 204 536 +3417 286 318 332 565 +3418 755 423 261 435 +3419 755 261 423 256 +3420 320 400 667 262 +3421 667 400 320 365 +3422 708 537 698 503 +3423 698 537 708 707 +3424 366 280 641 240 +3425 641 280 366 367 +3426 262 399 400 388 +3427 400 399 262 219 +3428 679 320 333 360 +3429 679 333 320 321 +3430 691 350 694 305 +3431 691 694 350 711 +3432 235 627 231 632 +3433 235 231 627 304 +3434 308 122 273 493 +3435 273 122 308 739 +3436 534 201 205 239 +3437 205 201 534 522 +3438 259 429 295 426 +3439 295 429 259 341 +3440 628 349 126 494 +3441 628 126 349 315 +3442 490 126 349 494 +3443 490 349 126 315 +3444 570 283 351 720 +3445 570 351 283 318 +3446 751 245 243 246 +3447 243 245 751 771 +3448 351 723 716 720 +3449 351 716 723 348 +3450 264 595 601 604 +3451 264 601 595 265 +3452 255 427 420 331 +3453 255 420 427 433 +3454 330 721 741 724 +3455 330 741 721 257 +3456 431 260 261 328 +3457 261 260 431 423 +3458 725 2 722 326 +3459 722 2 725 742 +3460 330 715 723 319 +3461 330 723 715 724 +3462 348 330 723 319 +3463 348 723 330 2 +$EndElements diff --git a/doc/source/tutorials/lbs/9pebs/xs_files/xs_50.xs b/doc/source/tutorials/lbs/9pebs/xs_files/xs_50.xs new file mode 100644 index 0000000000..53fd1f65b4 --- /dev/null +++ b/doc/source/tutorials/lbs/9pebs/xs_files/xs_50.xs @@ -0,0 +1,47 @@ +NUM_GROUPS 9 +NUM_MOMENTS 2 +SIGMA_T_BEGIN +0 3.810000E-04 +1 5.410000E-04 +2 1.230000E-03 +3 2.970000E-04 +4 1.860000E-04 +5 1.850000E-04 +6 1.850000E-04 +7 1.920000E-04 +8 2.310000E-04 +SIGMA_T_END +TRANSFER_MOMENTS_BEGIN +M_GPRIME_G_VAL 0 0 0 1.460000E-04 +M_GPRIME_G_VAL 0 0 1 1.690000E-04 +M_GPRIME_G_VAL 0 0 2 2.730000E-05 +M_GPRIME_G_VAL 0 1 1 2.840000E-04 +M_GPRIME_G_VAL 0 1 2 2.460000E-04 +M_GPRIME_G_VAL 0 2 2 8.140000E-04 +M_GPRIME_G_VAL 0 2 3 3.990000E-04 +M_GPRIME_G_VAL 0 3 3 1.940000E-04 +M_GPRIME_G_VAL 0 3 4 9.960000E-05 +M_GPRIME_G_VAL 0 4 4 1.690000E-04 +M_GPRIME_G_VAL 0 4 5 1.600000E-05 +M_GPRIME_G_VAL 0 5 5 1.680000E-04 +M_GPRIME_G_VAL 0 5 6 1.520000E-05 +M_GPRIME_G_VAL 0 6 6 1.390000E-04 +M_GPRIME_G_VAL 0 6 7 4.800000E-05 +M_GPRIME_G_VAL 0 7 6 3.630000E-09 +M_GPRIME_G_VAL 0 7 7 1.730000E-04 +M_GPRIME_G_VAL 0 7 8 1.920000E-05 +M_GPRIME_G_VAL 0 8 7 9.050000E-05 +M_GPRIME_G_VAL 0 8 8 1.380000E-04 +M_GPRIME_G_VAL 1 0 0 1.240000E-04 +M_GPRIME_G_VAL 1 0 1 5.580000E-05 +M_GPRIME_G_VAL 1 1 1 2.260000E-04 +M_GPRIME_G_VAL 1 1 2 4.060000E-05 +M_GPRIME_G_VAL 1 2 2 4.500000E-04 +M_GPRIME_G_VAL 1 4 4 2.510000E-05 +M_GPRIME_G_VAL 1 5 5 3.530000E-05 +M_GPRIME_G_VAL 1 6 6 4.220000E-05 +M_GPRIME_G_VAL 1 7 6 3.570000E-09 +M_GPRIME_G_VAL 1 7 7 3.270000E-05 +M_GPRIME_G_VAL 1 8 7 1.280000E-06 +M_GPRIME_G_VAL 1 8 8 2.690000E-05 +TRANSFER_MOMENTS_END diff --git a/doc/source/tutorials/lbs/9pebs/xs_files/xs_p_1.xs b/doc/source/tutorials/lbs/9pebs/xs_files/xs_p_1.xs new file mode 100644 index 0000000000..aa71bb4c5c --- /dev/null +++ b/doc/source/tutorials/lbs/9pebs/xs_files/xs_p_1.xs @@ -0,0 +1,90 @@ +NUM_GROUPS 9 +NUM_MOMENTS 2 +SIGMA_T_BEGIN +0 1.280000E-01 +1 1.390000E-01 +2 1.990000E-01 +3 3.230000E-01 +4 4.010000E-01 +5 4.170000E-01 +6 4.170000E-01 +7 4.200000E-01 +8 4.600000E-01 +SIGMA_T_END +SIGMA_F_BEGIN +0 1.950000E-04 +1 1.220000E-04 +2 6.630000E-05 +3 3.030000E-05 +4 7.270000E-05 +5 6.000000E-04 +6 9.500000E-04 +7 3.620000E-03 +8 1.020000E-02 +SIGMA_F_END +NU_BEGIN +0 3.620000E+00 +1 3.060000E+00 +2 2.610000E+00 +3 2.480000E+00 +4 2.420000E+00 +5 2.420000E+00 +6 2.420000E+00 +7 2.440000E+00 +8 2.430000E+00 +NU_END +CHI_BEGIN +0 9.5065957E-03 +1 1.2469302E-01 +2 6.7034967E-01 +3 1.6160215E-01 +4 3.3816747E-02 +5 3.1821658E-05 +6 0.0000000E+00 +7 0.0000000E+00 +8 0.0000000E+00 +CHI_END +TRANSFER_MOMENTS_BEGIN +M_GPRIME_G_VAL 0 0 0 5.740000E-02 +M_GPRIME_G_VAL 0 0 1 3.900000E-02 +M_GPRIME_G_VAL 0 0 2 2.030000E-02 +M_GPRIME_G_VAL 0 0 3 8.560000E-04 +M_GPRIME_G_VAL 0 0 4 1.420000E-04 +M_GPRIME_G_VAL 0 1 1 8.930000E-02 +M_GPRIME_G_VAL 0 1 2 4.780000E-02 +M_GPRIME_G_VAL 0 1 3 1.260000E-03 +M_GPRIME_G_VAL 0 1 4 1.950000E-04 +M_GPRIME_G_VAL 0 2 2 1.710000E-01 +M_GPRIME_G_VAL 0 2 3 2.650000E-02 +M_GPRIME_G_VAL 0 2 4 4.930000E-05 +M_GPRIME_G_VAL 0 3 3 2.830000E-01 +M_GPRIME_G_VAL 0 3 4 3.900000E-02 +M_GPRIME_G_VAL 0 4 4 3.870000E-01 +M_GPRIME_G_VAL 0 4 5 1.320000E-02 +M_GPRIME_G_VAL 0 5 5 4.000000E-01 +M_GPRIME_G_VAL 0 5 6 1.250000E-02 +M_GPRIME_G_VAL 0 6 5 3.720000E-05 +M_GPRIME_G_VAL 0 6 6 3.700000E-01 +M_GPRIME_G_VAL 0 6 7 3.990000E-02 +M_GPRIME_G_VAL 0 7 6 1.020000E-04 +M_GPRIME_G_VAL 0 7 7 3.980000E-01 +M_GPRIME_G_VAL 0 7 8 1.550000E-02 +M_GPRIME_G_VAL 0 8 7 1.000000E-01 +M_GPRIME_G_VAL 0 8 8 3.450000E-01 +M_GPRIME_G_VAL 1 0 0 4.500000E-02 +M_GPRIME_G_VAL 1 0 3 4.480000E-05 +M_GPRIME_G_VAL 1 0 4 5.590000E-05 +M_GPRIME_G_VAL 1 1 1 4.900000E-02 +M_GPRIME_G_VAL 1 1 3 1.140000E-04 +M_GPRIME_G_VAL 1 1 4 9.540000E-06 +M_GPRIME_G_VAL 1 2 2 2.990000E-02 +M_GPRIME_G_VAL 1 2 4 1.390000E-06 +M_GPRIME_G_VAL 1 3 3 4.450000E-02 +M_GPRIME_G_VAL 1 4 4 2.870000E-02 +M_GPRIME_G_VAL 1 5 5 2.670000E-02 +M_GPRIME_G_VAL 1 6 5 3.200000E-05 +M_GPRIME_G_VAL 1 6 6 3.390000E-02 +M_GPRIME_G_VAL 1 7 6 4.490000E-05 +M_GPRIME_G_VAL 1 7 7 2.230000E-02 +M_GPRIME_G_VAL 1 8 8 1.750000E-02 +TRANSFER_MOMENTS_END diff --git a/doc/source/tutorials/lbs/9pebs/xs_files/xs_p_5.xs b/doc/source/tutorials/lbs/9pebs/xs_files/xs_p_5.xs new file mode 100644 index 0000000000..a017d593d9 --- /dev/null +++ b/doc/source/tutorials/lbs/9pebs/xs_files/xs_p_5.xs @@ -0,0 +1,93 @@ +NUM_GROUPS 9 +NUM_MOMENTS 2 +SIGMA_T_BEGIN +0 1.280000E-01 +1 1.400000E-01 +2 1.990000E-01 +3 3.230000E-01 +4 4.010000E-01 +5 4.160000E-01 +6 4.170000E-01 +7 4.200000E-01 +8 4.600000E-01 +SIGMA_T_END +SIGMA_F_BEGIN +0 2.040000E-04 +1 1.200000E-04 +2 6.700000E-05 +3 3.030000E-05 +4 7.190000E-05 +5 5.980000E-04 +6 9.500000E-04 +7 3.580000E-03 +8 1.030000E-02 +SIGMA_F_END +NU_BEGIN +0 3.660000E+00 +1 3.060000E+00 +2 2.610000E+00 +3 2.470000E+00 +4 2.420000E+00 +5 2.410000E+00 +6 2.410000E+00 +7 2.430000E+00 +8 2.430000E+00 +NU_END +CHI_BEGIN +0 9.0013797E-03 +1 1.2474196E-01 +2 6.7061277E-01 +3 1.6166558E-01 +4 3.3929813E-02 +5 3.8520317E-05 +6 9.9793567E-06 +7 0.0000000E+00 +8 0.0000000E+00 +CHI_END +TRANSFER_MOMENTS_BEGIN +M_GPRIME_G_VAL 0 0 0 5.800000E-02 +M_GPRIME_G_VAL 0 0 1 3.850000E-02 +M_GPRIME_G_VAL 0 0 2 1.940000E-02 +M_GPRIME_G_VAL 0 0 3 8.230000E-04 +M_GPRIME_G_VAL 0 0 4 2.090000E-04 +M_GPRIME_G_VAL 0 1 1 8.950000E-02 +M_GPRIME_G_VAL 0 1 2 4.780000E-02 +M_GPRIME_G_VAL 0 1 3 1.220000E-03 +M_GPRIME_G_VAL 0 1 4 1.940000E-04 +M_GPRIME_G_VAL 0 1 5 1.920000E-07 +M_GPRIME_G_VAL 0 2 2 1.710000E-01 +M_GPRIME_G_VAL 0 2 3 2.660000E-02 +M_GPRIME_G_VAL 0 2 4 4.570000E-05 +M_GPRIME_G_VAL 0 2 5 8.490000E-08 +M_GPRIME_G_VAL 0 3 3 2.830000E-01 +M_GPRIME_G_VAL 0 3 4 3.910000E-02 +M_GPRIME_G_VAL 0 4 4 3.870000E-01 +M_GPRIME_G_VAL 0 4 5 1.310000E-02 +M_GPRIME_G_VAL 0 5 5 4.000000E-01 +M_GPRIME_G_VAL 0 5 6 1.250000E-02 +M_GPRIME_G_VAL 0 6 5 3.280000E-05 +M_GPRIME_G_VAL 0 6 6 3.700000E-01 +M_GPRIME_G_VAL 0 6 7 3.960000E-02 +M_GPRIME_G_VAL 0 7 6 1.000000E-04 +M_GPRIME_G_VAL 0 7 7 3.980000E-01 +M_GPRIME_G_VAL 0 7 8 1.550000E-02 +M_GPRIME_G_VAL 0 8 7 1.000000E-01 +M_GPRIME_G_VAL 0 8 8 3.450000E-01 +M_GPRIME_G_VAL 1 0 0 4.570000E-02 +M_GPRIME_G_VAL 1 0 3 4.550000E-06 +M_GPRIME_G_VAL 1 1 1 4.900000E-02 +M_GPRIME_G_VAL 1 1 3 8.200000E-05 +M_GPRIME_G_VAL 1 1 4 1.530000E-05 +M_GPRIME_G_VAL 1 1 5 1.450000E-07 +M_GPRIME_G_VAL 1 2 2 3.000000E-02 +M_GPRIME_G_VAL 1 2 4 7.210000E-07 +M_GPRIME_G_VAL 1 2 5 2.780000E-08 +M_GPRIME_G_VAL 1 3 3 4.430000E-02 +M_GPRIME_G_VAL 1 4 4 2.870000E-02 +M_GPRIME_G_VAL 1 5 5 2.660000E-02 +M_GPRIME_G_VAL 1 6 5 2.780000E-05 +M_GPRIME_G_VAL 1 6 6 3.400000E-02 +M_GPRIME_G_VAL 1 7 6 4.330000E-05 +M_GPRIME_G_VAL 1 7 7 2.230000E-02 +M_GPRIME_G_VAL 1 8 8 1.760000E-02 +TRANSFER_MOMENTS_END diff --git a/doc/source/tutorials/lbs/index.rst b/doc/source/tutorials/lbs/index.rst index d04dd593cb..b4ca5f3633 100644 --- a/doc/source/tutorials/lbs/index.rst +++ b/doc/source/tutorials/lbs/index.rst @@ -7,4 +7,5 @@ Linear Boltzmann Solver primer/index keigen/index + 9pebs/index src_driven/index