diff --git a/dune/mmdg/mmdg.hh b/dune/mmdg/mmdg.hh index 0983c462f89e97d8d17228071061f29e10611bd6..e575f73dc56e855ac29dbf67a89a2596e2c9f4c9 100644 --- a/dune/mmdg/mmdg.hh +++ b/dune/mmdg/mmdg.hh @@ -98,6 +98,8 @@ public: } } + std::cout << "ERROR bulk: " << bulkError << "\t ERROR interface: " << sqrt(interfaceError) << "\n"; + return sqrt(bulkError * bulkError + interfaceError); } @@ -232,10 +234,6 @@ private: } } - std::cout << iFrame[0] << "\t" << (*Base::A)[iElemIdxSLE][iElemIdxSLE] << "\t" << - (*Base::A)[iElemIdxSLE+1][iElemIdxSLE] << "\t" << (*Base::A)[iElemIdxSLE][iElemIdxSLE+1] << "\t" << - (*Base::A)[iElemIdxSLE+1][iElemIdxSLE+1] << "\t" << std::endl; - //get neighboring bulk grid elements of iELem const auto elemIn = intersct.inside(); const auto elemOut = intersct.outside(); @@ -619,10 +617,6 @@ private: interfaceUpdate5; } } - - std::cout << "neighborTerms: " << (*Base::A)[neighborIdxSLE][iElemIdxSLE] << "\t" << - (*Base::A)[neighborIdxSLE][iElemIdxSLE+1] << "\t" << (*Base::A)[neighborIdxSLE+1][iElemIdxSLE] << "\t" << - (*Base::A)[neighborIdxSLE+1][iElemIdxSLE+1] << "\t\n\n"; } else //boundary interface edge { @@ -678,9 +672,6 @@ private: } } } - std::cout << iFrame[0] << "\t" << (*Base::A)[iElemIdxSLE][iElemIdxSLE] << "\t" << - (*Base::A)[iElemIdxSLE+1][iElemIdxSLE] << "\t" << (*Base::A)[iElemIdxSLE][iElemIdxSLE+1] << "\t" << - (*Base::A)[iElemIdxSLE+1][iElemIdxSLE+1] << "\t" << std::endl; } } diff --git a/src/CMakeLists.txt b/src/CMakeLists.txt index 0733f2789cfe13704721b3793ae5163ddf4c2cc5..d497cb91e2723c1e2768cabe0cb75ff2e0a77421 100644 --- a/src/CMakeLists.txt +++ b/src/CMakeLists.txt @@ -15,4 +15,4 @@ target_compile_definitions("mmdg-2d" PRIVATE GRIDDIM=2) target_link_dune_default_libraries("mmdg-2d") dune_symlink_to_source_files(FILES grids parameterMMDG.ini parameterDG.ini - convergenceTest.py) + dgAnalysis.py mmdgAnalysis.py) diff --git a/src/convergenceTest.py b/src/dgAnalysis.py similarity index 97% rename from src/convergenceTest.py rename to src/dgAnalysis.py index 92ddfa4ad6d118cb0599ee057cfe85acbd450e79..cf7c9ba8a97de93b1d11882cb3c76703daa9bd02 100644 --- a/src/convergenceTest.py +++ b/src/dgAnalysis.py @@ -8,7 +8,7 @@ from math import log # writes a dgf file for a grid with N elements per directions and dimension dim def writeDGF(N, dim): # open dgf file - file = open(join("grids", "cube" + str(dim) + "d.dgf"), "w") + file = open(join("grids", "cube" + str(dim) + "d.dgf"), "rw") file.write("DGF\n\nInterval\n") diff --git a/src/grids/mmdg2_1.msh b/src/grids/mmdg2_1.msh new file mode 100644 index 0000000000000000000000000000000000000000..abbd088c8a54f5a618d99c0bc7c411627bbadd3a --- /dev/null +++ b/src/grids/mmdg2_1.msh @@ -0,0 +1,26 @@ +$MeshFormat +2.2 0 8 +$EndMeshFormat +$Nodes +8 +1 0 0 0 +2 1 0 0 +3 1 1 0 +4 0 1 0 +5 0.5 0 0 +6 0.5 1 0 +7 0.25 0.5 0 +8 0.75 0.5 0 +$EndNodes +$Elements +9 +1 1 2 2 7 5 6 +2 2 2 1 1 1 5 7 +3 2 2 1 1 5 2 8 +4 2 2 1 1 3 6 8 +5 2 2 1 1 6 4 7 +6 2 2 1 1 4 1 7 +7 2 2 1 1 2 3 8 +8 2 2 1 1 5 6 7 +9 2 2 1 1 6 5 8 +$EndElements diff --git a/src/grids/mmdg2_2e-1.msh b/src/grids/mmdg2_2e-1.msh new file mode 100644 index 0000000000000000000000000000000000000000..4f8dbeb6ac9bd68933693752b17eceb95346fd26 --- /dev/null +++ b/src/grids/mmdg2_2e-1.msh @@ -0,0 +1,146 @@ +$MeshFormat +2.2 0 8 +$EndMeshFormat +$Nodes +52 +1 0 0 0 +2 1 0 0 +3 1 1 0 +4 0 1 0 +5 0.5 0 0 +6 0.5 1 0 +7 0.1666666666662513 0 0 +8 0.3333333333328944 0 0 +9 0.6666666666666665 0 0 +10 0.8333333333333333 0 0 +11 1 0.1999999999995579 0 +12 1 0.3999999999989749 0 +13 1 0.5999999999989468 0 +14 1 0.7999999999994734 0 +15 0.8333333333333335 1 0 +16 0.6666666666666667 1 0 +17 0.3333333333337956 1 0 +18 0.1666666666673601 1 0 +19 0 0.7999999999999998 0 +20 0 0.6000000000013869 0 +21 0 0.4000000000016644 0 +22 0 0.2000000000008322 0 +23 0.5 0.1999999999995579 0 +24 0.5 0.3999999999989749 0 +25 0.5 0.5999999999989468 0 +26 0.5 0.7999999999994734 0 +27 0.2500000000005776 0.7611111111103405 0 +28 0.75 0.2388888888889752 0 +29 0.7499999999999999 0.7611111111112133 0 +30 0.2499999999995729 0.2388888888881981 0 +31 0.1870652435573258 0.4999999999998131 0 +32 0.8129347564428507 0.4999999999998247 0 +33 0.65925000365145 0.3725910018345227 0 +34 0.3407499963486918 0.3725910018334172 0 +35 0.32806328145738 0.6461149538633681 0 +36 0.6719367185428285 0.6461149538637824 0 +37 0.1517263981417085 0.6765826715524911 0 +38 0.1538246286620958 0.3373038746544156 0 +39 0.8461753713376103 0.3373038746545046 0 +40 0.8482736018586605 0.6765826715525562 0 +41 0.1363988441350373 0.8575080805862982 0 +42 0.6376527783863527 0.1506875316447077 0 +43 0.8632514507784905 0.1448063437813714 0 +44 0.3602327691320551 0.8524301276380233 0 +45 0.3623472216134552 0.1506875316444239 0 +46 0.8636011558653323 0.8575080805863023 0 +47 0.6397672308682492 0.8524301276381734 0 +48 0.1367485492212488 0.1448063437814694 0 +49 0.249907731717487 0.09990614287876229 0 +50 0.7503437129319981 0.900648151371187 0 +51 0.7500922682821269 0.09990614287877231 0 +52 0.2496562870685305 0.9006481513712714 0 +$EndNodes +$Elements +85 +1 1 2 2 7 5 23 +2 1 2 2 7 23 24 +3 1 2 2 7 24 25 +4 1 2 2 7 25 26 +5 1 2 2 7 26 6 +6 2 2 1 1 20 21 31 +7 2 2 1 1 12 13 32 +8 2 2 1 1 24 25 34 +9 2 2 1 1 34 25 35 +10 2 2 1 1 25 33 36 +11 2 2 1 1 25 24 33 +12 2 2 1 1 20 31 37 +13 2 2 1 1 12 32 39 +14 2 2 1 1 31 21 38 +15 2 2 1 1 32 13 40 +16 2 2 1 1 31 34 35 +17 2 2 1 1 33 32 36 +18 2 2 1 1 33 23 42 +19 2 2 1 1 35 26 44 +20 2 2 1 1 23 34 45 +21 2 2 1 1 26 36 47 +22 2 2 1 1 19 37 41 +23 2 2 1 1 11 39 43 +24 2 2 1 1 40 14 46 +25 2 2 1 1 38 22 48 +26 2 2 1 1 24 23 33 +27 2 2 1 1 23 24 34 +28 2 2 1 1 25 26 35 +29 2 2 1 1 26 25 36 +30 2 2 1 1 27 35 44 +31 2 2 1 1 28 33 42 +32 2 2 1 1 36 29 47 +33 2 2 1 1 34 30 45 +34 2 2 1 1 18 4 41 +35 2 2 1 1 6 17 44 +36 2 2 1 1 5 9 42 +37 2 2 1 1 10 2 43 +38 2 2 1 1 8 5 45 +39 2 2 1 1 16 6 47 +40 2 2 1 1 3 15 46 +41 2 2 1 1 1 7 48 +42 2 2 1 1 19 20 37 +43 2 2 1 1 11 12 39 +44 2 2 1 1 13 14 40 +45 2 2 1 1 21 22 38 +46 2 2 1 1 4 19 41 +47 2 2 1 1 26 6 44 +48 2 2 1 1 23 5 42 +49 2 2 1 1 2 11 43 +50 2 2 1 1 6 26 47 +51 2 2 1 1 5 23 45 +52 2 2 1 1 22 1 48 +53 2 2 1 1 14 3 46 +54 2 2 1 1 37 27 41 +55 2 2 1 1 39 28 43 +56 2 2 1 1 30 38 48 +57 2 2 1 1 29 40 46 +58 2 2 1 1 34 31 38 +59 2 2 1 1 30 34 38 +60 2 2 1 1 33 28 39 +61 2 2 1 1 35 27 37 +62 2 2 1 1 31 35 37 +63 2 2 1 1 36 32 40 +64 2 2 1 1 32 33 39 +65 2 2 1 1 29 36 40 +66 2 2 1 1 7 8 49 +67 2 2 1 1 17 18 52 +68 2 2 1 1 9 10 51 +69 2 2 1 1 15 16 50 +70 2 2 1 1 18 41 52 +71 2 2 1 1 44 17 52 +72 2 2 1 1 10 43 51 +73 2 2 1 1 42 9 51 +74 2 2 1 1 8 45 49 +75 2 2 1 1 48 7 49 +76 2 2 1 1 46 15 50 +77 2 2 1 1 16 47 50 +78 2 2 1 1 41 27 52 +79 2 2 1 1 27 44 52 +80 2 2 1 1 43 28 51 +81 2 2 1 1 28 42 51 +82 2 2 1 1 45 30 49 +83 2 2 1 1 29 46 50 +84 2 2 1 1 47 29 50 +85 2 2 1 1 30 48 49 +$EndElements diff --git a/src/grids/mmdg2_5e-1.msh b/src/grids/mmdg2_5e-1.msh new file mode 100644 index 0000000000000000000000000000000000000000..cb4bbe9bc94e2d8a76af6fe26bb654f6fdac3d01 --- /dev/null +++ b/src/grids/mmdg2_5e-1.msh @@ -0,0 +1,40 @@ +$MeshFormat +2.2 0 8 +$EndMeshFormat +$Nodes +13 +1 0 0 0 +2 1 0 0 +3 1 1 0 +4 0 1 0 +5 0.5 0 0 +6 0.5 1 0 +7 1 0.499999999998694 0 +8 0 0.5000000000020591 0 +9 0.5 0.499999999998694 0 +10 0.75 0.749999999999347 0 +11 0.7500000000000001 0.249999999999347 0 +12 0.2499999999983174 0.2500000000010296 0 +13 0.2499999999999999 0.7500000000010295 0 +$EndNodes +$Elements +18 +1 1 2 2 7 5 9 +2 1 2 2 7 9 6 +3 2 2 1 1 8 1 12 +4 2 2 1 1 9 6 13 +5 2 2 1 1 6 9 10 +6 2 2 1 1 7 3 10 +7 2 2 1 1 5 2 11 +8 2 2 1 1 8 12 13 +9 2 2 1 1 7 10 11 +10 2 2 1 1 10 9 11 +11 2 2 1 1 6 4 13 +12 2 2 1 1 12 9 13 +13 2 2 1 1 1 5 12 +14 2 2 1 1 3 6 10 +15 2 2 1 1 2 7 11 +16 2 2 1 1 9 5 11 +17 2 2 1 1 5 9 12 +18 2 2 1 1 4 8 13 +$EndElements diff --git a/src/grids/mmdg2_5e-2.msh b/src/grids/mmdg2_5e-2.msh new file mode 100644 index 0000000000000000000000000000000000000000..343b392ba70b95ac2d0dad7c7513764fa3bff4f5 --- /dev/null +++ b/src/grids/mmdg2_5e-2.msh @@ -0,0 +1,1618 @@ +$MeshFormat +2.2 0 8 +$EndMeshFormat +$Nodes +557 +1 0 0 0 +2 1 0 0 +3 1 1 0 +4 0 1 0 +5 0.5 0 0 +6 0.5 1 0 +7 0.04999999999990733 0 0 +8 0.09999999999977893 0 0 +9 0.1499999999996332 0 0 +10 0.1999999999994874 0 0 +11 0.2499999999993471 0 0 +12 0.2999999999994734 0 0 +13 0.349999999999605 0 0 +14 0.3999999999997367 0 0 +15 0.4499999999998684 0 0 +16 0.5499999999999999 0 0 +17 0.5999999999999998 0 0 +18 0.6499999999999998 0 0 +19 0.6999999999999998 0 0 +20 0.7499999999999999 0 0 +21 0.7999999999999998 0 0 +22 0.8499999999999999 0 0 +23 0.8999999999999999 0 0 +24 0.95 0 0 +25 1 0.04999999999989965 0 +26 1 0.09999999999981467 0 +27 1 0.1499999999997036 0 +28 1 0.1999999999995579 0 +29 1 0.2499999999994121 0 +30 1 0.2999999999992664 0 +31 1 0.3499999999991206 0 +32 1 0.3999999999989749 0 +33 1 0.4499999999988292 0 +34 1 0.4999999999986943 0 +35 1 0.5499999999988151 0 +36 1 0.5999999999989468 0 +37 1 0.6499999999990784 0 +38 1 0.69999999999921 0 +39 1 0.7499999999993417 0 +40 1 0.7999999999994734 0 +41 1 0.8499999999996051 0 +42 1 0.8999999999997368 0 +43 1 0.9499999999998684 0 +44 0.9500000000000001 1 0 +45 0.9000000000000002 1 0 +46 0.8500000000000002 1 0 +47 0.8000000000000002 1 0 +48 0.7500000000000001 1 0 +49 0.7000000000000002 1 0 +50 0.6500000000000001 1 0 +51 0.6000000000000001 1 0 +52 0.55 1 0 +53 0.4499999999997918 1 0 +54 0.3999999999999999 1 0 +55 0.3500000000003467 1 0 +56 0.3000000000006934 1 0 +57 0.2500000000010293 1 0 +58 0.2000000000008322 1 0 +59 0.1500000000006241 1 0 +60 0.1000000000004161 1 0 +61 0.05000000000020799 1 0 +62 0 0.9499999999997918 0 +63 0 0.8999999999995836 0 +64 0 0.8499999999996529 0 +65 0 0.7999999999999998 0 +66 0 0.7500000000003466 0 +67 0 0.7000000000006934 0 +68 0 0.6500000000010401 0 +69 0 0.6000000000013869 0 +70 0 0.5500000000017335 0 +71 0 0.5000000000020587 0 +72 0 0.4500000000018723 0 +73 0 0.4000000000016644 0 +74 0 0.3500000000014564 0 +75 0 0.3000000000012483 0 +76 0 0.2500000000010403 0 +77 0 0.2000000000008322 0 +78 0 0.1500000000006241 0 +79 0 0.100000000000416 0 +80 0 0.05000000000020799 0 +81 0.5 0.04999999999989965 0 +82 0.5 0.09999999999981467 0 +83 0.5 0.1499999999997036 0 +84 0.5 0.1999999999995579 0 +85 0.5 0.2499999999994121 0 +86 0.5 0.2999999999992664 0 +87 0.5 0.3499999999991206 0 +88 0.5 0.3999999999989749 0 +89 0.5 0.4499999999988292 0 +90 0.5 0.4999999999986943 0 +91 0.5 0.5499999999988151 0 +92 0.5 0.5999999999989468 0 +93 0.5 0.6499999999990784 0 +94 0.5 0.69999999999921 0 +95 0.5 0.7499999999993417 0 +96 0.5 0.7999999999994734 0 +97 0.5 0.8499999999996051 0 +98 0.5 0.8999999999997368 0 +99 0.5 0.9499999999998684 0 +100 0.2497306454724049 0.4751470136344981 0 +101 0.7504107934988311 0.3250644689190832 0 +102 0.7541557205800797 0.6233938201466901 0 +103 0.2520053594985455 0.7292589899853057 0 +104 0.258615605599164 0.2345278605513121 0 +105 0.6937473361112714 0.8080620191273222 0 +106 0.6688901146021897 0.1765251653364122 0 +107 0.8302549550514978 0.4749743355016683 0 +108 0.6749114451600884 0.4726801172811524 0 +109 0.8355907867833444 0.1693130895946516 0 +110 0.8369861072050986 0.7542802070227145 0 +111 0.1569169096430552 0.5990942195936173 0 +112 0.3434659090908446 0.6005681818169457 0 +113 0.1552325726926221 0.3519953277174926 0 +114 0.34518311292793 0.8472168725159825 0 +115 0.1538011190940286 0.8461988809059644 0 +116 0.3510035083859189 0.3558167174334695 0 +117 0.1404066520724719 0.1417731488421393 0 +118 0.3577499078835535 0.1391566449214899 0 +119 0.6338611036583141 0.6879789330485307 0 +120 0.6331882096429141 0.2969800656617549 0 +121 0.8685635713205109 0.3505171197307041 0 +122 0.8753075327751748 0.5919269339162788 0 +123 0.1237499999999057 0.4750000000019656 0 +124 0.8761159691153089 0.8731714181926674 0 +125 0.1238922307336762 0.7236642356702622 0 +126 0.3854478405545074 0.7279400211252394 0 +127 0.3747033499306782 0.4772888545945694 0 +128 0.3834066468414903 0.2481457194155206 0 +129 0.6046416858129871 0.5593425261612842 0 +130 0.6121475921977569 0.887852407802275 0 +131 0.2516318230640071 0.1099964826631142 0 +132 0.7648818524035651 0.8853333066624626 0 +133 0.1132308936221418 0.2482517516451791 0 +134 0.2542733046615823 0.8934549006270143 0 +135 0.7554467235587462 0.1062282479559618 0 +136 0.6032173568660381 0.3916602059288518 0 +137 0.898176202408011 0.2563364667948786 0 +138 0.6003587344026745 0.1003587344026751 0 +139 0.7535190281825558 0.717388810671866 0 +140 0.9083680237379925 0.09649025900273114 0 +141 0.2484769011598529 0.3753262844661949 0 +142 0.6068944408107303 0.7815291013992868 0 +143 0.754477429268941 0.2245341357583543 0 +144 0.7479658018869848 0.5268929343542361 0 +145 0.7491671873270114 0.4250422231600123 0 +146 0.2513437859144103 0.5740312296851994 0 +147 0.9090642022331659 0.6784748321151453 0 +148 0.4095916161658128 0.917717836917518 0 +149 0.09040838383437991 0.9177178369175267 0 +150 0.5844975177837035 0.2059734580928828 0 +151 0.9131103272335787 0.4352772128319431 0 +152 0.08553427419351038 0.08553427419360254 0 +153 0.07927189082016492 0.5569265851861156 0 +154 0.9184605241305219 0.7984414827595488 0 +155 0.2611834894349589 0.8145473487757714 0 +156 0.9234343491381443 0.5225989249732229 0 +157 0.8223685257281469 0.6689056586758331 0 +158 0.0844784483645685 0.6460282476767579 0 +159 0.4211370055192702 0.6440124169053252 0 +160 0.09206667148229644 0.4001830571881418 0 +161 0.421395702812658 0.558927526236394 0 +162 0.08748864264775096 0.8039120959232916 0 +163 0.4152979790023028 0.8015972779383835 0 +164 0.5815746904438992 0.4675279276198586 0 +165 0.7688082384350313 0.8037172572421479 0 +166 0.4163920455366587 0.4053636685265841 0 +167 0.8185859448541919 0.2783612597916292 0 +168 0.4186391752489541 0.08312516522499955 0 +169 0.1903683254077927 0.2777096138400372 0 +170 0.6916593118303604 0.9180394399796943 0 +171 0.6831292485556847 0.3726880921720838 0 +172 0.2191922357438176 0.6480049576410591 0 +173 0.4210342562668906 0.3250039274341257 0 +174 0.9214173100192421 0.1781554228486945 0 +175 0.6788361051646973 0.07863123694197484 0 +176 0.3120711491310889 0.6747691641710769 0 +177 0.1863704837547744 0.4258327005304101 0 +178 0.3249951816792428 0.9230499381266651 0 +179 0.1791167157972947 0.9201842488858732 0 +180 0.1863705369936639 0.5306191122804027 0 +181 0.1859949031340842 0.7734233755083383 0 +182 0.1819738879796866 0.07825704808046342 0 +183 0.3228675434921317 0.4353688146646459 0 +184 0.3094980846980015 0.5250765522666316 0 +185 0.07221025911054924 0.3207230126488837 0 +186 0.07587420536369546 0.1801742736528835 0 +187 0.8267394556388385 0.07449073900673024 0 +188 0.4244723158842222 0.1838214586671966 0 +189 0.2867665718781287 0.308120480682735 0 +190 0.6866132063586964 0.5739632648535631 0 +191 0.3198355568234174 0.07417028099120296 0 +192 0.2192012871602649 0.1751684591183832 0 +193 0.5745184861956141 0.6413856382677084 0 +194 0.9151708670395629 0.9294519573141526 0 +195 0.2894053806058324 0.1705983910855964 0 +196 0.8273319050366645 0.9303124821763326 0 +197 0.803006615505567 0.3751736859673994 0 +198 0.8181833053242898 0.5468078097497288 0 +199 0.6979194341411752 0.2742977285256413 0 +200 0.5653007550919654 0.3308178566829889 0 +201 0.5640712056187398 0.718452904181254 0 +202 0.3313373677454967 0.7776539256009961 0 +203 0.9468872555310363 0.3133513019424248 0 +204 0.6704807534502293 0.7418255873621956 0 +205 0.6977547586574115 0.6726268161460279 0 +206 0.9355560018417549 0.6228293243222048 0 +207 0.5562337871368216 0.8435858218472347 0 +208 0.5607645968288267 0.9318758135060774 0 +209 0.1566623215747155 0.6643181237126197 0 +210 0.8433241290929889 0.8290902083471572 0 +211 0.1588378913583081 0.2119811909766855 0 +212 0.04933843519678052 0.7244462727669088 0 +213 0.4425575972878264 0.725062338057693 0 +214 0.04536564655275577 0.4751695871402875 0 +215 0.4514957364418514 0.4763210104900831 0 +216 0.5608519475686982 0.2603547450560277 0 +217 0.3340724599855545 0.2258891964259196 0 +218 0.5609600551289829 0.5362714393363796 0 +219 0.6510238010745221 0.2294836820219968 0 +220 0.9378175077622266 0.7448776677516196 0 +221 0.6474698430657915 0.6301001483909576 0 +222 0.9353920420129255 0.8570750753207156 0 +223 0.7745540071974564 0.1663477983591544 0 +224 0.5646751937603898 0.1494703402464654 0 +225 0.8607415785507584 0.4299076817571233 0 +226 0.8063525421639948 0.5944565707773923 0 +227 0.9403596029486629 0.3884307061665472 0 +228 0.2058565395948245 0.8428032949413863 0 +229 0.6306556348157422 0.8239973333033894 0 +230 0.0616546975535686 0.8532528295898719 0 +231 0.4384661991929864 0.8531524019040699 0 +232 0.2314152364788986 0.3187675180228893 0 +233 0.2891436061167874 0.6210135961767875 0 +234 0.3417400155235063 0.3062164487851289 0 +235 0.4426660540837872 0.26781837813827 0 +236 0.6237810381395202 0.5197729849129894 0 +237 0.712297442367088 0.862389667948732 0 +238 0.6270992710387222 0.9597482476268537 0 +239 0.5585325353512823 0.05853253535129929 0 +240 0.1898656208185297 0.7117165058899525 0 +241 0.553444637647995 0.4280180192371889 0 +242 0.9441016056895579 0.0588023102531931 0 +243 0.7448666472342904 0.9370314503624821 0 +244 0.7229004743974459 0.474875925933374 0 +245 0.2319469489069827 0.05556118373229882 0 +246 0.271900287796512 0.9465072017744118 0 +247 0.05561459546658358 0.2693698325423718 0 +248 0.730900793656638 0.05423304948219992 0 +249 0.4511500630724099 0.1319740862371132 0 +250 0.9475510925818291 0.2284692868243311 0 +251 0.6512014235882909 0.4208814055994266 0 +252 0.8484893430458376 0.2305655228571079 0 +253 0.8230009985915636 0.8806356787612888 0 +254 0.1314077928305123 0.0498227148883904 0 +255 0.6313333675280879 0.05052191110669729 0 +256 0.05675199256871557 0.5997830267310952 0 +257 0.3714452381993278 0.05052586052915341 0 +258 0.9530590468386131 0.1335415697470696 0 +259 0.9567924536778196 0.5726512913155438 0 +260 0.04913825732135779 0.1302623570493534 0 +261 0.6188554879295255 0.1528267758846017 0 +262 0.9530840684299939 0.4775267145831908 0 +263 0.6345588225452584 0.3458893462719324 0 +264 0.1437949896359608 0.2993692691862573 0 +265 0.3764496648667942 0.1918487229288493 0 +266 0.3614967578717853 0.5544373666686648 0 +267 0.8705280936861197 0.04978540839375752 0 +268 0.7045563737321343 0.1273664152885849 0 +269 0.6958404588433797 0.5252827483568742 0 +270 0.3047404401122254 0.72759165235524 0 +271 0.5570181474625969 0.5911367564670716 0 +272 0.2918592572376946 0.8582558052204114 0 +273 0.3728691879742146 0.9502508144825167 0 +274 0.1271308120265834 0.9502508144822419 0 +275 0.7746780306288465 0.4722026519330647 0 +276 0.8076956880232988 0.1209444966855129 0 +277 0.8735012467982166 0.9560521076256783 0 +278 0.4433644185499273 0.6013268116506604 0 +279 0.8684419867894355 0.3036586211801991 0 +280 0.3646974312410845 0.6488031263485685 0 +281 0.8664328217949442 0.5173327399549288 0 +282 0.8028891901822308 0.7139281664893918 0 +283 0.7724164983817836 0.6733050390692457 0 +284 0.1338845309821435 0.4236335215835123 0 +285 0.7699966250111762 0.5755107587231617 0 +286 0.2070187465143608 0.1278586187353137 0 +287 0.1381201336883055 0.7751762333845719 0 +288 0.7234757732639632 0.1798984563893256 0 +289 0.8669334049452833 0.6592383970587211 0 +290 0.3655385635997133 0.4087577849125216 0 +291 0.1364365448651305 0.5522006464915769 0 +292 0.7327408400665278 0.7699426539516551 0 +293 0.4473890458786398 0.9473890458786357 0 +294 0.05261095412134759 0.9473890458786525 0 +295 0.2165346787111144 0.9581215291132464 0 +296 0.451499546149036 0.3707978437914529 0 +297 0.1985293500551384 0.4814858820670205 0 +298 0.8801115751415789 0.146820855884953 0 +299 0.5528331217280631 0.7686588399893091 0 +300 0.6571299162505164 0.8720757820692129 0 +301 0.7821857315431981 0.04843928563560247 0 +302 0.2667284131134753 0.4281028302297756 0 +303 0.3025726002501716 0.3828006248837429 0 +304 0.766476386806906 0.2762553836523597 0 +305 0.2059904527245422 0.224042609108708 0 +306 0.8619111066920923 0.7075540385346829 0 +307 0.6304890226706285 0.7316677580879795 0 +308 0.7341636136978019 0.3753807891887549 0 +309 0.8224363556898665 0.3363119720733953 0 +310 0.1995012979057304 0.3704690745184831 0 +311 0.4549990643429248 0.5241847626397635 0 +312 0.6791448609366528 0.3130226797776725 0 +313 0.2077488177976626 0.5990284148512793 0 +314 0.9634969791162666 0.6647962052270306 0 +315 0.2356443081730677 0.5227614203069652 0 +316 0.8012460268919619 0.2078933598082844 0 +317 0.8044692578215192 0.4214253599219586 0 +318 0.04916649878468431 0.3656315737839247 0 +319 0.3029298433386956 0.1262033428429042 0 +320 0.5573009699040377 0.8888385483067451 0 +321 0.04466111616963207 0.5230771881599143 0 +322 0.1065370159204032 0.5980215315196763 0 +323 0.04464169382010792 0.6770168329746329 0 +324 0.4553583061799669 0.6770168329737294 0 +325 0.392592870847595 0.6013459049377597 0 +326 0.05010477240069269 0.05010477240069206 0 +327 0.5503132618810201 0.3700436802476798 0 +328 0.3020086457419087 0.5727878615331763 0 +329 0.04635811509283821 0.422415696490968 0 +330 0.2638187168563809 0.672292138797551 0 +331 0.1312373625304227 0.8895794489140064 0 +332 0.368437549986012 0.890535538216415 0 +333 0.2351369964778813 0.7730849041288244 0 +334 0.8884049951290789 0.7532503648889161 0 +335 0.04183941607609104 0.215601557031923 0 +336 0.6570703759674044 0.78525496735917 0 +337 0.7017603008506905 0.4258849648128589 0 +338 0.6977740964025522 0.2298573501073247 0 +339 0.04387457313486425 0.7711595381912401 0 +340 0.4557589989037846 0.7731040863121632 0 +341 0.6967234899393037 0.622801232214066 0 +342 0.1014740114293527 0.1372149083101871 0 +343 0.4580002512139542 0.2168513736253956 0 +344 0.9493954216244296 0.2735749297949621 0 +345 0.6513468777520623 0.1208434500593076 0 +346 0.8769214150546155 0.4684605668655307 0 +347 0.5434529468419101 0.1056019694931136 0 +348 0.5405868722451829 0.2211317174071556 0 +349 0.2982405437994235 0.4773296079963544 0 +350 0.9613107997224472 0.4302469267158971 0 +351 0.1062171097439629 0.3482423182685905 0 +352 0.3971409068449803 0.3652173105832698 0 +353 0.5421336475271649 0.6775695489558734 0 +354 0.2846960337371849 0.04114211357170867 0 +355 0.1363743284537482 0.1003964977809092 0 +356 0.1090966515037232 0.8449669539768018 0 +357 0.3900644619796104 0.8460601295194985 0 +358 0.4529305037747518 0.04706949622523467 0 +359 0.7890100781512441 0.9578589281581202 0 +360 0.377076942948207 0.7700757918061019 0 +361 0.2074598605850783 0.8856946913951019 0 +362 0.6044545402295802 0.4372240179847279 0 +363 0.9640966721511109 0.7836044180624668 0 +364 0.4575273982374246 0.4252155215196717 0 +365 0.8790258024494598 0.1887231148714607 0 +366 0.7847689066005419 0.757755627744457 0 +367 0.1216799029621495 0.1874403012433654 0 +368 0.04568545438027648 0.898369782086087 0 +369 0.4543145456197192 0.8983697820860351 0 +370 0.6724760869515349 0.9599457906197852 0 +371 0.2384501358606186 0.2726149820829117 0 +372 0.2991207982579813 0.2634822598372923 0 +373 0.6520854438656802 0.5484072409368088 0 +374 0.09428340926402913 0.5056465562967709 0 +375 0.4028051033108814 0.5072941377884196 0 +376 0.08747757088592459 0.6976433908346764 0 +377 0.4031834944785209 0.6857396450527289 0 +378 0.9666621742491917 0.5245553861738934 0 +379 0.3750291935522773 0.09759458148340688 0 +380 0.6152441388758014 0.2503415086586079 0 +381 0.8548389371685288 0.1096408080947227 0 +382 0.7555871078973356 0.8423771030904679 0 +383 0.08928472889566144 0.450341403116941 0 +384 0.8895092164000596 0.3889652312655001 0 +385 0.8873175750483571 0.8306243261142839 0 +386 0.2827995691228924 0.7732969263311051 0 +387 0.790086180023164 0.5121282037081513 0 +388 0.4619046354198147 0.8155707532307389 0 +389 0.03860358266723676 0.8156648927408112 0 +390 0.184929282868171 0.3225018407369485 0 +391 0.09229025448440294 0.7546698607247982 0 +392 0.3125538980338607 0.8165617476481875 0 +393 0.5443393649072944 0.4881614546772152 0 +394 0.2536106626046513 0.1945349439849111 0 +395 0.9017514550032809 0.5726721693763399 0 +396 0.402400331351418 0.1416040445003074 0 +397 0.1292500059859802 0.6358915692859026 0 +398 0.963891307840411 0.8211933762441933 0 +399 0.6782140532698846 0.03667723950617441 0 +400 0.9646220060719894 0.1749999999996308 0 +401 0.4630399712999427 0.3227240298724472 0 +402 0.9583717868342486 0.8997737418632958 0 +403 0.1564121157439772 0.25223442347293 0 +404 0.6291949651270329 0.1965171877977409 0 +405 0.4064580054113563 0.4500520102717385 0 +406 0.8179793779032289 0.7956589828316165 0 +407 0.1790657259432605 0.03672818934023053 0 +408 0.4667245260341172 0.1765293837057934 0 +409 0.3908443588316791 0.2925015531732034 0 +410 0.5934802279678203 0.6937191864605945 0 +411 0.9573375754946977 0.9576741608523883 0 +412 0.3244657813728067 0.9649253417578207 0 +413 0.1745564413072897 0.9657113184962722 0 +414 0.03337037575222886 0.1752076375471232 0 +415 0.2252762589115147 0.4095432874794913 0 +416 0.172499125577336 0.8094370402800635 0 +417 0.7263918921769503 0.5717826285542447 0 +418 0.09110816821549891 0.04104102080637414 0 +419 0.03379279842942446 0.324783775174799 0 +420 0.6079729133336805 0.6110708596976053 0 +421 0.1478709450836389 0.5064670985501989 0 +422 0.5350088009142873 0.6292026721270922 0 +423 0.1809054269338307 0.5673228711726203 0 +424 0.3400381470732255 0.468174175485003 0 +425 0.1610208047914121 0.3919391291716187 0 +426 0.2539662905479687 0.1507267064050371 0 +427 0.8256361829811653 0.03412235512314868 0 +428 0.4603653277338254 0.08785192375097509 0 +429 0.4121480762490145 0.03963467226605669 0 +430 0.2923030575347897 0.9003800449534078 0 +431 0.3249999999995393 0.03305166159618313 0 +432 0.3519665027941716 0.7009164695391679 0 +433 0.8654121180848789 0.7935575953273727 0 +434 0.04019375432587061 0.090257343473235 0 +435 0.1807549093391457 0.6339255465172312 0 +436 0.1600168547851253 0.4592738226430081 0 +437 0.728426313458826 0.9006984662383427 0 +438 0.3656908561232857 0.8099350697968218 0 +439 0.2457396779346125 0.6156825657291778 0 +440 0.463951837141102 0.5668878201049159 0 +441 0.6259494802599692 0.4770436778294656 0 +442 0.6431671385988222 0.9198361339757751 0 +443 0.1655140737265952 0.7396810337643753 0 +444 0.03862406714556109 0.636569531136371 0 +445 0.4169989865780496 0.2216971305550464 0 +446 0.2736617344758935 0.5392798757365425 0 +447 0.4152278717393257 0.7595559030479162 0 +448 0.3304210237880728 0.1742320938753097 0 +449 0.5379519167578553 0.1853151031491531 0 +450 0.03836915922793568 0.5631443525379933 0 +451 0.2238385494575892 0.6897240455242099 0 +452 0.3111345127679581 0.3410155874505181 0 +453 0.7865447001122223 0.312757946018236 0 +454 0.853689251639804 0.5691187740050729 0 +455 0.5407519362896518 0.2980694096608765 0 +456 0.913415156155579 0.8913710594477399 0 +457 0.3262942568710334 0.6362552926042857 0 +458 0.7993481573901593 0.8380001243276335 0 +459 0.7080544959997271 0.7172560627711576 0 +460 0.349087106364916 0.5098474336177241 0 +461 0.08164780269813229 0.2201675432231446 0 +462 0.739508219438075 0.1449602294982567 0 +463 0.09929234430361432 0.2913260348380919 0 +464 0.1743320074185034 0.1694463455251626 0 +465 0.3791389406847964 0.3296480474351715 0 +466 0.9620569778868129 0.3468880834731886 0 +467 0.2487196454417913 0.8562473539884282 0 +468 0.347311731835843 0.2631044525758362 0 +469 0.3252641280299949 0.8822068751496304 0 +470 0.2782080411740575 0.08141468076024576 0 +471 0.8002850709191635 0.6357538508896741 0 +472 0.867273564373962 0.9117542456091825 0 +473 0.7321183475741471 0.8166805215764001 0 +474 0.08936745812991223 0.9619290107852583 0 +475 0.4106325418705449 0.9619290107855901 0 +476 0.911368084300317 0.3407290838525748 0 +477 0.5924814845045354 0.3052235527714857 0 +478 0.6700921027989531 0.7016499769131177 0 +479 0.2680733468150821 0.3406648079178646 0 +480 0.4623950320561547 0.638253455916426 0 +481 0.538351002488784 0.9616435507955764 0 +482 0.5924000895064813 0.03857094394166447 0 +483 0.961523706655898 0.09205846106363957 0 +484 0.8250000000000005 0.9697163215458744 0 +485 0.9041179323815716 0.2175817442266915 0 +486 0.5369170077790385 0.8133135869928283 0 +487 0.7930168996910204 0.08752569232095185 0 +488 0.2187428499099531 0.09118102345565593 0 +489 0.5301160374313432 0.7249999999992764 0 +490 0.7139957142357497 0.3403183483508636 0 +491 0.9204827990943742 0.9668192659370126 0 +492 0.7174349990280539 0.09161473741718035 0 +493 0.6623824235122697 0.2656638357921663 0 +494 0.841454432548047 0.3837168417860135 0 +495 0.6621819438670427 0.6631547325251469 0 +496 0.3388861253994859 0.109281212559751 0 +497 0.3501138188309216 0.740835572085349 0 +498 0.9637710557669383 0.7124730749275328 0 +499 0.2942076551402093 0.2105441242933902 0 +500 0.7955351283730758 0.2460611911973545 0 +501 0.643026712888818 0.3827797624930736 0 +502 0.5791175872800626 0.9624371246142227 0 +503 0.6023145198345515 0.8510392075423994 0 +504 0.778050656867791 0.3532549483869845 0 +505 0.7228576216102353 0.2973070747630232 0 +506 0.2244963716586193 0.4457928157605835 0 +507 0.593239762466594 0.3513644838226894 0 +508 0.1322011345022288 0.8159382408941387 0 +509 0.5294867498279504 0.9255362910062888 0 +510 0.6404687712118804 0.08758883312766372 0 +511 0.8954857705942025 0.6236806623649044 0 +512 0.911306083537109 0.03603188343090977 0 +513 0.7914592410140857 0.9177403467916418 0 +514 0.6177474469202523 0.6543631658887457 0 +515 0.9116835584435188 0.2963193499977084 0 +516 0.1585458675076888 0.6970958702026304 0 +517 0.6492510970626899 0.5909475453757141 0 +518 0.9045215032364184 0.7217084961382757 0 +519 0.7325032055590074 0.6651403966228191 0 +520 0.9725906148557412 0.6249999999990127 0 +521 0.2137316466502739 0.5578724247984167 0 +522 0.858423369274369 0.2672304676559537 0 +523 0.7118203828090949 0.9649893431552671 0 +524 0.7727773747134168 0.3923992838540097 0 +525 0.528524325710543 0.03548949374157363 0 +526 0.5270911050160706 0.5249999999987545 0 +527 0.670147250863278 0.8378118330095142 0 +528 0.7285365265746809 0.2577238879284484 0 +529 0.1225915595403771 0.6749583903538479 0 +530 0.9680327347565868 0.03196726524337641 0 +531 0.9113447992341193 0.4846025434218488 0 +532 0.8275144170389073 0.5114130073383303 0 +533 0.6937500268135806 0.7665890812886438 0 +534 0.2106198521839702 0.746429936301799 0 +535 0.5915406831840286 0.7446163567940703 0 +536 0.6574017626604933 0.5103158188384964 0 +537 0.212134210843817 0.8026591927268767 0 +538 0.08906577150172047 0.8788191549964326 0 +539 0.4107852174620618 0.8790975603048821 0 +540 0.526243823782489 0.07499999999985718 0 +541 0.2730928832126872 0.505556104817914 0 +542 0.5965995259747046 0.9250980459719914 0 +543 0.8461320677418123 0.618708165779209 0 +544 0.5345072264052757 0.5602149963532405 0 +545 0.5883628842632853 0.509477229934693 0 +546 0.5809747520174912 0.8136873151790744 0 +547 0.1287322232975094 0.384360729348862 0 +548 0.5312731906525275 0.3297861893179864 0 +549 0.361921121901421 0.4479283279856957 0 +550 0.1896952596580733 0.6731855129633265 0 +551 0.1754943195203298 0.8768921130084664 0 +552 0.9157389889343566 0.138752026870862 0 +553 0.4725124947351448 0.9725124947351792 0 +554 0.02681885849570535 0.9731811415042946 0 +555 0.1680211244877542 0.1235463317927976 0 +556 0.2250952835560232 0.9194542174261033 0 +557 0.02510455377172086 0.02510455377187056 0 +$EndNodes +$Elements +1052 +1 1 2 2 7 5 81 +2 1 2 2 7 81 82 +3 1 2 2 7 82 83 +4 1 2 2 7 83 84 +5 1 2 2 7 84 85 +6 1 2 2 7 85 86 +7 1 2 2 7 86 87 +8 1 2 2 7 87 88 +9 1 2 2 7 88 89 +10 1 2 2 7 89 90 +11 1 2 2 7 90 91 +12 1 2 2 7 91 92 +13 1 2 2 7 92 93 +14 1 2 2 7 93 94 +15 1 2 2 7 94 95 +16 1 2 2 7 95 96 +17 1 2 2 7 96 97 +18 1 2 2 7 97 98 +19 1 2 2 7 98 99 +20 1 2 2 7 99 6 +21 2 2 1 1 129 218 545 +22 2 2 1 1 189 371 372 +23 2 2 1 1 150 224 261 +24 2 2 1 1 171 308 337 +25 2 2 1 1 236 129 545 +26 2 2 1 1 203 30 466 +27 2 2 1 1 264 351 463 +28 2 2 1 1 410 307 535 +29 2 2 1 1 30 203 344 +30 2 2 1 1 102 341 417 +31 2 2 1 1 226 198 454 +32 2 2 1 1 211 133 367 +33 2 2 1 1 109 276 381 +34 2 2 1 1 97 207 320 +35 2 2 1 1 66 67 212 +36 2 2 1 1 94 95 213 +37 2 2 1 1 71 72 214 +38 2 2 1 1 10 245 407 +39 2 2 1 1 119 307 410 +40 2 2 1 1 298 109 381 +41 2 2 1 1 147 206 314 +42 2 2 1 1 402 43 411 +43 2 2 1 1 210 124 253 +44 2 2 1 1 89 90 215 +45 2 2 1 1 224 83 347 +46 2 2 1 1 226 454 543 +47 2 2 1 1 128 217 265 +48 2 2 1 1 99 293 369 +49 2 2 1 1 294 62 368 +50 2 2 1 1 290 183 303 +51 2 2 1 1 116 290 303 +52 2 2 1 1 189 232 371 +53 2 2 1 1 98 99 369 +54 2 2 1 1 62 63 368 +55 2 2 1 1 232 141 310 +56 2 2 1 1 218 129 271 +57 2 2 1 1 351 185 463 +58 2 2 1 1 64 230 368 +59 2 2 1 1 231 97 369 +60 2 2 1 1 10 11 245 +61 2 2 1 1 85 86 235 +62 2 2 1 1 206 147 511 +63 2 2 1 1 32 227 466 +64 2 2 1 1 42 43 402 +65 2 2 1 1 144 244 275 +66 2 2 1 1 223 135 276 +67 2 2 1 1 225 151 346 +68 2 2 1 1 107 225 346 +69 2 2 1 1 73 318 329 +70 2 2 1 1 477 263 507 +71 2 2 1 1 459 292 533 +72 2 2 1 1 120 263 477 +73 2 2 1 1 139 292 459 +74 2 2 1 1 476 279 515 +75 2 2 1 1 395 206 511 +76 2 2 1 1 121 279 476 +77 2 2 1 1 259 206 395 +78 2 2 1 1 270 103 330 +79 2 2 1 1 176 270 330 +80 2 2 1 1 235 86 401 +81 2 2 1 1 219 106 338 +82 2 2 1 1 85 216 455 +83 2 2 1 1 308 145 337 +84 2 2 1 1 31 32 466 +85 2 2 1 1 224 138 261 +86 2 2 1 1 168 249 396 +87 2 2 1 1 139 459 519 +88 2 2 1 1 456 194 472 +89 2 2 1 1 194 277 472 +90 2 2 1 1 379 168 396 +91 2 2 1 1 201 410 535 +92 2 2 1 1 222 41 402 +93 2 2 1 1 233 330 439 +94 2 2 1 1 218 393 545 +95 2 2 1 1 135 248 301 +96 2 2 1 1 135 301 487 +97 2 2 1 1 233 176 330 +98 2 2 1 1 198 285 387 +99 2 2 1 1 264 113 351 +100 2 2 1 1 75 76 247 +101 2 2 1 1 109 223 276 +102 2 2 1 1 130 320 503 +103 2 2 1 1 198 226 285 +104 2 2 1 1 393 218 526 +105 2 2 1 1 367 133 461 +106 2 2 1 1 283 139 519 +107 2 2 1 1 207 97 486 +108 2 2 1 1 217 128 468 +109 2 2 1 1 285 144 387 +110 2 2 1 1 5 81 358 +111 2 2 1 1 15 5 358 +112 2 2 1 1 225 107 317 +113 2 2 1 1 63 64 368 +114 2 2 1 1 97 98 369 +115 2 2 1 1 75 247 419 +116 2 2 1 1 234 189 372 +117 2 2 1 1 232 310 390 +118 2 2 1 1 184 266 328 +119 2 2 1 1 266 112 328 +120 2 2 1 1 222 124 385 +121 2 2 1 1 154 222 385 +122 2 2 1 1 82 83 249 +123 2 2 1 1 136 241 327 +124 2 2 1 1 245 182 407 +125 2 2 1 1 124 210 385 +126 2 2 1 1 19 20 248 +127 2 2 1 1 226 102 285 +128 2 2 1 1 223 143 288 +129 2 2 1 1 341 190 417 +130 2 2 1 1 305 211 464 +131 2 2 1 1 253 132 458 +132 2 2 1 1 112 266 325 +133 2 2 1 1 266 161 325 +134 2 2 1 1 83 82 347 +135 2 2 1 1 192 305 464 +136 2 2 1 1 137 250 344 +137 2 2 1 1 375 266 460 +138 2 2 1 1 50 51 238 +139 2 2 1 1 147 498 518 +140 2 2 1 1 132 382 458 +141 2 2 1 1 212 67 323 +142 2 2 1 1 94 213 324 +143 2 2 1 1 282 110 366 +144 2 2 1 1 215 90 311 +145 2 2 1 1 89 88 241 +146 2 2 1 1 71 214 321 +147 2 2 1 1 135 223 462 +148 2 2 1 1 224 150 449 +149 2 2 1 1 124 222 456 +150 2 2 1 1 143 223 316 +151 2 2 1 1 223 109 316 +152 2 2 1 1 56 57 246 +153 2 2 1 1 214 72 329 +154 2 2 1 1 8 9 254 +155 2 2 1 1 280 377 432 +156 2 2 1 1 19 248 399 +157 2 2 1 1 89 215 364 +158 2 2 1 1 35 36 259 +159 2 2 1 1 291 153 374 +160 2 2 1 1 371 104 372 +161 2 2 1 1 133 211 403 +162 2 2 1 1 242 140 512 +163 2 2 1 1 291 374 421 +164 2 2 1 1 454 122 543 +165 2 2 1 1 66 212 339 +166 2 2 1 1 213 95 340 +167 2 2 1 1 173 235 401 +168 2 2 1 1 251 171 337 +169 2 2 1 1 216 85 348 +170 2 2 1 1 108 244 269 +171 2 2 1 1 244 144 269 +172 2 2 1 1 318 160 329 +173 2 2 1 1 263 120 312 +174 2 2 1 1 220 39 363 +175 2 2 1 1 86 85 455 +176 2 2 1 1 320 130 542 +177 2 2 1 1 222 154 398 +178 2 2 1 1 138 224 347 +179 2 2 1 1 259 36 520 +180 2 2 1 1 366 110 406 +181 2 2 1 1 190 341 517 +182 2 2 1 1 341 221 517 +183 2 2 1 1 246 57 295 +184 2 2 1 1 110 282 306 +185 2 2 1 1 282 157 306 +186 2 2 1 1 320 207 503 +187 2 2 1 1 330 172 439 +188 2 2 1 1 248 20 301 +189 2 2 1 1 126 213 447 +190 2 2 1 1 220 154 334 +191 2 2 1 1 271 129 420 +192 2 2 1 1 28 29 250 +193 2 2 1 1 213 126 377 +194 2 2 1 1 280 159 377 +195 2 2 1 1 251 108 441 +196 2 2 1 1 244 145 275 +197 2 2 1 1 17 18 255 +198 2 2 1 1 150 216 348 +199 2 2 1 1 147 314 498 +200 2 2 1 1 140 267 512 +201 2 2 1 1 138 239 482 +202 2 2 1 1 144 275 387 +203 2 2 1 1 285 102 417 +204 2 2 1 1 78 79 260 +205 2 2 1 1 211 305 403 +206 2 2 1 1 106 268 288 +207 2 2 1 1 191 257 379 +208 2 2 1 1 26 27 258 +209 2 2 1 1 319 118 448 +210 2 2 1 1 13 14 257 +211 2 2 1 1 161 266 375 +212 2 2 1 1 229 142 336 +213 2 2 1 1 486 299 546 +214 2 2 1 1 287 181 416 +215 2 2 1 1 216 150 380 +216 2 2 1 1 297 100 315 +217 2 2 1 1 180 297 315 +218 2 2 1 1 267 140 381 +219 2 2 1 1 187 267 381 +220 2 2 1 1 154 220 363 +221 2 2 1 1 300 130 503 +222 2 2 1 1 229 300 503 +223 2 2 1 1 261 138 345 +224 2 2 1 1 106 219 404 +225 2 2 1 1 39 220 498 +226 2 2 1 1 151 225 384 +227 2 2 1 1 227 32 350 +228 2 2 1 1 151 227 350 +229 2 2 1 1 150 261 404 +230 2 2 1 1 115 228 551 +231 2 2 1 1 247 185 419 +232 2 2 1 1 289 147 306 +233 2 2 1 1 157 289 306 +234 2 2 1 1 48 243 359 +235 2 2 1 1 41 222 398 +236 2 2 1 1 305 169 403 +237 2 2 1 1 319 195 426 +238 2 2 1 1 106 288 338 +239 2 2 1 1 280 112 325 +240 2 2 1 1 159 280 325 +241 2 2 1 1 234 372 468 +242 2 2 1 1 41 42 402 +243 2 2 1 1 118 319 496 +244 2 2 1 1 140 242 483 +245 2 2 1 1 137 344 515 +246 2 2 1 1 123 374 383 +247 2 2 1 1 374 214 383 +248 2 2 1 1 98 97 320 +249 2 2 1 1 296 88 364 +250 2 2 1 1 137 252 485 +251 2 2 1 1 287 416 508 +252 2 2 1 1 255 138 482 +253 2 2 1 1 237 170 300 +254 2 2 1 1 191 379 496 +255 2 2 1 1 230 64 389 +256 2 2 1 1 97 231 388 +257 2 2 1 1 362 251 441 +258 2 2 1 1 142 229 546 +259 2 2 1 1 241 136 362 +260 2 2 1 1 164 241 362 +261 2 2 1 1 100 302 349 +262 2 2 1 1 302 183 349 +263 2 2 1 1 83 224 449 +264 2 2 1 1 299 201 535 +265 2 2 1 1 102 226 471 +266 2 2 1 1 136 327 507 +267 2 2 1 1 219 380 404 +268 2 2 1 1 393 164 545 +269 2 2 1 1 299 142 546 +270 2 2 1 1 233 146 328 +271 2 2 1 1 273 178 332 +272 2 2 1 1 148 273 332 +273 2 2 1 1 274 149 331 +274 2 2 1 1 179 274 331 +275 2 2 1 1 201 299 489 +276 2 2 1 1 193 271 420 +277 2 2 1 1 73 74 318 +278 2 2 1 1 244 108 337 +279 2 2 1 1 227 151 384 +280 2 2 1 1 297 421 436 +281 2 2 1 1 85 235 343 +282 2 2 1 1 56 246 412 +283 2 2 1 1 448 217 499 +284 2 2 1 1 380 150 404 +285 2 2 1 1 206 259 520 +286 2 2 1 1 248 175 399 +287 2 2 1 1 82 249 428 +288 2 2 1 1 45 46 277 +289 2 2 1 1 33 34 262 +290 2 2 1 1 247 76 335 +291 2 2 1 1 107 275 317 +292 2 2 1 1 386 202 392 +293 2 2 1 1 89 241 393 +294 2 2 1 1 282 139 283 +295 2 2 1 1 157 282 283 +296 2 2 1 1 235 173 409 +297 2 2 1 1 128 235 409 +298 2 2 1 1 258 140 483 +299 2 2 1 1 157 471 543 +300 2 2 1 1 453 309 504 +301 2 2 1 1 28 250 400 +302 2 2 1 1 112 233 328 +303 2 2 1 1 234 116 452 +304 2 2 1 1 189 234 452 +305 2 2 1 1 87 88 296 +306 2 2 1 1 230 162 356 +307 2 2 1 1 163 231 357 +308 2 2 1 1 309 197 504 +309 2 2 1 1 127 375 460 +310 2 2 1 1 297 180 421 +311 2 2 1 1 188 265 396 +312 2 2 1 1 265 118 396 +313 2 2 1 1 108 251 337 +314 2 2 1 1 22 23 267 +315 2 2 1 1 54 55 273 +316 2 2 1 1 59 60 274 +317 2 2 1 1 281 156 395 +318 2 2 1 1 252 365 485 +319 2 2 1 1 279 121 309 +320 2 2 1 1 239 138 347 +321 2 2 1 1 312 120 493 +322 2 2 1 1 245 11 354 +323 2 2 1 1 204 459 533 +324 2 2 1 1 228 115 416 +325 2 2 1 1 212 376 391 +326 2 2 1 1 376 125 391 +327 2 2 1 1 215 375 405 +328 2 2 1 1 196 253 472 +329 2 2 1 1 317 197 494 +330 2 2 1 1 225 317 494 +331 2 2 1 1 155 228 537 +332 2 2 1 1 194 402 411 +333 2 2 1 1 160 318 351 +334 2 2 1 1 318 185 351 +335 2 2 1 1 131 319 426 +336 2 2 1 1 228 155 467 +337 2 2 1 1 498 220 518 +338 2 2 1 1 110 306 334 +339 2 2 1 1 375 127 405 +340 2 2 1 1 241 88 327 +341 2 2 1 1 20 21 301 +342 2 2 1 1 128 265 445 +343 2 2 1 1 239 16 482 +344 2 2 1 1 158 256 322 +345 2 2 1 1 115 331 356 +346 2 2 1 1 332 114 357 +347 2 2 1 1 57 58 295 +348 2 2 1 1 183 302 303 +349 2 2 1 1 302 141 303 +350 2 2 1 1 153 291 322 +351 2 2 1 1 291 111 322 +352 2 2 1 1 162 230 389 +353 2 2 1 1 231 163 388 +354 2 2 1 1 50 238 370 +355 2 2 1 1 364 215 405 +356 2 2 1 1 147 289 511 +357 2 2 1 1 8 254 418 +358 2 2 1 1 345 138 510 +359 2 2 1 1 232 189 479 +360 2 2 1 1 141 232 479 +361 2 2 1 1 16 239 525 +362 2 2 1 1 241 164 393 +363 2 2 1 1 275 145 317 +364 2 2 1 1 246 178 412 +365 2 2 1 1 253 124 472 +366 2 2 1 1 242 25 483 +367 2 2 1 1 195 448 499 +368 2 2 1 1 186 260 342 +369 2 2 1 1 186 342 367 +370 2 2 1 1 270 176 432 +371 2 2 1 1 133 247 461 +372 2 2 1 1 155 386 392 +373 2 2 1 1 156 281 531 +374 2 2 1 1 250 174 400 +375 2 2 1 1 232 169 371 +376 2 2 1 1 122 395 511 +377 2 2 1 1 109 252 316 +378 2 2 1 1 238 51 502 +379 2 2 1 1 290 116 352 +380 2 2 1 1 166 290 352 +381 2 2 1 1 166 296 364 +382 2 2 1 1 490 312 505 +383 2 2 1 1 471 226 543 +384 2 2 1 1 277 46 484 +385 2 2 1 1 459 205 519 +386 2 2 1 1 33 262 350 +387 2 2 1 1 235 128 445 +388 2 2 1 1 252 167 500 +389 2 2 1 1 169 232 390 +390 2 2 1 1 265 217 448 +391 2 2 1 1 47 48 359 +392 2 2 1 1 171 263 312 +393 2 2 1 1 152 254 355 +394 2 2 1 1 30 31 466 +395 2 2 1 1 312 199 505 +396 2 2 1 1 25 242 530 +397 2 2 1 1 145 244 337 +398 2 2 1 1 243 132 513 +399 2 2 1 1 17 255 482 +400 2 2 1 1 256 153 322 +401 2 2 1 1 184 349 460 +402 2 2 1 1 134 246 556 +403 2 2 1 1 176 233 457 +404 2 2 1 1 281 395 454 +405 2 2 1 1 237 105 473 +406 2 2 1 1 233 112 457 +407 2 2 1 1 281 346 531 +408 2 2 1 1 167 279 309 +409 2 2 1 1 161 278 325 +410 2 2 1 1 278 159 325 +411 2 2 1 1 208 320 542 +412 2 2 1 1 106 261 345 +413 2 2 1 1 349 424 460 +414 2 2 1 1 129 236 373 +415 2 2 1 1 262 151 350 +416 2 2 1 1 254 182 355 +417 2 2 1 1 195 319 448 +418 2 2 1 1 116 234 465 +419 2 2 1 1 203 476 515 +420 2 2 1 1 275 107 387 +421 2 2 1 1 166 364 405 +422 2 2 1 1 250 29 344 +423 2 2 1 1 146 233 439 +424 2 2 1 1 356 331 538 +425 2 2 1 1 332 357 539 +426 2 2 1 1 268 106 345 +427 2 2 1 1 310 113 390 +428 2 2 1 1 150 348 449 +429 2 2 1 1 132 237 382 +430 2 2 1 1 288 143 338 +431 2 2 1 1 26 258 483 +432 2 2 1 1 260 79 434 +433 2 2 1 1 134 361 467 +434 2 2 1 1 361 228 467 +435 2 2 1 1 104 305 394 +436 2 2 1 1 305 192 394 +437 2 2 1 1 269 144 417 +438 2 2 1 1 190 269 417 +439 2 2 1 1 3 44 411 +440 2 2 1 1 43 3 411 +441 2 2 1 1 105 237 527 +442 2 2 1 1 196 277 484 +443 2 2 1 1 290 166 405 +444 2 2 1 1 76 77 335 +445 2 2 1 1 131 245 470 +446 2 2 1 1 170 237 437 +447 2 2 1 1 237 132 437 +448 2 2 1 1 338 199 493 +449 2 2 1 1 219 338 493 +450 2 2 1 1 228 361 551 +451 2 2 1 1 265 188 445 +452 2 2 1 1 200 477 507 +453 2 2 1 1 173 296 352 +454 2 2 1 1 296 166 352 +455 2 2 1 1 22 267 427 +456 2 2 1 1 326 7 418 +457 2 2 1 1 264 169 390 +458 2 2 1 1 113 264 390 +459 2 2 1 1 342 117 367 +460 2 2 1 1 289 157 543 +461 2 2 1 1 7 8 418 +462 2 2 1 1 270 202 386 +463 2 2 1 1 260 152 342 +464 2 2 1 1 269 190 373 +465 2 2 1 1 257 14 429 +466 2 2 1 1 377 126 432 +467 2 2 1 1 270 432 497 +468 2 2 1 1 96 95 299 +469 2 2 1 1 276 135 487 +470 2 2 1 1 11 12 354 +471 2 2 1 1 178 246 430 +472 2 2 1 1 246 134 430 +473 2 2 1 1 253 196 513 +474 2 2 1 1 316 252 500 +475 2 2 1 1 69 256 444 +476 2 2 1 1 175 268 345 +477 2 2 1 1 170 243 523 +478 2 2 1 1 256 69 450 +479 2 2 1 1 303 141 479 +480 2 2 1 1 267 187 427 +481 2 2 1 1 142 307 336 +482 2 2 1 1 243 48 523 +483 2 2 1 1 370 238 442 +484 2 2 1 1 90 91 311 +485 2 2 1 1 334 306 518 +486 2 2 1 1 307 204 336 +487 2 2 1 1 451 103 534 +488 2 2 1 1 240 451 534 +489 2 2 1 1 262 34 378 +490 2 2 1 1 156 262 378 +491 2 2 1 1 384 225 494 +492 2 2 1 1 70 71 321 +493 2 2 1 1 380 219 493 +494 2 2 1 1 67 68 323 +495 2 2 1 1 93 94 324 +496 2 2 1 1 247 133 463 +497 2 2 1 1 185 247 463 +498 2 2 1 1 136 251 362 +499 2 2 1 1 72 73 329 +500 2 2 1 1 80 326 434 +501 2 2 1 1 204 307 478 +502 2 2 1 1 307 119 478 +503 2 2 1 1 406 110 433 +504 2 2 1 1 79 80 434 +505 2 2 1 1 249 83 408 +506 2 2 1 1 188 249 408 +507 2 2 1 1 192 286 426 +508 2 2 1 1 286 131 426 +509 2 2 1 1 246 295 556 +510 2 2 1 1 512 24 530 +511 2 2 1 1 336 105 527 +512 2 2 1 1 229 336 527 +513 2 2 1 1 16 17 482 +514 2 2 1 1 52 481 502 +515 2 2 1 1 37 38 314 +516 2 2 1 1 252 109 365 +517 2 2 1 1 84 85 343 +518 2 2 1 1 288 268 462 +519 2 2 1 1 65 66 339 +520 2 2 1 1 95 96 340 +521 2 2 1 1 286 464 555 +522 2 2 1 1 309 121 494 +523 2 2 1 1 452 303 479 +524 2 2 1 1 247 335 461 +525 2 2 1 1 281 107 346 +526 2 2 1 1 25 26 483 +527 2 2 1 1 103 270 386 +528 2 2 1 1 88 87 327 +529 2 2 1 1 175 255 399 +530 2 2 1 1 255 18 399 +531 2 2 1 1 132 253 513 +532 2 2 1 1 258 27 400 +533 2 2 1 1 174 258 400 +534 2 2 1 1 278 92 480 +535 2 2 1 1 257 191 431 +536 2 2 1 1 13 257 431 +537 2 2 1 1 327 200 507 +538 2 2 1 1 182 245 488 +539 2 2 1 1 243 170 437 +540 2 2 1 1 132 243 437 +541 2 2 1 1 314 206 520 +542 2 2 1 1 245 131 488 +543 2 2 1 1 249 168 428 +544 2 2 1 1 249 188 396 +545 2 2 1 1 276 187 381 +546 2 2 1 1 39 40 363 +547 2 2 1 1 199 312 493 +548 2 2 1 1 121 384 494 +549 2 2 1 1 92 278 440 +550 2 2 1 1 120 380 493 +551 2 2 1 1 293 53 475 +552 2 2 1 1 61 294 474 +553 2 2 1 1 248 135 492 +554 2 2 1 1 96 97 388 +555 2 2 1 1 64 65 389 +556 2 2 1 1 242 512 530 +557 2 2 1 1 342 152 355 +558 2 2 1 1 254 9 407 +559 2 2 1 1 182 254 407 +560 2 2 1 1 257 168 379 +561 2 2 1 1 85 84 348 +562 2 2 1 1 273 55 412 +563 2 2 1 1 178 273 412 +564 2 2 1 1 274 179 413 +565 2 2 1 1 59 274 413 +566 2 2 1 1 259 156 378 +567 2 2 1 1 35 259 378 +568 2 2 1 1 258 174 552 +569 2 2 1 1 286 192 464 +570 2 2 1 1 260 186 414 +571 2 2 1 1 78 260 414 +572 2 2 1 1 153 256 450 +573 2 2 1 1 290 405 549 +574 2 2 1 1 175 248 492 +575 2 2 1 1 273 148 475 +576 2 2 1 1 149 274 474 +577 2 2 1 1 210 406 433 +578 2 2 1 1 169 264 403 +579 2 2 1 1 264 133 403 +580 2 2 1 1 256 158 444 +581 2 2 1 1 154 363 398 +582 2 2 1 1 221 341 495 +583 2 2 1 1 261 106 404 +584 2 2 1 1 81 82 428 +585 2 2 1 1 14 15 429 +586 2 2 1 1 344 203 515 +587 2 2 1 1 293 148 369 +588 2 2 1 1 149 294 368 +589 2 2 1 1 271 92 544 +590 2 2 1 1 139 282 366 +591 2 2 1 1 142 299 535 +592 2 2 1 1 29 30 344 +593 2 2 1 1 167 252 522 +594 2 2 1 1 170 370 442 +595 2 2 1 1 359 243 513 +596 2 2 1 1 118 379 396 +597 2 2 1 1 464 117 555 +598 2 2 1 1 254 152 418 +599 2 2 1 1 60 61 474 +600 2 2 1 1 53 54 475 +601 2 2 1 1 156 259 395 +602 2 2 1 1 51 52 502 +603 2 2 1 1 94 93 353 +604 2 2 1 1 171 251 501 +605 2 2 1 1 250 137 485 +606 2 2 1 1 174 250 485 +607 2 2 1 1 305 104 371 +608 2 2 1 1 169 305 371 +609 2 2 1 1 367 117 464 +610 2 2 1 1 211 367 464 +611 2 2 1 1 207 486 546 +612 2 2 1 1 272 155 392 +613 2 2 1 1 114 272 392 +614 2 2 1 1 54 273 475 +615 2 2 1 1 274 60 474 +616 2 2 1 1 292 139 366 +617 2 2 1 1 152 260 434 +618 2 2 1 1 183 290 549 +619 2 2 1 1 251 136 501 +620 2 2 1 1 392 202 438 +621 2 2 1 1 252 137 522 +622 2 2 1 1 284 123 383 +623 2 2 1 1 160 284 383 +624 2 2 1 1 109 298 365 +625 2 2 1 1 298 174 365 +626 2 2 1 1 341 205 495 +627 2 2 1 1 151 262 531 +628 2 2 1 1 358 81 428 +629 2 2 1 1 15 358 429 +630 2 2 1 1 138 255 510 +631 2 2 1 1 148 293 475 +632 2 2 1 1 294 149 474 +633 2 2 1 1 168 257 429 +634 2 2 1 1 32 33 350 +635 2 2 1 1 162 287 508 +636 2 2 1 1 210 253 458 +637 2 2 1 1 68 69 444 +638 2 2 1 1 409 234 468 +639 2 2 1 1 128 409 468 +640 2 2 1 1 374 123 421 +641 2 2 1 1 255 175 510 +642 2 2 1 1 69 70 450 +643 2 2 1 1 159 278 480 +644 2 2 1 1 92 93 480 +645 2 2 1 1 299 95 489 +646 2 2 1 1 165 292 366 +647 2 2 1 1 140 258 552 +648 2 2 1 1 262 156 531 +649 2 2 1 1 397 111 435 +650 2 2 1 1 153 321 374 +651 2 2 1 1 321 214 374 +652 2 2 1 1 278 161 440 +653 2 2 1 1 92 271 422 +654 2 2 1 1 155 272 467 +655 2 2 1 1 197 309 494 +656 2 2 1 1 125 287 391 +657 2 2 1 1 287 162 391 +658 2 2 1 1 387 107 532 +659 2 2 1 1 214 329 383 +660 2 2 1 1 329 160 383 +661 2 2 1 1 91 92 440 +662 2 2 1 1 339 212 391 +663 2 2 1 1 162 339 391 +664 2 2 1 1 212 323 376 +665 2 2 1 1 323 158 376 +666 2 2 1 1 159 324 377 +667 2 2 1 1 324 213 377 +668 2 2 1 1 136 263 501 +669 2 2 1 1 108 269 536 +670 2 2 1 1 6 99 481 +671 2 2 1 1 52 6 481 +672 2 2 1 1 263 136 507 +673 2 2 1 1 215 311 375 +674 2 2 1 1 263 171 501 +675 2 2 1 1 209 397 435 +676 2 2 1 1 118 265 448 +677 2 2 1 1 333 103 386 +678 2 2 1 1 37 314 520 +679 2 2 1 1 158 322 397 +680 2 2 1 1 218 271 544 +681 2 2 1 1 133 264 463 +682 2 2 1 1 117 342 355 +683 2 2 1 1 266 184 460 +684 2 2 1 1 45 277 491 +685 2 2 1 1 277 194 491 +686 2 2 1 1 187 276 487 +687 2 2 1 1 23 24 512 +688 2 2 1 1 271 193 422 +689 2 2 1 1 267 23 512 +690 2 2 1 1 97 96 486 +691 2 2 1 1 268 175 492 +692 2 2 1 1 134 272 430 +693 2 2 1 1 283 102 471 +694 2 2 1 1 179 295 413 +695 2 2 1 1 295 58 413 +696 2 2 1 1 268 135 462 +697 2 2 1 1 38 39 498 +698 2 2 1 1 124 456 472 +699 2 2 1 1 140 298 381 +700 2 2 1 1 111 291 423 +701 2 2 1 1 291 180 423 +702 2 2 1 1 135 268 492 +703 2 2 1 1 65 339 389 +704 2 2 1 1 340 96 388 +705 2 2 1 1 88 89 364 +706 2 2 1 1 137 279 522 +707 2 2 1 1 245 354 470 +708 2 2 1 1 202 270 497 +709 2 2 1 1 114 392 438 +710 2 2 1 1 421 123 436 +711 2 2 1 1 123 284 436 +712 2 2 1 1 284 177 436 +713 2 2 1 1 296 173 401 +714 2 2 1 1 87 296 401 +715 2 2 1 1 272 114 469 +716 2 2 1 1 92 91 544 +717 2 2 1 1 272 134 467 +718 2 2 1 1 311 161 375 +719 2 2 1 1 360 126 447 +720 2 2 1 1 343 235 445 +721 2 2 1 1 317 145 524 +722 2 2 1 1 287 125 443 +723 2 2 1 1 181 287 443 +724 2 2 1 1 277 196 472 +725 2 2 1 1 48 49 523 +726 2 2 1 1 279 167 522 +727 2 2 1 1 237 300 527 +728 2 2 1 1 176 280 432 +729 2 2 1 1 279 137 515 +730 2 2 1 1 336 204 533 +731 2 2 1 1 372 104 499 +732 2 2 1 1 217 372 499 +733 2 2 1 1 163 357 438 +734 2 2 1 1 353 93 422 +735 2 2 1 1 49 50 370 +736 2 2 1 1 144 285 417 +737 2 2 1 1 112 280 457 +738 2 2 1 1 177 284 425 +739 2 2 1 1 292 165 473 +740 2 2 1 1 34 35 378 +741 2 2 1 1 360 163 438 +742 2 2 1 1 321 153 450 +743 2 2 1 1 280 176 457 +744 2 2 1 1 295 179 556 +745 2 2 1 1 333 155 537 +746 2 2 1 1 353 193 410 +747 2 2 1 1 201 353 410 +748 2 2 1 1 322 111 397 +749 2 2 1 1 198 281 454 +750 2 2 1 1 187 301 427 +751 2 2 1 1 301 21 427 +752 2 2 1 1 356 162 508 +753 2 2 1 1 158 323 444 +754 2 2 1 1 84 83 449 +755 2 2 1 1 161 311 440 +756 2 2 1 1 223 288 462 +757 2 2 1 1 424 183 549 +758 2 2 1 1 107 281 532 +759 2 2 1 1 281 198 532 +760 2 2 1 1 324 159 480 +761 2 2 1 1 180 291 421 +762 2 2 1 1 193 353 422 +763 2 2 1 1 368 230 538 +764 2 2 1 1 369 148 539 +765 2 2 1 1 149 368 538 +766 2 2 1 1 231 369 539 +767 2 2 1 1 338 143 528 +768 2 2 1 1 157 283 471 +769 2 2 1 1 481 208 502 +770 2 2 1 1 339 162 389 +771 2 2 1 1 163 340 388 +772 2 2 1 1 110 334 433 +773 2 2 1 1 284 160 547 +774 2 2 1 1 382 237 473 +775 2 2 1 1 306 147 518 +776 2 2 1 1 328 146 446 +777 2 2 1 1 102 283 519 +778 2 2 1 1 182 286 555 +779 2 2 1 1 155 333 386 +780 2 2 1 1 302 100 506 +781 2 2 1 1 131 286 488 +782 2 2 1 1 90 89 393 +783 2 2 1 1 304 167 453 +784 2 2 1 1 286 182 488 +785 2 2 1 1 130 300 442 +786 2 2 1 1 84 343 408 +787 2 2 1 1 343 188 408 +788 2 2 1 1 353 201 489 +789 2 2 1 1 94 353 489 +790 2 2 1 1 304 143 500 +791 2 2 1 1 370 170 523 +792 2 2 1 1 167 309 453 +793 2 2 1 1 141 302 415 +794 2 2 1 1 326 152 434 +795 2 2 1 1 40 41 398 +796 2 2 1 1 516 125 529 +797 2 2 1 1 152 326 418 +798 2 2 1 1 298 140 552 +799 2 2 1 1 323 68 444 +800 2 2 1 1 177 310 415 +801 2 2 1 1 310 141 415 +802 2 2 1 1 177 297 436 +803 2 2 1 1 18 19 399 +804 2 2 1 1 27 28 400 +805 2 2 1 1 308 171 490 +806 2 2 1 1 105 292 473 +807 2 2 1 1 300 170 442 +808 2 2 1 1 186 335 414 +809 2 2 1 1 335 77 414 +810 2 2 1 1 93 92 422 +811 2 2 1 1 311 91 440 +812 2 2 1 1 313 172 435 +813 2 2 1 1 313 146 439 +814 2 2 1 1 70 321 450 +815 2 2 1 1 143 316 500 +816 2 2 1 1 53 293 553 +817 2 2 1 1 293 99 553 +818 2 2 1 1 294 61 554 +819 2 2 1 1 62 294 554 +820 2 2 1 1 171 312 490 +821 2 2 1 1 87 86 548 +822 2 2 1 1 86 87 401 +823 2 2 1 1 269 373 536 +824 2 2 1 1 105 336 533 +825 2 2 1 1 297 177 506 +826 2 2 1 1 125 376 529 +827 2 2 1 1 310 177 425 +828 2 2 1 1 113 310 425 +829 2 2 1 1 146 313 521 +830 2 2 1 1 292 105 533 +831 2 2 1 1 186 367 461 +832 2 2 1 1 9 10 407 +833 2 2 1 1 318 74 419 +834 2 2 1 1 313 111 423 +835 2 2 1 1 382 165 458 +836 2 2 1 1 83 84 408 +837 2 2 1 1 93 324 480 +838 2 2 1 1 185 318 419 +839 2 2 1 1 346 151 531 +840 2 2 1 1 195 394 426 +841 2 2 1 1 330 103 451 +842 2 2 1 1 111 313 435 +843 2 2 1 1 184 328 446 +844 2 2 1 1 116 303 452 +845 2 2 1 1 101 304 453 +846 2 2 1 1 355 182 555 +847 2 2 1 1 55 56 412 +848 2 2 1 1 58 59 413 +849 2 2 1 1 205 341 519 +850 2 2 1 1 341 102 519 +851 2 2 1 1 331 149 538 +852 2 2 1 1 148 332 539 +853 2 2 1 1 172 313 439 +854 2 2 1 1 77 78 414 +855 2 2 1 1 334 154 433 +856 2 2 1 1 143 304 528 +857 2 2 1 1 202 360 438 +858 2 2 1 1 74 75 419 +859 2 2 1 1 319 131 470 +860 2 2 1 1 191 319 470 +861 2 2 1 1 96 299 486 +862 2 2 1 1 44 45 491 +863 2 2 1 1 406 210 458 +864 2 2 1 1 165 406 458 +865 2 2 1 1 5 16 525 +866 2 2 1 1 81 5 525 +867 2 2 1 1 146 315 446 +868 2 2 1 1 90 393 526 +869 2 2 1 1 145 308 524 +870 2 2 1 1 101 490 505 +871 2 2 1 1 240 516 550 +872 2 2 1 1 301 187 487 +873 2 2 1 1 100 297 506 +874 2 2 1 1 180 315 521 +875 2 2 1 1 21 22 427 +876 2 2 1 1 12 13 431 +877 2 2 1 1 174 298 552 +878 2 2 1 1 24 2 530 +879 2 2 1 1 2 25 530 +880 2 2 1 1 197 317 524 +881 2 2 1 1 122 511 543 +882 2 2 1 1 511 289 543 +883 2 2 1 1 320 208 509 +884 2 2 1 1 98 320 509 +885 2 2 1 1 300 229 527 +886 2 2 1 1 101 308 490 +887 2 2 1 1 167 304 500 +888 2 2 1 1 304 101 505 +889 2 2 1 1 130 442 542 +890 2 2 1 1 165 366 406 +891 2 2 1 1 101 453 504 +892 2 2 1 1 179 331 551 +893 2 2 1 1 332 178 469 +894 2 2 1 1 315 100 541 +895 2 2 1 1 179 361 556 +896 2 2 1 1 372 217 468 +897 2 2 1 1 191 354 431 +898 2 2 1 1 354 12 431 +899 2 2 1 1 308 101 504 +900 2 2 1 1 363 40 398 +901 2 2 1 1 173 352 465 +902 2 2 1 1 360 202 497 +903 2 2 1 1 126 360 497 +904 2 2 1 1 307 142 535 +905 2 2 1 1 181 333 537 +906 2 2 1 1 103 333 534 +907 2 2 1 1 351 113 547 +908 2 2 1 1 172 330 451 +909 2 2 1 1 340 163 447 +910 2 2 1 1 213 340 447 +911 2 2 1 1 314 38 498 +912 2 2 1 1 380 120 477 +913 2 2 1 1 216 380 477 +914 2 2 1 1 442 238 542 +915 2 2 1 1 114 332 469 +916 2 2 1 1 432 126 497 +917 2 2 1 1 349 183 424 +918 2 2 1 1 319 191 496 +919 2 2 1 1 164 441 545 +920 2 2 1 1 46 47 484 +921 2 2 1 1 315 146 521 +922 2 2 1 1 175 345 510 +923 2 2 1 1 95 94 489 +924 2 2 1 1 188 343 445 +925 2 2 1 1 395 122 454 +926 2 2 1 1 335 186 461 +927 2 2 1 1 443 240 534 +928 2 2 1 1 181 443 534 +929 2 2 1 1 164 362 441 +930 2 2 1 1 352 116 465 +931 2 2 1 1 348 84 449 +932 2 2 1 1 416 115 508 +933 2 2 1 1 516 209 550 +934 2 2 1 1 357 114 438 +935 2 2 1 1 441 236 545 +936 2 2 1 1 165 382 473 +937 2 2 1 1 327 87 548 +938 2 2 1 1 200 327 548 +939 2 2 1 1 99 98 509 +940 2 2 1 1 199 338 528 +941 2 2 1 1 168 358 428 +942 2 2 1 1 358 168 429 +943 2 2 1 1 349 184 541 +944 2 2 1 1 6 53 553 +945 2 2 1 1 99 6 553 +946 2 2 1 1 4 62 554 +947 2 2 1 1 61 4 554 +948 2 2 1 1 154 385 433 +949 2 2 1 1 36 37 520 +950 2 2 1 1 91 90 526 +951 2 2 1 1 361 179 551 +952 2 2 1 1 7 326 557 +953 2 2 1 1 326 80 557 +954 2 2 1 1 82 81 540 +955 2 2 1 1 441 108 536 +956 2 2 1 1 236 441 536 +957 2 2 1 1 80 1 557 +958 2 2 1 1 1 7 557 +959 2 2 1 1 203 466 476 +960 2 2 1 1 466 227 476 +961 2 2 1 1 200 455 477 +962 2 2 1 1 455 216 477 +963 2 2 1 1 189 452 479 +964 2 2 1 1 160 351 547 +965 2 2 1 1 163 360 447 +966 2 2 1 1 384 121 476 +967 2 2 1 1 227 384 476 +968 2 2 1 1 220 334 518 +969 2 2 1 1 354 191 470 +970 2 2 1 1 331 115 551 +971 2 2 1 1 47 359 484 +972 2 2 1 1 359 196 484 +973 2 2 1 1 333 181 534 +974 2 2 1 1 365 174 485 +975 2 2 1 1 415 302 506 +976 2 2 1 1 373 190 517 +977 2 2 1 1 86 455 548 +978 2 2 1 1 385 210 433 +979 2 2 1 1 402 194 456 +980 2 2 1 1 222 402 456 +981 2 2 1 1 221 495 514 +982 2 2 1 1 376 158 529 +983 2 2 1 1 409 173 465 +984 2 2 1 1 420 129 517 +985 2 2 1 1 221 420 517 +986 2 2 1 1 394 192 426 +987 2 2 1 1 347 82 540 +988 2 2 1 1 239 347 540 +989 2 2 1 1 115 356 508 +990 2 2 1 1 113 425 547 +991 2 2 1 1 100 349 541 +992 2 2 1 1 411 44 491 +993 2 2 1 1 410 193 514 +994 2 2 1 1 313 423 521 +995 2 2 1 1 373 236 536 +996 2 2 1 1 424 127 460 +997 2 2 1 1 158 397 529 +998 2 2 1 1 193 420 514 +999 2 2 1 1 379 118 496 +1000 2 2 1 1 196 359 513 +1001 2 2 1 1 230 356 538 +1002 2 2 1 1 357 231 539 +1003 2 2 1 1 117 355 555 +1004 2 2 1 1 416 181 537 +1005 2 2 1 1 495 119 514 +1006 2 2 1 1 177 415 506 +1007 2 2 1 1 178 430 469 +1008 2 2 1 1 209 516 529 +1009 2 2 1 1 394 195 499 +1010 2 2 1 1 423 180 521 +1011 2 2 1 1 361 134 556 +1012 2 2 1 1 234 409 465 +1013 2 2 1 1 129 373 517 +1014 2 2 1 1 49 370 523 +1015 2 2 1 1 184 446 541 +1016 2 2 1 1 420 221 514 +1017 2 2 1 1 104 394 499 +1018 2 2 1 1 451 240 550 +1019 2 2 1 1 205 459 478 +1020 2 2 1 1 430 272 469 +1021 2 2 1 1 194 411 491 +1022 2 2 1 1 198 387 532 +1023 2 2 1 1 446 315 541 +1024 2 2 1 1 435 172 550 +1025 2 2 1 1 459 204 478 +1026 2 2 1 1 205 478 495 +1027 2 2 1 1 308 504 524 +1028 2 2 1 1 304 505 528 +1029 2 2 1 1 397 209 529 +1030 2 2 1 1 229 503 546 +1031 2 2 1 1 119 410 514 +1032 2 2 1 1 405 127 549 +1033 2 2 1 1 238 502 542 +1034 2 2 1 1 502 208 542 +1035 2 2 1 1 81 525 540 +1036 2 2 1 1 525 239 540 +1037 2 2 1 1 478 119 495 +1038 2 2 1 1 209 435 550 +1039 2 2 1 1 240 443 516 +1040 2 2 1 1 208 481 509 +1041 2 2 1 1 481 99 509 +1042 2 2 1 1 443 125 516 +1043 2 2 1 1 228 416 537 +1044 2 2 1 1 425 284 547 +1045 2 2 1 1 127 424 549 +1046 2 2 1 1 455 200 548 +1047 2 2 1 1 172 451 550 +1048 2 2 1 1 504 197 524 +1049 2 2 1 1 505 199 528 +1050 2 2 1 1 503 207 546 +1051 2 2 1 1 526 218 544 +1052 2 2 1 1 91 526 544 +$EndElements diff --git a/src/grids/mmdg2_8e-2.msh b/src/grids/mmdg2_8e-2.msh new file mode 100644 index 0000000000000000000000000000000000000000..9e623fd23091c0f4bafac4896cbb266c80eefce8 --- /dev/null +++ b/src/grids/mmdg2_8e-2.msh @@ -0,0 +1,728 @@ +$MeshFormat +2.2 0 8 +$EndMeshFormat +$Nodes +254 +1 0 0 0 +2 1 0 0 +3 1 1 0 +4 0 1 0 +5 0.5 0 0 +6 0.5 1 0 +7 0.07142857142843365 0 0 +8 0.1428571428567969 0 0 +9 0.2142857142851601 0 0 +10 0.2857142857137215 0 0 +11 0.357142857142481 0 0 +12 0.4285714285712405 0 0 +13 0.5714285714285713 0 0 +14 0.6428571428571426 0 0 +15 0.714285714285714 0 0 +16 0.7857142857142855 0 0 +17 0.857142857142857 0 0 +18 0.9285714285714286 0 0 +19 1 0.0769230769229268 0 +20 1 0.1538461538458463 0 +21 1 0.230769230768699 0 +22 1 0.3076923076915516 0 +23 1 0.3846153846144044 0 +24 1 0.4615384615372571 0 +25 1 0.5384615384603233 0 +26 1 0.6153846153836027 0 +27 1 0.6923076923068823 0 +28 1 0.7692307692301616 0 +29 1 0.8461538461534411 0 +30 1 0.9230769230767206 0 +31 0.9285714285714287 1 0 +32 0.8571428571428574 1 0 +33 0.785714285714286 1 0 +34 0.7142857142857145 1 0 +35 0.642857142857143 1 0 +36 0.5714285714285714 1 0 +37 0.4285714285712302 1 0 +38 0.3571428571431544 1 0 +39 0.2857142857150783 1 0 +40 0.214285714286606 1 0 +41 0.1428571428577374 1 0 +42 0.07142857142886866 1 0 +43 0 0.9230769230766027 0 +44 0 0.8461538461535256 0 +45 0 0.7692307692309824 0 +46 0 0.6923076923084392 0 +47 0 0.6153846153858957 0 +48 0 0.5384615384633524 0 +49 0 0.4615384615403821 0 +50 0 0.3846153846169851 0 +51 0 0.307692307693588 0 +52 0 0.230769230770191 0 +53 0 0.1538461538467939 0 +54 0 0.07692307692339706 0 +55 0.5 0.0769230769229268 0 +56 0.5 0.1538461538458463 0 +57 0.5 0.230769230768699 0 +58 0.5 0.3076923076915516 0 +59 0.5 0.3846153846144044 0 +60 0.5 0.4615384615372571 0 +61 0.5 0.5384615384603233 0 +62 0.5 0.6153846153836027 0 +63 0.5 0.6923076923068823 0 +64 0.5 0.7692307692301616 0 +65 0.5 0.8461538461534411 0 +66 0.5 0.9230769230767206 0 +67 0.2504748702923796 0.754092061387872 0 +68 0.2426000948918591 0.2354221082191981 0 +69 0.7494135811856533 0.2459173722334425 0 +70 0.7426000948920916 0.7645778917809913 0 +71 0.7490114449356149 0.5000000000000427 0 +72 0.2490114449349392 0.4999999999998578 0 +73 0.8423093187523188 0.3716717484253268 0 +74 0.3417883195505164 0.3844189718856555 0 +75 0.3423093187524725 0.6283282515735447 0 +76 0.8417883195504593 0.6155810281135683 0 +77 0.156318193618302 0.6263467758966023 0 +78 0.6563181936184592 0.3736532241034638 0 +79 0.6473847972135029 0.6257050520659939 0 +80 0.1473847972132165 0.3742949479343735 0 +81 0.1393850474076573 0.8627998836229871 0 +82 0.3604274020449172 0.1375146778630111 0 +83 0.3654903671653926 0.8597801846820881 0 +84 0.8605964193935006 0.1372401517338151 0 +85 0.6344909698341953 0.1401676749523836 0 +86 0.1344909698341386 0.1401676749523491 0 +87 0.8604274020452631 0.8624853221372143 0 +88 0.634490969834195 0.8598323250476363 0 +89 0.377452649550349 0.4999999999987901 0 +90 0.8774526495506947 0.4999999999987901 0 +91 0.6250802607616314 0.491506309169891 0 +92 0.1250802607613675 0.508493690830656 0 +93 0.3843286896479221 0.7381375954412405 0 +94 0.884328689648418 0.2618624045582963 0 +95 0.3850216999228802 0.2607508768400751 0 +96 0.8850216999226399 0.7392491231601106 0 +97 0.1154257652468023 0.7385593176500636 0 +98 0.6154257652469551 0.2614406823497976 0 +99 0.1154257652472399 0.2614406823497836 0 +100 0.6154257652471548 0.7385593176502698 0 +101 0.7499999999999997 0.1058422698804232 0 +102 0.7500000000000004 0.8941577301195426 0 +103 0.2500000000008422 0.8941577301201512 0 +104 0.2499999999994409 0.1058422698800476 0 +105 0.2485484593795464 0.6566789891863812 0 +106 0.7483968466500114 0.3433223584732912 0 +107 0.2595516656999467 0.4080519534778347 0 +108 0.7595516657002777 0.5919480465220666 0 +109 0.9205592886108833 0.574027370018245 0 +110 0.4183052437436979 0.5718245110795017 0 +111 0.4205592886111812 0.4259726299800154 0 +112 0.9183052437438046 0.4281754889185974 0 +113 0.5684026865309805 0.5827705148719613 0 +114 0.06840268653091772 0.4172294851286097 0 +115 0.5810882814062074 0.4209038847730853 0 +116 0.08108828140604736 0.5790961152274097 0 +117 0.08119467987291122 0.9135031520090845 0 +118 0.4187794843072439 0.08653482198035392 0 +119 0.4191485050887926 0.9132680105667792 0 +120 0.9187918033561278 0.08651481278710943 0 +121 0.5808379686792861 0.08671314234747157 0 +122 0.08083796867926865 0.0867131423474819 0 +123 0.5808379686792743 0.9132868576525188 0 +124 0.9187794843074339 0.913465178019806 0 +125 0.684420408706469 0.07145438561433357 0 +126 0.8111703117337525 0.06651156937069343 0 +127 0.3155775013038663 0.9285399155130027 0 +128 0.6844204087064627 0.9285456143856696 0 +129 0.811146166398291 0.9334492126108707 0 +130 0.1844204087060725 0.07145438561419015 0 +131 0.1888270406670318 0.9334941499660703 0 +132 0.3111461663978382 0.06655078738909841 0 +133 0.4213618882142314 0.6627187708561466 0 +134 0.919462161593025 0.6653829414008566 0 +135 0.9213618882150689 0.3372812291423883 0 +136 0.4194621615931646 0.3346170585981534 0 +137 0.07359172777164903 0.3329806586487578 0 +138 0.5735917277717208 0.6670193413515696 0 +139 0.5787083307188202 0.3369628301808118 0 +140 0.07870833071859898 0.6630371698193979 0 +141 0.4281932468805577 0.1959655046416312 0 +142 0.9262311748915383 0.1916980119124183 0 +143 0.4278750542577779 0.8060729680085905 0 +144 0.9281932468805347 0.8040344953581195 0 +145 0.073730616944245 0.8085917333575485 0 +146 0.5720735085265558 0.1937057005049465 0 +147 0.5720735085266573 0.8062942994947769 0 +148 0.07207350852669436 0.193705700505333 0 +149 0.7698607446500706 0.427573900330096 0 +150 0.7720333525642068 0.6819221590537349 0 +151 0.2720333525639895 0.318077840946372 0 +152 0.2698824036112326 0.5724262921922655 0 +153 0.6852734211013367 0.4495856653035472 0 +154 0.1852734211011287 0.5504143346966772 0 +155 0.6871086489905062 0.7011246837146126 0 +156 0.1871086489900948 0.2988753162859 0 +157 0.1870765457145895 0.4486755436390323 0 +158 0.6870765457149579 0.5513244563612582 0 +159 0.8157912565016698 0.5474017326058701 0 +160 0.3157912565010199 0.4525982673932686 0 +161 0.1869036659824387 0.7031828363386669 0 +162 0.6869036659823601 0.2968171636608022 0 +163 0.8128436484833228 0.2969323117701576 0 +164 0.3128436484835214 0.703067688228514 0 +165 0.210863314702639 0.8258551391417666 0 +166 0.2902187145495366 0.1731031917112337 0 +167 0.7891107930879283 0.174159909235389 0 +168 0.7902187145499709 0.8268968082888815 0 +169 0.3211341185431866 0.7836224726556783 0 +170 0.6788897976150619 0.2163531077759515 0 +171 0.6788897976151452 0.7836468922241852 0 +172 0.178889797615032 0.2163531077757877 0 +173 0.3163503529490885 0.2500521268966956 0 +174 0.8163503529491521 0.7499478731035168 0 +175 0.8212066729335721 0.2355775403625497 0 +176 0.1786453160591812 0.7645480380679014 0 +177 0.3424980418300198 0.5574680421915226 0 +178 0.8424969142286576 0.4425319678301068 0 +179 0.2875485164904966 0.8412382295813878 0 +180 0.7124351616769756 0.1587440991371499 0 +181 0.2104013803631071 0.1584958135695073 0 +182 0.7104013803630734 0.8415041864305377 0 +183 0.4442158189931258 0.4998394417671529 0 +184 0.9442158189932577 0.5001605582305148 0 +185 0.5593053529497274 0.5029892749213281 0 +186 0.05930535294957766 0.497010725079311 0 +187 0.3428702874309707 0.3096183916245994 0 +188 0.8428702874308257 0.6903816083755079 0 +189 0.448324027542816 0.7307692307685217 0 +190 0.9476171481727875 0.2680025868629772 0 +191 0.4478584387614118 0.2683790997376995 0 +192 0.9478584387612634 0.7316209002611042 0 +193 0.05259082607669224 0.7325937744960247 0 +194 0.5515677880504349 0.2692307692301254 0 +195 0.05156778805062174 0.2692307692318895 0 +196 0.5515677880505699 0.7307692307685217 0 +197 0.7495500735143081 0.04784766928378261 0 +198 0.7495476097045669 0.9521483288775889 0 +199 0.2504494430589252 0.9521523328045981 0 +200 0.2495476097040356 0.04785167112224068 0 +201 0.7149623015049694 0.3985337870525996 0 +202 0.2150056194275524 0.6014665979929816 0 +203 0.7114773151030253 0.6350010305203462 0 +204 0.2114773151026159 0.3649989694799308 0 +205 0.3087999984334385 0.5109007205495283 0 +206 0.808798812128134 0.4890992899945578 0 +207 0.3617455536232299 0.2034003767142119 0 +208 0.8617455536232643 0.7965996232860165 0 +209 0.4448853843756368 0.8654410722570247 0 +210 0.9448853843757364 0.1345589277426503 0 +211 0.4448853843757088 0.1345589277427527 0 +212 0.9448853843758317 0.8654410722572937 0 +213 0.05494769872229995 0.8656475139076739 0 +214 0.5549476987224454 0.134352486092255 0 +215 0.0549476987223394 0.1343524860923481 0 +216 0.5549476987223469 0.8656475139076906 0 +217 0.865210886274782 0.3169369234740422 0 +218 0.3652108862745368 0.6830630765248613 0 +219 0.1338605521777816 0.6829843628381995 0 +220 0.6338605521778853 0.3170156371616934 0 +221 0.1971059182710059 0.8769823169658498 0 +222 0.30295054816047 0.1229732894021746 0 +223 0.802891018357075 0.1230216399185951 0 +224 0.8029505481609682 0.8770267105978868 0 +225 0.6308777348057212 0.6831020986956114 0 +226 0.1308777348055501 0.3168979013047037 0 +227 0.1282128001129501 0.9481517829549954 0 +228 0.3717871998878636 0.9481517829553525 0 +229 0.371787199887214 0.05184821704495908 0 +230 0.871787199887476 0.05184821704486077 0 +231 0.6282128001125241 0.05184821704486055 0 +232 0.8717871998874611 0.9481517829551439 0 +233 0.6282128001125395 0.9481517829551438 0 +234 0.1282128001122552 0.0518482170447043 0 +235 0.8666814298886449 0.1992999344259772 0 +236 0.1332165671538149 0.8007957302783737 0 +237 0.4542381140797564 0.9522085470614905 0 +238 0.04576188592044404 0.9522085470612922 0 +239 0.4542381140794822 0.0477914529386344 0 +240 0.9542381140796515 0.04779145293854031 0 +241 0.5457618859203484 0.04779145293853972 0 +242 0.04576188592018556 0.04779145293861268 0 +243 0.545761885920217 0.9522085470614489 0 +244 0.9542381140797827 0.9522085470614482 0 +245 0.136681674516113 0.5686699633542317 0 +246 0.6366816745163169 0.4313300366460424 0 +247 0.365527713158599 0.442128725417602 0 +248 0.8655277131588617 0.5578712745805547 0 +249 0.625220010305692 0.2029167913957698 0 +250 0.625220010305788 0.797083208604217 0 +251 0.1252200103057762 0.2029167913958133 0 +252 0.3747070574035698 0.7969033051968992 0 +253 0.1251982739630294 0.4471225722197947 0 +254 0.6251982739632321 0.5528774277808349 0 +$EndNodes +$Elements +465 +1 1 2 2 7 5 55 +2 1 2 2 7 55 56 +3 1 2 2 7 56 57 +4 1 2 2 7 57 58 +5 1 2 2 7 58 59 +6 1 2 2 7 59 60 +7 1 2 2 7 60 61 +8 1 2 2 7 61 62 +9 1 2 2 7 62 63 +10 1 2 2 7 63 64 +11 1 2 2 7 64 65 +12 1 2 2 7 65 66 +13 1 2 2 7 66 6 +14 2 2 1 1 153 71 158 +15 2 2 1 1 72 154 157 +16 2 2 1 1 154 92 157 +17 2 2 1 1 91 153 158 +18 2 2 1 1 114 50 137 +19 2 2 1 1 62 113 138 +20 2 2 1 1 47 116 140 +21 2 2 1 1 115 59 139 +22 2 2 1 1 109 26 134 +23 2 2 1 1 110 62 133 +24 2 2 1 1 23 112 135 +25 2 2 1 1 59 111 136 +26 2 2 1 1 106 73 149 +27 2 2 1 1 75 105 152 +28 2 2 1 1 74 107 151 +29 2 2 1 1 108 76 150 +30 2 2 1 1 73 106 163 +31 2 2 1 1 105 75 164 +32 2 2 1 1 108 71 159 +33 2 2 1 1 72 107 160 +34 2 2 1 1 77 105 161 +35 2 2 1 1 106 78 162 +36 2 2 1 1 107 72 157 +37 2 2 1 1 71 108 158 +38 2 2 1 1 105 77 202 +39 2 2 1 1 78 106 201 +40 2 2 1 1 76 108 159 +41 2 2 1 1 107 74 160 +42 2 2 1 1 105 67 161 +43 2 2 1 1 69 106 162 +44 2 2 1 1 106 69 163 +45 2 2 1 1 67 105 164 +46 2 2 1 1 167 69 180 +47 2 2 1 1 166 68 181 +48 2 2 1 1 70 168 182 +49 2 2 1 1 72 152 154 +50 2 2 1 1 149 71 153 +51 2 2 1 1 165 67 179 +52 2 2 1 1 107 157 204 +53 2 2 1 1 157 80 204 +54 2 2 1 1 158 108 203 +55 2 2 1 1 79 158 203 +56 2 2 1 1 111 74 136 +57 2 2 1 1 112 73 135 +58 2 2 1 1 75 110 133 +59 2 2 1 1 76 109 134 +60 2 2 1 1 68 151 156 +61 2 2 1 1 150 70 155 +62 2 2 1 1 83 119 228 +63 2 2 1 1 117 81 227 +64 2 2 1 1 118 82 229 +65 2 2 1 1 86 122 234 +66 2 2 1 1 120 84 230 +67 2 2 1 1 85 121 231 +68 2 2 1 1 87 124 232 +69 2 2 1 1 123 88 233 +70 2 2 1 1 116 77 140 +71 2 2 1 1 78 115 139 +72 2 2 1 1 113 79 138 +73 2 2 1 1 80 114 137 +74 2 2 1 1 127 83 228 +75 2 2 1 1 81 131 227 +76 2 2 1 1 130 86 234 +77 2 2 1 1 82 132 229 +78 2 2 1 1 84 126 230 +79 2 2 1 1 125 85 231 +80 2 2 1 1 129 87 232 +81 2 2 1 1 88 128 233 +82 2 2 1 1 151 68 173 +83 2 2 1 1 70 150 174 +84 2 2 1 1 62 61 113 +85 2 2 1 1 49 50 114 +86 2 2 1 1 47 48 116 +87 2 2 1 1 60 59 115 +88 2 2 1 1 25 26 109 +89 2 2 1 1 23 24 112 +90 2 2 1 1 61 62 110 +91 2 2 1 1 59 60 111 +92 2 2 1 1 68 172 181 +93 2 2 1 1 69 170 180 +94 2 2 1 1 171 70 182 +95 2 2 1 1 67 169 179 +96 2 2 1 1 83 127 179 +97 2 2 1 1 85 125 180 +98 2 2 1 1 86 130 181 +99 2 2 1 1 128 88 182 +100 2 2 1 1 127 103 179 +101 2 2 1 1 125 101 180 +102 2 2 1 1 130 104 181 +103 2 2 1 1 102 128 182 +104 2 2 1 1 50 51 137 +105 2 2 1 1 63 62 138 +106 2 2 1 1 46 47 140 +107 2 2 1 1 59 58 139 +108 2 2 1 1 151 173 187 +109 2 2 1 1 174 150 188 +110 2 2 1 1 26 27 134 +111 2 2 1 1 62 63 133 +112 2 2 1 1 22 23 135 +113 2 2 1 1 58 59 136 +114 2 2 1 1 81 117 213 +115 2 2 1 1 121 85 214 +116 2 2 1 1 122 86 215 +117 2 2 1 1 88 123 216 +118 2 2 1 1 82 118 211 +119 2 2 1 1 84 120 210 +120 2 2 1 1 119 83 209 +121 2 2 1 1 124 87 212 +122 2 2 1 1 154 152 202 +123 2 2 1 1 149 153 201 +124 2 2 1 1 44 45 145 +125 2 2 1 1 57 56 146 +126 2 2 1 1 52 53 148 +127 2 2 1 1 65 64 147 +128 2 2 1 1 56 57 141 +129 2 2 1 1 20 21 142 +130 2 2 1 1 64 65 143 +131 2 2 1 1 28 29 144 +132 2 2 1 1 38 39 127 +133 2 2 1 1 16 17 126 +134 2 2 1 1 14 15 125 +135 2 2 1 1 34 35 128 +136 2 2 1 1 32 33 129 +137 2 2 1 1 8 9 130 +138 2 2 1 1 40 41 131 +139 2 2 1 1 10 11 132 +140 2 2 1 1 145 81 213 +141 2 2 1 1 85 146 214 +142 2 2 1 1 86 148 215 +143 2 2 1 1 147 88 216 +144 2 2 1 1 49 114 186 +145 2 2 1 1 113 61 185 +146 2 2 1 1 116 48 186 +147 2 2 1 1 60 115 185 +148 2 2 1 1 141 82 211 +149 2 2 1 1 142 84 210 +150 2 2 1 1 83 143 209 +151 2 2 1 1 87 144 212 +152 2 2 1 1 25 109 184 +153 2 2 1 1 61 110 183 +154 2 2 1 1 112 24 184 +155 2 2 1 1 111 60 183 +156 2 2 1 1 156 151 204 +157 2 2 1 1 150 155 203 +158 2 2 1 1 156 99 172 +159 2 2 1 1 100 155 171 +160 2 2 1 1 115 91 185 +161 2 2 1 1 92 116 186 +162 2 2 1 1 109 90 184 +163 2 2 1 1 110 89 183 +164 2 2 1 1 90 112 184 +165 2 2 1 1 89 111 183 +166 2 2 1 1 90 109 248 +167 2 2 1 1 111 89 247 +168 2 2 1 1 186 114 253 +169 2 2 1 1 113 185 254 +170 2 2 1 1 68 156 172 +171 2 2 1 1 155 70 171 +172 2 2 1 1 92 186 253 +173 2 2 1 1 185 91 254 +174 2 2 1 1 115 78 246 +175 2 2 1 1 77 116 245 +176 2 2 1 1 136 74 187 +177 2 2 1 1 76 134 188 +178 2 2 1 1 165 81 236 +179 2 2 1 1 167 84 235 +180 2 2 1 1 73 112 178 +181 2 2 1 1 112 90 178 +182 2 2 1 1 89 110 177 +183 2 2 1 1 110 75 177 +184 2 2 1 1 164 93 169 +185 2 2 1 1 162 98 170 +186 2 2 1 1 67 164 169 +187 2 2 1 1 69 162 170 +188 2 2 1 1 95 136 187 +189 2 2 1 1 134 96 188 +190 2 2 1 1 176 165 236 +191 2 2 1 1 166 82 207 +192 2 2 1 1 87 168 208 +193 2 2 1 1 175 167 235 +194 2 2 1 1 52 148 195 +195 2 2 1 1 145 45 193 +196 2 2 1 1 147 64 196 +197 2 2 1 1 57 146 194 +198 2 2 1 1 28 144 192 +199 2 2 1 1 64 143 189 +200 2 2 1 1 142 21 190 +201 2 2 1 1 141 57 191 +202 2 2 1 1 148 99 195 +203 2 2 1 1 97 145 193 +204 2 2 1 1 100 147 196 +205 2 2 1 1 146 98 194 +206 2 2 1 1 55 118 239 +207 2 2 1 1 43 117 238 +208 2 2 1 1 19 120 240 +209 2 2 1 1 121 55 241 +210 2 2 1 1 119 66 237 +211 2 2 1 1 122 54 242 +212 2 2 1 1 124 30 244 +213 2 2 1 1 66 123 243 +214 2 2 1 1 143 93 189 +215 2 2 1 1 144 96 192 +216 2 2 1 1 94 142 190 +217 2 2 1 1 95 141 191 +218 2 2 1 1 173 166 207 +219 2 2 1 1 168 174 208 +220 2 2 1 1 118 55 211 +221 2 2 1 1 120 19 210 +222 2 2 1 1 66 119 209 +223 2 2 1 1 30 124 212 +224 2 2 1 1 55 121 214 +225 2 2 1 1 117 43 213 +226 2 2 1 1 123 66 216 +227 2 2 1 1 54 122 215 +228 2 2 1 1 227 42 238 +229 2 2 1 1 13 231 241 +230 2 2 1 1 230 18 240 +231 2 2 1 1 233 36 243 +232 2 2 1 1 31 232 244 +233 2 2 1 1 37 228 237 +234 2 2 1 1 229 12 239 +235 2 2 1 1 7 234 242 +236 2 2 1 1 114 80 253 +237 2 2 1 1 79 113 254 +238 2 2 1 1 109 76 248 +239 2 2 1 1 74 111 247 +240 2 2 1 1 91 115 246 +241 2 2 1 1 116 92 245 +242 2 2 1 1 40 131 199 +243 2 2 1 1 16 126 197 +244 2 2 1 1 125 15 197 +245 2 2 1 1 127 39 199 +246 2 2 1 1 10 132 200 +247 2 2 1 1 34 128 198 +248 2 2 1 1 129 33 198 +249 2 2 1 1 130 9 200 +250 2 2 1 1 101 125 197 +251 2 2 1 1 126 101 197 +252 2 2 1 1 131 103 199 +253 2 2 1 1 132 104 200 +254 2 2 1 1 128 102 198 +255 2 2 1 1 102 129 198 +256 2 2 1 1 103 127 199 +257 2 2 1 1 104 130 200 +258 2 2 1 1 93 133 189 +259 2 2 1 1 96 134 192 +260 2 2 1 1 135 94 190 +261 2 2 1 1 136 95 191 +262 2 2 1 1 99 137 195 +263 2 2 1 1 138 100 196 +264 2 2 1 1 98 139 194 +265 2 2 1 1 140 97 193 +266 2 2 1 1 161 67 176 +267 2 2 1 1 133 63 189 +268 2 2 1 1 134 27 192 +269 2 2 1 1 22 135 190 +270 2 2 1 1 58 136 191 +271 2 2 1 1 163 69 175 +272 2 2 1 1 137 51 195 +273 2 2 1 1 63 138 196 +274 2 2 1 1 139 58 194 +275 2 2 1 1 46 140 193 +276 2 2 1 1 117 227 238 +277 2 2 1 1 231 121 241 +278 2 2 1 1 120 230 240 +279 2 2 1 1 123 233 243 +280 2 2 1 1 232 124 244 +281 2 2 1 1 118 229 239 +282 2 2 1 1 234 122 242 +283 2 2 1 1 228 119 237 +284 2 2 1 1 67 165 176 +285 2 2 1 1 69 167 175 +286 2 2 1 1 101 167 180 +287 2 2 1 1 168 102 182 +288 2 2 1 1 104 166 181 +289 2 2 1 1 103 165 179 +290 2 2 1 1 135 73 217 +291 2 2 1 1 75 133 218 +292 2 2 1 1 140 77 219 +293 2 2 1 1 78 139 220 +294 2 2 1 1 138 79 225 +295 2 2 1 1 80 137 226 +296 2 2 1 1 143 83 252 +297 2 2 1 1 13 14 231 +298 2 2 1 1 17 18 230 +299 2 2 1 1 11 12 229 +300 2 2 1 1 41 42 227 +301 2 2 1 1 35 36 233 +302 2 2 1 1 31 32 232 +303 2 2 1 1 7 8 234 +304 2 2 1 1 37 38 228 +305 2 2 1 1 146 85 249 +306 2 2 1 1 148 86 251 +307 2 2 1 1 88 147 250 +308 2 2 1 1 71 149 206 +309 2 2 1 1 152 72 205 +310 2 2 1 1 149 73 178 +311 2 2 1 1 75 152 177 +312 2 2 1 1 68 166 173 +313 2 2 1 1 168 70 174 +314 2 2 1 1 55 56 211 +315 2 2 1 1 19 20 210 +316 2 2 1 1 65 66 209 +317 2 2 1 1 29 30 212 +318 2 2 1 1 84 142 235 +319 2 2 1 1 131 81 221 +320 2 2 1 1 132 82 222 +321 2 2 1 1 126 84 223 +322 2 2 1 1 87 129 224 +323 2 2 1 1 56 55 214 +324 2 2 1 1 43 44 213 +325 2 2 1 1 66 65 216 +326 2 2 1 1 53 54 215 +327 2 2 1 1 81 145 236 +328 2 2 1 1 137 99 226 +329 2 2 1 1 100 138 225 +330 2 2 1 1 133 93 218 +331 2 2 1 1 94 135 217 +332 2 2 1 1 14 125 231 +333 2 2 1 1 126 17 230 +334 2 2 1 1 38 127 228 +335 2 2 1 1 128 35 233 +336 2 2 1 1 32 129 232 +337 2 2 1 1 8 130 234 +338 2 2 1 1 132 11 229 +339 2 2 1 1 131 41 227 +340 2 2 1 1 139 98 220 +341 2 2 1 1 97 140 219 +342 2 2 1 1 153 91 246 +343 2 2 1 1 92 154 245 +344 2 2 1 1 97 161 176 +345 2 2 1 1 94 163 175 +346 2 2 1 1 80 157 253 +347 2 2 1 1 158 79 254 +348 2 2 1 1 82 141 207 +349 2 2 1 1 144 87 208 +350 2 2 1 1 151 107 204 +351 2 2 1 1 108 150 203 +352 2 2 1 1 152 105 202 +353 2 2 1 1 106 149 201 +354 2 2 1 1 142 94 235 +355 2 2 1 1 93 143 252 +356 2 2 1 1 141 95 207 +357 2 2 1 1 96 144 208 +358 2 2 1 1 74 151 187 +359 2 2 1 1 150 76 188 +360 2 2 1 1 76 159 248 +361 2 2 1 1 160 74 247 +362 2 2 1 1 145 97 236 +363 2 2 1 1 98 146 249 +364 2 2 1 1 99 148 251 +365 2 2 1 1 147 100 250 +366 2 2 1 1 56 141 211 +367 2 2 1 1 20 142 210 +368 2 2 1 1 143 65 209 +369 2 2 1 1 144 29 212 +370 2 2 1 1 44 145 213 +371 2 2 1 1 146 56 214 +372 2 2 1 1 65 147 216 +373 2 2 1 1 148 53 215 +374 2 2 1 1 73 163 217 +375 2 2 1 1 164 75 218 +376 2 2 1 1 169 83 179 +377 2 2 1 1 170 85 180 +378 2 2 1 1 172 86 181 +379 2 2 1 1 88 171 182 +380 2 2 1 1 162 78 220 +381 2 2 1 1 77 161 219 +382 2 2 1 1 72 160 205 +383 2 2 1 1 159 71 206 +384 2 2 1 1 90 159 206 +385 2 2 1 1 160 89 205 +386 2 2 1 1 101 126 223 +387 2 2 1 1 103 131 221 +388 2 2 1 1 104 132 222 +389 2 2 1 1 129 102 224 +390 2 2 1 1 173 95 187 +391 2 2 1 1 96 174 188 +392 2 2 1 1 60 61 183 +393 2 2 1 1 24 25 184 +394 2 2 1 1 153 78 201 +395 2 2 1 1 77 154 202 +396 2 2 1 1 48 49 186 +397 2 2 1 1 61 60 185 +398 2 2 1 1 78 153 246 +399 2 2 1 1 154 77 245 +400 2 2 1 1 155 79 203 +401 2 2 1 1 80 156 204 +402 2 2 1 1 4 43 238 +403 2 2 1 1 54 1 242 +404 2 2 1 1 66 6 237 +405 2 2 1 1 30 3 244 +406 2 2 1 1 6 66 243 +407 2 2 1 1 5 55 239 +408 2 2 1 1 2 19 240 +409 2 2 1 1 55 5 241 +410 2 2 1 1 79 155 225 +411 2 2 1 1 156 80 226 +412 2 2 1 1 155 100 225 +413 2 2 1 1 99 156 226 +414 2 2 1 1 157 92 253 +415 2 2 1 1 91 158 254 +416 2 2 1 1 159 90 248 +417 2 2 1 1 89 160 247 +418 2 2 1 1 165 103 221 +419 2 2 1 1 166 104 222 +420 2 2 1 1 167 101 223 +421 2 2 1 1 102 168 224 +422 2 2 1 1 63 64 189 +423 2 2 1 1 21 22 190 +424 2 2 1 1 57 58 191 +425 2 2 1 1 27 28 192 +426 2 2 1 1 45 46 193 +427 2 2 1 1 58 57 194 +428 2 2 1 1 51 52 195 +429 2 2 1 1 64 63 196 +430 2 2 1 1 161 97 219 +431 2 2 1 1 98 162 220 +432 2 2 1 1 163 94 217 +433 2 2 1 1 93 164 218 +434 2 2 1 1 81 165 221 +435 2 2 1 1 82 166 222 +436 2 2 1 1 84 167 223 +437 2 2 1 1 168 87 224 +438 2 2 1 1 6 37 237 +439 2 2 1 1 12 5 239 +440 2 2 1 1 42 4 238 +441 2 2 1 1 18 2 240 +442 2 2 1 1 5 13 241 +443 2 2 1 1 1 7 242 +444 2 2 1 1 3 31 244 +445 2 2 1 1 36 6 243 +446 2 2 1 1 169 93 252 +447 2 2 1 1 170 98 249 +448 2 2 1 1 100 171 250 +449 2 2 1 1 172 99 251 +450 2 2 1 1 39 40 199 +451 2 2 1 1 15 16 197 +452 2 2 1 1 9 10 200 +453 2 2 1 1 33 34 198 +454 2 2 1 1 83 169 252 +455 2 2 1 1 85 170 249 +456 2 2 1 1 171 88 250 +457 2 2 1 1 86 172 251 +458 2 2 1 1 149 178 206 +459 2 2 1 1 177 152 205 +460 2 2 1 1 89 177 205 +461 2 2 1 1 178 90 206 +462 2 2 1 1 95 173 207 +463 2 2 1 1 174 96 208 +464 2 2 1 1 94 175 235 +465 2 2 1 1 97 176 236 +$EndElements diff --git a/src/mmdg.cc b/src/mmdg.cc index fc9d19b5e0de7b0157c381b5f60be5fe254321bd..503e4ebce7c4b368e2010fe15c9a98b483af26ef 100644 --- a/src/mmdg.cc +++ b/src/mmdg.cc @@ -101,10 +101,33 @@ int main(int argc, char** argv) MMDG mmdg(gridView, mapper, iGridView, iMapper, *problem); mmdg(mu, xi); - //print total run time + //error and total run time + const double error = mmdg.computeL2error(); const double timeTotal = ClockHelper::elapsedRuntime(timeInitial); + + //determine maximum edge length + double largestEdgeLength = 0.; + for (const auto& edge : edges(gridView)) + { + largestEdgeLength = std::max(largestEdgeLength, edge.geometry().volume()); + } + + //write error, runtime and maximum edge length + std::ofstream errorFile("convergenceData", std::ios::out | std::ios::trunc); + + if (errorFile.is_open()) + { + errorFile << error << "\n" << timeTotal << "\n" << largestEdgeLength; + errorFile.close(); + } + else + { + std::cout << "Unable to write error data to file.\n"; + } + + //print total run time and error std::cout << "Finished with a total run time of " << timeTotal << - " Seconds.\nL2-error: " << mmdg.computeL2error() << ".\n"; + " Seconds.\nL2-error: " << error << ".\n"; return 0; } diff --git a/src/mmdgAnalysis.py b/src/mmdgAnalysis.py new file mode 100644 index 0000000000000000000000000000000000000000..ff746bd0a7b2e3e58779824cbe25e1e07af73f2b --- /dev/null +++ b/src/mmdgAnalysis.py @@ -0,0 +1,114 @@ +import numpy as np +from os import mkdir, remove +from os.path import join, exists +from subprocess import call +from math import log + + +# changes the parameter file such that gridfile is used +def writeINI(gridfile): + # storage for file content + lines = None + + # new line to write + newline = "gridfile = " + gridfile + + # indicates whether the parameter gridfile already exists + hasGridlineParameter = False + + # read file content and replace gridfile parameter if existent + with open("parameterMMDG.ini", "r") as file: + lines = file.read().splitlines() + + for i in range(len(lines)): + if (lines[i].startswith("gridfile")): + lines[i] = newline + hasGridlineParameter = True + break + + # add parameter gridfile if not already existent + if (not hasGridlineParameter): + lines.append(newline) + + # write parameters to file + with open("parameterMMDG.ini", "w") as file: + file.write("\n".join(lines)) + + +# reads error, run time and maximum edge length, i.e. the output from the +# excuted discontinuous Galerkin scheme +def readData(): + file = open("convergenceData", "r") + + # read data + error = float(file.readline()) + timeTotal = float(file.readline()) + numberOfElements = float(file.readline()) + + file.close() + + # return data + return np.array([error, timeTotal, numberOfElements]) + + +# === main === + +# create directory for the plots if necessary +if not exists("plots"): + mkdir("plots") + +# number of runs from which the run times are taken as mean value +numberOfMeans = 1 + +# names of the grid files sorted by refinement +gridfiles = ["mmdg2_1", "mmdg2_5e-1", "mmdg2_2e-1", "mmdg2_8e-2", "mmdg2_5e-2"] + +dim = 2 + +# storage for errors, run times and largestEdgeLength +# column 0: error +# column 1: total run time +# column 2: largest edge length +# different rows correspond to different numbers of grid elements +data = np.zeros((len(gridfiles), 3)) + +for i in range(len(gridfiles)): + # write dgf file + writeINI(gridfiles[i]) + + # take mean value over numberOfMeans runs + for j in range(numberOfMeans): + # run executable + call(join(".", "mmdg-" + str(dim) + "d")) + + # store data + data[i, :] += readData() + # end for j +# end for i + +# divide by numberOfMeans +data /= numberOfMeans + +# determine experimental order of convergence (EOC) +EOC = np.empty(len(gridfiles) - 1) + +for i in range(len(gridfiles) - 1): + if data[i, 0] != 0: + EOC[i] = log(data[i, 0] / data[i+1, 0]) / \ + log(data[i, 2] / data[i+1, 2]) + else: + EOC[i] = -1.0*float('inf') +# end for i + +print("\n\n-------------------------------------------------------------\n") +# print EOCs +print("EOCs (" + str(dim) + "d):\t") +print(EOC) +print("\n\n") + +# write data to file +np.savetxt(join("plots", "mmdgAnalysis_" + str(dim) + "d"), data) +np.savetxt(join("plots", "mmdgEOC_" + str(dim) + "d"), EOC) + +# remove output files +remove("convergenceData") diff --git a/src/parameterMMDG.ini b/src/parameterMMDG.ini index 4abb5eff094bb8d4b40b7563301f02ff7a9b0f1b..fb6290810e44ee0860375165a8e320bef7eea5d2 100644 --- a/src/parameterMMDG.ini +++ b/src/parameterMMDG.ini @@ -1,3 +1,4 @@ mu = 1000 xi = 0.75 problem = mmdg1 +gridfile = mmdg2_5e-2 \ No newline at end of file