diff --git a/test/blob2d-order1-h4e-2.msh b/test/blob2d-order1-h4e-2.msh new file mode 100644 index 0000000000000000000000000000000000000000..23bde7e921ecc3fcdbe19679f95f04eea6b606ca --- /dev/null +++ b/test/blob2d-order1-h4e-2.msh @@ -0,0 +1,659 @@ +$MeshFormat +2.2 0 8 +$EndMeshFormat +$Nodes +217 +1 -0.4725430905820004 0.1646517366169588 0 +2 -0.4762661263952118 0.1250877152932441 0 +3 -0.4773599478381062 0.0853681458554567 0 +4 -0.4748733565534289 0.04572043802992175 0 +5 -0.467300319710205 0.006745496374242575 0 +6 -0.4523421022370778 -0.02997415663449131 0 +7 -0.4274865474217647 -0.06073263012049521 0 +8 -0.3931711449409717 -0.0803035602652334 0 +9 -0.3541954021335921 -0.08739963522756089 0 +10 -0.3144827270400216 -0.08674056357209309 0 +11 -0.2749104845112884 -0.0830682885299959 0 +12 -0.2353497013168931 -0.07926456065881671 0 +13 -0.1956604144229968 -0.07733118395197389 0 +14 -0.155966711578534 -0.07891235244169709 0 +15 -0.1167806337081991 -0.08536004992487563 0 +16 -0.07870433262703255 -0.09668867888745789 0 +17 -0.04169550447504759 -0.1111600376162981 0 +18 -0.005394997590983816 -0.1273390031579322 0 +19 0.03068646220171178 -0.1440042339504319 0 +20 0.06720514627158052 -0.1596756388422167 0 +21 0.1051153372298121 -0.1714461165998769 0 +22 0.1442704922483436 -0.1699860411042714 0 +23 0.1708662293012103 -0.1417406133594973 0 +24 0.1814413120367157 -0.1035684055214138 0 +25 0.1840528408594294 -0.06394601525820293 0 +26 0.1818595180123657 -0.02427921318389751 0 +27 0.176079759682738 0.01503105664964789 0 +28 0.1671905436836027 0.05375823464281488 0 +29 0.1552733991155098 0.0916622810283626 0 +30 0.1401003964773046 0.1283802621586579 0 +31 0.1210947092852233 0.1632575488526254 0 +32 0.09719691393341853 0.1949492174977272 0 +33 0.06689182935603925 0.22046921999074 0 +34 0.03077656443659951 0.2368039039685216 0 +35 -0.0079326157896339 0.2456641997307331 0 +36 -0.04731354998143419 0.2510024761943727 0 +37 -0.08672919137787519 0.2560933870385579 0 +38 -0.1258389919785983 0.2631224885473071 0 +39 -0.1644749323606459 0.2724307460902149 0 +40 -0.2028773114587941 0.2826705652598434 0 +41 -0.2413464056950951 0.292654208924524 0 +42 -0.2801629829280154 0.3011670139090882 0 +43 -0.3195170438536954 0.3065894611756452 0 +44 -0.3592029068447401 0.3063536016624879 0 +45 -0.3975412661328094 0.2964955501668801 0 +46 -0.4298890146845241 0.2738428388971937 0 +47 -0.452149762197038 0.2411179486306202 0 +48 -0.4654203077373789 0.203725528435598 0 +49 -0.283131791830473 0.1087235266165297 0 +50 -0.02839939422862658 0.07140366512189922 0 +51 -0.151735664363485 0.1314140059362366 0 +52 0.05655933132983804 -0.01946510083338654 0 +53 -0.3470636551304878 0.197019004755758 0 +54 -0.2051005191441966 0.03263461012333729 0 +55 -0.3563436913434244 0.02294108079688978 0 +56 -0.1121571577412265 0.01431095053540177 0 +57 -0.2384564299642795 0.1914091686131413 0 +58 -0.05984938805875462 0.1504387161188163 0 +59 0.06316507564539464 0.07365027818586706 0 +60 -0.3869380201876764 0.1096564244386179 0 +61 -0.02638064310117438 -0.02463176520181481 0 +62 0.02395518807678415 0.1552318657628732 0 +63 -0.2809915859417534 0.0006398812766845174 0 +64 0.100906320596154 -0.08852032556538908 0 +65 -0.2885003909752318 0.2311010134564113 0 +66 -0.1600429100264193 0.2003272736756557 0 +67 -0.09620089592479175 0.0818486298422228 0 +68 -0.2141127352736424 0.09357637153817633 0 +69 0.1081787751645238 0.02273100140656015 0 +70 0.04599623070750147 -0.08034882337394095 0 +71 -0.4090895905975009 0.1695297890348602 0 +72 -0.4142723799902692 0.04967448893706428 0 +73 -0.1629536830801517 -0.01729366921131759 0 +74 0.01558072583762482 0.03218137936855958 0 +75 -0.2964244213754144 0.1736203291923601 0 +76 -0.157700881930831 0.06443607283687824 0 +77 -0.2592200366139046 0.05454278586311223 0 +78 -0.3305449028939836 -0.03004170946876822 0 +79 -0.3367911021531139 0.07792976185740536 0 +80 -0.2244017953170984 -0.02218591142233278 0 +81 0.1266383405526554 -0.03690388273696893 0 +82 -0.1051706534055906 0.1973367970054319 0 +83 -0.340908838683926 0.252381795956068 0 +84 0.07584917044075079 0.1287383674485626 0 +85 -0.08508645748375987 -0.03812688613477899 0 +86 -0.01555755782618679 0.1896396709319706 0 +87 -0.3956368108035659 0.2266984763870937 0 +88 -0.3993502917439344 -0.01257408304061092 0 +89 -0.05281099781371601 0.01689997609577099 0 +90 -0.3349406093364702 0.1342262246720746 0 +91 -0.2392482962867322 0.2504618354595648 0 +92 -0.03053813900737234 0.1220310625162465 0 +93 0.02344264906675819 0.08824154092149872 0 +94 -0.1894111155648078 0.1453925890103812 0 +95 -0.258427141331394 0.1320968035124849 0 +96 -0.01123016340549049 -0.07127328113766834 0 +97 -0.3076037086651191 0.03544789526366843 0 +98 -0.1040163047466756 0.1304754153685062 0 +99 -0.4329286868645948 0.1040194168872023 0 +100 0.1094420290950928 0.07318464516789902 0 +101 0.1264891722530739 -0.1265993368540926 0 +102 -0.1886341923184308 0.2323564930415953 0 +103 0.03207949271685234 0.1915926509427627 0 +104 -0.05521311059825917 0.2122124554902238 0 +105 0.01395870358797302 -0.03397780012163154 0 +106 0.06724155442814206 -0.1147575718078307 0 +107 0.06289328399286487 0.02733737800192848 0 +108 -0.126964017775521 -0.04452941207594663 0 +109 -0.1511198661007538 0.01978735835094075 0 +110 -0.1310653884299348 0.1676884636616916 0 +111 -0.2627650775744838 -0.0417517812155524 0 +112 -0.06618856343121793 0.1071821956711284 0 +113 0.1430213068477506 -0.08005961766992309 0 +114 -0.3696021480816857 -0.04506372985463836 0 +115 -0.3694101369137729 0.1518170412197574 0 +116 -0.05100206464336021 -0.06709826621028275 0 +117 0.07252434527135243 -0.05764049432709661 0 +118 0.06151435952557038 0.1662525388166215 0 +119 -0.2428869357254709 0.01739802249928257 0 +120 -0.306373827643624 0.2674034534056846 0 +121 -0.4279422495728668 0.01601876977174471 0 +122 -0.3730357452811778 0.06513547388541019 0 +123 -0.439462148508511 0.1445812575817263 0 +124 -0.01676477891108569 0.01258219889845025 0 +125 -0.122360923310017 0.05220975115837262 0 +126 -0.3810929422026038 0.2639177563028782 0 +127 -0.2708885073658144 0.1969324122857966 0 +128 -0.05758994484851993 0.06496099024058383 0 +129 -0.4247168108523073 0.2080652479324283 0 +130 -0.1795720388483829 0.1037201982774379 0 +131 0.1394636811589191 -0.005188381019605448 0 +132 -0.1243666951778905 0.2255354550176586 0 +133 -0.3031627190525917 -0.05419727567397337 0 +134 0.1084843101897359 0.1112760665121997 0 +135 -0.3102521030984001 0.2003083022489857 0 +136 -0.2693471878464848 0.2611282798242717 0 +137 0.09566516590775487 -0.01781769087268814 0 +138 -0.2241016445718222 0.06507828843913877 0 +139 -0.2510594542295104 0.09293247606623428 0 +140 -0.1326326237226652 0.09735093119210697 0 +141 -0.296644290952609 0.06885958769946859 0 +142 -0.193513089863136 -0.03637987962585265 0 +143 -0.1958746684515904 0.186376933880656 0 +144 -0.4437609320493616 0.0677758771303823 0 +145 -0.3913889667186382 0.0225276075015296 0 +146 -0.3249359242683438 0.006172068462894031 0 +147 -0.08045196232044804 -0.002744496137503105 0 +148 -0.3267988519083628 0.1655251150933931 0 +149 -0.1656259493672474 0.1662398532329242 0 +150 0.04608126662182122 0.118935126183763 0 +151 -0.1950140634179905 0.001675369537059934 0 +152 0.09347189477240914 -0.1346674636885036 0 +153 -0.3854805080499073 0.1960982466360818 0 +154 -0.359618974663202 -0.01142621294596191 0 +155 -0.09648805801448768 0.1664668849502332 0 +156 0.0009878483777253644 0.064302521998855 0 +157 -0.01141902210062656 0.1554277831617881 0 +158 0.1343672406392891 0.04455413160320354 0 +159 0.02410275875804385 -0.0007919850213387036 0 +160 -0.1928733591024794 0.06617982189798599 0 +161 -0.3614687926351708 0.2280042756094886 0 +162 0.1557441729285529 -0.04561726578018516 0 +163 0.01832872536601649 -0.1000299397565083 0 +164 -0.4075657166144112 0.1359207858324329 0 +165 0.004019370419764801 0.211368537863055 0 +166 -0.2976861648532819 0.1407106768733521 0 +167 -0.1313488370318916 -0.008752696458298126 0 +168 -0.4081645667904638 0.08223954000343248 0 +169 -0.027648982594903 0.2185979500896451 0 +170 0.08718909523129281 0.05012446709599111 0 +171 -0.3175199124632994 0.1055547209600123 0 +172 -0.08692864699311986 0.03791430028914144 0 +173 -0.3345004020628536 -0.06102316736083863 0 +174 -0.3531059210559184 0.1073866078388797 0 +175 0.03651295551895878 0.0570276192232001 0 +176 -0.1546324507311705 0.2393093088376199 0 +177 -0.4422463896555397 0.1781107119203143 0 +178 0.04193459853285303 -0.04904381497245158 0 +179 -0.1628348170724352 -0.04849270426675628 0 +180 -0.22241261148305 0.2232634167635006 0 +181 -0.08272767168861184 0.2240191219031123 0 +182 -0.2783739167048174 0.03157729585418878 0 +183 -0.4017297805469668 -0.04348235281031765 0 +184 0.1033984038841308 -0.05559435660024836 0 +185 -0.3400218831942955 0.2816878880988483 0 +186 -0.02487931034065821 0.0428559587586631 0 +187 -0.09107135410490336 -0.06860224081273397 0 +188 -0.335191174049523 0.04670999777378404 0 +189 -0.2223380156989216 -0.0513826633749056 0 +190 -0.2655576580603503 0.1662459123908351 0 +191 0.1544864872571956 -0.112295226022302 0 +192 -0.3204863282747889 0.231249000047193 0 +193 0.006579076210037116 0.1226894869667563 0 +194 -0.3004800419462312 -0.02383576332374304 0 +195 0.09832712775733066 0.1502878747897323 0 +196 -0.4183342993592698 0.2433298618087805 0 +197 0.08060437039855906 0.1011568966996583 0 +198 -0.01188089962277111 0.09575903105093664 0 +199 -0.04363883084956099 0.1808902524118093 0 +200 -0.4285150932232066 -0.01314601379098435 0 +201 -0.2206092787351329 0.1272876404597573 0 +202 -0.2526947661092331 -0.01231936604174039 0 +203 -0.446680732028625 0.03833876942034325 0 +204 -0.1054128410785217 -0.01600289262712828 0 +205 -0.2280560486845925 0.1581348413112093 0 +206 -0.211810858327057 0.2549939157956889 0 +207 -0.01580691290897474 -0.09879631300501679 0 +208 0.07245910723005045 -0.08590299335515234 0 +209 -0.1785448874445977 0.03728096995677022 0 +210 -0.1309247474375324 0.1965380710306298 0 +211 -0.05914642507249168 -0.02314028751772174 0 +212 -0.2564209952907257 0.2252206304166123 0 +213 0.02179761895777071 -0.0669347318724402 0 +214 -0.03937591083020087 0.09367152920489702 0 +215 0.007730177925584392 0.1780017352959937 0 +216 -0.07384795210254413 0.1885607046466045 0 +217 -0.002222742105310149 -0.01202024886506414 0 +$EndNodes +$Elements +433 +1 15 2 0 1 1 +2 1 2 0 1 1 2 +3 1 2 0 1 2 3 +4 1 2 0 1 3 4 +5 1 2 0 1 4 5 +6 1 2 0 1 5 6 +7 1 2 0 1 6 7 +8 1 2 0 1 7 8 +9 1 2 0 1 8 9 +10 1 2 0 1 9 10 +11 1 2 0 1 10 11 +12 1 2 0 1 11 12 +13 1 2 0 1 12 13 +14 1 2 0 1 13 14 +15 1 2 0 1 14 15 +16 1 2 0 1 15 16 +17 1 2 0 1 16 17 +18 1 2 0 1 17 18 +19 1 2 0 1 18 19 +20 1 2 0 1 19 20 +21 1 2 0 1 20 21 +22 1 2 0 1 21 22 +23 1 2 0 1 22 23 +24 1 2 0 1 23 24 +25 1 2 0 1 24 25 +26 1 2 0 1 25 26 +27 1 2 0 1 26 27 +28 1 2 0 1 27 28 +29 1 2 0 1 28 29 +30 1 2 0 1 29 30 +31 1 2 0 1 30 31 +32 1 2 0 1 31 32 +33 1 2 0 1 32 33 +34 1 2 0 1 33 34 +35 1 2 0 1 34 35 +36 1 2 0 1 35 36 +37 1 2 0 1 36 37 +38 1 2 0 1 37 38 +39 1 2 0 1 38 39 +40 1 2 0 1 39 40 +41 1 2 0 1 40 41 +42 1 2 0 1 41 42 +43 1 2 0 1 42 43 +44 1 2 0 1 43 44 +45 1 2 0 1 44 45 +46 1 2 0 1 45 46 +47 1 2 0 1 46 47 +48 1 2 0 1 47 48 +49 1 2 0 1 48 1 +50 2 2 0 1 118 150 62 +51 2 2 0 1 115 174 60 +52 2 2 0 1 134 195 30 +53 2 2 0 1 49 166 95 +54 2 2 0 1 92 112 58 +55 2 2 0 1 51 130 94 +56 2 2 0 1 30 195 31 +57 2 2 0 1 84 150 118 +58 2 2 0 1 90 174 115 +59 2 2 0 1 58 157 92 +60 2 2 0 1 61 211 89 +61 2 2 0 1 37 181 132 +62 2 2 0 1 166 190 95 +63 2 2 0 1 89 124 61 +64 2 2 0 1 130 201 94 +65 2 2 0 1 93 150 59 +66 2 2 0 1 95 139 49 +67 2 2 0 1 61 105 96 +68 2 2 0 1 33 118 103 +69 2 2 0 1 37 132 38 +70 2 2 0 1 60 174 122 +71 2 2 0 1 64 113 101 +72 2 2 0 1 58 112 98 +73 2 2 0 1 21 101 22 +74 2 2 0 1 96 213 163 +75 2 2 0 1 163 207 96 +76 2 2 0 1 73 151 109 +77 2 2 0 1 113 184 81 +78 2 2 0 1 96 116 61 +79 2 2 0 1 100 134 29 +80 2 2 0 1 51 110 98 +81 2 2 0 1 101 152 64 +82 2 2 0 1 151 209 109 +83 2 2 0 1 29 158 100 +84 2 2 0 1 68 201 130 +85 2 2 0 1 126 196 46 +86 2 2 0 1 45 126 46 +87 2 2 0 1 107 137 52 +88 2 2 0 1 41 206 91 +89 2 2 0 1 2 99 3 +90 2 2 0 1 32 118 33 +91 2 2 0 1 26 162 131 +92 2 2 0 1 81 162 113 +93 2 2 0 1 122 168 60 +94 2 2 0 1 22 101 23 +95 2 2 0 1 95 201 139 +96 2 2 0 1 91 136 41 +97 2 2 0 1 64 152 106 +98 2 2 0 1 116 207 17 +99 2 2 0 1 103 118 62 +100 2 2 0 1 131 158 27 +101 2 2 0 1 98 112 67 +102 2 2 0 1 63 146 97 +103 2 2 0 1 69 137 107 +104 2 2 0 1 77 138 119 +105 2 2 0 1 119 138 54 +106 2 2 0 1 113 191 101 +107 2 2 0 1 112 128 67 +108 2 2 0 1 66 143 102 +109 2 2 0 1 96 207 116 +110 2 2 0 1 28 158 29 +111 2 2 0 1 58 199 157 +112 2 2 0 1 89 211 147 +113 2 2 0 1 59 175 93 +114 2 2 0 1 52 137 117 +115 2 2 0 1 119 151 80 +116 2 2 0 1 67 140 98 +117 2 2 0 1 98 140 51 +118 2 2 0 1 52 159 107 +119 2 2 0 1 27 158 28 +120 2 2 0 1 102 206 39 +121 2 2 0 1 60 164 115 +122 2 2 0 1 137 184 117 +123 2 2 0 1 97 182 63 +124 2 2 0 1 19 106 20 +125 2 2 0 1 29 134 30 +126 2 2 0 1 139 201 68 +127 2 2 0 1 150 197 59 +128 2 2 0 1 163 213 70 +129 2 2 0 1 141 188 79 +130 2 2 0 1 54 151 119 +131 2 2 0 1 94 149 51 +132 2 2 0 1 97 188 141 +133 2 2 0 1 161 192 53 +134 2 2 0 1 39 206 40 +135 2 2 0 1 2 123 99 +136 2 2 0 1 33 103 34 +137 2 2 0 1 16 116 17 +138 2 2 0 1 26 131 27 +139 2 2 0 1 21 152 101 +140 2 2 0 1 112 214 128 +141 2 2 0 1 11 111 12 +142 2 2 0 1 44 126 45 +143 2 2 0 1 131 162 81 +144 2 2 0 1 8 114 9 +145 2 2 0 1 76 140 125 +146 2 2 0 1 107 159 74 +147 2 2 0 1 14 108 15 +148 2 2 0 1 44 185 126 +149 2 2 0 1 122 174 79 +150 2 2 0 1 81 137 131 +151 2 2 0 1 131 137 69 +152 2 2 0 1 74 186 156 +153 2 2 0 1 125 140 67 +154 2 2 0 1 64 184 113 +155 2 2 0 1 99 164 60 +156 2 2 0 1 124 186 74 +157 2 2 0 1 99 144 3 +158 2 2 0 1 143 180 102 +159 2 2 0 1 130 160 68 +160 2 2 0 1 59 197 100 +161 2 2 0 1 145 154 88 +162 2 2 0 1 19 163 106 +163 2 2 0 1 98 155 58 +164 2 2 0 1 57 212 180 +165 2 2 0 1 24 113 25 +166 2 2 0 1 53 192 135 +167 2 2 0 1 60 168 99 +168 2 2 0 1 132 181 82 +169 2 2 0 1 84 195 134 +170 2 2 0 1 7 200 183 +171 2 2 0 1 75 190 166 +172 2 2 0 1 100 170 59 +173 2 2 0 1 76 125 109 +174 2 2 0 1 108 187 15 +175 2 2 0 1 85 187 108 +176 2 2 0 1 111 189 12 +177 2 2 0 1 47 129 48 +178 2 2 0 1 40 206 41 +179 2 2 0 1 6 200 7 +180 2 2 0 1 55 154 145 +181 2 2 0 1 106 163 70 +182 2 2 0 1 127 135 65 +183 2 2 0 1 75 135 127 +184 2 2 0 1 104 181 36 +185 2 2 0 1 101 191 23 +186 2 2 0 1 69 170 158 +187 2 2 0 1 158 170 100 +188 2 2 0 1 74 159 124 +189 2 2 0 1 109 125 56 +190 2 2 0 1 87 196 126 +191 2 2 0 1 130 140 76 +192 2 2 0 1 51 140 130 +193 2 2 0 1 114 173 9 +194 2 2 0 1 68 160 138 +195 2 2 0 1 77 141 139 +196 2 2 0 1 139 141 49 +197 2 2 0 1 102 176 66 +198 2 2 0 1 42 120 43 +199 2 2 0 1 63 194 146 +200 2 2 0 1 61 217 105 +201 2 2 0 1 159 217 124 +202 2 2 0 1 147 172 89 +203 2 2 0 1 47 196 129 +204 2 2 0 1 127 212 57 +205 2 2 0 1 11 133 111 +206 2 2 0 1 1 123 2 +207 2 2 0 1 115 164 71 +208 2 2 0 1 79 171 141 +209 2 2 0 1 13 179 14 +210 2 2 0 1 80 189 111 +211 2 2 0 1 105 213 96 +212 2 2 0 1 111 194 63 +213 2 2 0 1 39 176 102 +214 2 2 0 1 110 155 98 +215 2 2 0 1 107 175 59 +216 2 2 0 1 74 175 107 +217 2 2 0 1 106 208 64 +218 2 2 0 1 62 215 103 +219 2 2 0 1 120 136 65 +220 2 2 0 1 115 153 53 +221 2 2 0 1 56 167 109 +222 2 2 0 1 109 167 73 +223 2 2 0 1 71 153 115 +224 2 2 0 1 36 169 104 +225 2 2 0 1 142 179 13 +226 2 2 0 1 116 211 61 +227 2 2 0 1 55 145 122 +228 2 2 0 1 122 145 72 +229 2 2 0 1 67 172 125 +230 2 2 0 1 103 165 34 +231 2 2 0 1 126 185 83 +232 2 2 0 1 128 172 67 +233 2 2 0 1 106 152 20 +234 2 2 0 1 88 154 114 +235 2 2 0 1 83 192 161 +236 2 2 0 1 56 172 147 +237 2 2 0 1 36 181 37 +238 2 2 0 1 81 184 137 +239 2 2 0 1 51 149 110 +240 2 2 0 1 42 136 120 +241 2 2 0 1 110 149 66 +242 2 2 0 1 105 159 52 +243 2 2 0 1 10 133 11 +244 2 2 0 1 107 170 69 +245 2 2 0 1 121 200 5 +246 2 2 0 1 99 168 144 +247 2 2 0 1 69 158 131 +248 2 2 0 1 78 173 114 +249 2 2 0 1 53 148 115 +250 2 2 0 1 115 148 90 +251 2 2 0 1 14 179 108 +252 2 2 0 1 52 178 105 +253 2 2 0 1 72 168 122 +254 2 2 0 1 41 136 42 +255 2 2 0 1 70 208 106 +256 2 2 0 1 59 170 107 +257 2 2 0 1 5 200 6 +258 2 2 0 1 128 214 50 +259 2 2 0 1 123 164 99 +260 2 2 0 1 87 153 129 +261 2 2 0 1 141 182 97 +262 2 2 0 1 129 153 71 +263 2 2 0 1 77 182 141 +264 2 2 0 1 108 204 85 +265 2 2 0 1 113 162 25 +266 2 2 0 1 180 212 91 +267 2 2 0 1 82 155 110 +268 2 2 0 1 73 167 108 +269 2 2 0 1 133 194 111 +270 2 2 0 1 144 168 72 +271 2 2 0 1 108 179 73 +272 2 2 0 1 100 197 134 +273 2 2 0 1 55 188 146 +274 2 2 0 1 109 209 76 +275 2 2 0 1 146 188 97 +276 2 2 0 1 32 195 118 +277 2 2 0 1 8 183 114 +278 2 2 0 1 110 210 82 +279 2 2 0 1 50 186 128 +280 2 2 0 1 128 186 89 +281 2 2 0 1 66 210 110 +282 2 2 0 1 114 154 78 +283 2 2 0 1 120 185 43 +284 2 2 0 1 57 205 190 +285 2 2 0 1 38 176 39 +286 2 2 0 1 77 139 138 +287 2 2 0 1 138 139 68 +288 2 2 0 1 116 187 85 +289 2 2 0 1 16 187 116 +290 2 2 0 1 54 209 151 +291 2 2 0 1 31 195 32 +292 2 2 0 1 121 145 88 +293 2 2 0 1 80 202 119 +294 2 2 0 1 141 171 49 +295 2 2 0 1 72 145 121 +296 2 2 0 1 148 166 90 +297 2 2 0 1 75 166 148 +298 2 2 0 1 57 190 127 +299 2 2 0 1 3 144 4 +300 2 2 0 1 136 212 65 +301 2 2 0 1 76 160 130 +302 2 2 0 1 129 196 87 +303 2 2 0 1 135 148 53 +304 2 2 0 1 75 148 135 +305 2 2 0 1 63 202 111 +306 2 2 0 1 111 202 80 +307 2 2 0 1 85 211 116 +308 2 2 0 1 83 185 120 +309 2 2 0 1 120 192 83 +310 2 2 0 1 65 192 120 +311 2 2 0 1 132 176 38 +312 2 2 0 1 24 191 113 +313 2 2 0 1 92 214 112 +314 2 2 0 1 66 176 132 +315 2 2 0 1 142 151 73 +316 2 2 0 1 80 151 142 +317 2 2 0 1 20 152 21 +318 2 2 0 1 63 182 119 +319 2 2 0 1 119 182 77 +320 2 2 0 1 126 161 87 +321 2 2 0 1 83 161 126 +322 2 2 0 1 114 183 88 +323 2 2 0 1 84 197 150 +324 2 2 0 1 89 172 128 +325 2 2 0 1 93 175 156 +326 2 2 0 1 7 183 8 +327 2 2 0 1 35 169 36 +328 2 2 0 1 94 205 143 +329 2 2 0 1 143 205 57 +330 2 2 0 1 5 203 121 +331 2 2 0 1 121 203 72 +332 2 2 0 1 48 177 1 +333 2 2 0 1 86 199 169 +334 2 2 0 1 91 212 136 +335 2 2 0 1 169 199 104 +336 2 2 0 1 117 208 70 +337 2 2 0 1 117 178 52 +338 2 2 0 1 122 188 55 +339 2 2 0 1 66 149 143 +340 2 2 0 1 79 188 122 +341 2 2 0 1 23 191 24 +342 2 2 0 1 181 216 82 +343 2 2 0 1 146 154 55 +344 2 2 0 1 78 154 146 +345 2 2 0 1 4 203 5 +346 2 2 0 1 70 178 117 +347 2 2 0 1 87 161 153 +348 2 2 0 1 153 161 53 +349 2 2 0 1 157 199 86 +350 2 2 0 1 25 162 26 +351 2 2 0 1 18 163 19 +352 2 2 0 1 117 184 64 +353 2 2 0 1 104 216 181 +354 2 2 0 1 129 177 48 +355 2 2 0 1 71 164 123 +356 2 2 0 1 34 165 35 +357 2 2 0 1 133 173 78 +358 2 2 0 1 10 173 133 +359 2 2 0 1 64 208 117 +360 2 2 0 1 143 149 94 +361 2 2 0 1 57 180 143 +362 2 2 0 1 156 186 50 +363 2 2 0 1 103 215 165 +364 2 2 0 1 123 177 71 +365 2 2 0 1 12 189 13 +366 2 2 0 1 119 202 63 +367 2 2 0 1 118 195 84 +368 2 2 0 1 18 207 163 +369 2 2 0 1 9 173 10 +370 2 2 0 1 1 177 123 +371 2 2 0 1 183 200 88 +372 2 2 0 1 165 215 86 +373 2 2 0 1 138 160 54 +374 2 2 0 1 190 205 95 +375 2 2 0 1 125 172 56 +376 2 2 0 1 88 200 121 +377 2 2 0 1 71 177 129 +378 2 2 0 1 72 203 144 +379 2 2 0 1 43 185 44 +380 2 2 0 1 15 187 16 +381 2 2 0 1 89 186 124 +382 2 2 0 1 50 198 156 +383 2 2 0 1 156 198 93 +384 2 2 0 1 46 196 47 +385 2 2 0 1 62 193 157 +386 2 2 0 1 157 193 92 +387 2 2 0 1 13 189 142 +388 2 2 0 1 17 207 18 +389 2 2 0 1 127 190 75 +390 2 2 0 1 124 217 61 +391 2 2 0 1 146 194 78 +392 2 2 0 1 180 206 102 +393 2 2 0 1 91 206 180 +394 2 2 0 1 150 193 62 +395 2 2 0 1 93 193 150 +396 2 2 0 1 65 212 127 +397 2 2 0 1 73 179 142 +398 2 2 0 1 76 209 160 +399 2 2 0 1 166 171 90 +400 2 2 0 1 78 194 133 +401 2 2 0 1 49 171 166 +402 2 2 0 1 135 192 65 +403 2 2 0 1 134 197 84 +404 2 2 0 1 167 204 108 +405 2 2 0 1 132 210 66 +406 2 2 0 1 82 210 132 +407 2 2 0 1 142 189 80 +408 2 2 0 1 79 174 171 +409 2 2 0 1 171 174 90 +410 2 2 0 1 165 169 35 +411 2 2 0 1 86 169 165 +412 2 2 0 1 156 175 74 +413 2 2 0 1 144 203 4 +414 2 2 0 1 105 217 159 +415 2 2 0 1 82 216 155 +416 2 2 0 1 147 204 56 +417 2 2 0 1 85 204 147 +418 2 2 0 1 147 211 85 +419 2 2 0 1 157 215 62 +420 2 2 0 1 155 216 58 +421 2 2 0 1 86 215 157 +422 2 2 0 1 160 209 54 +423 2 2 0 1 56 204 167 +424 2 2 0 1 193 198 92 +425 2 2 0 1 93 198 193 +426 2 2 0 1 178 213 105 +427 2 2 0 1 70 213 178 +428 2 2 0 1 198 214 92 +429 2 2 0 1 199 216 104 +430 2 2 0 1 58 216 199 +431 2 2 0 1 95 205 201 +432 2 2 0 1 201 205 94 +433 2 2 0 1 50 214 198 +$EndElements diff --git a/test/blob2d-order1-h6e-2.msh b/test/blob2d-order1-h6e-2.msh new file mode 100644 index 0000000000000000000000000000000000000000..52441bc224b76f18a77a8690ba919f575acb3d64 --- /dev/null +++ b/test/blob2d-order1-h6e-2.msh @@ -0,0 +1,302 @@ +$MeshFormat +2.2 0 8 +$EndMeshFormat +$Nodes +98 +1 -0.4725430905820004 0.1646517366169588 0 +2 -0.4771912568175954 0.105238158239418 0 +3 -0.4748733565263847 0.0457204378037239 0 +4 -0.46091775595722 -0.01206375870819189 0 +5 -0.4274865470764456 -0.06073263041779322 0 +6 -0.3739345359090098 -0.08520373639895074 0 +7 -0.3144827265548338 -0.08674056353825642 0 +8 -0.2551415644056517 -0.08104586276635461 0 +9 -0.1956604111535049 -0.07733118391659669 0 +10 -0.1362623741781192 -0.08146152721542142 0 +11 -0.07870432944187021 -0.09668868000667932 0 +12 -0.02348447239394116 -0.1191129146124495 0 +13 0.0306864654811038 -0.1440042354437775 0 +14 0.08591464382523883 -0.1663634011715317 0 +15 0.1442704952911986 -0.1699860398610476 0 +16 0.1774994439552023 -0.1230346364318028 0 +17 0.1840528408528492 -0.06394601234996893 0 +18 0.1793767008813017 -0.004564314530393861 0 +19 0.1671905428638292 0.05375823766942801 0 +20 0.1481205846334989 0.1102006893092006 0 +21 0.121094707488876 0.1632575516962436 0 +22 0.08292490390047319 0.2087582151136879 0 +23 0.03077656225280126 0.236803904658542 0 +24 -0.02759450399315626 0.2485413612966029 0 +25 -0.08672919351570241 0.2560933873527119 0 +26 -0.1452047948085679 0.267577388703267 0 +27 -0.2028773133348656 0.2826705657632713 0 +28 -0.2606956396597746 0.2971809588711993 0 +29 -0.3195170527135355 0.3065894618766041 0 +30 -0.3787582918960273 0.3029095267506502 0 +31 -0.4298890155459034 0.2738428380008051 0 +32 -0.4597214774834101 0.2227562121673124 0 +33 -0.285951277396536 0.1025804175065633 0 +34 -0.03566195777522857 0.07603087041697076 0 +35 -0.1534282128340498 0.1328789925358654 0 +36 0.05949382795053582 -0.01331339875324145 0 +37 -0.2068101715503699 0.03118820980709802 0 +38 -0.3503400972958695 0.196378982307797 0 +39 -0.3609209447920571 0.02515429249668077 0 +40 -0.2322764890432984 0.1914258886900886 0 +41 -0.1083325253627921 0.01178473395175389 0 +42 -0.06381055113230943 0.1616613026345616 0 +43 0.0613343502421656 0.07435392698956694 0 +44 -0.391906607520493 0.1083274349043393 0 +45 0.01736429545361701 0.1466086210621322 0 +46 -0.03018292959165764 -0.01915496433045011 0 +47 0.1055129788042974 -0.0851443775833976 0 +48 -0.2779539212651522 -0.00164895976424223 0 +49 -0.2902784466659754 0.2322059943599127 0 +50 -0.1039063241015989 0.07790445943622171 0 +51 -0.1230107110696168 0.1965360201666355 0 +52 -0.2162816393385686 0.1022291168391369 0 +53 0.1064059580073422 0.02367588665400916 0 +54 -0.164782789880846 -0.02204630445542383 0 +55 0.02037955864676078 -0.09242368002757761 0 +56 -0.4061332215724401 0.168052911787027 0 +57 0.01207112800027765 0.02907437843154281 0 +58 -0.4110672680056128 0.04957439898889293 0 +59 -0.2949531282594263 0.1696964684194627 0 +60 -0.3326337315558012 -0.02394056681171675 0 +61 -0.2591888307180861 0.05430298180422782 0 +62 -0.1844155319547547 0.225434844249134 0 +63 0.1241016241831161 -0.03647468372279483 0 +64 -0.2278481569914171 -0.02945516640958801 0 +65 -0.3478585245793106 0.2526913215845158 0 +66 -0.161205844314998 0.07463023671902519 0 +67 -0.3490702123863808 0.07598286426060075 0 +68 -0.3944934638826911 -0.01786866680851323 0 +69 0.07465432859153553 0.1328055091229944 0 +70 -0.3968639988992036 0.2330388255838553 0 +71 -0.01308622101161622 0.1889124611262808 0 +72 -0.08611705524327407 -0.04043037809088007 0 +73 -0.3356280325637821 0.1330945170345451 0 +74 -0.05610664395771314 0.02319076719143104 0 +75 -0.2341086841317337 0.2457836503867212 0 +76 -0.02572339439088428 0.1251759535595107 0 +77 0.01526807763460925 0.09482594390179022 0 +78 0.1121925646221863 0.07420405321443024 0 +79 -0.2032101632206738 0.1500732738934704 0 +80 -0.06276692144716189 0.2115099827024199 0 +81 -0.2503567830652025 0.1427843969385303 0 +82 -0.3094693246362213 0.03757542812143383 0 +83 -0.1033873234905965 0.1331281950774517 0 +84 -0.4350852003635158 0.1259755856984184 0 +85 0.1220010138429677 -0.1250919361177174 0 +86 -0.1514120427139066 0.03097273353090529 0 +87 0.07205497072710532 -0.1147701856174109 0 +88 -0.03390390062410364 -0.07019243891382632 0 +89 -0.2865556811816471 -0.04680729956271823 0 +90 0.03852677383736225 0.1827777422167275 0 +91 0.06083120697957822 0.02852472247025517 0 +92 0.01590516910001094 -0.03779922032616123 0 +93 -0.06649791017812356 0.1147801562249433 0 +94 -0.1711667442448505 0.1799979597303938 0 +95 0.06624135490197103 -0.06332092433843051 0 +96 0.1445980894625798 -0.08592748991927873 0 +97 -0.1260864413260195 -0.03641057853630991 0 +98 -0.4380247378467204 0.08696720312695849 0 +$EndNodes +$Elements +195 +1 15 2 0 1 1 +2 1 2 0 1 1 2 +3 1 2 0 1 2 3 +4 1 2 0 1 3 4 +5 1 2 0 1 4 5 +6 1 2 0 1 5 6 +7 1 2 0 1 6 7 +8 1 2 0 1 7 8 +9 1 2 0 1 8 9 +10 1 2 0 1 9 10 +11 1 2 0 1 10 11 +12 1 2 0 1 11 12 +13 1 2 0 1 12 13 +14 1 2 0 1 13 14 +15 1 2 0 1 14 15 +16 1 2 0 1 15 16 +17 1 2 0 1 16 17 +18 1 2 0 1 17 18 +19 1 2 0 1 18 19 +20 1 2 0 1 19 20 +21 1 2 0 1 20 21 +22 1 2 0 1 21 22 +23 1 2 0 1 22 23 +24 1 2 0 1 23 24 +25 1 2 0 1 24 25 +26 1 2 0 1 25 26 +27 1 2 0 1 26 27 +28 1 2 0 1 27 28 +29 1 2 0 1 28 29 +30 1 2 0 1 29 30 +31 1 2 0 1 30 31 +32 1 2 0 1 31 32 +33 1 2 0 1 32 1 +34 2 2 0 1 38 73 56 +35 2 2 0 1 6 68 60 +36 2 2 0 1 28 49 29 +37 2 2 0 1 5 68 6 +38 2 2 0 1 37 64 48 +39 2 2 0 1 48 61 37 +40 2 2 0 1 3 58 4 +41 2 2 0 1 56 73 44 +42 2 2 0 1 58 68 4 +43 2 2 0 1 65 70 30 +44 2 2 0 1 30 70 31 +45 2 2 0 1 60 82 48 +46 2 2 0 1 18 53 19 +47 2 2 0 1 58 67 39 +48 2 2 0 1 44 67 58 +49 2 2 0 1 49 65 29 +50 2 2 0 1 69 90 22 +51 2 2 0 1 21 69 22 +52 2 2 0 1 39 82 60 +53 2 2 0 1 34 74 50 +54 2 2 0 1 50 74 41 +55 2 2 0 1 18 63 53 +56 2 2 0 1 53 63 36 +57 2 2 0 1 25 51 26 +58 2 2 0 1 41 86 50 +59 2 2 0 1 46 74 57 +60 2 2 0 1 20 78 69 +61 2 2 0 1 42 83 51 +62 2 2 0 1 40 59 49 +63 2 2 0 1 20 69 21 +64 2 2 0 1 32 56 1 +65 2 2 0 1 49 59 38 +66 2 2 0 1 51 80 42 +67 2 2 0 1 51 83 35 +68 2 2 0 1 50 86 66 +69 2 2 0 1 32 70 56 +70 2 2 0 1 63 95 36 +71 2 2 0 1 37 61 52 +72 2 2 0 1 52 61 33 +73 2 2 0 1 61 82 33 +74 2 2 0 1 57 92 46 +75 2 2 0 1 9 54 10 +76 2 2 0 1 38 65 49 +77 2 2 0 1 33 82 67 +78 2 2 0 1 25 80 51 +79 2 2 0 1 35 66 52 +80 2 2 0 1 52 66 37 +81 2 2 0 1 57 74 34 +82 2 2 0 1 50 66 35 +83 2 2 0 1 51 62 26 +84 2 2 0 1 12 55 13 +85 2 2 0 1 48 89 60 +86 2 2 0 1 28 75 49 +87 2 2 0 1 69 77 45 +88 2 2 0 1 41 74 72 +89 2 2 0 1 72 74 46 +90 2 2 0 1 54 64 37 +91 2 2 0 1 50 93 34 +92 2 2 0 1 54 97 10 +93 2 2 0 1 49 75 40 +94 2 2 0 1 43 77 69 +95 2 2 0 1 17 63 18 +96 2 2 0 1 35 94 51 +97 2 2 0 1 60 68 39 +98 2 2 0 1 35 83 50 +99 2 2 0 1 45 90 69 +100 2 2 0 1 39 68 58 +101 2 2 0 1 59 73 38 +102 2 2 0 1 9 64 54 +103 2 2 0 1 8 89 64 +104 2 2 0 1 26 62 27 +105 2 2 0 1 45 76 71 +106 2 2 0 1 8 64 9 +107 2 2 0 1 71 76 42 +108 2 2 0 1 33 81 52 +109 2 2 0 1 10 72 11 +110 2 2 0 1 23 71 24 +111 2 2 0 1 36 92 57 +112 2 2 0 1 52 79 35 +113 2 2 0 1 12 88 55 +114 2 2 0 1 53 78 19 +115 2 2 0 1 64 89 48 +116 2 2 0 1 69 78 43 +117 2 2 0 1 6 60 7 +118 2 2 0 1 34 77 57 +119 2 2 0 1 57 77 43 +120 2 2 0 1 56 70 38 +121 2 2 0 1 38 70 65 +122 2 2 0 1 55 87 13 +123 2 2 0 1 29 65 30 +124 2 2 0 1 62 75 27 +125 2 2 0 1 43 78 53 +126 2 2 0 1 47 95 63 +127 2 2 0 1 53 91 43 +128 2 2 0 1 31 70 32 +129 2 2 0 1 36 91 53 +130 2 2 0 1 88 92 55 +131 2 2 0 1 4 68 5 +132 2 2 0 1 33 73 59 +133 2 2 0 1 62 94 40 +134 2 2 0 1 17 96 63 +135 2 2 0 1 14 85 15 +136 2 2 0 1 27 75 28 +137 2 2 0 1 40 81 59 +138 2 2 0 1 23 90 71 +139 2 2 0 1 58 98 44 +140 2 2 0 1 40 75 62 +141 2 2 0 1 56 84 1 +142 2 2 0 1 14 87 85 +143 2 2 0 1 10 97 72 +144 2 2 0 1 41 97 54 +145 2 2 0 1 37 86 54 +146 2 2 0 1 19 78 20 +147 2 2 0 1 55 95 87 +148 2 2 0 1 54 86 41 +149 2 2 0 1 40 94 79 +150 2 2 0 1 67 73 33 +151 2 2 0 1 44 73 67 +152 2 2 0 1 59 81 33 +153 2 2 0 1 44 84 56 +154 2 2 0 1 1 84 2 +155 2 2 0 1 57 91 36 +156 2 2 0 1 46 88 72 +157 2 2 0 1 72 88 11 +158 2 2 0 1 3 98 58 +159 2 2 0 1 43 91 57 +160 2 2 0 1 60 89 7 +161 2 2 0 1 63 96 47 +162 2 2 0 1 48 82 61 +163 2 2 0 1 24 80 25 +164 2 2 0 1 71 90 45 +165 2 2 0 1 45 77 76 +166 2 2 0 1 46 92 88 +167 2 2 0 1 67 82 39 +168 2 2 0 1 15 85 16 +169 2 2 0 1 13 87 14 +170 2 2 0 1 44 98 84 +171 2 2 0 1 87 95 47 +172 2 2 0 1 16 96 17 +173 2 2 0 1 66 86 37 +174 2 2 0 1 11 88 12 +175 2 2 0 1 71 80 24 +176 2 2 0 1 85 87 47 +177 2 2 0 1 42 80 71 +178 2 2 0 1 2 98 3 +179 2 2 0 1 7 89 8 +180 2 2 0 1 36 95 92 +181 2 2 0 1 92 95 55 +182 2 2 0 1 72 97 41 +183 2 2 0 1 22 90 23 +184 2 2 0 1 76 77 34 +185 2 2 0 1 51 94 62 +186 2 2 0 1 52 81 79 +187 2 2 0 1 76 93 42 +188 2 2 0 1 79 81 40 +189 2 2 0 1 42 93 83 +190 2 2 0 1 34 93 76 +191 2 2 0 1 79 94 35 +192 2 2 0 1 85 96 16 +193 2 2 0 1 83 93 50 +194 2 2 0 1 47 96 85 +195 2 2 0 1 84 98 2 +$EndElements diff --git a/test/blob2d-order1-h8e-2.msh b/test/blob2d-order1-h8e-2.msh new file mode 100644 index 0000000000000000000000000000000000000000..458dd8080af9883b5a6ef6ad25d0298569c1fb5c --- /dev/null +++ b/test/blob2d-order1-h8e-2.msh @@ -0,0 +1,191 @@ +$MeshFormat +2.2 0 8 +$EndMeshFormat +$Nodes +61 +1 -0.4725430905820004 0.1646517366169588 0 +2 -0.4773599478310085 0.08536814529218768 0 +3 -0.4673003194035515 0.006745495289718781 0 +4 -0.4274865461408011 -0.06073263122332467 0 +5 -0.3541954026674466 -0.0873996351975409 0 +6 -0.2749104861843841 -0.08306828870430079 0 +7 -0.1956604127888046 -0.07733118393429128 0 +8 -0.1167806313540748 -0.08536005048181899 0 +9 -0.04169549862147468 -0.1111600401018526 0 +10 0.03068646850135696 -0.1440042368191183 0 +11 0.105115344896881 -0.1714461181374037 0 +12 0.1708662314746014 -0.1417406086771747 0 +13 0.1840528408459943 -0.06394600932027859 0 +14 0.1760797584763959 0.01503106303733582 0 +15 0.155273396719484 0.09166228767345191 0 +16 0.1210947052116375 0.1632575553011128 0 +17 0.06689182154281072 0.220469224870406 0 +18 -0.007932619693392521 0.2456642003733992 0 +19 -0.08672919877000161 0.2560933881248318 0 +20 -0.1644749394707306 0.2724307479360701 0 +21 -0.2413464133764506 0.2926542108090772 0 +22 -0.3195170592956217 0.3065894623973552 0 +23 -0.3975412741958552 0.2964955465946852 0 +24 -0.4521497617068335 0.2411179496461016 0 +25 -0.2877066332189855 0.1074352417349247 0 +26 -0.02967539994291918 0.06940577338656762 0 +27 -0.1519591500011385 0.1340964507511098 0 +28 0.05636333569205443 -0.02150015435707153 0 +29 -0.2100008319785572 0.03288947751568416 0 +30 -0.3488151832811237 0.1899639632680639 0 +31 -0.3533596425773649 0.01755290244392255 0 +32 -0.2237006036557083 0.1975718606623247 0 +33 -0.05870269643742826 0.1624567848912205 0 +34 -0.09496176491301939 0.02396963659710527 0 +35 -0.4006908906205165 0.127544374001792 0 +36 0.06155424142750854 0.09087780220354552 0 +37 -0.0215522133838158 -0.01709743190359331 0 +38 0.09887082457626195 -0.09266716974157363 0 +39 -0.2839774703351805 0.006076845187414082 0 +40 0.02127917470037199 0.1578123568797543 0 +41 -0.2951064441886064 0.2365401984896358 0 +42 0.1053025359059351 0.03254813698369635 0 +43 -0.1004918428432163 0.08346343156095137 0 +44 -0.1371133176670015 0.2045298464731114 0 +45 -0.2236092738718777 0.1086994325157413 0 +46 0.02247586384968048 -0.07727887490083041 0 +47 -0.4098421275061732 0.05099858116136879 0 +48 0.1241338590993283 -0.02610682667957831 0 +49 -0.1550374982025248 -0.01372171564088193 0 +50 0.03207952495314315 0.03237579220853598 0 +51 -0.4053840647762471 0.1959854965170016 0 +52 -0.2936743443918267 0.1707302751834731 0 +53 -0.328891299386235 -0.03231707848678983 0 +54 -0.3697522979073813 0.2444487694854738 0 +55 -0.2657865018278374 0.0634449307228146 0 +56 -0.1576748114705166 0.06817034238378727 0 +57 -0.3340130632489786 0.07254214496019823 0 +58 -0.2239426006322122 -0.02593572831442603 0 +59 -0.3901792229469287 -0.01752539433544087 0 +60 -0.3340509240107165 0.133957182583822 0 +61 -0.08450731847966086 -0.03933390273190407 0 +$EndNodes +$Elements +121 +1 15 2 0 1 1 +2 1 2 0 1 1 2 +3 1 2 0 1 2 3 +4 1 2 0 1 3 4 +5 1 2 0 1 4 5 +6 1 2 0 1 5 6 +7 1 2 0 1 6 7 +8 1 2 0 1 7 8 +9 1 2 0 1 8 9 +10 1 2 0 1 9 10 +11 1 2 0 1 10 11 +12 1 2 0 1 11 12 +13 1 2 0 1 12 13 +14 1 2 0 1 13 14 +15 1 2 0 1 14 15 +16 1 2 0 1 15 16 +17 1 2 0 1 16 17 +18 1 2 0 1 17 18 +19 1 2 0 1 18 19 +20 1 2 0 1 19 20 +21 1 2 0 1 20 21 +22 1 2 0 1 21 22 +23 1 2 0 1 22 23 +24 1 2 0 1 23 24 +25 1 2 0 1 24 1 +26 2 2 0 1 20 32 21 +27 2 2 0 1 20 44 32 +28 2 2 0 1 33 40 26 +29 2 2 0 1 26 37 34 +30 2 2 0 1 36 40 16 +31 2 2 0 1 1 35 2 +32 2 2 0 1 18 33 19 +33 2 2 0 1 26 40 36 +34 2 2 0 1 32 44 27 +35 2 2 0 1 18 40 33 +36 2 2 0 1 24 54 51 +37 2 2 0 1 15 36 16 +38 2 2 0 1 16 40 17 +39 2 2 0 1 27 45 32 +40 2 2 0 1 10 38 11 +41 2 2 0 1 45 52 32 +42 2 2 0 1 25 52 45 +43 2 2 0 1 23 54 24 +44 2 2 0 1 1 51 35 +45 2 2 0 1 33 44 19 +46 2 2 0 1 27 44 33 +47 2 2 0 1 26 43 33 +48 2 2 0 1 33 43 27 +49 2 2 0 1 10 46 38 +50 2 2 0 1 32 41 21 +51 2 2 0 1 31 57 39 +52 2 2 0 1 39 57 55 +53 2 2 0 1 38 46 28 +54 2 2 0 1 34 43 26 +55 2 2 0 1 17 40 18 +56 2 2 0 1 5 59 53 +57 2 2 0 1 4 59 5 +58 2 2 0 1 26 50 37 +59 2 2 0 1 11 38 12 +60 2 2 0 1 35 51 30 +61 2 2 0 1 12 38 13 +62 2 2 0 1 30 60 35 +63 2 2 0 1 49 56 34 +64 2 2 0 1 15 42 36 +65 2 2 0 1 29 56 49 +66 2 2 0 1 35 47 2 +67 2 2 0 1 36 50 26 +68 2 2 0 1 37 50 28 +69 2 2 0 1 51 54 30 +70 2 2 0 1 28 48 38 +71 2 2 0 1 9 61 37 +72 2 2 0 1 28 46 37 +73 2 2 0 1 37 46 9 +74 2 2 0 1 14 42 15 +75 2 2 0 1 7 49 8 +76 2 2 0 1 49 61 8 +77 2 2 0 1 37 61 34 +78 2 2 0 1 21 41 22 +79 2 2 0 1 8 61 9 +80 2 2 0 1 19 44 20 +81 2 2 0 1 39 58 6 +82 2 2 0 1 6 53 39 +83 2 2 0 1 38 48 13 +84 2 2 0 1 2 47 3 +85 2 2 0 1 42 50 36 +86 2 2 0 1 5 53 6 +87 2 2 0 1 28 50 42 +88 2 2 0 1 47 57 31 +89 2 2 0 1 6 58 7 +90 2 2 0 1 35 57 47 +91 2 2 0 1 34 56 43 +92 2 2 0 1 9 46 10 +93 2 2 0 1 34 61 49 +94 2 2 0 1 13 48 14 +95 2 2 0 1 39 53 31 +96 2 2 0 1 55 57 25 +97 2 2 0 1 29 58 39 +98 2 2 0 1 39 55 29 +99 2 2 0 1 53 59 31 +100 2 2 0 1 24 51 1 +101 2 2 0 1 42 48 28 +102 2 2 0 1 22 54 23 +103 2 2 0 1 3 59 4 +104 2 2 0 1 41 52 30 +105 2 2 0 1 32 52 41 +106 2 2 0 1 14 48 42 +107 2 2 0 1 30 54 41 +108 2 2 0 1 41 54 22 +109 2 2 0 1 27 56 45 +110 2 2 0 1 45 56 29 +111 2 2 0 1 29 55 45 +112 2 2 0 1 45 55 25 +113 2 2 0 1 43 56 27 +114 2 2 0 1 47 59 3 +115 2 2 0 1 49 58 29 +116 2 2 0 1 31 59 47 +117 2 2 0 1 52 60 30 +118 2 2 0 1 7 58 49 +119 2 2 0 1 25 60 52 +120 2 2 0 1 35 60 57 +121 2 2 0 1 57 60 25 +$EndElements diff --git a/test/make-blob-meshes.sh b/test/make-blob-meshes.sh index eb9c446531fdf145640ccb988b0fd9549a190621..c5e633dc5d6c9f24ec6634b2c8d40ef3cb6b2dce 100755 --- a/test/make-blob-meshes.sh +++ b/test/make-blob-meshes.sh @@ -1,10 +1,11 @@ #! /bin/bash -order=4 -for i in 8e-2 6e-2 4e-2 ; do - gmsh -2 -order $order \ - -string "Mesh.CharacteristicLengthMax = $i;" \ - -string "Geometry.OCCTargetUnit='MM';" \ - -format msh2 \ - -o blob2d-order$order-h$i.msh \ - blob-2d.step +for order in 1 4; do + for i in 8e-2 6e-2 4e-2 ; do + gmsh -2 -order $order \ + -string "Mesh.CharacteristicLengthMax = $i;" \ + -string "Geometry.OCCTargetUnit='MM';" \ + -format msh2 \ + -o blob2d-order$order-h$i.msh \ + blob-2d.step + done done diff --git a/test/test_refinement.py b/test/test_refinement.py index 8655ed774852e3d54b0713b5ecc95c80f372c771..8914bc5716315b68624f5a4c541e315e26a59214 100644 --- a/test/test_refinement.py +++ b/test/test_refinement.py @@ -22,6 +22,9 @@ OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE. """ +import logging +from functools import partial + import pytest import pyopencl as cl import pyopencl.clmath # noqa @@ -41,20 +44,22 @@ from meshmode.discretization.poly_element import ( PolynomialEquidistantSimplexGroupFactory, ) -import logging logger = logging.getLogger(__name__) -from functools import partial +def get_blob_mesh(mesh_par, order=4): + # from meshmode.mesh.io import generate_gmsh, FileSource + # return generate_gmsh( + # FileSource("blob-2d.step"), 2, order=order, + # force_ambient_dim=2, + # other_options=[ + # "-string", "Mesh.CharacteristicLengthMax = %s;" % mesh_par] + # ) -def gen_blob_mesh(h=0.2, order=1): - from meshmode.mesh.io import generate_gmsh, FileSource - return generate_gmsh( - FileSource("blob-2d.step"), 2, order=order, - force_ambient_dim=2, - other_options=[ - "-string", "Mesh.CharacteristicLengthMax = %s;" % h] - ) + from meshmode.mesh.io import read_gmsh + return read_gmsh( + "blob2d-order%d-h%s.msh" % (order, mesh_par), + force_ambient_dim=2) def random_refine_flags(fract, mesh): @@ -117,7 +122,7 @@ def uniform_refine_flags(mesh): 3), ("blob2d_rand", - gen_blob_mesh, + partial(get_blob_mesh, "6e-2", order=1), partial(random_refine_flags, 0.4), 4), @@ -145,7 +150,7 @@ def test_refinement(case_name, mesh_gen, flag_gen, num_generations): r = Refiner(mesh) - for igen in range(num_generations): + for _ in range(num_generations): flags = flag_gen(mesh) mesh = r.refine(flags) @@ -163,11 +168,11 @@ def test_refinement(case_name, mesh_gen, flag_gen, num_generations): ]) @pytest.mark.parametrize(("mesh_name", "dim", "mesh_pars"), [ ("circle", 1, [20, 30, 40]), - ("blob", 2, [8e-2, 4e-2, 2e-2]), + ("blob", 2, ["8e-2", "6e-2", "4e-2"]), ("warp", 2, [4, 5, 6]), ("warp", 3, [4, 5, 6]), ]) -@pytest.mark.parametrize("mesh_order", [1, 5]) +@pytest.mark.parametrize("mesh_order", [1, 4, 5]) @pytest.mark.parametrize("refine_flags", [ # FIXME: slow #uniform_refine_flags, @@ -194,10 +199,6 @@ def test_refinement_connection( from pytools.convergence import EOCRecorder eoc_rec = EOCRecorder() - def f(x): - from six.moves import reduce - return 0.1 * reduce(lambda x, y: x * cl.clmath.sin(5 * y), x) - for mesh_par in mesh_pars: # {{{ get mesh @@ -211,8 +212,8 @@ def test_refinement_connection( if mesh_order == 5: pytest.xfail("https://gitlab.tiker.net/inducer/meshmode/issues/2") assert dim == 2 - h = mesh_par - mesh = gen_blob_mesh(h, mesh_order) + mesh = get_blob_mesh(mesh_par, mesh_order) + h = float(mesh_par) elif mesh_name == "warp": from meshmode.mesh.generation import generate_warped_rect_mesh mesh = generate_warped_rect_mesh(dim, order=mesh_order, n=mesh_par) @@ -222,6 +223,22 @@ def test_refinement_connection( # }}} + from meshmode.mesh.processing import find_bounding_box + mesh_bbox_low, mesh_bbox_high = find_bounding_box(mesh) + mesh_ext = mesh_bbox_high-mesh_bbox_low + + def f(x): + result = 1 + if mesh_name == "blob": + factor = 15 + else: + factor = 9 + + for iaxis in range(len(x)): + result = result * cl.clmath.sin(factor * (x[iaxis]/mesh_ext[iaxis])) + + return result + discr = Discretization(cl_ctx, mesh, group_factory(order)) refiner = refiner_cls(mesh) @@ -266,9 +283,13 @@ def test_refinement_connection( err = la.norm((f_interp - f_true).get(queue), np.inf) eoc_rec.add_data_point(h, err) + order_slack = 0.5 + if mesh_name == "blob" and order > 1: + order_slack = 1 + print(eoc_rec) assert ( - eoc_rec.order_estimate() >= order-0.5 + eoc_rec.order_estimate() >= order-order_slack or eoc_rec.max_error() < 1e-14)