From 6fd1abebde73dca7f6ae7ed3d8474492fe8707fa Mon Sep 17 00:00:00 2001 From: Yoann Audouin Date: Wed, 26 May 2021 16:30:10 +0200 Subject: [PATCH] [test] Adding test on geometry for VoroGaus plugin --- src/Plugins/VoroGauss/Test/test_VoroGauss.py | 103 +++++++++++++++---- 1 file changed, 84 insertions(+), 19 deletions(-) diff --git a/src/Plugins/VoroGauss/Test/test_VoroGauss.py b/src/Plugins/VoroGauss/Test/test_VoroGauss.py index caee1ee4..85a75e5a 100644 --- a/src/Plugins/VoroGauss/Test/test_VoroGauss.py +++ b/src/Plugins/VoroGauss/Test/test_VoroGauss.py @@ -18,12 +18,14 @@ # # Author : Yoann Audouin (EDF R&D) +from os import path, remove #### import the simple module from the paraview from paraview.simple import * from vtk.util import numpy_support import numpy as np from medcoupling import * +import MEDLoader as ml from MEDLoader import * #### disable automatic camera reset on 'Show' @@ -34,9 +36,9 @@ def MyAssert(clue): raise RuntimeError("Assertion failed !") -def test_VoroGauss(result, ref, field_name): +def test_fields(result, ref, field_name): """ - Test result of VoroGauss filter + Test fields content of result from VoroGauss filter """ ds0 = servermanager.Fetch(result) block = ds0.GetBlock(0) @@ -45,35 +47,68 @@ def test_VoroGauss(result, ref, field_name): MyAssert(np.allclose(data, ref)) +def test_geom(result, ref_mass_center, ref_volumes): + """ + Test center of mass and volumes of result from VoroGauss filter + """ + tmp_file = 'tmp_mesh.med' + if path.exists(tmp_file): + remove(tmp_file) + + # If we write the time steps It will move the time to the last time step + SaveData(tmp_file, proxy=result, WriteAllTimeSteps=0) + + data = ml.MEDFileData(tmp_file) + mesh = data.getMeshes()[0] + + mesh0 = mesh.getLevel0Mesh() + + centerOfMass = mesh0.computeCellCenterOfMass() + MyAssert(np.allclose(np.array(centerOfMass.getValuesAsTuple()), ref_mass_center)) + + volume = mesh0.getMeasureField(True).getArray() + MyAssert(np.allclose(np.array(volume.getValues()), ref_volumes)) + + ### # Test of testMEDReader14 ### # create a new 'MED Reader' +file_name = 'testMEDReader14.med' testMEDReader14med = MEDReader(registrationName='testMEDReader14.med', - FileName='testMEDReader14.med') + FileName=file_name) testMEDReader14med.AllTimeSteps = ['0000', '0001', '0002', '0003', '0004'] +print("Testing {}".format(file_name)) fields = [(['zeField0'], ['TS0/Mesh/ComSup0/zeField0_MM0@@][@@GAUSS'], - [[[100.,101.], [102.,103.], [104.,105.], [106.,107.], [108.,109.], [110.,111.], [112.,113.], [114.,115.], [ 0., 1.], [ 2., 3.], [ 4., 5.], [ 6., 7.], [ 8., 9.], [ 10., 11.], [ 12., 13.], [ 14., 15.], [ 16., 17.], [ 18., 19.], [ 20., 21.], [ 22., 23.], [ 24., 25.], [ 26., 27.], [ 28., 29.], [ 30., 31.], [ 32., 33.], [ 34., 35.], [ 36., 37.], [ 38., 39.], [ 40., 41.], [ 42., 43.], [ 44., 45.], [ 46., 47.], [ 48., 49.], [ 50., 51.], [ 52., 53.], [ 54., 55.], [ 56., 57.], [ 58., 59.], [ 60., 61.], [ 62., 63.], [ 64., 65.], [ 66., 67.], [ 68., 69.], [ 70., 71.], [ 72., 73.], [ 74., 75.], [ 76., 77.], [ 78., 79.], [ 80., 81.], [ 82., 83.], [ 84., 85.], [ 86., 87.], [ 88., 89.], [ 90., 91.], [ 92., 93.], [ 94., 95.], [ 96., 97.], [ 98., 99.], [100.,101.], [102.,103.], [104.,105.], [106.,107.], [108.,109.], [110.,111.], [112.,113.], [114.,115.], [116.,117.], [118.,119.], [120.,121.], [122.,123.], [124.,125.], [126.,127.], [128.,129.], [130.,131.], [132.,133.], [134.,135.], [136.,137.], [138.,139.], [140.,141.], [142.,143.], [144.,145.], [146.,147.], [148.,149.], [150.,151.], [152.,153.], [154.,155.], [156.,157.], [158.,159.], [160.,161.], [162.,163.], [164.,165.], [166.,167.]]] + [[[100.,101.], [102.,103.], [104.,105.], [106.,107.], [108.,109.], [110.,111.], [112.,113.], [114.,115.], [ 0., 1.], [ 2., 3.], [ 4., 5.], [ 6., 7.], [ 8., 9.], [ 10., 11.], [ 12., 13.], [ 14., 15.], [ 16., 17.], [ 18., 19.], [ 20., 21.], [ 22., 23.], [ 24., 25.], [ 26., 27.], [ 28., 29.], [ 30., 31.], [ 32., 33.], [ 34., 35.], [ 36., 37.], [ 38., 39.], [ 40., 41.], [ 42., 43.], [ 44., 45.], [ 46., 47.], [ 48., 49.], [ 50., 51.], [ 52., 53.], [ 54., 55.], [ 56., 57.], [ 58., 59.], [ 60., 61.], [ 62., 63.], [ 64., 65.], [ 66., 67.], [ 68., 69.], [ 70., 71.], [ 72., 73.], [ 74., 75.], [ 76., 77.], [ 78., 79.], [ 80., 81.], [ 82., 83.], [ 84., 85.], [ 86., 87.], [ 88., 89.], [ 90., 91.], [ 92., 93.], [ 94., 95.], [ 96., 97.], [ 98., 99.], [100.,101.], [102.,103.], [104.,105.], [106.,107.], [108.,109.], [110.,111.], [112.,113.], [114.,115.], [116.,117.], [118.,119.], [120.,121.], [122.,123.], [124.,125.], [126.,127.], [128.,129.], [130.,131.], [132.,133.], [134.,135.], [136.,137.], [138.,139.], [140.,141.], [142.,143.], [144.,145.], [146.,147.], [148.,149.], [150.,151.], [152.,153.], [154.,155.], [156.,157.], [158.,159.], [160.,161.], [162.,163.], [164.,165.], [166.,167.]]], + [ [5.5499999999999989, 0.22500000000000014, 0], [5.1858543417366949, 0.1780112044817927, 0], [5.1858543417366949, 0.63613445378151257, 0], [5.7750000000000004, 0.45000000000000001, 0], [5.3638655462184888, 0.81414565826330532, 0], [5.8219887955182079, 0.81414565826330532, 0], [4.5500000000000007, 1.2250000000000001, 0], [4.1858543417366949, 1.178011204481793, 0], [4.1858543417366958, 1.6361344537815128, 0], [4.7750000000000004, 1.45, 0], [4.3638655462184879, 1.8141456582633053, 0], [4.8219887955182088, 1.8141456582633049, 0], [5.5499999999999989, 1.2250000000000001, 0], [5.1858543417366949, 1.1780112044817923, 0], [5.1858543417366949, 1.6361344537815128, 0], [5.7750000000000004, 1.4499999999999997, 0], [5.3638655462184897, 1.8141456582633051, 0], [5.821988795518207, 1.8141456582633055, 0], [4.5500000000000007, 2.2250000000000005, 0], [4.1858543417366949, 2.1780112044817934, 0], [4.1858543417366967, 2.6361344537815126, 0], [4.7750000000000004, 2.4500000000000002, 0], [4.3638655462184879, 2.8141456582633051, 0], [4.8219887955182088, 2.8141456582633051, 0], [5.5499999999999998, 2.2250000000000001, 0], [5.1858543417366949, 2.1780112044817934, 0], [5.1858543417366949, 2.6361344537815126, 0], [5.7750000000000004, 2.4500000000000002, 0], [5.3638655462184879, 2.8141456582633051, 0], [5.821988795518207, 2.8141456582633064, 0], [5.5499999999999998, 3.2250000000000005, 0], [5.1858543417366949, 3.1780112044817921, 0], [5.1858543417366949, 3.6361344537815126, 0], [5.7750000000000004, 3.4500000000000006, 0], [5.363865546218487, 3.8141456582633051, 0], [5.821988795518207, 3.8141456582633064, 0], [1.75, 0.25, 0], [1.25, 0.24999999999999994, 0], [1.25, 0.75, 0], [1.75, 0.75, 0], [2.75, 0.25, 0], [2.25, 0.24999999999999994, 0], [2.25, 0.75, 0], [2.75, 0.75, 0], [3.75, 0.25, 0], [3.25, 0.25, 0], [3.25, 0.75, 0], [3.75, 0.75, 0], [4.75, 0.25, 0], [4.25, 0.24999999999999989, 0], [4.25, 0.75, 0], [4.75, 0.75, 0], [2.75, 1.25, 0], [2.25, 1.25, 0], [2.25, 1.75, 0], [2.75, 1.75, 0], [3.75, 1.25, 0], [3.25, 1.25, 0], [3.25, 1.75, 0], [3.75, 1.75, 0], [2.75, 2.25, 0], [2.25, 2.25, 0], [2.25, 2.75, 0], [2.75, 2.75, 0], [3.75, 2.25, 0], [3.25, 2.25, 0], [3.25, 2.75, 0], [3.75, 2.75, 0], [1.75, 3.25, 0], [1.25, 3.25, 0], [1.25, 3.75, 0], [1.75, 3.75, 0], [2.75, 3.25, 0], [2.25, 3.25, 0], [2.25, 3.75, 0], [2.75, 3.75, 0], [3.75, 3.25, 0], [3.25, 3.25, 0], [3.25, 3.75, 0], [3.75, 3.75, 0], [4.75, 3.25, 0], [4.25, 3.25, 0], [4.25, 3.75, 0], [4.75, 3.75, 0] ], +[0.2025, 0.14875, 0.14875, 0.2025, 0.14875, 0.14875, 0.2025, 0.14875, 0.14875, 0.2025, 0.14875, 0.14875, 0.2025, 0.14875, 0.14875, 0.2025, 0.14875, 0.14875, 0.2025, 0.14875, 0.14875, 0.2025, 0.14875, 0.14875, 0.2025, 0.14875, 0.14875, 0.2025, 0.14875, 0.14875, 0.2025, 0.14875, 0.14875, 0.2025, 0.14875, 0.14875, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, ], ), (['zeField0'], ['TS0/Mesh/ComSup1/zeField0_MM1@@][@@GAUSS'], - [[[116.,117.], [118.,119.], [120.,121.], [122.,123.], [124.,125.], [126.,127.], [128.,129.], [130.,131.], [132.,133.], [134.,135.], [136.,137.], [138.,139.], [140.,141.], [142.,143.], [144.,145.], [146.,147.], [148.,149.], [150.,151.]]] + [[[116.,117.], [118.,119.], [120.,121.], [122.,123.], [124.,125.], [126.,127.], [128.,129.], [130.,131.], [132.,133.], [134.,135.], [136.,137.], [138.,139.], [140.,141.], [142.,143.], [144.,145.], [146.,147.], [148.,149.], [150.,151.]]], + [ [0.1875, 0, 0], [0.5, 0, 0], [0.8125, 0, 0], [1.1875, 0, 0], [1.5, 0, 0], [1.8125, 0, 0], [2.1875, 0, 0], [2.5, 0, 0], [2.8125, 0, 0], [3.1875, 0, 0], [3.5, 0, 0], [3.8125, 0, 0], [4.1875, 0, 0], [4.5, 0, 0], [4.8125, 0, 0], [5.1875, 0, 0], [5.5, 0, 0], [5.8125, 0, 0]], + [0.375, 0.25, 0.375, 0.375, 0.25, 0.375, 0.375, 0.25, 0.375, 0.375, 0.25, 0.375, 0.375, 0.25, 0.375, 0.375, 0.25, 0.375] ), (['zeField0', 'zeField1'], ['TS0/Mesh/ComSup2/zeField0_MM2@@][@@GAUSS', 'TS0/Mesh/ComSup2/zeField1_MM0@@][@@GAUSS'], [[[152.,153.], [154.,155.], [156.,157.], [158.,159.], [160.,161.], [162.,163.], [164.,165.], [166.,167.], [168.,169.], [170.,171.], [172.,173.], [174.,175.], [176.,177.], [178.,179.], [180.,181.], [182.,183.]], - [[500.,501.], [502.,503.], [504.,505.], [506.,507.], [508.,509.], [510.,511.], [512.,513.], [514.,515.], [516.,517.], [518.,519.], [520.,521.], [522.,523.], [524.,525.], [526.,527.], [528.,529.], [530.,531.]]] + [[500.,501.], [502.,503.], [504.,505.], [506.,507.], [508.,509.], [510.,511.], [512.,513.], [514.,515.], [516.,517.], [518.,519.], [520.,521.], [522.,523.], [524.,525.], [526.,527.], [528.,529.], [530.,531.]]], + [ [6, 0.14999999999999999, 0], [6, 0.40000000000000002, 0], [6, 0.59999999999999998, 0], [6, 0.84999999999999998, 0], [6, 1.1499999999999999, 0], [6, 1.3999999999999999, 0], [6, 1.6000000000000001, 0], [6, 1.8500000000000001, 0], [6, 2.1499999999999999, 0], [6, 2.3999999999999999, 0], [6, 2.5999999999999996, 0], [6, 2.8499999999999996, 0], [6, 3.1499999999999999, 0], [6, 3.3999999999999999, 0], [6, 3.5999999999999996, 0], [6, 3.8499999999999996, 0]], + [0.3, 0.2, 0.2, 0.3, 0.3, 0.2, 0.2, 0.3, 0.3, 0.2, 0.2, 0.3, 0.3, 0.2, 0.2, 0.3] ), (['zeField1'], ['TS0/Mesh/ComSup3/zeField1_MM1@@][@@GAUSS'], - [[[532.,533.], [534.,535.], [536.,537.], [538.,539.], [540.,541.], [542.,543.], [544.,545.], [546.,547.], [548.,549.], [550.,551.], [552.,553.], [554.,555.], [556.,557.], [558.,559.], [560.,561.], [562.,563.], [564.,565.], [566.,567.], [568.,569.], [570.,571.], [572.,573.], [574.,575.], [576.,577.], [578.,579.], [580.,581.], [582.,583.], [584.,585.], [586.,587.], [588.,589.], [590.,591.]]] + [[[532.,533.], [534.,535.], [536.,537.], [538.,539.], [540.,541.], [542.,543.], [544.,545.], [546.,547.], [548.,549.], [550.,551.], [552.,553.], [554.,555.], [556.,557.], [558.,559.], [560.,561.], [562.,563.], [564.,565.], [566.,567.], [568.,569.], [570.,571.], [572.,573.], [574.,575.], [576.,577.], [578.,579.], [580.,581.], [582.,583.], [584.,585.], [586.,587.], [588.,589.], [590.,591.]]], + [ [5.875, 4, 0], [5.6666666666666661, 4, 0], [5.5, 4, 0], [5.3333333333333339, 4, 0], [5.125, 4, 0], [4.875, 4, 0], [4.666666666666667, 4, 0], [4.5, 4, 0], [4.333333333333333, 4, 0], [4.125, 4, 0], [3.875, 4, 0], [3.666666666666667, 4, 0], [3.5, 4, 0], [3.333333333333333, 4, 0], [3.125, 4, 0], [2.875, 4, 0], [2.6666666666666665, 4, 0], [2.5, 4, 0], [2.3333333333333335, 4, 0], [2.125, 4, 0], [1.875, 4, 0], [1.6666666666666667, 4, 0], [1.5, 4, 0], [1.3333333333333333, 4, 0], [1.125, 4, 0], [0.875, 4, 0], [0.66666666666666674, 4, 0], [0.5, 4, 0], [0.33333333333333331, 4, 0], [0.125, 4, 0]], # Center mass +[0.25, 0.16666667, 0.16666667, 0.16666667, 0.25, 0.25, 0.16666667, 0.16666667, 0.16666667, 0.25, 0.25, 0.16666667, 0.16666667, 0.16666667, 0.25, 0.25, 0.16666667, 0.16666667, 0.16666667, 0.25, 0.25, 0.16666667, 0.16666667, 0.16666667, 0.25, 0.25, 0.16666667, 0.16666667, 0.16666667, 0.25, ], ) ] -for field_names, field, refs in fields: +for field_names, field, refs, center_mass, volume in fields: testMEDReader14med.AllArrays = field testMEDReader14med.UpdatePipeline() @@ -81,18 +116,25 @@ for field_names, field, refs in fields: eLGAfieldToSurface1 = ELGAfieldToSurface(registrationName='ELGAfieldToSurface1', Input=testMEDReader14med) eLGAfieldToSurface1.UpdatePipeline() + + print(" ~> Geometry check") + test_geom(eLGAfieldToSurface1, np.array(center_mass), np.array(volume)) + eLGAfieldToSurface1.UpdatePipeline() + for field_name, ref in zip(field_names, refs): - print("Testing for leaf {} and field {}".format(field, field_name)) - test_VoroGauss(eLGAfieldToSurface1, np.array(ref), field_name) + print(" ~> Field check for leaf {} and field {}".format(field, field_name)) + test_fields(eLGAfieldToSurface1, np.array(ref), field_name) ### # Test of PG_3D.med ### # create a new 'MED Reader' -pG_3Dmed = MEDReader(registrationName='PG_3D.med', FileName='PG_3D.med') +file_name = 'PG_3D.med' +pG_3Dmed = MEDReader(registrationName='PG_3D.med', FileName=file_name) pG_3Dmed.AllArrays = ['TS0/Extruded/ComSup0/Extruded@@][@@P0', 'TS0/Extruded/ComSup0/MyFieldPG@@][@@GAUSS'] +print("Testing {}".format(file_name)) pG_3Dmed.UpdatePipeline() fields = [('Extruded', @@ -103,13 +145,19 @@ fields = [('Extruded', ), ] +ref_mass_center = [ [0.94762333243336117, 0.39671811506554749, 0.50000000000000022], [1.2856123659358007, 1.2937176028321535, 0.5], [0.60317166845782599, 1.1866543795235607, 0.49999999999999983], [0.31754942022615412, 0.3498246714004748, 0.49999999999999983], [0.67005418565845731, 0.66874178318868727, 0.49999999999999989], [1.0693170841349982, 0.90625807035941752, 0.49999999999999961], [0.92174856006299788, 1.1991729071340309, 0.49999999999999989], [1.9612936341755762, 0.50348318537092351, 0.50000000000000022], [2.2064075305015223, 1.382444756245496, 0.50000000000000033], [1.566693611827279, 1.2761686507020438, 0.4999999999999985], [1.3634692328068108, 0.45627567992173185, 0.49999999999999956], [1.6733772009850372, 0.76935226549794333, 0.49999999999999967], [2.0337072751412824, 1.0027753030374018, 0.5], [1.8684620659172329, 1.2892131924280825, 0.50000000000000122], [2.9189931607763295, 0.61023957197171785, 0.50000000000000011], [3.0712285806234125, 1.4711617857383916, 0.50000000000000011], [2.4742622100419083, 1.3656631479714718, 0.49999999999999778], [2.353412093878398, 0.56271918249158237, 0.49999999999999994], [2.6207114821730149, 0.86995893152151471, 0.50000000000000011], [2.9420980471289018, 1.0992924611855683, 0.50000000000000056], [2.7591912351279548, 1.3792470911050265, 0.50000000000000155], [3.8207265974811908, 0.71698620423188686, 0.5], [3.8800804506420783, 1.5598670948805409, 0.49999999999999983], [3.3258864500883116, 1.4551347200409823, 0.49999999999999617], [3.2873817694613265, 0.66915422974951189, 0.49999999999999917], [3.5120589475853734, 0.97056125117737468, 0.50000000000000011], [3.7944894614288533, 1.195809536997656, 0.50000000000000033], [3.5939390686114931, 1.4692736164053282, 0.50000000000000289], [1.395304554372264, 1.729157639348992, 0.50000000000000033], [1.7268056439430948, 2.4443749198861071, 0.49999999999999994], [1.139810258382987, 2.3653479708368494, 0.500000000000001], [0.85011258317843807, 1.6959098655987328, 0.5], [1.1699597425216433, 1.9493227901497743, 0.50000000000000067], [1.5237848261919678, 2.1363367169617011, 0.49999999999999933], [1.4119777502635156, 2.3723592442384547, 0.50000000000000122], [2.2689962780471022, 1.8079150043709842, 0.49999999999999917], [2.5076153246062343, 2.5050911293932603, 0.50000000000000056], [1.9633506244485432, 2.4268414277162744, 0.49999999999999983], [1.7560450450810097, 1.774354325997584, 0.49999999999999895], [2.0332907472607107, 2.0219293493545893, 0.50000000000000133], [2.3481769148963831, 2.2048536494459738, 0.49999999999999978], [2.2186967567517417, 2.4343932432074089, 0.49999999999999761], [3.0867278716275299, 1.8866611503992379, 0.50000000000000089], [3.2324609533010822, 2.5657933194298672, 0.49999999999999745], [2.7309545026101438, 2.4883076033333404, 0.50000000000000067], [2.6060084315793559, 1.8527890730110954, 0.50000000000000033], [2.8406372064211651, 2.0945307999080098, 0.49999999999999989], [3.1165699125765931, 2.2733704576614162, 0.49999999999999956], [2.969437069956574, 2.49641854082433, 0.49999999999999784], [3.8485072219503031, 1.9653944053160137, 0.49999999999999956], [3.9013512429912391, 2.6264786853937379, 0.49999999999999778], [3.4426374332033656, 2.5497409925362771, 0.49999999999999883], [3.4000089860484537, 1.9312126164197227, 0.49999999999999928], [3.592002405616304, 2.167126263118837, 0.500000000000002], [3.8289639444548826, 2.341887124172374, 0.500000000000002], [3.6642038641453327, 2.5584334431518974, 0.50000000000000233], [1.8429344844616757, 2.7816785778249042, 0.50000000000000189], [2.1679659779622691, 3.315102842720786, 0.49999999999999806], [1.6764000665080054, 3.264177701250552, 0.50000000000000133], [1.382638564442336, 2.7620638276480811, 0.49999999999999833], [1.6698418778859232, 2.9499351625225807, 0.499999999999995], [1.9782478727650958, 3.0864157767117208, 0.499999999999999], [1.9021965887377699, 3.2655916574943942, 0.49999999999999778], [2.5766471835242943, 2.8324298764678009, 0.500000000000001], [2.8087867799030333, 3.3478080591320842, 0.49999999999999872], [2.3599467247399768, 3.2976509153453302, 0.4999999999999985], [2.1485800293339876, 2.8125032620844279, 0.49999999999999944], [2.3931802866873291, 2.994538072875764, 0.49999999999999778], [2.6626428805340403, 3.1269322229991268, 0.50000000000000389], [2.5689167885216895, 3.2996199266994686, 0.50000000000000289], [3.2544173077164329, 2.883166276531282, 0.49999999999999772], [3.3936602075182747, 3.3804928180330172, 0.500000000000002], [2.9875831612767874, 3.3310847042152356, 0.50000000000000111], [2.8585647653384383, 2.8629298029738561, 0.50000000000000155], [3.0605410224700473, 3.039133864938464, 0.49999999999998429], [3.291039491843855, 3.1674484336067099, 0.49999999999999845], [3.1796671276405557, 3.3336358308716743, 0.49999999999999467], [3.8762597976675011, 2.9338849468785351, 0.50000000000000044], [3.9226037904554003, 3.4131515175119778, 0.50000000000000933], [3.5593394534481599, 3.3644682457320285, 0.49999999999999983], [3.5126042694442949, 2.913340920443205, 0.499999999999999], [3.6719304269254969, 3.0837209228569815, 0.50000000000000178], [3.863438019509482, 3.2079643608852657, 0.49999999999999906], [3.7344576117904769, 3.367636113884751, 0.4999999999999985], [2.2904814857498028, 3.5543397387612528, 0.50000000000000033], [2.6090903782288888, 3.9059611268460475, 0.50000000000000011], [2.2129758537694886, 3.8832556414138799, 0.4999999999999955], [1.9151197415276484, 3.5483358259005939, 0.499999999999999], [2.1696872340043147, 3.6706038571015411, 0.50000000000000089], [2.4326996779192753, 3.7564960681507791, 0.50000000000000544], [2.3924218055122894, 3.8789065824914792, 0.50000000000000289], [2.8841998984957486, 3.5770929042548816, 0.5], [3.1099010903342279, 3.9106625899725782, 0.49999999999999734], [2.7564812412331463, 3.8887241240831254, 0.50000000000000044], [2.5410518329576623, 3.5707774093800757, 0.49999999999999778], [2.7530259835441089, 3.6872060967081546, 0.50000000000000733], [2.9770995873755957, 3.7690116443984794, 0.50000000000000022], [2.919126559752371, 3.8849347437228317, 0.49999999999999017], [3.4220066794895612, 3.5998259701187703, 0.50000000000000044], [3.5547923202814071, 3.9153327789611954, 0.49999999999999695], [3.2441142769906111, 3.8941339824469603, 0.49999999999999833], [3.1110465701398939, 3.5932017735474724, 0.50000000000000067], [3.2803987405891997, 3.7037980999659528, 0.49999999999999811], [3.4655029620569326, 3.7815266686093754, 0.50000000000000155], [3.3898743690402782, 3.8909449368802123, 0.50000000000000033], [3.9039362336420136, 3.6225335025360939, 0.50000000000000155], [3.9438071133822659, 3.9199580829825589, 0.49999999999999684], [3.6759426175928152, 3.8994602051249383, 0.50000000000000133], [3.6251287520188078, 3.6156041296616945, 0.50000000000000044], [3.7518202398148928, 3.7203763939539267, 0.50000000000000144], [3.8979109204672659, 3.7940409553624206, 0.499999999999999], [3.8046879269395544, 3.896929832233293, 0.50000000000000044], [0.94762333243336094, 0.39671811506554761, 1.5000000000000004], [1.2856123659358001, 1.2937176028321535, 1.5], [0.60317166845782566, 1.1866543795235622, 1.5000000000000031], [0.31754942022615407, 0.34982467140047474, 1.4999999999999991], [0.67005418565845787, 0.66874178318868716, 1.5], [1.0693170841349993, 0.90625807035941641, 1.4999999999999987], [0.92174856006299843, 1.1991729071340302, 1.5], [1.9612936341755707, 0.50348318537092407, 1.5000000000000007], [2.2064075305015241, 1.3824447562454969, 1.5000000000000022], [1.5666936118272796, 1.2761686507020471, 1.4999999999999949], [1.3634692328068088, 0.45627567992173196, 1.4999999999999991], [1.6733772009850372, 0.76935226549794167, 1.4999999999999976], [2.0337072751412859, 1.0027753030374011, 1.5], [1.8684620659172388, 1.2892131924280823, 1.5000000000000036], [2.9189931607763513, 0.61023957197172063, 1.5000000000000024], [3.0712285806234059, 1.4711617857383894, 1.5000000000000013], [2.474262210041883, 1.3656631479714754, 1.4999999999999902], [2.3534120938784038, 0.5627191824915817, 1.5000000000000004], [2.6207114821730166, 0.86995893152151482, 1.5000000000000009], [2.9420980471288853, 1.0992924611855623, 1.5000000000000016], [2.7591912351279668, 1.3792470911050214, 1.5000000000000056], [3.8207265974811802, 0.71698620423188375, 1.499999999999998], [3.8800804506420734, 1.559867094880542, 1.499999999999998], [3.3258864500882677, 1.4551347200409841, 1.4999999999999918], [3.2873817694613083, 0.66915422974951322, 1.4999999999999973], [3.5120589475853627, 0.97056125117738057, 1.4999999999999998], [3.7944894614287961, 1.1958095369976693, 1.5000000000000029], [3.59393906861151, 1.4692736164053313, 1.5000000000000155], [1.3953045543722722, 1.7291576393489876, 1.5000000000000004], [1.7268056439431132, 2.444374919886104, 1.5000000000000022], [1.1398102583829886, 2.3653479708368526, 1.5000000000000042], [0.85011258317844074, 1.695909865598733, 1.5000000000000009], [1.1699597425216419, 1.9493227901497767, 1.5000000000000018], [1.5237848261919607, 2.1363367169617011, 1.4999999999999978], [1.4119777502635074, 2.3723592442384609, 1.5000000000000042], [2.2689962780471076, 1.8079150043709773, 1.4999999999999978], [2.5076153246062627, 2.5050911293932585, 1.5000000000000011], [1.963350624448543, 2.4268414277162651, 1.4999999999999953], [1.756045045081005, 1.7743543259975834, 1.4999999999999969], [2.0332907472607049, 2.0219293493545956, 1.5000000000000069], [2.348176914896404, 2.2048536494459752, 1.5000000000000007], [2.2186967567517373, 2.4343932432074138, 1.4999999999999873], [3.0867278716275921, 1.8866611503992379, 1.5000000000000016], [3.232460953301, 2.5657933194298574, 1.4999999999999858], [2.7309545026101589, 2.4883076033333489, 1.500000000000006], [2.6060084315793866, 1.8527890730110959, 1.5000000000000016], [2.8406372064211922, 2.0945307999080107, 1.5], [3.1165699125766215, 2.2733704576614162, 1.4999999999999987], [2.9694370699565842, 2.4964185408243402, 1.499999999999994], [3.8485072219503316, 1.9653944053159991, 1.5], [3.9013512429912356, 2.6264786853937485, 1.4999999999999938], [3.4426374332034193, 2.5497409925362899, 1.4999999999999925], [3.4000089860485074, 1.9312126164197219, 1.4999999999999967], [3.5920024056163022, 2.1671262631188304, 1.5000000000000051], [3.8289639444548964, 2.3418871241724011, 1.5000000000000062], [3.6642038641453101, 2.558433443151872, 1.5000000000000124], [1.8429344844616906, 2.7816785778248603, 1.5000000000000049], [2.1679659779622833, 3.315102842720806, 1.4999999999999867], [1.676400066507993, 3.2641777012505768, 1.5000000000000027], [1.38263856444233, 2.7620638276480851, 1.4999999999999969], [1.6698418778859192, 2.9499351625225949, 1.4999999999999818], [1.9782478727651263, 3.0864157767117475, 1.4999999999999991], [1.9021965887377461, 3.2655916574944626, 1.4999999999999905], [2.5766471835242943, 2.832429876467796, 1.5000000000000067], [2.8087867799030342, 3.3478080591321104, 1.4999999999999962], [2.359946724739971, 3.2976509153453177, 1.4999999999999871], [2.1485800293339685, 2.8125032620844159, 1.4999999999999989], [2.3931802866872931, 2.9945380728757933, 1.4999999999999871], [2.6626428805340914, 3.1269322229991352, 1.5000000000000107], [2.5689167885216131, 3.2996199266994837, 1.5000000000000246], [3.2544173077166034, 2.8831662765312651, 1.4999999999999958], [3.3936602075180655, 3.3804928180329901, 1.5000000000000049], [2.9875831612767874, 3.3310847042152365, 1.499999999999996], [2.8585647653384383, 2.862929802973853, 1.5000000000000051], [3.0605410224700447, 3.0391338649384441, 1.499999999999946], [3.2910394918437409, 3.1674484336067428, 1.4999999999999929], [3.1796671276407276, 3.3336358308716401, 1.4999999999999598], [3.8762597976675579, 2.9338849468785986, 1.5000000000000031], [3.9226037904554705, 3.4131515175118317, 1.5000000000000311], [3.5593394534480507, 3.3644682457320649, 1.4999999999999811], [3.5126042694442998, 2.9133409204432201, 1.5000000000000016], [3.6719304269255759, 3.0837209228569984, 1.5000000000000089], [3.8634380195095148, 3.2079643608852968, 1.4999999999999891], [3.7344576117906332, 3.3676361138847315, 1.4999999999999951], [2.290481485749821, 3.5543397387612243, 1.4999999999999969], [2.6090903782288115, 3.9059611268460737, 1.4999999999999982], [2.2129758537694983, 3.8832556414138604, 1.5000000000000036], [1.9151197415276635, 3.5483358259005606, 1.499999999999994], [2.1696872340043738, 3.6706038571015149, 1.5000000000000064], [2.432699677919318, 3.7564960681509265, 1.5000000000000053], [2.3924218055122042, 3.8789065824914539, 1.5000000000000093], [2.8841998984957997, 3.5770929042549278, 1.5], [3.1099010903337425, 3.9106625899726088, 1.4999999999999944], [2.7564812412330895, 3.8887241240829891, 1.5000000000000042], [2.5410518329576646, 3.5707774093800544, 1.4999999999999956], [2.7530259835441089, 3.6872060967081288, 1.500000000000034], [2.9770995873754651, 3.7690116443983488, 1.500000000000002], [2.9191265597521787, 3.8849347437227726, 1.4999999999999547], [3.4220066794895558, 3.59982597011878, 1.5], [3.5547923202815519, 3.9153327789611412, 1.4999999999999887], [3.2441142769906923, 3.8941339824468657, 1.5000000000000067], [3.1110465701399188, 3.5932017735474604, 1.5000000000000007], [3.28039874058928, 3.7037980999658422, 1.4999999999999925], [3.4655029620569708, 3.7815266686090645, 1.499999999999988], [3.3898743690403488, 3.8909449368805227, 1.499999999999988], [3.9039362336421823, 3.6225335025360001, 1.4999999999999873], [3.9438071133820896, 3.9199580829825589, 1.4999999999999947], [3.6759426175929462, 3.8994602051250937, 1.500000000000004], [3.6251287520186279, 3.6156041296616728, 1.5000000000000029], [3.7518202398147968, 3.7203763939539147, 1.500000000000006], [3.8979109204674489, 3.7940409553624894, 1.4999999999999942], [3.8046879269396392, 3.8969298322331238, 1.4999999999999987], [0.94762333243336028, 0.39671811506554844, 2.4999999999999991], [1.2856123659357976, 1.29371760283215, 2.4999999999999996], [0.60317166845782633, 1.1866543795235613, 2.4999999999999987], [0.3175494202261539, 0.34982467140047502, 2.5000000000000009], [0.67005418565845565, 0.66874178318868638, 2.4999999999999991], [1.069317084135001, 0.90625807035941719, 2.4999999999999991], [0.92174856006299577, 1.199172907134034, 2.5000000000000009], [1.9612936341755758, 0.50348318537092374, 2.4999999999999982], [2.2064075305015249, 1.3824447562454893, 2.5000000000000027], [1.5666936118272798, 1.2761686507020435, 2.5000000000000018], [1.3634692328068097, 0.45627567992173157, 2.4999999999999987], [1.6733772009850341, 0.76935226549794677, 2.5000000000000004], [2.0337072751412739, 1.0027753030374067, 2.5000000000000018], [1.8684620659172364, 1.2892131924280832, 2.5000000000000027], [2.918993160776314, 0.6102395719717193, 2.4999999999999973], [3.0712285806234219, 1.4711617857383916, 2.5000000000000053], [2.4742622100418856, 1.365663147971484, 2.4999999999999964], [2.3534120938783807, 0.56271918249158026, 2.4999999999999996], [2.6207114821730402, 0.86995893152152015, 2.5000000000000004], [2.9420980471289147, 1.0992924611855712, 2.5000000000000044], [2.7591912351279779, 1.3792470911050136, 2.4999999999999978], [3.8207265974811753, 0.7169862042318863, 2.5000000000000022], [3.880080450642049, 1.5598670948805389, 2.5000000000000018], [3.325886450088325, 1.4551347200409837, 2.5000000000000022], [3.2873817694613523, 0.66915422974950733, 2.5000000000000004], [3.5120589475853272, 0.97056125117737002, 2.4999999999999942], [3.7944894614287863, 1.1958095369976756, 2.5000000000000027], [3.5939390686115966, 1.4692736164053464, 2.5000000000000164], [1.3953045543722791, 1.7291576393490014, 2.4999999999999991], [1.7268056439431096, 2.4443749198861191, 2.499999999999992], [1.1398102583829821, 2.365347970836853, 2.500000000000008], [0.85011258317843852, 1.6959098655987432, 2.5000000000000022], [1.1699597425216424, 1.9493227901497725, 2.4999999999999973], [1.5237848261919529, 2.1363367169616891, 2.5000000000000013], [1.4119777502635102, 2.3723592442384933, 2.4999999999999987], [2.2689962780470916, 1.8079150043709833, 2.5000000000000009], [2.5076153246063635, 2.5050911293932909, 2.4999999999999964], [1.9633506244485393, 2.426841427716294, 2.5000000000000018], [1.7560450450810232, 1.7743543259975736, 2.5000000000000013], [2.0332907472606956, 2.0219293493545774, 2.5000000000000107], [2.3481769148963778, 2.2048536494459801, 2.4999999999999925], [2.2186967567517168, 2.4343932432074098, 2.499999999999964], [3.0867278716275921, 1.8866611503992439, 2.4999999999999991], [3.2324609533010351, 2.5657933194299463, 2.4999999999999725], [2.7309545026101425, 2.4883076033333391, 2.5000000000000409], [2.606008431579395, 1.8527890730110934, 2.5000000000000053], [2.8406372064211518, 2.0945307999080045, 2.4999999999999964], [3.1165699125766579, 2.2733704576613998, 2.5000000000000044], [2.9694370699565367, 2.4964185408243704, 2.4999999999999769], [3.8485072219502339, 1.9653944053159877, 2.5000000000000067], [3.9013512429911867, 2.6264786853937472, 2.4999999999999964], [3.4426374332034739, 2.5497409925362451, 2.4999999999999711], [3.4000089860484564, 1.9312126164197168, 2.5000000000000018], [3.5920024056163009, 2.1671262631188624, 2.500000000000008], [3.8289639444550061, 2.3418871241724233, 2.5000000000000102], [3.6642038641454362, 2.5584334431518601, 2.5000000000000284], [1.8429344844616891, 2.7816785778248896, 2.5000000000000031], [2.1679659779624019, 3.3151028427208762, 2.4999999999999667], [1.6764000665079744, 3.2641777012506101, 2.5000000000000013], [1.3826385644423203, 2.7620638276480682, 2.4999999999999911], [1.6698418778859365, 2.9499351625226171, 2.5000000000000178], [1.9782478727651034, 3.0864157767115725, 2.4999999999999938], [1.9021965887377579, 3.2655916574944355, 2.49999999999998], [2.5766471835243205, 2.8324298764677471, 2.5000000000000151], [2.8087867799029373, 3.3478080591320865, 2.5000000000000004], [2.359946724739999, 3.2976509153452707, 2.4999999999999516], [2.1485800293339685, 2.8125032620844257, 2.4999999999999987], [2.3931802866872425, 2.9945380728757951, 2.4999999999999671], [2.6626428805339533, 3.1269322229990482, 2.5000000000000187], [2.5689167885217552, 3.2996199266995361, 2.5000000000000688], [3.2544173077165501, 2.8831662765313384, 2.4999999999999907], [3.3936602075181073, 3.3804928180328342, 2.5000000000000231], [2.9875831612768309, 3.3310847042152769, 2.4999999999999951], [2.8585647653383583, 2.8629298029738339, 2.5000000000000075], [3.0605410224700171, 3.0391338649384565, 2.5000000000000351], [3.2910394918442449, 3.1674484336069049, 2.4999999999999929], [3.1796671276404629, 3.3336358308715983, 2.4999999999998925], [3.8762597976674127, 2.933884946878635, 2.5000000000000089], [3.9226037904552831, 3.413151517511932, 2.5000000000000573], [3.559339453447997, 3.3644682457321391, 2.4999999999999698], [3.5126042694444024, 2.913340920443261, 2.4999999999999853], [3.6719304269255235, 3.0837209228569811, 2.5000000000000044], [3.8634380195097635, 3.2079643608849144, 2.5000000000000098], [3.7344576117904378, 3.3676361138850051, 2.4999999999999916], [2.2904814857498299, 3.5543397387612425, 2.4999999999999876], [2.6090903782287569, 3.9059611268460439, 2.4999999999999853], [2.2129758537695241, 3.8832556414135975, 2.4999999999999885], [1.9151197415276122, 3.548335825900649, 2.499999999999972], [2.1696872340044644, 3.6706038571015505, 2.5000000000000244], [2.4326996779196191, 3.7564960681506614, 2.5000000000000053], [2.392421805512154, 3.8789065824915041, 2.5000000000000169], [2.8841998984957526, 3.5770929042548167, 2.5000000000000302], [3.1099010903336737, 3.9106625899727501, 2.4999999999999685], [2.75648124123326, 3.8887241240830686, 2.4999999999999996], [2.5410518329576428, 3.5707774093800695, 2.4999999999999627], [2.7530259835441471, 3.6872060967082798, 2.5000000000000644], [2.9770995873754078, 3.7690116443978852, 2.5], [2.9191265597522729, 3.8849347437230159, 2.499999999999909], [3.4220066794897224, 3.5998259701187529, 2.4999999999999756], [3.5547923202814182, 3.9153327789611172, 2.4999999999999996], [3.2441142769908078, 3.8941339824473071, 2.5000000000000462], [3.1110465701397647, 3.5932017735476247, 2.4999999999999951], [3.2803987405892912, 3.7037980999653715, 2.5000000000000031], [3.4655029620563949, 3.7815266686093367, 2.5000000000000031], [3.3898743690402782, 3.8909449368807021, 2.5000000000000169], [3.9039362336419212, 3.6225335025360135, 2.5000000000000244], [3.9438071133824417, 3.9199580829825589, 2.5000000000000187], [3.6759426175928054, 3.8994602051250835, 2.500000000000008], [3.625128752018655, 3.6156041296617243, 2.4999999999999947], [3.7518202398148373, 3.7203763939540027, 2.5000000000000009], [3.89791092046748, 3.7940409553618566, 2.4999999999999716], [3.804687926939653, 3.896929832233476, 2.4999999999999925], [0.94762333243336416, 0.39671811506554755, 3.5000000000000018], [1.2856123659357994, 1.2937176028321593, 3.5000000000000031], [0.6031716684578271, 1.1866543795235631, 3.4999999999999876], [0.31754942022615379, 0.34982467140047518, 3.499999999999996], [0.67005418565845754, 0.66874178318868649, 3.4999999999999987], [1.0693170841349975, 0.90625807035941519, 3.4999999999999987], [0.92174856006300498, 1.1991729071340285, 3.499999999999996], [1.9612936341755924, 0.50348318537092507, 3.5000000000000067], [2.2064075305015378, 1.38244475624551, 3.5000000000000178], [1.5666936118272796, 1.2761686507020438, 3.4999999999999742], [1.3634692328068054, 0.45627567992173107, 3.4999999999999925], [1.6733772009850356, 0.76935226549793778, 3.4999999999999876], [2.0337072751412704, 1.0027753030373998, 3.4999999999999982], [1.8684620659172446, 1.2892131924280887, 3.5000000000000187], [2.9189931607763624, 0.61023957197171852, 3.5000000000000107], [3.0712285806234285, 1.4711617857383898, 3.4999999999999902], [2.4742622100418905, 1.3656631479714536, 3.4999999999999507], [2.3534120938783971, 0.5627191824915827, 3.4999999999999964], [2.6207114821730113, 0.86995893152151282, 3.5000000000000036], [2.9420980471288529, 1.099292461185561, 3.5000000000000075], [2.7591912351280268, 1.3792470911050205, 3.5000000000000351], [3.8207265974811992, 0.7169862042318853, 3.4999999999999951], [3.8800804506420516, 1.5598670948805367, 3.4999999999999964], [3.325886450088313, 1.4551347200409834, 3.4999999999999605], [3.287381769461339, 0.66915422974951422, 3.4999999999999929], [3.5120589475853419, 0.97056125117737757, 3.4999999999999925], [3.794489461428828, 1.19580953699767, 3.4999999999999991], [3.5939390686114643, 1.4692736164053213, 3.5000000000000693], [1.3953045543722569, 1.7291576393489811, 3.5000000000000022], [1.7268056439431108, 2.444374919886092, 3.5000000000000093], [1.1398102583829848, 2.3653479708368441, 3.5000000000000089], [0.85011258317843674, 1.6959098655987417, 3.5000000000000009], [1.1699597425216446, 1.9493227901497889, 3.5000000000000084], [1.5237848261919511, 2.1363367169617513, 3.4999999999999996], [1.4119777502635158, 2.3723592442384533, 3.5000000000000271], [2.2689962780470849, 1.8079150043709571, 3.4999999999999933], [2.5076153246062267, 2.5050911293932079, 3.5000000000000058], [1.9633506244485581, 2.4268414277162913, 3.5000000000000009], [1.756045045081011, 1.7743543259976018, 3.4999999999999902], [2.0332907472607191, 2.0219293493545858, 3.5000000000000284], [2.348176914896388, 2.2048536494459801, 3.5000000000000075], [2.2186967567517799, 2.4343932432074165, 3.4999999999999303], [3.086727871627553, 1.8866611503992281, 3.5000000000000053], [3.2324609533012052, 2.5657933194298144, 3.4999999999999152], [2.7309545026101065, 2.4883076033333822, 3.5000000000000582], [2.6060084315793923, 1.8527890730110881, 3.5000000000000018], [2.8406372064212224, 2.0945307999080307, 3.5000000000000071], [3.1165699125764452, 2.2733704576614358, 3.4999999999999858], [2.9694370699566295, 2.4964185408242807, 3.4999999999999845], [3.8485072219503054, 1.9653944053160162, 3.4999999999999964], [3.9013512429911938, 2.6264786853937419, 3.4999999999999725], [3.4426374332033562, 2.5497409925362442, 3.4999999999999694], [3.4000089860483369, 1.9312126164197634, 3.4999999999999969], [3.5920024056163591, 2.1671262631188344, 3.500000000000016], [3.8289639444548156, 2.3418871241723713, 3.5000000000000173], [3.6642038641453181, 2.5584334431518885, 3.5000000000000688], [1.8429344844616937, 2.7816785778248478, 3.5000000000000213], [2.1679659779623264, 3.3151028427205889, 3.4999999999999303], [1.6764000665079972, 3.26417770125052, 3.4999999999999849], [1.382638564442314, 2.7620638276481086, 3.4999999999999942], [1.6698418778859234, 2.9499351625225456, 3.4999999999999032], [1.9782478727652961, 3.0864157767118998, 3.5000000000000049], [1.902196588737683, 3.2655916574944408, 3.4999999999999454], [2.5766471835242455, 2.8324298764677316, 3.500000000000004], [2.8087867799030972, 3.3478080591320873, 3.4999999999999818], [2.3599467247399155, 3.2976509153453497, 3.4999999999999041], [2.1485800293339761, 2.8125032620843697, 3.4999999999999907], [2.3931802866873486, 2.9945380728758337, 3.4999999999999569], [2.6626428805340234, 3.1269322229993306, 3.5000000000000262], [2.5689167885215776, 3.2996199266991875, 3.5000000000001332], [3.254417307716384, 2.8831662765312891, 3.4999999999999916], [3.3936602075183631, 3.3804928180330749, 3.5000000000000266], [2.987583161276731, 3.3310847042152267, 3.4999999999999836], [2.8585647653383965, 2.8629298029739085, 3.5000000000000049], [3.0605410224700829, 3.0391338649384618, 3.4999999999997975], [3.2910394918434371, 3.1674484336066007, 3.499999999999996], [3.1796671276408031, 3.333635830871716, 3.4999999999997913], [3.8762597976675131, 2.9338849468786012, 3.4999999999999942], [3.9226037904556139, 3.4131515175118459, 3.5000000000001279], [3.5593394534482914, 3.3644682457318242, 3.4999999999999316], [3.5126042694442918, 2.9133409204431922, 3.5000000000000044], [3.6719304269257864, 3.0837209228569757, 3.5000000000000293], [3.8634380195096143, 3.2079643608854114, 3.5000000000000253], [3.7344576117906843, 3.367636113884684, 3.4999999999999956], [2.2904814857497899, 3.5543397387612901, 3.4999999999999951], [2.6090903782285717, 3.9059611268463406, 3.500000000000016], [2.2129758537694992, 3.8832556414138666, 3.500000000000028], [1.9151197415276859, 3.5483358259005695, 3.4999999999999858], [2.1696872340042499, 3.6706038571015545, 3.5000000000000018], [2.4326996779195516, 3.7564960681506863, 3.5000000000000382], [2.3924218055123876, 3.8789065824917568, 3.5000000000000271], [2.8841998984958397, 3.577092904254763, 3.4999999999999862], [3.1099010903341222, 3.9106625899725977, 3.4999999999999973], [2.7564812412332174, 3.8887241240826294, 3.4999999999999982], [2.5410518329576406, 3.5707774093801148, 3.5000000000000191], [2.7530259835440725, 3.6872060967078593, 3.5000000000001483], [2.9770995873757795, 3.7690116443986925, 3.4999999999999845], [2.9191265597523062, 3.8849347437229427, 3.4999999999997664], [3.4220066794894444, 3.5998259701187663, 3.4999999999999654], [3.5547923202818414, 3.9153327789612136, 3.4999999999999356], [3.2441142769903375, 3.8941339824466588, 3.5000000000000098], [3.1110465701400205, 3.5932017735473507, 3.4999999999999996], [3.2803987405888613, 3.7037980999664271, 3.4999999999999574], [3.4655029620570916, 3.7815266686086901, 3.5000000000000218], [3.3898743690404456, 3.8909449368801337, 3.4999999999999796], [3.9039362336417995, 3.6225335025364989, 3.4999999999999574], [3.943807113381737, 3.9199580829826028, 3.4999999999999596], [3.6759426175933441, 3.899460205124559, 3.5000000000000635], [3.6251287520184414, 3.6156041296619916, 3.5000000000000409], [3.7518202398145579, 3.7203763939539862, 3.500000000000036], [3.8979109204676017, 3.7940409553624357, 3.5000000000000213], [3.8046879269392848, 3.8969298322328534, 3.4999999999999694]] +ref_volume = [0.2307479, 0.15204144, 0.19605815, 0.24402906, 0.25322044, 0.13539794, 0.15197308, 0.2160346, 0.14138917, 0.18234545, 0.22845154, 0.23639212, 0.12632931, 0.14146982, 0.2018309, 0.13110621, 0.16910217, 0.21340714, 0.22014299, 0.11757428, 0.13132831, 0.1881368, 0.12119257, 0.15632831, 0.19889586, 0.20447304, 0.10913285, 0.12154857, 0.1593864, 0.10358834, 0.13421564, 0.16917946, 0.17428482, 0.0929404, 0.10399294, 0.1472211, 0.09478267, 0.12285004, 0.15626754, 0.1603524, 0.08543977, 0.09529848, 0.1355654, 0.08634631, 0.11195386, 0.14388874, 0.14699917, 0.07825274, 0.08696578, 0.1244193, 0.07827927, 0.1015271, 0.13204306, 0.13422512, 0.07137931, 0.07899484, 0.1007649, 0.06436824, 0.08410863, 0.10765786, 0.1098287, 0.05832286, 0.06505681, 0.0911476, 0.05740917, 0.07509013, 0.09741154, 0.09879218, 0.05239023, 0.05817114, 0.0820399, 0.05081941, 0.06654105, 0.08769834, 0.08833485, 0.0467712, 0.05164724, 0.0734418, 0.04459897, 0.05846139, 0.07851826, 0.0784567, 0.04146578, 0.0454851, 0.0548834, 0.03438114, 0.04573713, 0.05946426, 0.05985208, 0.03154533, 0.03516467, 0.0478141, 0.02926867, 0.03906573, 0.05188354, 0.05171146, 0.0271807, 0.03008781, 0.0412544, 0.02452551, 0.03286375, 0.04483594, 0.04415003, 0.02312967, 0.02537271, 0.0352043, 0.02015167, 0.02713119, 0.03832146, 0.03716778, 0.01939224, 0.02101937, 0.2307479, 0.15204144, 0.19605815, 0.24402906, 0.25322044, 0.13539794, 0.15197308, 0.2160346, 0.14138917, 0.18234545, 0.22845154, 0.23639212, 0.12632931, 0.14146982, 0.2018309, 0.13110621, 0.16910217, 0.21340714, 0.22014299, 0.11757428, 0.13132831, 0.1881368, 0.12119257, 0.15632831, 0.19889586, 0.20447304, 0.10913285, 0.12154857, 0.1593864, 0.10358834, 0.13421564, 0.16917946, 0.17428482, 0.0929404, 0.10399294, 0.1472211, 0.09478267, 0.12285004, 0.15626754, 0.1603524, 0.08543977, 0.09529848, 0.1355654, 0.08634631, 0.11195386, 0.14388874, 0.14699917, 0.07825274, 0.08696578, 0.1244193, 0.07827927, 0.1015271, 0.13204306, 0.13422512, 0.07137931, 0.07899484, 0.1007649, 0.06436824, 0.08410863, 0.10765786, 0.1098287, 0.05832286, 0.06505681, 0.0911476, 0.05740917, 0.07509013, 0.09741154, 0.09879218, 0.05239023, 0.05817114, 0.0820399, 0.05081941, 0.06654105, 0.08769834, 0.08833485, 0.0467712, 0.05164724, 0.0734418, 0.04459897, 0.05846139, 0.07851826, 0.0784567, 0.04146578, 0.0454851, 0.0548834, 0.03438114, 0.04573713, 0.05946426, 0.05985208, 0.03154533, 0.03516467, 0.0478141, 0.02926867, 0.03906573, 0.05188354, 0.05171146, 0.0271807, 0.03008781, 0.0412544, 0.02452551, 0.03286375, 0.04483594, 0.04415003, 0.02312967, 0.02537271, 0.0352043, 0.02015167, 0.02713119, 0.03832146, 0.03716778, 0.01939224, 0.02101937, 0.2307479, 0.15204144, 0.19605815, 0.24402906, 0.25322044, 0.13539794, 0.15197308, 0.2160346, 0.14138917, 0.18234545, 0.22845154, 0.23639212, 0.12632931, 0.14146982, 0.2018309, 0.13110621, 0.16910217, 0.21340714, 0.22014299, 0.11757428, 0.13132831, 0.1881368, 0.12119257, 0.15632831, 0.19889586, 0.20447304, 0.10913285, 0.12154857, 0.1593864, 0.10358834, 0.13421564, 0.16917946, 0.17428482, 0.0929404, 0.10399294, 0.1472211, 0.09478267, 0.12285004, 0.15626754, 0.1603524, 0.08543977, 0.09529848, 0.1355654, 0.08634631, 0.11195386, 0.14388874, 0.14699917, 0.07825274, 0.08696578, 0.1244193, 0.07827927, 0.1015271, 0.13204306, 0.13422512, 0.07137931, 0.07899484, 0.1007649, 0.06436824, 0.08410863, 0.10765786, 0.1098287, 0.05832286, 0.06505681, 0.0911476, 0.05740917, 0.07509013, 0.09741154, 0.09879218, 0.05239023, 0.05817114, 0.0820399, 0.05081941, 0.06654105, 0.08769834, 0.08833485, 0.0467712, 0.05164724, 0.0734418, 0.04459897, 0.05846139, 0.07851826, 0.0784567, 0.04146578, 0.0454851, 0.0548834, 0.03438114, 0.04573713, 0.05946426, 0.05985208, 0.03154533, 0.03516467, 0.0478141, 0.02926867, 0.03906573, 0.05188354, 0.05171146, 0.0271807, 0.03008781, 0.0412544, 0.02452551, 0.03286375, 0.04483594, 0.04415003, 0.02312967, 0.02537271, 0.0352043, 0.02015167, 0.02713119, 0.03832146, 0.03716778, 0.01939224, 0.02101937, 0.2307479, 0.15204144, 0.19605815, 0.24402906, 0.25322044, 0.13539794, 0.15197308, 0.2160346, 0.14138917, 0.18234545, 0.22845154, 0.23639212, 0.12632931, 0.14146982, 0.2018309, 0.13110621, 0.16910217, 0.21340714, 0.22014299, 0.11757428, 0.13132831, 0.1881368, 0.12119257, 0.15632831, 0.19889586, 0.20447304, 0.10913285, 0.12154857, 0.1593864, 0.10358834, 0.13421564, 0.16917946, 0.17428482, 0.0929404, 0.10399294, 0.1472211, 0.09478267, 0.12285004, 0.15626754, 0.1603524, 0.08543977, 0.09529848, 0.1355654, 0.08634631, 0.11195386, 0.14388874, 0.14699917, 0.07825274, 0.08696578, 0.1244193, 0.07827927, 0.1015271, 0.13204306, 0.13422512, 0.07137931, 0.07899484, 0.1007649, 0.06436824, 0.08410863, 0.10765786, 0.1098287, 0.05832286, 0.06505681, 0.0911476, 0.05740917, 0.07509013, 0.09741154, 0.09879218, 0.05239023, 0.05817114, 0.0820399, 0.05081941, 0.06654105, 0.08769834, 0.08833485, 0.0467712, 0.05164724, 0.0734418, 0.04459897, 0.05846139, 0.07851826, 0.0784567, 0.04146578, 0.0454851, 0.0548834, 0.03438114, 0.04573713, 0.05946426, 0.05985208, 0.03154533, 0.03516467, 0.0478141, 0.02926867, 0.03906573, 0.05188354, 0.05171146, 0.0271807, 0.03008781, 0.0412544, 0.02452551, 0.03286375, 0.04483594, 0.04415003, 0.02312967, 0.02537271, 0.0352043, 0.02015167, 0.02713119, 0.03832146, 0.03716778, 0.01939224, 0.02101937] + # create a new 'ELGA field To Surface' eLGAfieldToSurface2 = ELGAfieldToSurface(registrationName='ELGAfieldToSurface2', Input=pG_3Dmed) eLGAfieldToSurface2.UpdatePipeline() +print(" ~> Geometry check") +test_geom(eLGAfieldToSurface2, np.array(ref_mass_center), np.array(ref_volume)) + for field_name, ref in fields: - print("Testing PG_3D.med for field {}".format(field_name)) - test_VoroGauss(eLGAfieldToSurface2, np.array(ref), field_name) + print(" ~> Field check for field {}".format(field_name)) + test_fields(eLGAfieldToSurface2, np.array(ref), field_name) ### # Test of VoroGauss1.med @@ -168,16 +216,24 @@ f1.setArray(arr) WriteFieldUsingAlreadyWrittenMesh(fname,f1) # create a new 'MED Reader' -voroGauss1med = MEDReader(registrationName='VoroGauss1.med', FileName='VoroGauss1.med') +voroGauss1med = MEDReader(registrationName='VoroGauss1.med', FileName=fname) voroGauss1med.AllArrays = ['TS0/mesh/ComSup0/MyFieldCell@@][@@P0', 'TS0/mesh/ComSup0/MyFieldPG@@][@@GAUSS'] voroGauss1med.AllTimeSteps = ['0000'] +print("Testing {}".format(fname)) voroGauss1med.UpdatePipeline() # create a new 'ELGA field To Surface' eLGAfieldToSurface3 = ELGAfieldToSurface(registrationName='ELGAfieldToSurface3', Input=voroGauss1med) eLGAfieldToSurface3.UpdatePipeline() +ref_mass_center = [ [2.2250000000000001, 0.22499999999999998, 0], [2.6361344537815126, 0.1780112044817927, 0], [2.1780112044817925, 0.63613445378151257, 0], [3.2250000000000001, 0.22499999999999998, 0], [3.636134453781513, 0.17801120448179258, 0], [3.178011204481793, 0.63613445378151268, 0], [4.2249999999999996, 0.22499999999999998, 0], [4.6361344537815121, 0.1780112044817927, 0], [4.1780112044817912, 0.63613445378151257, 0], [0.22499999999999998, 1.2250000000000001, 0], [0.63613445378151257, 1.1780112044817928, 0], [0.1780112044817927, 1.6361344537815126, 0], [1.2250000000000001, 1.2250000000000001, 0], [1.6361344537815126, 1.1780112044817928, 0], [1.1780112044817919, 1.6361344537815115, 0], [0.25, 0.25, 0], [0.75, 0.24999999999999994, 0], [0.75, 0.75, 0], [0.24999999999999997, 0.75, 0], [1.25, 0.25, 0], [1.75, 0.25, 0], [1.75, 0.75, 0], [1.25, 0.75, 0]] + +ref_volume = [0.2025, 0.14875, 0.14875, 0.2025, 0.14875, 0.14875, 0.2025, 0.14875, 0.14875, 0.2025, 0.14875, 0.14875, 0.2025, 0.14875, 0.14875, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, 0.25, ] + +print(" ~> Geometry check") +test_geom(eLGAfieldToSurface3, np.array(ref_mass_center), np.array(ref_volume)) + fields = [('MyFieldCell', [[0., 0.], [0., 0.], [0., 0.], [1., 1.], [1., 1.], [1., 1.], [2., 2.], [2., 2.], [2., 2.], [3., 3.], [3., 3.], [3., 3.], [4., 4.], [4., 4.], [4., 4.], [5., 5.], [5., 5.], [5., 5.], [5., 5.], [6., 6.], [6., 6.], [6., 6.], [6., 6.]] ), @@ -187,8 +243,8 @@ fields = [('MyFieldCell', ] for field_name, ref in fields: - print("Testing VoroGauss.med for field {}".format(field_name)) - test_VoroGauss(eLGAfieldToSurface3, np.array(ref), field_name) + print(" ~> Field check for field {}".format(field_name)) + test_fields(eLGAfieldToSurface3, np.array(ref), field_name) ### # Test for hexa element @@ -297,16 +353,18 @@ fieldGauss2.checkConsistencyLight(); WriteFieldUsingAlreadyWrittenMesh(file_name, fieldGauss2) # create a new 'MED Reader' -simple_meshmed = MEDReader(registrationName='simple_mesh.med', FileName='simple_mesh.med') +simple_meshmed = MEDReader(registrationName='simple_mesh.med', FileName=file_name) simple_meshmed.AllArrays = ['TS0/MESH/ComSup0/RESU____EPSI_NOEU@@][@@GAUSS', 'TS0/MESH/ComSup0/RESU____SIGM_NOEU@@][@@GAUSS'] simple_meshmed.AllTimeSteps = ['0000'] +print("Testing {}".format(file_name)) simple_meshmed.UpdatePipeline() # create a new 'ELGA field To Surface' eLGAfieldToSurface4 = ELGAfieldToSurface(registrationName='ELGAfieldToSurface4', Input=simple_meshmed) eLGAfieldToSurface4.UpdatePipeline() + fields = [('RESU____EPSI_NOEU', [ [-5.65022e-05,-5.6733e-05 , 0.000179768 , -3.56665e-07, -7.22919e-07, -0.00020561], @@ -332,6 +390,13 @@ fields = [('RESU____EPSI_NOEU', ]) ] +ref_mass_center = [ [0.030000000000000041, 0.029999999999999947, 1.3000000000000034], [0.029999999999999968, 0.042000000000000086, 1.2999999999999967], [0.04200000000000001, 0.030000000000000013, 1.3000000000000036], [0.042000000000000051, 0.042000000000000107, 1.2999999999999954], [0.030000000000000058, 0.030000000000000103, 1.499999999999992], [0.029999999999999912, 0.041999999999999808, 1.5000000000000067], [0.041999999999999836, 0.030000000000000041, 1.4999999999999989], [0.042000000000000134, 0.041999999999999961, 1.5000000000000013]] + +ref_volume = [2.88e-05, 2.88e-05, 2.88e-05, 2.88e-05, 2.88e-05, 2.88e-05, 2.88e-05, 2.88e-05] + +print(" ~> Geometry check") +test_geom(eLGAfieldToSurface4, np.array(ref_mass_center), np.array(ref_volume)) + for field_name, ref in fields: - print("Testing simple_mesh.med for field {}".format(field_name)) - test_VoroGauss(eLGAfieldToSurface4, np.array(ref), field_name) + print(" ~> Field check for field {}".format(field_name)) + test_fields(eLGAfieldToSurface4, np.array(ref), field_name) -- 2.39.2