From 4d6c06d4ebae1b342861d2cf937f0ce3bb331d5b Mon Sep 17 00:00:00 2001 From: Jean-Paul Pelteret Date: Thu, 20 Sep 2018 10:11:36 +0200 Subject: [PATCH] Fix GridIn::read_abaqus() for codimension 1 case. --- .../20180920YuxiangWangJean-PaulPelteret | 4 + source/grid/grid_in.cc | 41 +- tests/grid/grid_in_abaqus_02.cc | 115 +++ tests/grid/grid_in_abaqus_02.output | 7 + .../abaqus/3d/codim_1-jagged_surface.inp | 838 ++++++++++++++++++ tests/grid/grids/abaqus/3d/codim_1-tube.inp | 221 +++++ 6 files changed, 1209 insertions(+), 17 deletions(-) create mode 100644 doc/news/changes/minor/20180920YuxiangWangJean-PaulPelteret create mode 100644 tests/grid/grid_in_abaqus_02.cc create mode 100644 tests/grid/grid_in_abaqus_02.output create mode 100644 tests/grid/grids/abaqus/3d/codim_1-jagged_surface.inp create mode 100644 tests/grid/grids/abaqus/3d/codim_1-tube.inp diff --git a/doc/news/changes/minor/20180920YuxiangWangJean-PaulPelteret b/doc/news/changes/minor/20180920YuxiangWangJean-PaulPelteret new file mode 100644 index 0000000000..e2028eeb00 --- /dev/null +++ b/doc/news/changes/minor/20180920YuxiangWangJean-PaulPelteret @@ -0,0 +1,4 @@ +Fixed: The GridIn::read_abaqus() function now correctly reads in grids designed +for the co-dimension 1 case (specifically with dim == 2 and spacedim == 3). +
+(Yuxiang Wang, Jean-Paul Pelteret, 2018/09/20) diff --git a/source/grid/grid_in.cc b/source/grid/grid_in.cc index fbc76f5745..a24174e4ac 100644 --- a/source/grid/grid_in.cc +++ b/source/grid/grid_in.cc @@ -916,7 +916,7 @@ GridIn::read_ucd(std::istream &in, namespace { - template + template class Abaqus_to_UCD { public: @@ -953,12 +953,18 @@ GridIn::read_abaqus(std::istream &in, const bool apply_all_indicators_to_manifolds) { Assert(tria != nullptr, ExcNoTriangulationSelected()); - Assert(dim == 2 || dim == 3, ExcNotImplemented()); + // This implementation has only been verified for: + // - 2d grids with codimension 0 + // - 3d grids with codimension 0 + // - 3d grids with codimension 1 + Assert((spacedim == 2 && dim == spacedim) || + (spacedim == 3 && (dim == spacedim || dim == spacedim - 1)), + ExcNotImplemented()); AssertThrow(in, ExcIO()); // Read in the Abaqus file into an intermediate object // that is to be passed along to the UCD reader - Abaqus_to_UCD abaqus_to_ucd; + Abaqus_to_UCD abaqus_to_ucd; abaqus_to_ucd.read_in_abaqus(in); std::stringstream in_ucd; @@ -3220,12 +3226,12 @@ GridIn::get_format_names() namespace { - template - Abaqus_to_UCD::Abaqus_to_UCD() + template + Abaqus_to_UCD::Abaqus_to_UCD() : tolerance(5e-16) // Used to offset Cubit tolerance error when outputting // value close to zero { - AssertThrow(dim == 2 || dim == 3, ExcNotImplemented()); + AssertThrow(spacedim == 2 || spacedim == 3, ExcNotImplemented()); } // Convert from a string to some other data type @@ -3256,9 +3262,9 @@ namespace return number; } - template + template void - Abaqus_to_UCD::read_in_abaqus(std::istream &input_stream) + Abaqus_to_UCD::read_in_abaqus(std::istream &input_stream) { // References: // http://www.egr.msu.edu/software/abaqus/Documentation/docs/v6.7/books/usb/default.htm?startat=pt01ch02.html @@ -3298,11 +3304,11 @@ namespace if (line[0] == '*') goto cont; - std::vector node(dim + 1); + std::vector node(spacedim + 1); std::istringstream iss(line); char comma; - for (unsigned int i = 0; i < dim + 1; ++i) + for (unsigned int i = 0; i < spacedim + 1; ++i) iss >> node[i] >> comma; node_list.push_back(node); @@ -3572,10 +3578,11 @@ namespace } } - template + template std::vector - Abaqus_to_UCD::get_global_node_numbers(const int face_cell_no, - const int face_cell_face_no) const + Abaqus_to_UCD::get_global_node_numbers( + const int face_cell_no, + const int face_cell_face_no) const { std::vector quad_node_list(GeometryInfo::vertices_per_face); @@ -3669,9 +3676,9 @@ namespace return quad_node_list; } - template + template void - Abaqus_to_UCD::write_out_avs_ucd(std::ostream &output) const + Abaqus_to_UCD::write_out_avs_ucd(std::ostream &output) const { // References: // http://www.dealii.org/developer/doxygen/deal.II/structGeometryInfo.html @@ -3725,7 +3732,7 @@ namespace // Node coordinates output.setf(std::ios::scientific, std::ios::floatfield); - for (unsigned int jj = 1; jj < dim + 1; ++jj) + for (unsigned int jj = 1; jj < spacedim + 1; ++jj) { // invoke tolerance -> set points close to zero equal to zero if (std::abs(node_list[ii][jj]) > tolerance) @@ -3733,7 +3740,7 @@ namespace else output << 0.0 << "\t"; } - if (dim == 2) + if (spacedim == 2) output << 0.0 << "\t"; output << std::endl; diff --git a/tests/grid/grid_in_abaqus_02.cc b/tests/grid/grid_in_abaqus_02.cc new file mode 100644 index 0000000000..df1984c0c8 --- /dev/null +++ b/tests/grid/grid_in_abaqus_02.cc @@ -0,0 +1,115 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2004 - 2017 by the deal.II authors +// +// This file is part of the deal.II library. +// +// The deal.II library is free software; you can use it, redistribute +// it, and/or modify it under the terms of the GNU Lesser General +// Public License as published by the Free Software Foundation; either +// version 2.1 of the License, or (at your option) any later version. +// The full text of the license can be found in the file LICENSE.md at +// the top level directory of deal.II. +// +// --------------------------------------------------------------------- + + + +// Check whether we can read in with the abaqus format +// Specifically, this tests for the codimension 1 case; previously the +// third dimension was being ignored, effectively projecting the mesh to +// the XY plane. + +#include +#include +#include +#include +#include + +#include "../tests.h" + + + +template +void +abaqus_grid(const std::string path_and_name, + const bool output_grid = false, + const std::string gridout_name = "") +{ + Triangulation tria; + GridIn grid_in; + grid_in.attach_triangulation(tria); + std::ifstream input_file(path_and_name); + grid_in.read_abaqus(input_file); + + deallog << " " << tria.n_active_cells() << " active cells" << std::endl; + + int hash = 0; + int index = 0; + for (typename Triangulation::active_cell_iterator c = + tria.begin_active(); + c != tria.end(); + ++c, ++index) + for (unsigned int i = 0; i < GeometryInfo::vertices_per_cell; ++i) + hash += (index * i * c->vertex_index(i)) % (tria.n_active_cells() + 1); + deallog << " hash=" << hash << std::endl; + + if (output_grid) + { + // Write in uncompressed VTK so that its possible to verify by eye + // that the last coordinate of all nodes is non-zero. + std::ofstream fout_vtk(gridout_name + ".vtk"); + GridOutFlags::Vtk flags_vtk; + flags_vtk.compression_level = DataOutBase::VtkFlags::no_compression; + GridOut gridout; + gridout.set_flags(flags_vtk); + gridout.write_vtk(tria, fout_vtk); + } +} + + +int +main() +{ + initlog(); + + try + { + deallog << "codim_1-tube" << std::endl; + abaqus_grid<2, 3>(SOURCE_DIR "/grids/abaqus/3d/codim_1-tube.inp", + false, + "codim_1-tube"); + deallog << "codim_1-jagged_surface" << std::endl; + abaqus_grid<2, 3>(SOURCE_DIR + "/grids/abaqus/3d/codim_1-jagged_surface.inp", + false, + "codim_1-jagged_surface"); + } + catch (std::exception &exc) + { + deallog << std::endl + << std::endl + << "----------------------------------------------------" + << std::endl; + deallog << "Exception on processing: " << std::endl + << exc.what() << std::endl + << "Aborting!" << std::endl + << "----------------------------------------------------" + << std::endl; + return 1; + } + catch (...) + { + deallog << std::endl + << std::endl + << "----------------------------------------------------" + << std::endl; + deallog << "Unknown exception!" << std::endl + << "Aborting!" << std::endl + << "----------------------------------------------------" + << std::endl; + return 1; + }; + + return 0; +} diff --git a/tests/grid/grid_in_abaqus_02.output b/tests/grid/grid_in_abaqus_02.output new file mode 100644 index 0000000000..4e455ff909 --- /dev/null +++ b/tests/grid/grid_in_abaqus_02.output @@ -0,0 +1,7 @@ + +DEAL::codim_1-tube +DEAL:: 96 active cells +DEAL:: hash=13187 +DEAL::codim_1-jagged_surface +DEAL:: 394 active cells +DEAL:: hash=233401 diff --git a/tests/grid/grids/abaqus/3d/codim_1-jagged_surface.inp b/tests/grid/grids/abaqus/3d/codim_1-jagged_surface.inp new file mode 100644 index 0000000000..e9fdd3189d --- /dev/null +++ b/tests/grid/grids/abaqus/3d/codim_1-jagged_surface.inp @@ -0,0 +1,838 @@ +*Part, name=MESH +*Node + 141, 34.5765038, 6.54757595, -5.44943762 + 165, 34.3235245, 6.43613148, -5.61718273 + 166, 34.9154816, 6.35121059, -5.30172157 + 167, 34.7116585, 6.31439734, -5.58330107 + 185, 35.180088, 6.41555166, -5.08476877 + 198, 34.4002953, 6.2120142, -5.70087242 + 199, 34.7464294, 6.00455332, -5.57756519 + 200, 34.9904747, 6.0333662, -5.32775497 + 218, 35.22155, 6.0728693, -5.06722021 + 219, 35.4351501, 6.47886086, -4.85016203 + 232, 34.477066, 5.9878974, -5.78456163 + 233, 34.9950294, 5.6949029, -5.30289698 + 234, 34.746685, 5.67930794, -5.55249643 + 252, 35.2266884, 5.71758509, -5.03206491 + 253, 35.4555016, 6.11660862, -4.8017168 + 254, 35.6650658, 6.53899097, -4.5805521 + 267, 34.4732056, 5.67101574, -5.7636404 + 268, 34.7338715, 5.34977245, -5.53212214 + 269, 34.9816628, 5.34871244, -5.27756548 + 289, 35.241272, 5.35557222, -5.02212858 + 290, 35.4600716, 5.74758863, -4.75630093 + 291, 35.700901, 6.16409349, -4.54051018 + 292, 35.903286, 6.5995903, -4.311306 + 305, 34.4693489, 5.35413456, -5.7427187 + 306, 34.986496, 4.99560022, -5.29194689 + 307, 34.7216263, 5.01674557, -5.55119133 + 328, 35.2729683, 4.99107933, -5.04971123 + 329, 35.5003395, 5.37437344, -4.75965834 + 330, 35.7083778, 5.78489876, -4.48894978 + 331, 35.9484138, 6.21357775, -4.2731595 + 332, 36.146183, 6.65823841, -4.03994989 + 345, 34.4425049, 5.04855824, -5.78111649 + 346, 34.7061348, 4.67484045, -5.58186007 + 347, 35.0015259, 4.63518953, -5.33591938 + 370, 35.3113785, 4.6221118, -5.09872723 + 371, 35.5561638, 5.00211811, -4.79659891 + 372, 35.746273, 5.40205145, -4.47890472 + 373, 35.932663, 5.82893372, -4.19792223 + 374, 36.1417389, 6.26879072, -3.96145558 + 375, 36.3566895, 6.71575975, -3.73165178 + 388, 34.4156609, 4.74298239, -5.8195138 + 389, 35.0138702, 4.2704587, -5.37992287 + 390, 34.6870193, 4.30920506, -5.62281132 + 415, 35.3444977, 4.2395339, -5.15238905 + 416, 35.6212006, 4.62921762, -4.85947514 + 417, 35.8199425, 5.02292728, -4.51857281 + 418, 35.9836121, 5.44080448, -4.18642664 + 419, 36.1577606, 5.87845945, -3.90479636 + 420, 36.3443298, 6.31824398, -3.6511693 + 421, 36.5316086, 6.76344204, -3.39801455 + 434, 34.3702545, 4.36980915, -5.88573503 + 435, 34.6878433, 3.95363903, -5.67292118 + 436, 35.028595, 3.89573145, -5.42297792 + 464, 35.3796692, 3.85138202, -5.21094894 + 465, 35.6896629, 4.25227213, -4.93575239 + 466, 35.9194298, 4.6511445, -4.60269547 + 467, 36.0938072, 5.06366873, -4.24720383 + 468, 36.2464943, 5.49538517, -3.91719222 + 469, 36.373867, 5.92680502, -3.60162854 + 470, 36.4696655, 6.36991167, -3.29443622 + 471, 36.6892776, 6.81544399, -3.05822086 + 484, 34.3248444, 3.99663591, -5.95195627 + 485, 35.0291061, 3.56047845, -5.47468472 + 486, 34.720314, 3.60122609, -5.70528698 + 515, 35.3850822, 3.39882612, -5.30107594 + 516, 35.7402229, 3.85671353, -4.99818563 + 517, 36.0191917, 4.27339792, -4.69596338 + 518, 36.2312851, 4.69091272, -4.34339571 + 519, 36.3846817, 5.12650394, -3.98640037 + 520, 36.4762917, 5.54776764, -3.6181972 + 521, 36.5693092, 5.97475815, -3.28405333 + 522, 36.6937294, 6.41352081, -2.98536015 + 523, 36.8814278, 6.86592293, -2.73688602 + 536, 34.3827209, 3.61679554, -5.99705887 + 537, 34.7449837, 3.2943542, -5.74536037 + 538, 34.9893074, 3.33446336, -5.54736662 + 569, 35.1252975, 3.14373803, -5.5130043 + 570, 35.7939034, 3.45260572, -5.04895115 + 571, 36.083065, 3.88125682, -4.76618481 + 572, 36.3558464, 4.30701685, -4.44664526 + 573, 36.5504608, 4.77336311, -4.09842014 + 574, 36.6598396, 5.22774124, -3.73534727 + 575, 36.724247, 5.61880398, -3.33609796 + 576, 36.7747307, 6.01956177, -2.96302605 + 577, 36.8734856, 6.46036959, -2.64798117 + 578, 37.0786819, 6.92040348, -2.42234993 + 591, 34.4406013, 3.2369554, -6.04216194 + 592, 34.8397789, 3.01253581, -5.72251272 + 629, 35.2497749, 2.84872413, -5.50760794 + 630, 35.545311, 2.98046494, -5.29183149 + 631, 36.1363335, 3.48582363, -4.82289553 + 632, 35.8629608, 3.05648327, -5.06669617 + 633, 36.4344826, 3.90957379, -4.5480299 + 634, 36.7736053, 4.37481451, -4.19605064 + 635, 36.8443375, 4.91314363, -3.8680141 + 636, 36.9311333, 5.32356024, -3.47314715 + 637, 36.9570274, 5.67590475, -3.03753662 + 638, 36.991333, 6.06486702, -2.64738846 + 639, 37.0949135, 6.5158267, -2.33106112 + 640, 37.3335152, 6.99437094, -2.15413857 + 653, 34.5573349, 2.90545177, -6.03548717 + 654, 34.9530602, 2.70113969, -5.73330927 + 692, 35.3650398, 2.52397847, -5.51677465 + 693, 35.6366577, 2.62303782, -5.29159117 + 694, 36.1890297, 3.0937736, -4.85311174 + 695, 36.4783707, 3.50933146, -4.61870956 + 696, 35.9351158, 2.67981935, -5.07601404 + 697, 36.8171387, 3.90876985, -4.34551334 + 698, 37.0175285, 4.63633251, -3.95978522 + 699, 37.1253395, 5.05874252, -3.63011384 + 700, 37.1936493, 5.41773033, -3.20634007 + 701, 37.2099609, 5.74505711, -2.75648355 + 702, 37.2371368, 6.11988306, -2.34562087 + 703, 37.3573723, 6.58751583, -2.04276562 + 704, 37.5768127, 7.08219481, -1.87247705 + 717, 34.6740646, 2.57394791, -6.02881193 + 718, 35.089962, 2.41131186, -5.74389219 + 761, 35.4761391, 2.21255016, -5.52871466 + 762, 35.7350197, 2.28087854, -5.3086853 + 763, 36.2259293, 2.71025538, -4.85390854 + 764, 36.5094032, 3.11261463, -4.65081882 + 765, 36.8263512, 3.52085137, -4.45366144 + 766, 35.9919205, 2.31923842, -5.07716703 + 767, 37.2521553, 3.83627629, -4.18078041 + 768, 37.1707916, 4.22939253, -4.07349014 + 769, 37.2920532, 4.78788233, -3.76453328 + 770, 37.178524, 4.46524954, -3.97418857 + 771, 37.4169388, 5.1857996, -3.39668798 + 772, 37.4605446, 5.50271082, -2.94503236 + 773, 37.4850502, 5.82906246, -2.50148129 + 774, 37.5306206, 6.20384026, -2.08922219 + 775, 37.650753, 6.67487335, -1.78100455 + 776, 37.8848801, 7.1963582, -1.65547621 + 789, 34.8374519, 2.28534865, -5.99589634 + 790, 35.2265015, 2.1140008, -5.74349022 + 836, 35.6044197, 1.90943956, -5.56960297 + 837, 35.8325119, 1.95230496, -5.33725214 + 838, 36.2636452, 2.33501697, -4.85247946 + 839, 36.5321617, 2.71727633, -4.65594292 + 840, 36.8391457, 3.11704063, -4.48010588 + 841, 37.1595764, 3.47337365, -4.31132078 + 842, 36.0500679, 1.96722186, -5.08422899 + 843, 37.621624, 3.91630292, -3.95865488 + 844, 37.5112915, 4.26583338, -3.9156878 + 845, 37.4308815, 4.53658724, -3.83182549 + 846, 37.5926971, 4.90774775, -3.52944446 + 847, 37.7009277, 5.28630114, -3.15523481 + 848, 37.7426262, 5.60184526, -2.71218085 + 849, 37.7780457, 5.92261362, -2.27101445 + 850, 37.8513489, 6.31130648, -1.88278317 + 851, 37.9927864, 6.78119802, -1.5959146 + 865, 35.0008392, 1.99674952, -5.96298122 + 866, 35.3468246, 1.82584417, -5.75279713 + 913, 35.7143097, 1.60310769, -5.60822201 + 914, 35.9337959, 1.62868679, -5.37826395 + 915, 36.2843628, 1.96771622, -4.83602571 + 916, 36.537899, 2.33079958, -4.62564087 + 917, 36.8438683, 2.70518184, -4.47126245 + 918, 37.1558189, 3.09200311, -4.34926939 + 919, 37.51091, 3.4841485, -4.1175456 + 920, 37.369278, 3.3643198, -4.22176075 + 921, 36.1313095, 1.62562501, -5.11918068 + 922, 37.7778206, 3.56715727, -3.95137763 + 923, 37.9509315, 3.97199941, -3.73400187 + 924, 37.8472137, 4.30767679, -3.69717598 + 925, 37.7375221, 4.6215744, -3.63634205 + 926, 37.8818016, 4.98903751, -3.28093553 + 927, 37.9831467, 5.36536551, -2.91464448 + 928, 38.0354309, 5.70540905, -2.4973805 + 929, 38.0847321, 6.01916552, -2.06219244 + 930, 38.1844406, 6.41830492, -1.70326161 + 931, 38.3447723, 6.87557697, -1.42408693 + 933, 38.6390648, 7.3538866, -1.28658402 + 947, 35.1119537, 1.73318994, -5.94817924 + 948, 35.4571152, 1.53777182, -5.77167177 + 1000, 35.8267517, 1.30435264, -5.6855092 + 1001, 36.0759239, 1.32852221, -5.46658325 + 1002, 36.3279381, 1.60805666, -4.84873247 + 1003, 36.526207, 1.95326662, -4.5838418 + 1004, 36.8109894, 2.30953383, -4.3911109 + 1005, 37.1614037, 2.67580652, -4.30546904 + 1006, 37.4291382, 3.06546402, -4.22431183 + 1007, 37.6616325, 3.10243225, -4.09213114 + 1008, 36.2423668, 1.30070066, -5.18671703 + 1009, 37.8943672, 3.17556763, -3.95647097 + 1010, 38.047451, 3.61669111, -3.75140667 + 1011, 38.2383347, 4.00293636, -3.49104285 + 1012, 38.1471252, 4.34767199, -3.44838238 + 1013, 38.0364609, 4.68286371, -3.38976073 + 1014, 38.1881943, 5.07971907, -3.07408667 + 1015, 38.2644882, 5.43390322, -2.67786574 + 1016, 38.3242226, 5.76703691, -2.26464057 + 1017, 38.4008713, 6.11766863, -1.87589931 + 1018, 38.5127258, 6.50654936, -1.52560496 + 1019, 38.6987991, 6.94282913, -1.25349402 + 1020, 39.0165253, 7.39226723, -1.10065258 + 1035, 35.2230682, 1.46963036, -5.93337679 + 1036, 35.5355568, 1.24705398, -5.78820992 + 1089, 35.880085, 0.996387362, -5.67628765 + 1090, 36.1286926, 1.00038087, -5.4910841 + 1091, 36.426754, 1.26673782, -4.91132975 + 1092, 36.5263023, 1.57961416, -4.56797171 + 1093, 36.766983, 1.92784715, -4.32043457 + 1094, 37.1065865, 2.26808643, -4.18352222 + 1095, 37.4459534, 2.64212823, -4.13335705 + 1096, 37.7119064, 2.66522646, -3.97510219 + 1097, 36.3607788, 0.990840793, -5.26596689 + 1098, 37.9738808, 2.77294612, -3.89774895 + 1099, 38.128891, 3.242419, -3.7609911 + 1100, 38.3114548, 3.64429331, -3.51428986 + 1101, 38.4979858, 4.02422762, -3.22951412 + 1102, 38.4249878, 4.38240194, -3.19871616 + 1103, 38.3152847, 4.73014545, -3.13928914 + 1104, 38.4420547, 5.10968876, -2.79672313 + 1105, 38.5453186, 5.48452187, -2.43667531 + 1106, 38.6195793, 5.82714462, -2.04414749 + 1107, 38.7001343, 6.16900063, -1.65780556 + 1108, 38.8269844, 6.55417919, -1.3132273 + 1109, 39.0222549, 6.97055864, -1.02525425 + 1111, 39.3741455, 7.3789835, -0.890217125 + 1126, 35.2592125, 1.17815495, -5.93566895 + 1127, 35.6025009, 0.969341457, -5.80589485 + 1178, 35.8974953, 0.722779572, -5.64702606 + 1179, 36.1671791, 0.69369626, -5.49908257 + 1180, 36.5723038, 0.960073173, -5.01517963 + 1181, 36.6009979, 1.21914911, -4.62321377 + 1182, 36.7478752, 1.54790223, -4.29366446 + 1183, 37.024128, 1.89204621, -4.06279755 + 1184, 37.3987503, 2.23462963, -3.94906616 + 1185, 37.746315, 2.21009731, -3.72494936 + 1186, 36.4521675, 0.684450686, -5.34502697 + 1187, 38.0128746, 2.4533453, -3.72122693 + 1188, 38.1613388, 2.88462591, -3.75062585 + 1189, 38.3665962, 3.28634071, -3.54112267 + 1190, 38.5601921, 3.66566515, -3.26107287 + 1191, 38.7385254, 4.03975344, -2.95715952 + 1192, 38.6644058, 4.39950371, -2.91793323 + 1193, 38.5688972, 4.75799036, -2.86864066 + 1194, 38.699337, 5.13756609, -2.53095317 + 1195, 38.8189583, 5.51429462, -2.18298125 + 1196, 38.879303, 5.83942842, -1.78290939 + 1197, 38.9939003, 6.20494032, -1.43286812 + 1198, 39.1351776, 6.58199644, -1.09190547 + 1199, 39.3576355, 6.97725153, -0.815669775 + 1201, 39.737175, 7.34147787, -0.766311586 + 1217, 35.2953568, 0.886679411, -5.93796158 + 1218, 35.7016754, 0.754672885, -5.74502373 + 1268, 35.8795471, 0.45454669, -5.64280605 + 1269, 36.1664085, 0.394487768, -5.49358749 + 1270, 36.7126389, 0.664177179, -5.13165569 + 1271, 36.7483711, 0.894825816, -4.73348904 + 1272, 36.7897186, 1.1710701, -4.33271408 + 1273, 36.9684906, 1.51088679, -4.01028776 + 1274, 37.2695351, 1.85896921, -3.78462553 + 1275, 37.5268822, 1.83979237, -3.50138974 + 1276, 36.4832306, 0.362274528, -5.36618042 + 1277, 38.004528, 2.05499434, -3.3974278 + 1278, 38.1992912, 2.33019757, -3.48403955 + 1279, 38.1693878, 2.62975645, -3.69874501 + 1280, 38.372345, 2.9296999, -3.55148673 + 1281, 38.6101532, 3.30860591, -3.29880476 + 1282, 38.7982712, 3.6793766, -2.99378705 + 1283, 38.9621048, 4.04635382, -2.67024183 + 1284, 38.8896065, 4.40876198, -2.62923455 + 1285, 38.8054733, 4.77319384, -2.58607316 + 1286, 38.9424438, 5.1479125, -2.25238371 + 1287, 39.0345306, 5.49643087, -1.88202631 + 1288, 39.1430588, 5.85049486, -1.52922773 + 1289, 39.2696571, 6.21334457, -1.18781364 + 1290, 39.4366379, 6.5844841, -0.867383659 + 1291, 39.6989441, 6.97010756, -0.638158381 + 1312, 35.5997391, 0.597805858, -5.78629255 + 1359, 35.8523941, 0.156282887, -5.65318108 + 1361, 36.7997894, 0.345469117, -5.22199011 + 1362, 36.9534874, 0.625772536, -4.88649368 + 1363, 36.9340172, 0.830581963, -4.44659519 + 1364, 36.9824524, 1.12269199, -4.03544521 + 1365, 37.1911392, 1.47659326, -3.7179873 + 1366, 37.7881966, 1.77826571, -3.20693064 + 1367, 37.4182205, 1.44436514, -3.42119527 + 1369, 38.2738266, 1.99300706, -3.12676525 + 1370, 38.4452744, 2.28422928, -3.24255443 + 1371, 38.3351822, 2.60533571, -3.53812885 + 1372, 38.612854, 2.9507401, -3.3136723 + 1373, 38.8474541, 3.32029319, -3.03482223 + 1374, 39.0207405, 3.68501592, -2.70964026 + 1375, 39.1591682, 4.04334497, -2.36510444 + 1376, 39.1019249, 4.40904284, -2.33188581 + 1377, 39.0229912, 4.77511406, -2.29051781 + 1378, 39.1584625, 5.13963223, -1.95433879 + 1379, 39.2734108, 5.4944582, -1.60743642 + 1380, 39.3961945, 5.84967709, -1.26860929 + 1381, 39.5568581, 6.21578455, -0.957187057 + 1382, 39.7540054, 6.57990217, -0.665589452 + 1383, 40.0492744, 6.96437645, -0.485698581 + 1405, 35.262394, 0.467456043, -5.97783756 + 1406, 35.5527573, 0.280540854, -5.80499315 + 1457, 37.0712318, 0.324649334, -5.01549721 + 1458, 37.1468697, 0.558611333, -4.60628891 + 1459, 37.114418, 0.765027821, -4.14856815 + 1460, 37.2010803, 1.08491921, -3.74724293 + 1461, 37.6603966, 1.41256773, -3.1326344 + 1462, 38.0748444, 1.72208464, -2.96610689 + 1463, 37.4089737, 1.04710495, -3.44614577 + 1465, 38.545723, 1.96261883, -2.88215423 + 1466, 38.7117195, 2.2614882, -2.99730253 + 1467, 38.5758247, 2.60041976, -3.31599164 + 1468, 38.8651962, 2.95906973, -3.06504536 + 1469, 39.0618782, 3.32376242, -2.74606228 + 1470, 39.2148018, 3.68168616, -2.40384316 + 1471, 39.346981, 4.03551531, -2.05412602 + 1472, 39.292408, 4.40021515, -2.02129745 + 1473, 39.2322121, 4.76847458, -1.99044192 + 1474, 39.369976, 5.12825155, -1.65812361 + 1475, 39.4982414, 5.48270321, -1.32477629 + 1476, 39.6601181, 5.84450436, -1.02088237 + 1477, 39.8445282, 6.20570803, -0.729069471 + 1478, 40.0696754, 6.57108831, -0.468850106 + 1551, 37.2919006, 0.2451341, -4.76742363 + 1552, 37.3510323, 0.457334876, -4.32699156 + 1553, 37.3074188, 0.712967753, -3.8486836 + 1554, 37.6203232, 1.01254523, -3.13962841 + 1555, 37.9218178, 1.38121653, -2.85790539 + 1556, 38.3680038, 1.66903615, -2.7425611 + 1557, 37.4857407, 0.66005218, -3.53552461 + 1559, 38.8301086, 1.92110622, -2.65684175 + 1560, 38.9568977, 2.25331378, -2.73103714 + 1561, 38.8193588, 2.6003027, -3.05332208 + 1562, 39.0758057, 2.96184754, -2.77152491 + 1563, 39.2610703, 3.32198429, -2.44353008 + 1564, 39.3997803, 3.67467904, -2.09175706 + 1565, 39.5321846, 4.02509594, -1.74178171 + 1566, 39.4830093, 4.38924837, -1.71216524 + 1567, 39.4332733, 4.75768185, -1.68764246 + 1568, 39.5769424, 5.1144948, -1.36307263 + 1569, 39.7483673, 5.47611761, -1.06774068 + 1570, 39.9399567, 5.83794308, -0.791532278 + 1571, 40.1443787, 6.19492579, -0.519479513 + 1647, 37.5047646, 0.384817958, -3.99143839 + 1648, 37.6720352, 0.614769876, -3.21648645 + 1649, 37.8794327, 0.984944522, -2.86715889 + 1650, 38.2192116, 1.33880746, -2.63161421 + 1651, 38.6646347, 1.62029588, -2.53288889 + 1652, 37.6850319, 0.330207616, -3.67911935 + 1654, 39.0745621, 1.9146539, -2.3861146 + 1655, 39.1894722, 2.2490797, -2.45103216 + 1656, 39.0410995, 2.60247302, -2.76822734 + 1657, 39.2737732, 2.96080661, -2.46789193 + 1658, 39.4431458, 3.31578708, -2.12934041 + 1659, 39.5803642, 3.66476822, -1.77671051 + 1660, 39.7222786, 4.01363516, -1.43325055 + 1661, 39.6855545, 4.37908554, -1.41282105 + 1662, 39.6370697, 4.74586773, -1.38966477 + 1663, 39.8157082, 5.10709906, -1.0967443 + 1664, 40.0204506, 5.47163963, -0.836458027 + 1665, 40.2241745, 5.82797527, -0.573889971 + 1744, 37.8905411, 0.583898187, -2.912328 + 1745, 38.1514053, 0.960013509, -2.61052155 + 1746, 38.5199509, 1.30187786, -2.42012954 + 1747, 38.9415894, 1.5943414, -2.30015802 + 1750, 39.3148727, 1.9062103, -2.1113286 + 1751, 39.4092445, 2.24537969, -2.1591289 + 1752, 39.2544823, 2.60119367, -2.47464705 + 1753, 39.4648666, 2.95698547, -2.15824366 + 1754, 39.6189232, 3.30669427, -1.81096876 + 1755, 39.7668686, 3.65320134, -1.46536517 + 1756, 39.9521713, 4.00698328, -1.15496588 + 1757, 39.9105492, 4.37161541, -1.1329577 + 1758, 39.8650627, 4.73842669, -1.11387825 + 1759, 40.0791397, 5.10392237, -0.859588265 + 1760, 40.2980728, 5.46595955, -0.616647363 + 1842, 38.1513252, 0.569698453, -2.64153671 + 1843, 38.452549, 0.938202679, -2.3962028 + 1844, 38.8237343, 1.26773632, -2.22184086 + 1845, 39.1968422, 1.58201575, -2.03972363 + 1848, 39.5263634, 1.90506661, -1.81289625 + 1849, 39.598484, 2.24503565, -1.84586108 + 1850, 39.4484634, 2.59942174, -2.16636395 + 1851, 39.6384964, 2.94984818, -1.83804047 + 1852, 39.8046761, 3.29555964, -1.49859297 + 1853, 39.9956856, 3.64578843, -1.18518102 + 1854, 40.2111816, 4.00406599, -0.909250736 + 1855, 40.1727295, 4.37039471, -0.894539535 + 1856, 40.1325493, 4.73739243, -0.882226348 + 1857, 40.3570099, 5.1029439, -0.647261977 + 1941, 38.4439201, 0.567996681, -2.41639638 + 1942, 38.7537537, 0.920588076, -2.19042993 + 1943, 39.0936012, 1.25266182, -1.97781503 + 1944, 39.4449158, 1.57014728, -1.77122462 + 1947, 39.7666435, 1.89068842, -1.53466594 + 1948, 39.8104362, 2.23568869, -1.54706144 + 1949, 39.6326027, 2.59518933, -1.8508817 + 1950, 39.8325996, 2.94030285, -1.52971184 + 1951, 40.0287666, 3.28647208, -1.21261775 + 1952, 40.2563095, 3.64149046, -0.938407958 + 1953, 40.5071106, 4.00797653, -0.724397957 + 1954, 40.4549828, 4.37174845, -0.693157494 + 1955, 40.4049911, 4.73714399, -0.666203797 + 2043, 38.7391014, 0.568356454, -2.20667505 + 2044, 39.0391769, 0.90763098, -1.96500075 + 2045, 39.3736458, 1.2334075, -1.74055278 + 2046, 39.7132416, 1.54981828, -1.51775217 + 2049, 40.019062, 1.87138748, -1.26409769 + 2050, 40.0460205, 2.22041607, -1.26536763 + 2051, 39.830677, 2.58631206, -1.54393411 + 2052, 40.0495834, 2.9287591, -1.23793912 + 2053, 40.2875175, 3.27896333, -0.961250722 + 2054, 40.5424309, 3.64160109, -0.733088076 + 2145, 39.0277443, 0.562029481, -1.99374688 + 2146, 39.3270035, 0.890259266, -1.73961377 + 2147, 39.6626396, 1.2094965, -1.51105785 + 2148, 39.9904709, 1.52511632, -1.26878631 + 2151, 40.3108025, 1.8470192, -1.03659248 + 2152, 40.3171082, 2.20144224, -1.02117872 + 2153, 40.0564919, 2.57313728, -1.25708199 + 2154, 40.3117714, 2.91811395, -0.988921463 + 2155, 40.5718002, 3.27402759, -0.746706367 + 2248, 39.3146095, 0.546521485, -1.77641821 + 2249, 39.6237259, 0.867994487, -1.52436471 + 2250, 39.9531631, 1.18217003, -1.27692759 + 2251, 40.287014, 1.49581075, -1.03883135 + 2254, 40.6198349, 1.82273686, -0.841227829 + 2255, 40.6205292, 2.18397474, -0.830564678 + 2256, 40.3206215, 2.55871749, -1.00868785 + 2257, 40.5993996, 2.9091177, -0.781215608 + 2354, 39.6090279, 0.526046753, -1.56946611 + 2355, 39.9107552, 0.838355184, -1.29499745 + 2356, 40.2641106, 1.14648616, -1.06629324 + 2357, 40.6162758, 1.46274948, -0.858284652 + 2360, 40.9543037, 1.80483508, -0.699903309 + 2362, 40.6129074, 2.54585147, -0.807650089 + 2454, 39.9090996, 0.499893278, -1.37884927 + 2455, 40.2230873, 0.796499074, -1.09488297 + 2456, 40.5981674, 1.10136247, -0.884527266 + 2457, 40.9625359, 1.42922926, -0.702434063 + 2554, 40.2192726, 0.460563034, -1.21433496 + 2555, 40.5606117, 0.739465535, -0.932800949 + 2556, 40.9582825, 1.05046546, -0.753893912 + 2651, 40.9295082, 0.679916978, -0.860399663 +*Element, type=S4 + 160, 141, 165, 198, 167 + 161, 199, 200, 166, 167 + 178, 166, 200, 218, 185 + 192, 198, 232, 199, 167 + 193, 233, 200, 199, 234 + 210, 218, 200, 233, 252 + 211, 218, 253, 219, 185 + 225, 199, 232, 267, 234 + 226, 268, 269, 233, 234 + 244, 233, 269, 289, 252 + 245, 290, 253, 218, 252 + 246, 219, 253, 291, 254 + 260, 267, 305, 268, 234 + 261, 306, 269, 268, 307 + 280, 289, 269, 306, 328 + 281, 289, 329, 290, 252 + 282, 291, 253, 290, 330 + 283, 291, 331, 292, 254 + 297, 268, 305, 345, 307 + 298, 346, 347, 306, 307 + 318, 306, 347, 370, 328 + 319, 371, 329, 289, 328 + 320, 290, 329, 372, 330 + 321, 373, 331, 291, 330 + 322, 292, 331, 374, 332 + 336, 345, 388, 346, 307 + 337, 389, 347, 346, 390 + 359, 370, 347, 389, 415 + 360, 370, 416, 371, 328 + 361, 372, 329, 371, 417 + 362, 372, 418, 373, 330 + 363, 374, 331, 373, 419 + 364, 374, 420, 375, 332 + 378, 346, 388, 434, 390 + 379, 435, 436, 389, 390 + 404, 389, 436, 464, 415 + 405, 465, 416, 370, 415 + 406, 371, 416, 466, 417 + 407, 467, 418, 372, 417 + 408, 373, 418, 468, 419 + 409, 469, 420, 374, 419 + 410, 375, 420, 470, 421 + 424, 434, 484, 435, 390 + 425, 485, 436, 435, 486 + 452, 464, 436, 485, 515 + 453, 464, 516, 465, 415 + 454, 466, 416, 465, 517 + 455, 466, 518, 467, 417 + 456, 468, 418, 467, 519 + 457, 468, 520, 469, 419 + 458, 470, 420, 469, 521 + 459, 470, 522, 471, 421 + 473, 435, 484, 536, 486 + 474, 537, 538, 485, 486 + 502, 485, 538, 569, 515 + 503, 570, 516, 464, 515 + 504, 465, 516, 571, 517 + 505, 572, 518, 466, 517 + 506, 467, 518, 573, 519 + 507, 574, 520, 468, 519 + 508, 469, 520, 575, 521 + 509, 576, 522, 470, 521 + 510, 471, 522, 577, 523 + 524, 536, 591, 537, 486 + 525, 569, 538, 537, 592 + 556, 569, 629, 630, 515 + 557, 571, 516, 570, 631 + 558, 630, 632, 570, 515 + 559, 571, 633, 572, 517 + 560, 573, 518, 572, 634 + 561, 573, 635, 574, 519 + 562, 575, 520, 574, 636 + 563, 575, 637, 576, 521 + 564, 577, 522, 576, 638 + 565, 577, 639, 578, 523 + 579, 537, 591, 653, 592 + 580, 654, 629, 569, 592 + 615, 630, 629, 692, 693 + 616, 570, 632, 694, 631 + 617, 695, 633, 571, 631 + 618, 696, 632, 630, 693 + 619, 572, 633, 697, 634 + 620, 698, 635, 573, 634 + 621, 574, 635, 699, 636 + 622, 700, 637, 575, 636 + 623, 576, 637, 701, 638 + 624, 702, 639, 577, 638 + 625, 578, 639, 703, 640 + 639, 653, 717, 654, 592 + 640, 692, 629, 654, 718 + 679, 692, 761, 762, 693 + 680, 694, 632, 696, 763 + 681, 694, 764, 695, 631 + 682, 697, 633, 695, 765 + 683, 762, 766, 696, 693 + 684, 697, 767, 768, 634 + 685, 699, 635, 698, 769 + 686, 768, 770, 698, 634 + 687, 699, 771, 700, 636 + 688, 701, 637, 700, 772 + 689, 701, 773, 702, 638 + 690, 703, 639, 702, 774 + 691, 703, 775, 704, 640 + 705, 654, 717, 789, 718 + 706, 790, 761, 692, 718 + 748, 762, 761, 836, 837 + 749, 696, 766, 838, 763 + 750, 839, 764, 694, 763 + 751, 695, 764, 840, 765 + 752, 841, 767, 697, 765 + 753, 842, 766, 762, 837 + 754, 768, 767, 843, 844 + 755, 698, 770, 845, 769 + 756, 846, 771, 699, 769 + 757, 845, 770, 768, 844 + 758, 700, 771, 847, 772 + 759, 848, 773, 701, 772 + 760, 702, 773, 849, 774 + 761, 850, 775, 703, 774 + 762, 704, 775, 851, 776 + 776, 789, 865, 790, 718 + 777, 836, 761, 790, 866 + 821, 836, 913, 914, 837 + 822, 838, 766, 842, 915 + 823, 838, 916, 839, 763 + 824, 840, 764, 839, 917 + 825, 840, 918, 841, 765 + 826, 919, 767, 841, 920 + 827, 914, 921, 842, 837 + 828, 843, 767, 919, 922 + 829, 843, 923, 924, 844 + 830, 845, 925, 846, 769 + 831, 847, 771, 846, 926 + 832, 924, 925, 845, 844 + 833, 847, 927, 848, 772 + 834, 849, 773, 848, 928 + 835, 849, 929, 850, 774 + 836, 851, 775, 850, 930 + 852, 790, 865, 947, 866 + 853, 948, 913, 836, 866 + 901, 914, 913, 1000, 1001 + 902, 842, 921, 1002, 915 + 903, 1003, 916, 838, 915 + 904, 839, 916, 1004, 917 + 905, 1005, 918, 840, 917 + 906, 841, 918, 1006, 920 + 907, 1006, 1007, 919, 920 + 908, 1008, 921, 914, 1001 + 909, 919, 1007, 1009, 922 + 910, 1010, 923, 843, 922 + 911, 924, 923, 1011, 1012 + 912, 846, 925, 1013, 926 + 913, 1014, 927, 847, 926 + 914, 1013, 925, 924, 1012 + 915, 848, 927, 1015, 928 + 916, 1016, 929, 849, 928 + 917, 850, 929, 1017, 930 + 918, 1018, 931, 851, 930 + 935, 947, 1035, 948, 866 + 936, 1000, 913, 948, 1036 + 987, 1000, 1089, 1090, 1001 + 988, 1002, 921, 1008, 1091 + 989, 1002, 1092, 1003, 915 + 990, 1004, 916, 1003, 1093 + 991, 1004, 1094, 1005, 917 + 992, 1006, 918, 1005, 1095 + 993, 1096, 1007, 1006, 1095 + 994, 1090, 1097, 1008, 1001 + 995, 1009, 1007, 1096, 1098 + 996, 1009, 1099, 1010, 922 + 997, 1011, 923, 1010, 1100 + 998, 1011, 1101, 1102, 1012 + 999, 1013, 1103, 1014, 926 +1000, 1015, 927, 1014, 1104 +1001, 1102, 1103, 1013, 1012 +1002, 1015, 1105, 1016, 928 +1003, 1017, 929, 1016, 1106 +1004, 1017, 1107, 1018, 930 +1005, 1019, 931, 1018, 1108 +1006, 1019, 1109, 1020, 933 +1022, 948, 1035, 1126, 1036 +1023, 1127, 1089, 1000, 1036 +1073, 1090, 1089, 1178, 1179 +1074, 1008, 1097, 1180, 1091 +1075, 1181, 1092, 1002, 1091 +1076, 1003, 1092, 1182, 1093 +1077, 1183, 1094, 1004, 1093 +1078, 1005, 1094, 1184, 1095 +1079, 1184, 1185, 1096, 1095 +1080, 1186, 1097, 1090, 1179 +1081, 1096, 1185, 1187, 1098 +1082, 1188, 1099, 1009, 1098 +1083, 1010, 1099, 1189, 1100 +1084, 1190, 1101, 1011, 1100 +1085, 1102, 1101, 1191, 1192 +1086, 1014, 1103, 1193, 1104 +1087, 1194, 1105, 1015, 1104 +1088, 1193, 1103, 1102, 1192 +1089, 1016, 1105, 1195, 1106 +1090, 1196, 1107, 1017, 1106 +1091, 1018, 1107, 1197, 1108 +1092, 1198, 1109, 1019, 1108 +1093, 1020, 1109, 1199, 1111 +1110, 1126, 1217, 1127, 1036 +1111, 1178, 1089, 1127, 1218 +1160, 1178, 1268, 1269, 1179 +1161, 1180, 1097, 1186, 1270 +1162, 1180, 1271, 1181, 1091 +1163, 1182, 1092, 1181, 1272 +1164, 1182, 1273, 1183, 1093 +1165, 1184, 1094, 1183, 1274 +1166, 1275, 1185, 1184, 1274 +1167, 1269, 1276, 1186, 1179 +1168, 1187, 1185, 1277, 1278 +1169, 1187, 1279, 1188, 1098 +1170, 1189, 1099, 1188, 1280 +1171, 1189, 1281, 1190, 1100 +1172, 1191, 1101, 1190, 1282 +1173, 1191, 1283, 1284, 1192 +1174, 1193, 1285, 1194, 1104 +1175, 1195, 1105, 1194, 1286 +1176, 1284, 1285, 1193, 1192 +1177, 1195, 1287, 1196, 1106 +1178, 1197, 1107, 1196, 1288 +1179, 1197, 1289, 1198, 1108 +1180, 1199, 1109, 1198, 1290 +1181, 1199, 1291, 1201, 1111 +1200, 1127, 1217, 1312, 1218 +1201, 1312, 1268, 1178, 1218 +1249, 1186, 1276, 1361, 1270 +1250, 1362, 1271, 1180, 1270 +1251, 1181, 1271, 1363, 1272 +1252, 1364, 1273, 1182, 1272 +1253, 1183, 1273, 1365, 1274 +1254, 1277, 1185, 1275, 1366 +1255, 1365, 1367, 1275, 1274 +1257, 1277, 1369, 1370, 1278 +1258, 1371, 1279, 1187, 1278 +1259, 1188, 1279, 1371, 1280 +1260, 1372, 1281, 1189, 1280 +1261, 1190, 1281, 1373, 1282 +1262, 1374, 1283, 1191, 1282 +1263, 1284, 1283, 1375, 1376 +1264, 1194, 1285, 1377, 1286 +1265, 1378, 1287, 1195, 1286 +1266, 1377, 1285, 1284, 1376 +1267, 1196, 1287, 1379, 1288 +1268, 1380, 1289, 1197, 1288 +1269, 1198, 1289, 1381, 1290 +1270, 1382, 1291, 1199, 1290 +1292, 1312, 1217, 1405, 1406 +1293, 1359, 1268, 1312, 1406 +1341, 1361, 1457, 1362, 1270 +1342, 1363, 1271, 1362, 1458 +1343, 1363, 1459, 1364, 1272 +1344, 1365, 1273, 1364, 1460 +1345, 1275, 1367, 1461, 1366 +1346, 1462, 1369, 1277, 1366 +1347, 1463, 1367, 1365, 1460 +1349, 1370, 1369, 1465, 1466 +1350, 1370, 1467, 1371, 1278 +1351, 1371, 1467, 1372, 1280 +1352, 1373, 1281, 1372, 1468 +1353, 1373, 1469, 1374, 1282 +1354, 1375, 1283, 1374, 1470 +1355, 1375, 1471, 1472, 1376 +1356, 1377, 1473, 1378, 1286 +1357, 1379, 1287, 1378, 1474 +1358, 1472, 1473, 1377, 1376 +1359, 1379, 1475, 1380, 1288 +1360, 1381, 1289, 1380, 1476 +1361, 1381, 1477, 1382, 1290 +1362, 1383, 1291, 1382, 1478 +1434, 1362, 1457, 1551, 1458 +1435, 1552, 1459, 1363, 1458 +1436, 1364, 1459, 1553, 1460 +1437, 1461, 1367, 1463, 1554 +1438, 1461, 1555, 1462, 1366 +1439, 1465, 1369, 1462, 1556 +1440, 1553, 1557, 1463, 1460 +1442, 1465, 1559, 1560, 1466 +1443, 1561, 1467, 1370, 1466 +1444, 1372, 1467, 1561, 1468 +1445, 1562, 1469, 1373, 1468 +1446, 1374, 1469, 1563, 1470 +1447, 1564, 1471, 1375, 1470 +1448, 1472, 1471, 1565, 1566 +1449, 1378, 1473, 1567, 1474 +1450, 1568, 1475, 1379, 1474 +1451, 1567, 1473, 1472, 1566 +1452, 1380, 1475, 1569, 1476 +1453, 1570, 1477, 1381, 1476 +1454, 1382, 1477, 1571, 1478 +1527, 1553, 1459, 1552, 1647 +1528, 1463, 1557, 1648, 1554 +1529, 1649, 1555, 1461, 1554 +1530, 1462, 1555, 1650, 1556 +1531, 1651, 1559, 1465, 1556 +1532, 1652, 1557, 1553, 1647 +1534, 1560, 1559, 1654, 1655 +1535, 1560, 1656, 1561, 1466 +1536, 1561, 1656, 1562, 1468 +1537, 1563, 1469, 1562, 1657 +1538, 1563, 1658, 1564, 1470 +1539, 1565, 1471, 1564, 1659 +1540, 1565, 1660, 1661, 1566 +1541, 1567, 1662, 1568, 1474 +1542, 1569, 1475, 1568, 1663 +1543, 1661, 1662, 1567, 1566 +1544, 1569, 1664, 1570, 1476 +1545, 1571, 1477, 1570, 1665 +1621, 1648, 1744, 1649, 1554 +1622, 1650, 1555, 1649, 1745 +1623, 1650, 1746, 1651, 1556 +1624, 1654, 1559, 1651, 1747 +1627, 1654, 1750, 1751, 1655 +1628, 1752, 1656, 1560, 1655 +1629, 1562, 1656, 1752, 1657 +1630, 1753, 1658, 1563, 1657 +1631, 1564, 1658, 1754, 1659 +1632, 1755, 1660, 1565, 1659 +1633, 1661, 1660, 1756, 1757 +1634, 1568, 1662, 1758, 1663 +1635, 1759, 1664, 1569, 1663 +1636, 1758, 1662, 1661, 1757 +1637, 1570, 1664, 1760, 1665 +1716, 1649, 1744, 1842, 1745 +1717, 1843, 1746, 1650, 1745 +1718, 1651, 1746, 1844, 1747 +1719, 1845, 1750, 1654, 1747 +1722, 1751, 1750, 1848, 1849 +1723, 1751, 1850, 1752, 1655 +1724, 1752, 1850, 1753, 1657 +1725, 1754, 1658, 1753, 1851 +1726, 1754, 1852, 1755, 1659 +1727, 1756, 1660, 1755, 1853 +1728, 1756, 1854, 1855, 1757 +1729, 1758, 1856, 1759, 1663 +1730, 1760, 1664, 1759, 1857 +1731, 1855, 1856, 1758, 1757 +1813, 1842, 1941, 1843, 1745 +1814, 1844, 1746, 1843, 1942 +1815, 1844, 1943, 1845, 1747 +1816, 1848, 1750, 1845, 1944 +1819, 1848, 1947, 1948, 1849 +1820, 1949, 1850, 1751, 1849 +1821, 1753, 1850, 1949, 1851 +1822, 1950, 1852, 1754, 1851 +1823, 1755, 1852, 1951, 1853 +1824, 1952, 1854, 1756, 1853 +1825, 1855, 1854, 1953, 1954 +1826, 1759, 1856, 1955, 1857 +1828, 1955, 1856, 1855, 1954 +1912, 1843, 1941, 2043, 1942 +1913, 2044, 1943, 1844, 1942 +1914, 1845, 1943, 2045, 1944 +1915, 2046, 1947, 1848, 1944 +1918, 1948, 1947, 2049, 2050 +1919, 1948, 2051, 1949, 1849 +1920, 1949, 2051, 1950, 1851 +1921, 1951, 1852, 1950, 2052 +1922, 1951, 2053, 1952, 1853 +1923, 1953, 1854, 1952, 2054 +2011, 2043, 2145, 2044, 1942 +2012, 2045, 1943, 2044, 2146 +2013, 2045, 2147, 2046, 1944 +2014, 2049, 1947, 2046, 2148 +2017, 2049, 2151, 2152, 2050 +2018, 2153, 2051, 1948, 2050 +2019, 1950, 2051, 2153, 2052 +2020, 2154, 2053, 1951, 2052 +2021, 1952, 2053, 2155, 2054 +2111, 2044, 2145, 2248, 2146 +2112, 2249, 2147, 2045, 2146 +2113, 2046, 2147, 2250, 2148 +2114, 2251, 2151, 2049, 2148 +2117, 2152, 2151, 2254, 2255 +2118, 2152, 2256, 2153, 2050 +2119, 2153, 2256, 2154, 2052 +2120, 2155, 2053, 2154, 2257 +2213, 2248, 2354, 2249, 2146 +2214, 2250, 2147, 2249, 2355 +2215, 2250, 2356, 2251, 2148 +2216, 2254, 2151, 2251, 2357 +2220, 2362, 2256, 2152, 2255 +2221, 2154, 2256, 2362, 2257 +2313, 2249, 2354, 2454, 2355 +2314, 2455, 2356, 2250, 2355 +2315, 2251, 2356, 2456, 2357 +2316, 2457, 2360, 2254, 2357 +2410, 2454, 2554, 2455, 2355 +2411, 2456, 2356, 2455, 2555 +2412, 2456, 2556, 2457, 2357 +2509, 2651, 2556, 2456, 2555 +*End Part + diff --git a/tests/grid/grids/abaqus/3d/codim_1-tube.inp b/tests/grid/grids/abaqus/3d/codim_1-tube.inp new file mode 100644 index 0000000000..e41fa8f09a --- /dev/null +++ b/tests/grid/grids/abaqus/3d/codim_1-tube.inp @@ -0,0 +1,221 @@ +*Part, name=Part-1 +*Node + 1, -8.66025448, 5., 10. + 2, -8.66025448, 5., 0. + 3, 0., -10., 0. + 4, 10., 8.13201453e-08, 0. + 5, 0., -10., 10. + 6, 10., -8.13201453e-08, 10. + 7, -8.66025352, 5.00000048, 7.5 + 8, -8.66025352, 5.00000048, 5. + 9, -8.66025352, 5.00000048, 2.5 + 10, -9.65925789, 2.58819032, 0. + 11, -10., -8.13201453e-08, 0. + 12, -9.65925789, -2.58819032, 0. + 13, -8.66025448, -5., 0. + 14, -7.07106781, -7.07106781, 0. + 15, -5., -8.66025448, 0. + 16, -2.58819032, -9.65925789, 0. + 17, 2.58819032, -9.65925789, 0. + 18, 5., -8.66025448, 0. + 19, 7.07106781, -7.07106781, 0. + 20, 8.66025448, -5., 0. + 21, 9.65925789, -2.58819032, 0. + 22, 10., 0., 2.5 + 23, 10., 0., 5. + 24, 10., 0., 7.5 + 25, -2.58819032, -9.65925789, 10. + 26, -5., -8.66025448, 10. + 27, -7.07106781, -7.07106781, 10. + 28, -8.66025448, -5., 10. + 29, -9.65925789, -2.58819032, 10. + 30, -10., 8.13201453e-08, 10. + 31, -9.65925789, 2.58819032, 10. + 32, 9.65925789, -2.58819032, 10. + 33, 8.66025448, -5., 10. + 34, 7.07106781, -7.07106781, 10. + 35, 5., -8.66025448, 10. + 36, 2.58819032, -9.65925789, 10. + 37, -7.07106781, 7.07106781, 10. + 38, -5., 8.66025448, 10. + 39, -2.58819032, 9.65925789, 10. + 40, 8.13201453e-08, 10., 10. + 41, 2.58819032, 9.65925789, 10. + 42, 5., 8.66025448, 10. + 43, 7.07106781, 7.07106781, 10. + 44, 8.66025448, 5., 10. + 45, 9.65925789, 2.58819032, 10. + 46, 9.65925789, 2.58819032, 0. + 47, 8.66025448, 5., 0. + 48, 7.07106781, 7.07106781, 0. + 49, 5., 8.66025448, 0. + 50, 2.58819032, 9.65925789, 0. + 51, -8.13201453e-08, 10., 0. + 52, -2.58819032, 9.65925789, 0. + 53, -5., 8.66025448, 0. + 54, -7.07106781, 7.07106781, 0. + 55, -9.65925884, 2.58818913, 7.5 + 56, -9.65925884, 2.58818913, 5. + 57, -9.65925884, 2.58818913, 2.5 + 58, -10., -1.509958e-06, 7.5 + 59, -10., -1.509958e-06, 5. + 60, -10., -1.509958e-06, 2.5 + 61, -9.65925789, -2.58819151, 7.5 + 62, -9.65925789, -2.58819151, 5. + 63, -9.65925789, -2.58819151, 2.5 + 64, -8.66025352, -5.00000048, 7.5 + 65, -8.66025352, -5.00000048, 5. + 66, -8.66025352, -5.00000048, 2.5 + 67, -7.07106781, -7.07106781, 7.5 + 68, -7.07106781, -7.07106781, 5. + 69, -7.07106781, -7.07106781, 2.5 + 70, -5.00000048, -8.66025352, 7.5 + 71, -5.00000048, -8.66025352, 5. + 72, -5.00000048, -8.66025352, 2.5 + 73, -2.58819151, -9.65925789, 7.5 + 74, -2.58819151, -9.65925789, 5. + 75, -2.58819151, -9.65925789, 2.5 + 76, -1.62920685e-06, -10., 7.5 + 77, -1.62920685e-06, -10., 5. + 78, -1.62920685e-06, -10., 2.5 + 79, 2.58819079, -9.65925789, 7.5 + 80, 2.58819079, -9.65925789, 5. + 81, 2.58819079, -9.65925789, 2.5 + 82, 4.99999952, -8.66025448, 7.5 + 83, 4.99999952, -8.66025448, 5. + 84, 4.99999952, -8.66025448, 2.5 + 85, 7.07106781, -7.07106781, 7.5 + 86, 7.07106781, -7.07106781, 5. + 87, 7.07106781, -7.07106781, 2.5 + 88, 8.66025448, -4.99999952, 7.5 + 89, 8.66025448, -4.99999952, 5. + 90, 8.66025352, -5., 2.5 + 91, 9.65925789, -2.58819056, 7.5 + 92, 9.65925789, -2.58819056, 5. + 93, 9.65925789, -2.58819056, 2.5 + 94, -7.07106781, 7.07106781, 2.5 + 95, -7.07106781, 7.07106781, 5. + 96, -7.07106781, 7.07106781, 7.5 + 97, -5.00000048, 8.66025352, 2.5 + 98, -5.00000048, 8.66025352, 5. + 99, -5.00000048, 8.66025352, 7.5 + 100, -2.58819151, 9.65925789, 2.5 + 101, -2.58819151, 9.65925789, 5. + 102, -2.58819151, 9.65925789, 7.5 + 103, -1.62920685e-06, 10., 2.5 + 104, -1.62920685e-06, 10., 5. + 105, -1.62920685e-06, 10., 7.5 + 106, 2.58819079, 9.65925789, 2.5 + 107, 2.58819079, 9.65925789, 5. + 108, 2.58819079, 9.65925789, 7.5 + 109, 4.99999952, 8.66025448, 2.5 + 110, 4.99999952, 8.66025448, 5. + 111, 4.99999952, 8.66025448, 7.5 + 112, 7.07106781, 7.07106781, 2.5 + 113, 7.07106781, 7.07106781, 5. + 114, 7.07106781, 7.07106781, 7.5 + 115, 8.66025448, 4.99999952, 2.5 + 116, 8.66025448, 4.99999952, 5. + 117, 8.66025352, 5., 7.5 + 118, 9.65925789, 2.58819056, 2.5 + 119, 9.65925789, 2.58819056, 5. + 120, 9.65925789, 2.58819056, 7.5 +*Element, type=S4R + 1, 1, 7, 55, 31 + 2, 7, 8, 56, 55 + 3, 8, 9, 57, 56 + 4, 9, 2, 10, 57 + 5, 31, 55, 58, 30 + 6, 55, 56, 59, 58 + 7, 56, 57, 60, 59 + 8, 57, 10, 11, 60 + 9, 30, 58, 61, 29 +10, 58, 59, 62, 61 +11, 59, 60, 63, 62 +12, 60, 11, 12, 63 +13, 29, 61, 64, 28 +14, 61, 62, 65, 64 +15, 62, 63, 66, 65 +16, 63, 12, 13, 66 +17, 28, 64, 67, 27 +18, 64, 65, 68, 67 +19, 65, 66, 69, 68 +20, 66, 13, 14, 69 +21, 27, 67, 70, 26 +22, 67, 68, 71, 70 +23, 68, 69, 72, 71 +24, 69, 14, 15, 72 +25, 26, 70, 73, 25 +26, 70, 71, 74, 73 +27, 71, 72, 75, 74 +28, 72, 15, 16, 75 +29, 25, 73, 76, 5 +30, 73, 74, 77, 76 +31, 74, 75, 78, 77 +32, 75, 16, 3, 78 +33, 5, 76, 79, 36 +34, 76, 77, 80, 79 +35, 77, 78, 81, 80 +36, 78, 3, 17, 81 +37, 36, 79, 82, 35 +38, 79, 80, 83, 82 +39, 80, 81, 84, 83 +40, 81, 17, 18, 84 +41, 35, 82, 85, 34 +42, 82, 83, 86, 85 +43, 83, 84, 87, 86 +44, 84, 18, 19, 87 +45, 34, 85, 88, 33 +46, 85, 86, 89, 88 +47, 86, 87, 90, 89 +48, 87, 19, 20, 90 +49, 33, 88, 91, 32 +50, 88, 89, 92, 91 +51, 89, 90, 93, 92 +52, 90, 20, 21, 93 +53, 32, 91, 24, 6 +54, 91, 92, 23, 24 +55, 92, 93, 22, 23 +56, 93, 21, 4, 22 +57, 2, 9, 94, 54 +58, 9, 8, 95, 94 +59, 8, 7, 96, 95 +60, 7, 1, 37, 96 +61, 54, 94, 97, 53 +62, 94, 95, 98, 97 +63, 95, 96, 99, 98 +64, 96, 37, 38, 99 +65, 53, 97, 100, 52 +66, 97, 98, 101, 100 +67, 98, 99, 102, 101 +68, 99, 38, 39, 102 +69, 52, 100, 103, 51 +70, 100, 101, 104, 103 +71, 101, 102, 105, 104 +72, 102, 39, 40, 105 +73, 51, 103, 106, 50 +74, 103, 104, 107, 106 +75, 104, 105, 108, 107 +76, 105, 40, 41, 108 +77, 50, 106, 109, 49 +78, 106, 107, 110, 109 +79, 107, 108, 111, 110 +80, 108, 41, 42, 111 +81, 49, 109, 112, 48 +82, 109, 110, 113, 112 +83, 110, 111, 114, 113 +84, 111, 42, 43, 114 +85, 48, 112, 115, 47 +86, 112, 113, 116, 115 +87, 113, 114, 117, 116 +88, 114, 43, 44, 117 +89, 47, 115, 118, 46 +90, 115, 116, 119, 118 +91, 116, 117, 120, 119 +92, 117, 44, 45, 120 +93, 46, 118, 22, 4 +94, 118, 119, 23, 22 +95, 119, 120, 24, 23 +96, 120, 45, 6, 24 +*End Part + -- 2.39.5