From 42fb81d04300abcde628831af2262956cf4cd0f1 Mon Sep 17 00:00:00 2001 From: Roy Stogner Date: Thu, 23 Mar 2023 16:53:33 -0500 Subject: [PATCH 01/17] discontinuous_spline_extraction file option This gets us a lot more flexibility with what IGA meshes we can support. Refs #18768 --- framework/src/meshgenerators/FileMeshGenerator.C | 11 +++++++++++ 1 file changed, 11 insertions(+) diff --git a/framework/src/meshgenerators/FileMeshGenerator.C b/framework/src/meshgenerators/FileMeshGenerator.C index 1141ec97b522..c3f44099f670 100644 --- a/framework/src/meshgenerators/FileMeshGenerator.C +++ b/framework/src/meshgenerators/FileMeshGenerator.C @@ -45,6 +45,14 @@ FileMeshGenerator::validParams() "If clear_spline_nodes=true, IsoGeometric Analyis spline nodes " "and constraints are removed from an IGA mesh, after which only " "C^0 Rational-Bernstein-Bezier elements will remain."); + params.addParam("discontinuous_spline_extraction", + false, + "If discontinuous_spline_extraction=true, " + "Rational-Bernstein-Bezier elements extracted from a spline mesh " + "will be disconnected from neighboring elements, coupled only via " + "their extraction operators. This may be less efficient than the " + "default C^0 extracted mesh, but may be necessary if the extracted " + "mesh is non-conforming."); params.addClassDescription("Read a mesh from a file."); return params; } @@ -87,6 +95,9 @@ FileMeshGenerator::generate() { if (mesh->processor_id() == 0) { + if (getParam("discontinuous_spline_extraction")) + exreader->set_discontinuous_bex(true); + exreader->read(_file_name); if (getParam("clear_spline_nodes")) From 3ea748ee3027cd916fc3abbd9fcdfd6ddb041812 Mon Sep 17 00:00:00 2001 From: Roy Stogner Date: Thu, 23 Mar 2023 17:10:46 -0500 Subject: [PATCH 02/17] Discontinuous spline extraction regression test --- .../2d_discontinuous_iga.i | 53 ++++++++++++++++++ .../gold/2d_discontinuous_iga_out.e | Bin 0 -> 35480 bytes .../meshgenerators/file_mesh_generator/tests | 15 +++++ 3 files changed, 68 insertions(+) create mode 100644 test/tests/meshgenerators/file_mesh_generator/2d_discontinuous_iga.i create mode 100644 test/tests/meshgenerators/file_mesh_generator/gold/2d_discontinuous_iga_out.e diff --git a/test/tests/meshgenerators/file_mesh_generator/2d_discontinuous_iga.i b/test/tests/meshgenerators/file_mesh_generator/2d_discontinuous_iga.i new file mode 100644 index 000000000000..de02a4ba24f8 --- /dev/null +++ b/test/tests/meshgenerators/file_mesh_generator/2d_discontinuous_iga.i @@ -0,0 +1,53 @@ +[Mesh] + [cyl2d_iga] + type = FileMeshGenerator + file = PressurizedCyl_Patch6_4Elem.e + discontinuous_spline_extraction = true + [] + allow_renumbering = false + parallel_type = replicated +[] + +[Variables] + [u] + order = SECOND # Must match mesh order + family = RATIONAL_BERNSTEIN + [] +[] + +[Kernels] + [time] + type = TimeDerivative + variable = u + block = 0 # Avoid direct calculations on spline nodes + [] + [diff] + type = Diffusion + variable = u + block = 0 # Avoid direct calculations on spline nodes + [] + [null] + type = NullKernel + variable = u + block = 1 # Keep kernel coverage check happy + [] +[] + +[ICs] + [u] + type = FunctionIC + variable = u + function = 'sin(x)' + [] +[] + +[Executioner] + type = Transient + num_steps = 2 + solve_type = NEWTON + dtmin = 1 +[] + +[Outputs] + exodus = true +[] diff --git a/test/tests/meshgenerators/file_mesh_generator/gold/2d_discontinuous_iga_out.e b/test/tests/meshgenerators/file_mesh_generator/gold/2d_discontinuous_iga_out.e new file mode 100644 index 0000000000000000000000000000000000000000..c6d6edce714398cb36699fa70c0b0a265860412b GIT binary patch literal 35480 zcmeHQ3zQsHnVv{u9)yqx$V;FQBq+{IGMPL{TZEZOm^CkwNr(n0s=KRuigZ`CRo#;r zbOm)0R33o<;bDYb7GZhl0*WhVS=sp5BglffEGCl~CJ=?VyUNJ|a+GEE`)}2)s(w^= z-b*roHSO5M0d+-0>|K7UQ(YgG6O-)UQAWlO(3~@?EGZaJ3YE4ZOaN0CA zuV*#I&TBb3bHZe!kX4d;md+^uR5~}zq-HzU+i{pe|7jVGZaNZinsctiGh5hAly8Rn zpKj3KbJ_nB@q4(7r{pvX80^B-visZE{WHt&XBUY0Lpe^LYSGQ~AEQe*QYKCdbWY{X z=myH$bS~m_#L0+H} z92?1|NqIe=p$8%mkB@IlLR5CbG(cFvEZB;kRu!|5Z*guq6!*+PoW|t5wbyhr6D<^Z zcs~+ToRiD%z;UV;_X)r`-Ag>;foO8?pR|e6xIbr^IT9wnmv~L-iB5fE^d8dPiGQ3j z@%$Iz{QU_N-0Y`IEDzOL9zffKJDihA#4P9D0rn5!huk}v@QE&ZJn@nE#mm5Ho2D$+ zGntKQrr_|4?jcMHNYg2WDfi%K8qw&XIT!g2f_T#iZ#uOpg(`7YrRG;!7UJ8F2u(wzH7 zV_*5tAK>~Q%AWt5JHC;+bYNB4bN64c`F#1y`{V<+Aly39JuZ}AIPREpM{l@7=J`tL z=lM$M=J}qT-*V)WCwx10-Qpj}qoY3y$1f@8^AFs6?WLWY%gQ%6`_xNDFZ{GeS1CTv zSIQ@zuY^BK9|D&~M_*@hdGo8?^6%;NpLV)@o3!)+nf}vvGnig z^5`v0E_43vhmM@)AM48N$fb9)>#s4nv`#b*5I=x;1R{+`jzat(;?ao5ARdc&9OChaA3{6<@kGRv5GBNu5kHK03gW4V zA3;0~aW>-Vh#y7#7~&a-A4hCPY(Z>AoP$W_2H8>b5NUiw^bjrdT%v(!BYKDyP9M=j zw9s=UM52x8AzCEaXdv2%9-@WQHv#Ff@OP~Fi{wKC zv-fg1r9x8{7+eWwkd zyyMnuynN4%K5>8k+6|tK`r^bd&FQ`Ja{1YZbCdU0An=r?M~^ykh3?sSzX`+NncW@R zvrB$)3ZICb@7kDdiUsl&`>RU;{_#Bf|%W~B>R&B)slE;WI!U!G6+7kl_l zz8&(8dBSwwt3BS<~6rbS5O{S>hV>Y-OG#QL#E>GYTLX8Z47K~eE$5lWSh2FU9dR0 zxV>FXEL7VT;o{6Rm0+_zMxc`GN*NSdx^2p?{{JDQbe6G4lfX4N#I|(z0Z7 zot9yZ-0szLt5icxYgQN?5*nZdkycKjnzm{|-zx2FjF=dAaT{mheU`|h{m_v(%#n|B1 zLp9l_8i`7sLUbW%|ssgD{2PESSkg~CIWD6-Z(zwUb7|D?hk+jlj22Y9v4zlmCG(pVNfkqXwXEVSX(Xv%%n|A&-Z5iLdPVh zg7T(S>nKSJVcFbEMSE`j0r`k`_j|knMPCl zI@NGn0VAPR^Cnl*VqL5xWt0k)B`9UlYLme_jxtdc0nS0FRa?C4L}8_#-$QhGIN!@q zD*-@wUR{0cp#%~!$#oCZLUOMc0V+9oN#XJc;GrJow7i|D6OL2>aj-vyWh{kU5)SPk z7^P!Ud)wmn#S0d;Ev~EF08FLOe}#@9W75`K-Q?j$Nivv}4bvcIs+N+2iI_7zGZ8`& zkg)>@r=gWEM8QZY%-<}!Lh&o^s@O!$l=Qz#CFTlAcw^^<3prZOETjgKkLL)*~oYzl3e zSXoKRqH5G}Ux!-lF==(zc|B`Z^CHzcN|Fr@<>}Q8znH;P1fkdt-zwR|Oeh#U3fVR` z9FbLAK@o()8eXX9CJLaYCI{GCAsXh9^Kd4j`gm*{BcRhy1j!vEAO;thp$K}i#dQji zLvJbRc&9)!%FIOZY7+!xB`&s@#spXl<;W!(ARhvh@Q)xF59iiz#Yz$uP8qg>>RsB; za0*37QVbB5O_t3?*?GKnl+bM3bm!N_bxIJhtiO(;q#WOc=m!r z)URf-Y(rP!y)QQ)!~wCyEsoW4jtI(sMFHOhQ|jFvc0;DpYw4QE*ClClOKYKlG?G}W zu3BPimb#rqJ{^1--oX6eO80qDR& zormB@Oj-_)4k^O!j(S8-JY=hjC{Zay;bGh?7)jL{sDuX!MaRpMRCN+53J=NcRb(v$ zNsURiZk*m<$C@;U!<5&Z4Jgoxwj5NmW+eRLi*{cGNH4D{80TPFgJKs_m|&rQ-4LZ7 z@^L_=!6>QGdTNVz57EI@6DO+Z7?T1G)yo#1px&mXwSm3OyFdkeQnL&#^A1)6Oa+if zvr%%PDpFeoIDGLt{w<78F1K2V9v=U+Tx?GK)D*PMnkF~ zBQc6lYz?nTg@q`V0}al2I5 zWO|t_jdzk>F0y>6u8r&E(x%q~k#%ActdG5AS*tQp?mxF)9v|M;kO~RqJvhA06tc8D zmr7qAqUe+wdU>`pV3w#;p`Bx@2O?>Fy!`GPFBhqd>LSL~hL4xu1)g@=Risf$o)}!i zi%m`p=l$hQ*U6Y%YngaSBi9*rg zt!^YT1-$H!-4d{75Q~%3y{7H2vJ&rLoyEx|Fw{*tK^%~|#T7lw;U_K>2e5QaGAA&C zl=<9k2$2yebSHz$wkY8|6s?Cai0!;pz^kKdXNQb#V?8Hr3BtCorBf8UehW;JSS`Wg z9k=MZR=jdayvXH%sEjFRFBX(q9I)ydf+#J;X8RiU6z#HAFECa}5}=1tW;E3c#JkxV%O;42w0Tf@&&&ZzkCKE9cxcrZ zFPhip+=MuVw8+IS#FKd4pWb*ZUEU{NZJh0Bu)kCx8^?I5)bNZzVJ6qBDwc_zz=KOI zdC|owjBEKScoAAk)p^G#O)p;1GfCPk!0m0!tr~WS@McxSd6>ZV8J4MI-wNJD&^(-X z;fbf>JoNiRuLE@YL$7-TX#J{`OcCtCZJnci*5WayQlifi&> z(z>p-E4w?eEm>!k`s`i8@Z0Yb0)%-<@fV^i3pW7RF33X%&|OgsI}9&Hr~Q>~=xGcM zi+Q~w)k;W`b&2A1IQS4n1<}N$Pu0RGy6##Ij=iC5?ZxLdXq>@bxM0N=?;P;jjNd?uGcmjfUMZ>D4-+;V4l^J5kcHVu)4U7npVgwp+=VBujyI{bkmI z-O6^ir2l;rqIE4^UbG7q^mMXgAfv1$wY@-T?}A&&{kTqrTGK>SRn>b@r0n=bDi)Y* z8zGg#C}a!!wqjeGN?rnS-xriXFE-+&@2FHFG7^Q=#-mXHNh@Leh>FPSqx)v6ORhLz zgBL-)2c5c=TF;pNdmHrfxRl2vk+&oBxRl4v)y-?1T|gROyE??wuwO|@(eD4I;XK%} zUR@$876-Ohj^680D(wPAowj(79RMUnfT`RfT_tFVG{KmojXl|?0E>Y5i8zG`GgG>4 z1qY*}kSzzUvM~VOGYGYoA_0cK{znvw&YHQSl%N@0Tx3z(CSmsl0ZP;2LP?q_e4eiq z$2%7X)YumTutm22?!1SFL`<69y}G-nd)Z2p~$!ic9%o zJER6WPY~xJwy$8HHl_6u7z>n=D-Kw792G^%&NxcZQ*4(K_`$GEvF`6$B`#&`4UxhQ zVrnk<&WA4{@gn6B7HCT$-P`_5HV=wW}- z5-8LawpBD-BM=GPnQ6=pba}OkXRfDJt4s|0)nNjHP zAz4;xF&a4!eZZ1Zi@`A|GYXlkGiLb$I`JZtgTQYXWhwEd=WPMlixmiF*~{+{C>$^n z&R!j5)X^W}whea3X{?g!%WtVjn)|1Ck^7TsqXC1vag||5pfDSz@EEu#;Jd)hvXO6a z+_1q}Sl1ize>MRB^KVqb`QnTn`B~8hCA>gXHQa$k*?YWhwGr|SZgHv%A)?tGGkQ%6Vf4)$%>N~P+z zhGNqCu8uXU>zcVGNmHr3jZo>NG(;(K<%wk%IVP>5Zw&IfcQ*Djx?$w!w&5$NWz@l=F|MPN7^4 zpEXNW8!$(q5c^N#<%EzjV--P0#L+@9&R8l0F-kpT%adWLSENzA$nPT8z0AMu%8rK} zCGXz#f$0an@P_wo)&sK+x%q;7Hsf2bbo{LR+?ovse{p?tY`Z%1($UeE<^4Miefe=q z+@mW0Vx@fDtUtVfZ>A37JMxPAP1$JQmOWE@RDS7^ukY%(H!kyU#D;$R)%_-I)o;r# zJpj5VTsu0t-}@%5^X=ED-T7&rysoZn~!alM=tpHMVp?OERUqGOWd^eWn5o8{mo}P+;E&cqTY1c_1j*R2lpI)=(*pK*e7q z-uKE?ldqRYzV!4=bLDnk<>3czKYsPe z56VL?_8;@?#dpXMk0Gk-ntEnMHZYW)-Yj8QAoUh$F>nEqN! ze-3o-@0fDUPrQ6XZ$8lc*e7?%LqiWw*nHypaQwTU9(nk! Date: Tue, 27 Jun 2023 14:46:30 -0500 Subject: [PATCH 03/17] Mark *why* we turn off parallelism and renumbering --- .../meshgenerators/file_mesh_generator/2d_diffusion_iga.i | 4 ++-- .../file_mesh_generator/2d_diffusion_iga_nosplines.i | 4 ++-- .../file_mesh_generator/3d_steady_diffusion_iga.i | 4 ++-- test/tests/meshgenerators/file_mesh_generator/tests | 3 +++ 4 files changed, 9 insertions(+), 6 deletions(-) diff --git a/test/tests/meshgenerators/file_mesh_generator/2d_diffusion_iga.i b/test/tests/meshgenerators/file_mesh_generator/2d_diffusion_iga.i index 21156bf482a1..090fba318e84 100644 --- a/test/tests/meshgenerators/file_mesh_generator/2d_diffusion_iga.i +++ b/test/tests/meshgenerators/file_mesh_generator/2d_diffusion_iga.i @@ -3,8 +3,8 @@ type = FileMeshGenerator file = PressurizedCyl_Patch6_4Elem.e [] - allow_renumbering = false - parallel_type = replicated + allow_renumbering = false # VTK diffs via XMLDiff are + parallel_type = replicated # really fragile [] [Variables] diff --git a/test/tests/meshgenerators/file_mesh_generator/2d_diffusion_iga_nosplines.i b/test/tests/meshgenerators/file_mesh_generator/2d_diffusion_iga_nosplines.i index 87d991293583..52bfa132a5a4 100644 --- a/test/tests/meshgenerators/file_mesh_generator/2d_diffusion_iga_nosplines.i +++ b/test/tests/meshgenerators/file_mesh_generator/2d_diffusion_iga_nosplines.i @@ -4,8 +4,8 @@ file = PressurizedCyl_Patch6_4Elem.e clear_spline_nodes = true [] - allow_renumbering = false - parallel_type = replicated + allow_renumbering = false # VTK diffs via XMLDiff are + parallel_type = replicated # really fragile [] [Variables] diff --git a/test/tests/meshgenerators/file_mesh_generator/3d_steady_diffusion_iga.i b/test/tests/meshgenerators/file_mesh_generator/3d_steady_diffusion_iga.i index 27fb0501c7b0..bced72c9273d 100644 --- a/test/tests/meshgenerators/file_mesh_generator/3d_steady_diffusion_iga.i +++ b/test/tests/meshgenerators/file_mesh_generator/3d_steady_diffusion_iga.i @@ -3,8 +3,8 @@ type = FileMeshGenerator file = Cube_With_Sidesets.e [] - allow_renumbering = false - parallel_type = replicated + allow_renumbering = false # VTK diffs via XMLDiff are + parallel_type = replicated # really fragile [] [Variables] diff --git a/test/tests/meshgenerators/file_mesh_generator/tests b/test/tests/meshgenerators/file_mesh_generator/tests index 78e99198edf9..85304abb7b78 100644 --- a/test/tests/meshgenerators/file_mesh_generator/tests +++ b/test/tests/meshgenerators/file_mesh_generator/tests @@ -35,6 +35,7 @@ dof_id_bytes = 8 # Set to >=9.0 once #21449 is resolved vtk_version = '>=9.1' + # VTK diffs via XMLDiff are really fragile max_parallel = 1 vtk = true xmldiff = '2d_diffusion_iga_out_002.pvtu 2d_diffusion_iga_out_002_0.vtu' @@ -57,6 +58,7 @@ dof_id_bytes = 8 # Set to >=9.0 once #21449 is resolved vtk_version = '>=9.1' + # VTK diffs via XMLDiff are really fragile max_parallel = 1 vtk = true xmldiff = '2d_diffusion_iga_nosplines_out_002.pvtu 2d_diffusion_iga_nosplines_out_002_0.vtu' @@ -78,6 +80,7 @@ dof_id_bytes = 8 # Set to >=9.0 once #21449 is resolved vtk_version = '>=9.1' + # VTK diffs via XMLDiff are really fragile max_parallel = 1 vtk = true xmldiff = '3d_steady_diffusion_iga_out_001.pvtu 3d_steady_diffusion_iga_out_001_0.vtu' From 821ced87649fca82f6a3b5fc17bec34bd939887f Mon Sep 17 00:00:00 2001 From: Roy Stogner Date: Tue, 27 Jun 2023 14:47:31 -0500 Subject: [PATCH 04/17] Rename test I like being able to use `--re _iga` with `run_tests` --- test/tests/meshgenerators/file_mesh_generator/tests | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/test/tests/meshgenerators/file_mesh_generator/tests b/test/tests/meshgenerators/file_mesh_generator/tests index 85304abb7b78..6035ad4ff90a 100644 --- a/test/tests/meshgenerators/file_mesh_generator/tests +++ b/test/tests/meshgenerators/file_mesh_generator/tests @@ -90,7 +90,7 @@ issues = '#18768' [] - [./exodus_discontinuous_spline_2d_test] + [./exodus_discontinuous_iga_2d_test] # Make this an XMLDiff once #21449 is resolved type = Exodiff input = '2d_discontinuous_iga.i' From 2bcb2cb76625339fd3e95176ccaa174662e61c9f Mon Sep 17 00:00:00 2001 From: Roy Stogner Date: Thu, 5 Oct 2023 17:19:44 -0500 Subject: [PATCH 05/17] Add some exodiff IGA tests The discontinuous IGA option is giving me problems in parallel, and we don't have enough test coverage (at the MOOSE level, vs the libMesh level) of parallel IGA, so let's fix that. These ExodusII outputs are kind of garbage, since they're turning properly weighted Rational-Bernstein-Bezier elements as plain Lagrange elements, but they're good enough to tell whether something is changing in a regression test. --- .../gold/2d_diffusion_iga_nosplines_out.e | Bin 0 -> 32484 bytes .../gold/2d_diffusion_iga_out.e | Bin 0 -> 35640 bytes .../meshgenerators/file_mesh_generator/tests | 35 +++++++++++++++++- 3 files changed, 34 insertions(+), 1 deletion(-) create mode 100644 test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_nosplines_out.e create mode 100644 test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_out.e diff --git a/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_nosplines_out.e b/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_nosplines_out.e new file mode 100644 index 0000000000000000000000000000000000000000..c15f6d9f39d2f6eb71fd373f6e89e0b69c622951 GIT binary patch literal 32484 zcmeHQ36LDsd0w!D4#X8MIkw>$1ej2}t6gaq2{hF#tA(u5g|rBRAnl%+o|%?&w1Ht8laFk;R#vtu_Wr;(s#Bma1gTZXR|MfXD zJ;!Se2^8wu-hO?3^S}T9-~ayizoQ42pL1YGN5{vo&&7T)_Sw2>$c9o>J30=)X~!H# zE2^^Xs3kgcw^>HHC}*`IozeSq=-e=~s_owYI1aPvnyRbx&<23) zO&(P(TQd!ucFe?nz`y}dZozR@o$K>+x~Vvtk(W!RW;pU@)5^-Wb~#}ao&5K4oRiaZ z@W{OpxkvYpEGlDKk@y3i%3yBZDd2MOK0|sQ!21`2$%<1Z7W&vocFc4%N2edeM!f!X zZx*c2<;s*1S<5S6yj`LU+p<}9y1YlZEDyo+bGiIB7fcU?XrkTU4nktKdvb+?;~X^| z5`uGjmU!iFL>ujNn8Lr=c)nzrB@&YVEb*MuGh995bboAkIj*?x`tR?-`9T>ITpi>0 z=pE8UePZiwYI1>bwwtve>cBq3363Ov6w!4)K_P zh%cnm{%_Brawa?DB~SSLne#6gSiTewVI#aR^Hj&2ojA_p&sY9FJ>zebKHNQ-tK&G2 z>n@0v-GWE>p3Vd`w8+o7SK!@4Mit#AgF$$_bi4Vg=0imOSI~PXCplBL)C@Gnmn+!? zvR(cL+DdP^CkN9zT#(Xl+$+*+U%M-J)%Z&3mw!6{{OPUNNPF+VcIOQ5{sL*wr5i4h zzO>7`e!BGU&(ijL_xRjtslU1Qr0MDDH&Zt+eO#KJepY(rug35AyH-HT{jBK_jq-`aB3@(tek#nNjzZ27~z z<6!E#pFIxz-=}U`umeX~`pHjzeEh-(rlsd^x^`{P_JzwH!PV*M-|>80`Q9!!AHBCN zyZ?B)zMa$a?)AfNKIT8n>3Mqg)@9$8=>995uFP=O%kNAtp9e|*Nw*){=H>S-_xk7e zxaT)K#PjjqGo#D?mY$n_ljq~g+yB7n>z;Svvaj>&cW^ov{_~H#{PlVHI?6qt{dxDi zkLrj0+kOswlsQ0liR>uVcVus=e)8=)*;e1Kkxxc;p6c>Lu^)z=Z1WM={{;Jy*pI?~ zH1=b#AB%lHb}D17{07nB`+h`A^tXr}`aPnL?)kq5t|*sR{4UtV8#VVqGVj5>OVp5ZPuBLdTox#El>|_~83y18| zp2gX|Os0P^GY!_euy;}C!lj+P3uB>IXz7(;zn|L>vadI{v{&t0vM85L_XAR0+^_Ve zm-H#hqW=CRi+fM+Q~P?BEcTq%5PQwg$hu4|9Ak^RdX{#@6A`{0lUlc8XiYu9C0&bR z-e_`ZuxwdSV>W`ea+qJRa->jZt1}K`eXQs7GkSW@=`~BaRk^0&6&ysmqR%x1&g2I>m-ptshw7wZ-Eg zDz7o78ARKt8Usd`QtUH^DQDH3Qr4YQD_@Q~BS6j36kX0)W>MyW6N*q!fQkP^imezJ z8Co|I_j25|I6$;2LvoX3*`+XxavQUW1M22x1hUd*8TDSwCJw0RfRqh255*W&<;ym! zTE2i<)Z0?bYgETEEan>2ri?nYh$x(y*%ZL|a7m=>0ivOfIxhtBQ{q{skOfSv%n^N< z@Q9Lf!N}^Eyw$nHSREiAQXJ4m=0x(3Ow~i?hMgLnS+^7@0UDA#Pwvv8gqY zk?o9v{)DJjgPG!hnSz?xBrD~ys3^j;0QZ2>gbXJk!))RTXb+tKi z!*%DX0d>dpwL5`O+=^aF+j0mK6-tlV+1IjOcLbh8$@ShsBrlbiLQ91c8)P)Zg;3Pa zg-aq_YMf*&tHD}3a;b#g2K4rXP;N?G53`z0Ewc)<5a^;|T1A(Mz;=dE)*g2rRBjhm3m4w(r{Y_BT=h5re95qnOMej zIYwEcQWjMk9o^V9qNo5|h%MBMCMlvY&hiI{E)V;fz=XE}7}m!DN+gkDUH~I}r(voB zP%a_Fi=a#d4|-Tq9XrzyBzluE;^0^gz4+x)7BT!N7^PEeaqrT_OZ!jnUD{H<0hr6d ze${P3rr73El!9hy#tbfH!!(GQiY1reB9_d+O@vScWK6Xlymd#hoN|eeE((l-pLKEL zV7Qu{)kb6Wht%dYA>YOWq}Ws(VqyyGX*XAFQ3~s2*Hs&2xksevA%C2?Dl?}M%wHsF(0mTUbW61-4iU}~#r;rSFTcsRh`GdT3DGOBrRB56~b~)M< zre5ZeVj2R*D^oc=4u=BSZaw4&9%ATG9-UQj65%{-9n6z50EvNh0b6SC5EkDMokRAr znvypS8Y+Z}na5rHDldT^Vr(3B>IG7$N{CgWq}bfT%Vvqh5d*uAfXFiX`JzPbD_~9I zCjp=vwy5Vmh?_63!Q%ysYub>BDU9(oIV;7@mq%ee2S=pXk!C|}Tprdl5aPxGNXi(G zNFnZ(anK`Ov58a20THkWJkAO{6g|-%A2%#hw@oT86y89VM{Rkd5=5)uhM{MilbR5a zB1ssTK%iz?Qw&u3dDSVXRsw4S3dy%|`?t~xh7eS2NqAVcd4wy7JAZF0r-gJ^lngiBYWJmO{9Miiye1cv!VD#f%S#yw{=r&BHK_4`D5h2m1U zNh`_f7)5rev(*jKgqY$yv@z(qn8ScnVr86^MXOOad|lUTPq9@)7mlo6 z?$8MKFq_d-0-@M0-{SmXE)>z|oQ*(agHTWcq40uNH*%8%P*t+y)L236IjcGdCPIAz zHtr*!YbT=Q?js-?7wDmgMzY206e5QvUuYPjtr|6MqPW^b0eOxqJ*IsEtb}sqk_^yq z0#xyjAUO|vTK8fZ!>BLAmZ9F&cE+brv?av>;n@^6<2Jst3?bf(G*1QYaIx2K#b8#V}PJ5g4+M0fFFV` z_2EIgEnVr&>YBv&8Qaj+U2Y?dB-V!Gs`#j?VRsnPpx=$}c5W+K1U@q0S$AZ3G67Q0 zUfpq3;?}0gxGM3Y2f#+X^x~}m7_hp*L+~TT&OtnUPUBrKy1=LiX3oAHI6I zqRMDN${)QJJ!F{u!&>;hJgeF{iS07$cL5hJccq_`gXLFbyu+AfjefN+RytSEsgDtrd1uGttj zLX)HUDq(agBpAlHmY|c+qf}EXdq~A8tObNJ=GwNymJsh^1n55fK zkLrOA;)?A~%*R?F~xVpeLIHxI!k%?eY+T65AcIkJTbhRd_0$`8fB;{4C z9iY^Vf=!}_c?By~wX^8?-W&AL!Ng!%FE+|aaVfj2$x=sG?0ipgKpN((t4gkhwm`~4 z4k9KRpxeXND*ca0p2f>Djq=rXnj0KDIw8f@4y{=+JcuP^mN(Gn;1&V|p-qSo79_=| zhQ$jv0=QU^hc2KuwhEpih9Ty|K7UvmL&dPBAWO{>j9HhnA3`}XccG>jW{58=UgXH-dwv?2IFd28NBPyXh)zD*CK4c!4HVcrOEf9C>mA>2ZzJz!$U zu$UCR+0cg-S8NCcRSE_d@kPpcqc&Rnk^mwUM1t0gvRTFuTNjx3V!MN_gsK(F^83&i zKnq1;vH|t1`%710m3I}3x6!-`p`mF+AuW_f$Ji~^!aA6=0;b(vRB2I4yAZfZb~j&w z4YhP3^qsobVE&3-w&0*-2gh}J4f!km^5P*x{DKAMquN7MHLU7|Bm!gtbc0otirJC? z@=+a|v8F*>=tjJ~j3#&itFY0UJbr} z1^8Yi8q_VvEdfGCTijwt0I4EC|GH=V-GiVf(gA0VR(0c_!>IuBPsAzAn0ob`Wgm== z;{G8;=P(B@3KC5x0<1sxgD4d3%Jf<(Q8&1fLQ&5lV~<7ws?*{^$(lKQnx+~jI2Q+0 z`GVG1A2pmf8K5DNVy6tR8Xg%QSiweyS7Q0!p*8Z*s^xXZc)(@n9!(UZA}uau?g!US z;`LF*OG>`jW})Pba!EHU3$RZfj^aWo&~o+lH}NdDy-S^%Ocut9*CA3|${X!E8tDFl zI1e#*nSb7k2CbkiP))8lV8eD)5-GdwC|S$#De4IQVU5a22uxhcm>8VH@-a#&I^Z{y zka(6uK&+fuyVSV2l&b(SIJpR9`jP?P04T8J5MQ+;$TFvpPi>7LV$q@43ccP%8kYuL z0Rmr23nD=TH;*oRe87V%Xo4(90r>(MSkEG?VUhy!z(A&;7(`vPshdE`oWgJcv6@_q z(dhFq1gy%n7#x=}r;sa;Zq|-!6wfjR2*QF~mAX&QR3p@%4LZ3xi)SO>E1uE+0FM@Ml3`##L0Z}_2!6k9hm zxO!Dfoz*01%2lum9-UNgNwOT@m+m2_*h>1QmaltfV=*gVtNhS08##+Huy@7v^8`vb6KF z*ZlB?uWprgbRGJ{SMFJZ^Y0&a{Gk`(c=jVpd+(Aa?QcA&&-k7+bIL;)uqYr95?jnbw86PSAJ5TKek?)Tru+e{R=)JO+JwN@z?ra zLH>>#d-ATAaDDfzy>FMK$(!drkp8YFP5#}1+aLbbNz&x4sfkmU-zQD{eETOJT6>x_ zvGe-7cb|H*H1XD>y}w>@i!|}_^Kaj2-78JJbIg^~-<70^x3i10>%S;Xo>#nd^QuRs z$@_8(zP@V&=ilEoe(DL*|lrkxU>ELRsBi-BYBys;#P? z%*0ht)S&VdBmqQS6+r|QQ504d7uWR@#JIXFCX*Qkh(dH%{j%~{Wa<6Stw&dPS3hn| zqwq2Jo0_`y=)cZA=iGbFJ@?kF?w;juX=!OW2-_jp4#hSjqw2D*WYw0IDL8DIX=_wkmY3_7P~RJ!RXY==0+ zbb3BNKNbHEck$$$Y662_m|1au2fu$#<^B8u5kJer8B|QVna**#G(BbD@Pl+r<;`e1 z%G+`-w%ORGW1FSqG1jD0Z*L^fopXA$gDD?Wr?4Vfehl(_5e=%zh}tL2Kb_d5Q<{wiRB(P zvMtkW&CbvRk%-5~b|fLHR4Pvek+rlU8+p6UxnUNrABOD^F5|5OhWkUJgElW6jLZz@ zVEZ}P&s5{ZhvS&;CH{Dc=yC6#)=zodpEHdd3DUcl_)Ph!j;?W959;f|8K*qn?=Qmf z!3hJL9OBpLH>9s#0&P=>CZuUe{B!Od<|k_0d+=SjgXp5i6Ay`3UKu>v4w08@nap`g zChzcz?jcM{$kQRE>3_n%Lx@Hn&AI5^Ac!}U@Mcqel1YKH3RS)GvJj8Ij16fP-9IB? z=(?J)iO2L@k{Ri-_w1Qe7iD`m=M=@c=iCj;db)5E65(CTb1k!XV(;NmWgi&gQw6Lr=$kHh+F;dFYv5qWoIt;Lri=XLJ1r6J)3I5)ITwQSxZ0 zBOgyxhkg&gIdo9TSky=n9&b^ca+UKnpikS09@MpzA)9If8tCZ>)f=M0OQ3DY<=tWQ zEW-32K1KSYbDO(ARoe5B;hX+%-)*taT=A|#q42!)Uq7YvqO*VO4C#rHUp|ria9iwp zWzMCA!gJDZQb_4J?!H_4#no4=oBzPP*wv?w;PbAng7gm(3tQ*Z)bUJ*0 z-#(ZB2UpNdeEiVPmn-NcKEC?sV0pR9%i$^Q|tor>*U*iOSX7aREs)PJC75I$iMEZ8h*pmlx|e91f^8x^o9ORLc%BC_Z{~DLNNRcdh^vmh zZP{kY*$+#*{$c*7lK(Ns_k(`>-~HFjoc()g|ISf8{pfAdzMuRu{ioxvk@g?)mi?!& zZc|lkHOb<`Y}J}I>wDQLEq}V<;fD&xO258!-6hZ7_53;ZYQ?4l#_DR3I=z*h;3L)^6mgn{!z3KzfSn8?=A2@r5OXp^3eBZST z4!m-sv}fbC6U_{h{(Z`{=l|UOEd09{-n;$#KbIamYyUrg{!sw3ZX=sj^dwt}?kZbm zrt?`YT+Hb;wBOt(a}*Ro2?7nxI;( zC*!SgC9%aem4tc<0iwj`IUR7;2PZsYIRmU&Zlbe?PWcv&XJ*k?H}oZe8$`2?rnp2d zVZ{V?vIL}sU9xNb!enP6v1lPPbk;GiV?pb@uGWruq0d^W#aDs-Ub%~qogJyJ4z+Xf zf>biT2$1T+MM_6}ai^jzShQ&I!j3aL)z0~g7j~7SuI=q#pQwamY(d-nuC{O@$}vig z*RSkdTMck=+k%pN8#?IDn4g_q9oZI2Mz)*TCxCt9Ss52d*EOHIyDY`dM1T55vSCrYy1(Do* zjPa0WxOd4oCubLF(q#s7zjWtRZpUsFmw!+ z1cVcz>b?G3Y->&Y5`^U)&IAhUWW|eISF!s|MYl9m_aN`9<{^~0VWMf!Wy{E*JxnQD zCP-$Ci8+Y?OeEACN5oJTtkjU2bjQOqTQ#{kMulz|wHz>~B@_&+MZ)52aP?404k~&g ztdS6JU`es5H<6I7go44upjLyK;(&<(HL*oj@(nhL7~i7TRR)p>lWX0$T{9y5y{O(rqER3 z#QFs-;zB5D>%7HvdFS20SW-i4ZOgeFM$Iq^T7+_P;(D0WEb5t+Gz-2i>V}z>x#~Bh zgdmeBbO7CEVFFmZ)Bp_kT)m=;D|XUGJLvtqa!*1zV}Xd;08rB_fxI)~6skJ!Eqegv z`zU}(nR+JWNba$Ck$K>Vj8u5p#VHJ$rH~3u6pGd562MH_Fg1H1TP%N!IST4kt$KUL zf`DkE6POHu>f!eX#H(HqkSw!S?DPmYa;ThvZ%_0hDtH_#!A1z=IyI5?ET?3R2ki5PwmjPfzIu%l~X z*P=5!x*94s08=TnUsZdMF}5|=FnPF0JmANe00TV=siAJFltWs6zwBJf zLMDKWCc0$X(1sbAJdYI9bVzPY<@hih3KYBXkRN!6DS2rOTZKu4^RRI+SIPh+237@Z zsJ}y4JRv%V;w3dDZRj-R4ThP=T|AQ)K@Txw6rIN^^pV<#2F7^{FB>@$M+oc$0a0X( z^Z7c%g)$)uD?2|CD0IRW)!Ya1=F9EyxWVFzK4fGHLwrrnO7Z5)rLdZVBU0^1vzhBc zV=`vAlOx_Z07)6*5h=vI5;n%BOOiN+91sBu-{UO7LroDK@o~a3Rr{pkLg5``dC-zC zQT%8X+%Sxcb5bJ$QX~m86A09ZYl;q&pH}Sw)r`>EfI{+Zobj!=g2@CWDGCqs7LRcG zac9s~Oi_wl3CKM>*S#bAiZE3}-X>RUFNl^l4Y+g_$|GKvEksdDo4{guk&7`co6Bek z&2Clmnd)Pr+@LrVZqahGIz*9O8frso-qp?vgvJ(Vljic2tu(OzJ>Y2TqvT^ISYZv8lj*FLg5v!YUCygpsFN? z`O6|IM!IbT6QMpn8z%_pv=c#cCkTkn1x6@>k!*3BLgdhkOB!A^Q1uEoQCw|;fV{+| z9@B&XOQ9UOL<97?0A>6mNY2CgjeD_-VGfvX$x!d|cE+Pnv?s*@;n`&QSd?YsZKQ;1 zS%y2lE^bqTfEDdE6i6<2I(SF@N{qA7^pS2_^lfS9cp z#{@w~1eO0q0p9>a>V4f-Q-;#(*_z1L8Qav>o^K+JB-Wbqs`#v`K4%(Iw>OOM3~nn) z1U?erS$kl5G6B-aUe$S3;;l`cc~#<74}i5s>BV~iXuzrl55bQZTaJhh8N$BqMnq3M z8r z3slqw6p>n+w)Q5rHg5nE@Ltu_)y(VI3~&=b5zS`Fg;b=z32^vQZ1}fuKDpd#CwhE^ zBEU%#KwK!M80+!g$@lqGD?`6T=#V1VG@Oa47+9Nt0#@E+p(lXYk8Ns0o-BMiBT>i!ONBcjKrhnc zg9I~CY)X~FI?q>L&2 z0|E3XfK4{^@El5-mjHE?jWLX6Q(VFmW1W4|J%>8a6Ti&OsGiScUQcF}=<$m}e4jL< z#D!u)P!}w6ybhmHH8S+ezqA<@ciJb-XzBY`QCe)rE>~&9q5wACjQT@mc#q4Qlk96Y z>@sHRR=mwwRki7dlXlr>AJ0TFqmy>Is*3$}E^xa{p2q9RE|+*dq-&FQxxDK2Lu5S| z1RLaUS=Q4X`t@JhE>CuE>qv!!@+R!wX7X8@o=c^#Oi^q~9lJc&=`b^DRA|we+KEV- z>@UCZ`pYFIqqc}~wdwukH$bGF);(#KktYV%iDFX_!+C$X#|<(j*O~@i63IG2Z6c>s zPIb&Y0g4X#dm*6E_$C304#PEl#fh2@F;OVmz18$YCXbi>v04J=3}SL}dcd%Xv#i7$ zSZ8u_84NX(P7(*?VR2cj7Vr}niUU|aC)pz~f|R|n+bTqg>!A-APNX;>9AXP=KH=3- zzOq9`voN2N7A@hc+R`D4-M9zF7}pX^-f=0_)Y6$t;uKZ_LK#!eYB3OvxRh(oA&62_ zEWS=c)m1S>_~*2T7R}*IVlfath0|0q>lcef({xTHMMK!8nxcip8U@A-NdojyDx9W9 zfxKqcopQ-1{`$2em=Tw9&3PCk6k0E-vZ?ndD&I4PHHT;#na8Ujb(9w`GBvP#rfE5p zb1V)>({QYOJc$Jf=$*{+35eo~?JSssHM;WIIC^yDu5knkGr0jpHVrIa9-Iv7F+iNc zxN6tHi%^$qj*LcW9O8K`lcePg-1d+bH42RTi1RRk1wKqe!^#?7B}N5y-qt6ciu178 zV*5?B(_;G#I6xaW_T(&Kl^|zTBu}!zeH01_uxgEtv5gm8SYBak9Mf=h;05g0el@J z611l0jXYj`ae#UF-I)X}gh@qy4W_!$QqFh@kY>DjOG@HPL@E|lrS~0$VpT*TEnrIf z&=hN6eQH{n+G@|Lv{d@f^@R$z}j?8SvL zKp&K;p2UmX^dWWXGMRXXs|u0gQoh8hqJhp6#CeD{MffLZY0?V%0_Eh21J>+EMUk@8 zkCL?%UpfY%Kdc=c2!V-984Gr#u&kSs3%-a`l#qCli-1_sy7Chs;!-XH#7je2Ak!Cy zcmklnQh?MOEeGV%E+S$vz)*u;?Ig`hgRb!PyAc|JND#qIW5^yK*5V2pp~yi%zHk-R znJTuAQ2}{iATglmL|w3}8$rsP!s5c76{QxV(c@teu&mT#a9qlqLN4o!Q8}woyvP(F zC{~P$oOspq*JT{Qv=F0WMS=(v4(kbL0hJ2s=!AGw1j`mR*GSMh@it!(i*O{=PM zAJIb5o$us$)X`qNfqf0ux2gF7s2E${+r4IWLxa^MX)2Y!IxQWPrzk}ZA4_+UV{8?D z>(bM`_h6ApPvt`Q^fON=q2FE;|MFTG;9&9FYf*}H81<=ya?Xg_6#7-uN6k{TI?Pch#G>GMW28!%F>xXz5@;b9r!Q3nF^V4Y zIpgruLunK*@*600FL~b>=KDjBmUch<*4g_%{F48z;QdD&v~%+voAAwH+J8WLV$FsF zzxZ@3z8yU0(n8@mY2U57w)|~N+`S_GVx@H55idTCZvl@K3I*B!W^j~m9G|5=AU%8E zXCLmqGcG+d_40WGSA0zR%|yO&yyBABbzKihg~CrkZ+3TV&%@F)=S}1r$20tIHUH6R zhNbaO8qeQ#OHvxY`N-Vz|M{Xcvir%FPtLtx8aZvxwZ~k^r12{rKJxS<|5qB{dDi`R zT(MIc`PGs?zgwOvjhueq(yu&tu{6H+f*AYxXOTuaHJI_q=`9 z+{dNyjhFrKt2dk^jVm{v_UUV1lt%U(KI`0X_evw`;te}4+E#r2+1vFk8+J%z+u!oX zE2rNgjbFZNjguGS~Y3x>yO?QyY5A4?C#H>u=?chN~6yV9rySpcSxhV zcV%8UtPj_Bu&wO-xSsm?2d7QN^;=i1f8;r7bkFqF&uz!`+u!+~Pncho#%{a*`dikX zAdT)%PdWU~%cRk#F247%?`)CAzVeQhOMkLW8vF9cy1sPf1JdYkd#Cl>ce^zDwD#C( z&pujw{tdV7{@uU)K^i;fQ#T*+Z|!*AcW(dMbl_kAugBg$^JSb*>}GMCzjNPx@9=}bI#g?f^G@;%b% fs~@{pF?UL%=9.0 once #21449 is resolved vtk_version = '>=9.1' - # VTK diffs via XMLDiff are really fragile + # VTK diffs via XMLDiff are really fragile ... and we can't add an + # exodiff version of this test, because exodiff gets confused + # when a spline control node is at a finite element's centroid. max_parallel = 1 vtk = true xmldiff = '3d_steady_diffusion_iga_out_001.pvtu 3d_steady_diffusion_iga_out_001_0.vtu' From d6cd5009c60c98e0ae2b747b31d74290c8a21462 Mon Sep 17 00:00:00 2001 From: Roy Stogner Date: Tue, 10 Oct 2023 15:49:27 -0500 Subject: [PATCH 06/17] Add a 1D discontinuous IGA test This seems to be enough to reproduce the problems we're having in parallel, and it should be much easier to debug. --- .../1d_discontinuous_iga.i | 52 ++++++++++++++++++ .../gold/1d_discontinuous_iga_out.e | Bin 0 -> 34352 bytes .../file_mesh_generator/test_2edge.e | Bin 0 -> 3172 bytes .../meshgenerators/file_mesh_generator/tests | 14 +++++ 4 files changed, 66 insertions(+) create mode 100644 test/tests/meshgenerators/file_mesh_generator/1d_discontinuous_iga.i create mode 100644 test/tests/meshgenerators/file_mesh_generator/gold/1d_discontinuous_iga_out.e create mode 100644 test/tests/meshgenerators/file_mesh_generator/test_2edge.e diff --git a/test/tests/meshgenerators/file_mesh_generator/1d_discontinuous_iga.i b/test/tests/meshgenerators/file_mesh_generator/1d_discontinuous_iga.i new file mode 100644 index 000000000000..b3fcf6a83548 --- /dev/null +++ b/test/tests/meshgenerators/file_mesh_generator/1d_discontinuous_iga.i @@ -0,0 +1,52 @@ +[Mesh] + [cyl2d_iga] + type = FileMeshGenerator + file = test_2edge.e + discontinuous_spline_extraction = true + [] + allow_renumbering = false +[] + +[Variables] + [u] + order = SECOND # Must match mesh order + family = RATIONAL_BERNSTEIN + [] +[] + +[Kernels] + [time] + type = TimeDerivative + variable = u + block = 0 # Avoid direct calculations on spline nodes + [] + [diff] + type = Diffusion + variable = u + block = 0 # Avoid direct calculations on spline nodes + [] + [null] + type = NullKernel + variable = u + block = 1 # Keep kernel coverage check happy + [] +[] + +[ICs] + [u] + type = FunctionIC + variable = u + function = 'sin(x)' + [] +[] + +[Executioner] + type = Transient + num_steps = 2 + solve_type = NEWTON + dtmin = 1 +[] + +[Outputs] + exodus = true +[] diff --git a/test/tests/meshgenerators/file_mesh_generator/gold/1d_discontinuous_iga_out.e b/test/tests/meshgenerators/file_mesh_generator/gold/1d_discontinuous_iga_out.e new file mode 100644 index 0000000000000000000000000000000000000000..56bb10bb4718f35bcf4b3d6c47f3b17bfb34126c GIT binary patch literal 34352 zcmeHQU5F&tbsoi%_E)iS5Ev(f8!tgJW`1{OwX>w@jCW?`LAx_+&1fBY6IZIMZg-cq zx~g42voneUd5{Q!7(+kDx&pG#;bI-kX>+0sUUzwYmdk_B3;_os1op!lrc;lot;S+ZoUiJ$r_ZJ>TY`yiQoB=^S@C-SjB_&dOt>j$xXYXkGp84y1q0|UhMqCQTi&oh+I1N_n| zgknZ9v0P!JIQK~8L@xaxHsbNGFWV4RX!#xpi~Tq>oStd;akLckd6aZs=lfoXN_x)>)tZmSn`(cq2H$sWjdtoQR2s0s!xh0FRB}+#8WIdO?c<2{TZb2 z1CyHSXuTiCIg%T;b9Db{XoSbDh~z=POFSgKSHFFR+O82D4P{m2JzJYkZ(YPq*a+_f zqBM8zHS850H9o>4*`2}PIJ{GPH$S%Z>=y2lc<&RXxpQB~{=E48Vtwd0)t3MLZ{g&5 z(O$p=LpFAzf$R`%iiQdDQKCBc6Fl?Sf$4^1)Cf;lAX%?bzUFlAl_mNuvMolyEoiIK z6KZ2bgW5pb;qS`SXAEiOp;6AB+9_+doy0%YAzNKfd}87as~g z*jD?R`-@MX9G!6BgaaoWIN`v3zyY%P!Yr*hbKc z2i%K7hM#IcD_iXQ1HKqWelO@R_MB*O;QJw;?*=^VFGf5}P8d|2w>v#H@35*4$b!S> z7Zw)g=6MkK0h^cmc^@6Cz_B?_c>a9P?#y>g>rNDy76-xXmh6tPZG>1gHMkZaK1aXn zj<_ippgZP!R~l4RJt760{PGo+ZKyr7z+CM3`pPBKww_vFX1>Q(mR2s#FE7n6ttpJc z{<_n-0dA?vWRPns-Sri|cKKr0?p(6hm-*@?bER{6%``7wx^#JU<*7BkwtRVYy?};; zn_K&PRvC`5iwn!^3lx)TZhP0aZq@-_UP$799n7n75J2nLK2ix;eZk88n6q6gV#`Zx zdF{&TMf?-&n*_df$gMasJ1!TP=m-yn0yq!2dSz+t%F0ty0Z#hCcm6u*|F$=bBZfXE zuhIT@l+4?IFyo(R%!obb(=j(pi<+_MEoTjou@>8~?Ps|w8hO}vhBr*l>~WFnsiK9U z!Fo~PyGCU8!aPz^Zn5p1TN~Gt6hsw7@f>50N34_6%# zpc6T!YjgvDV2C^sC1xo9CVryDu5avb?d?}RTy?Gr5FOf#+$L3)6xvaSk}K&0y8d0{ z`8s|KH&wCe1DY-%4UhMr7zf<=^j+@s`fz&5)?!J2zKUTL=b$z*T;2dlO5eJZ!YH^V zQceJoBY;cGfK*Ap%9JgEi8TbGg30!1DbMZiR@ANL6k|0&3Q`|1gkP^C52*Tpl$%oI zRP0;q>5b?1_G&6vK0R~uRsM-%p1CdE9xNdbW@2lnFt9Sj&HeC(FvYnO`ARWn#*vRa zFVbJ04Kv|;$cE>H)1c|y{}wyk?Yw|Ys^mDADCX~pt)?i;!WpqXm9&$rP z5a^RRc8pM>!@GXKsa1wNGCEc`m24MJ%yUkI1PaGMcwGqC!G@YPTs z*Bo?g(?D2HQrIEu>MsFEX zQff?Lz!Aiz70qxV6h6Onc@pc=>x|hBtaW4zhZyLstv&E$W>KcU$N~7c2GTD-jh&1a#2TZ0O)vSAazEcLhgy) zvH}?0HvzN_>Y0opCCBfd5gLDRN1 zp}<%k5GLBf_$``8+8@wwdU-(d%ud$nAv%x8Lulk|0V`~n5Q>)ha(dtyr$JhWG(0rG zNO;xAS8iI*#4_d@70QyAvaZ^=7?CfDq5!yrZGy67mx;nE%by}j9#$tGGnlLdK;(uK zf1IM^3N4l-FtTUHmI{Ez5E)z~Wpa4X!y%7CYb;3qB4fnC!!9Ny#6ue?_dFPtTWob@ zeRciPQ!DFJ)f<4_F6>ur53Podim6vsz z9;z%1HFZu(A|XAT9*k59QkK@C8+!=L%bQRe4aG?TE6D?Xi)k=WQAmb*rc$o3{At~V zl(k#{IZbrQUW5(v-B@IbZKvB6&eZ4}A)r9Dr%(BThnPjx!#GrxLF#%C! zjPs>BBZP7+3d^0J1{BJ)MP2xyZ@!WaPm(My^&uNmSP^RqR_dFtgu=Q2N2l45X0s%O z#{9YJAV=Rg0LeTQ87btwtO$e6g-xGA0f>Y}8gUlz&{A|qd}Latwoj@r6wyI84npGv zGtE{Z48zE{AhjVNWs)$FfI$6@V|sA;Js$OW(7@V&LW*tV_*TaZ3oV(fjXn?KP-M8$ zyff&^wkB0B1BwVQ3Gc|iB26__v?&x@CDD3;kC3jUJo;@pL>8rJ0<-#})MCzH=sK1Y z&GXo;9}^W0MN)Xj85;bMGP^X`Dw8xBralitkB5UU+A^`SO3J!!RHmxFb}mFBa!e|yflx%ltIgar0pzBAM4c7n zo^2i>nF#es*?3GqX(#gJ9up9q3ye_YGuiq+g~*|IQyk0==UypH)K{B4pr~;%VmcOJ z5lYIX8KCL{jPQ>j1rL{}k760aq%JQspx&c)Mp3BSlM;Z4YzAU1Dva>5iN(XvPsZ2v zeM%m%WUrezDVe*Fc?>fW*;B#N4)0O2S6{Kot`0D<#4(Y*F9#6%fS7*Q!F)a`g7W{m zfDgcs`u5duCPV3yd`;u)jNM*XjAuwAiFNF}Dlw~STTVl|s)q69;I?Ta@nL~yi*Z&V z0;G|>+Idy_)@H)ID*d($z;UDW`n>=cu-fDy_|an5kkKJWxP3L_RDFqfC}tN?p=xWw zLouVu434VsK%s7bnX$2hNKJT1X|F=GkSDdpLdSdZaEe)TlCSh^G8z!moVy`hGxw~l z`nu+e0O^%Q1N|IKYcRsNivboon`}+=Q1t_vg2!3ellllOn3HZ?RlMtpFH3 zMb>6vafa6B0dN7&@WA8lee4E=3!seVtmHx}GT8-4zGQ~~j^IzXP zLK6xxxntY63WOC&drZ6F}}K3Y`K@M`gd6o@*S(+IDAAKvhxmTtMu`jAhNLc7 zB%_X$QDX}9AG@E(_%lbmWc{4!?hmfq&9 zRc(6k#4o4p6NM&j^u#aMs@OYoLHK2gH0~$AT!?%~*G~NMXw%z=$Tl$ub|BueoTNMS z^xN&1PrA1gq(VY@7j|#mctF#0srKcqiBFl}m!FUwW=2Nki3uX1Nl+*K{#Y#ma|SUvx!3o@Y?hV2ft8b!M_|ZI%1s|oq{R)V zF5#yy6bU#wCs`5LK+0D=S|S!Ar9)JE z`Us3Mp(U8SlThmL;>;y|3d?{{#;#l~2BOiI^4K{9OjBerPP0ox*a7d-I$+ZRV}>LFrYNPLX<8s+WewGc(SYUK`ovRx9%dFhHPO;yrv@CL(>r?#mas~YTop-4 zHotE|ApzF5=ob6E=RY~xV84y=*QXERDx8-FQsvEnZ?Qc+rps2OmhN5KW6cUwmT1Q4Ub}f)b>@9;Otg zR=NOs3_z?kfjQy@G>AqYP=s?b2sG~Y@#TwRY@%-u@h($Z*;HAwf-Z#S+-#@r6j8f0iFyl5Ers73))jAW(=z+(UW6iSbfEY zP~cK9=@+kvVRa0OnDWg}-vx+uj7ZSQXIK(ok$#txptXpkD!&C&-FzEgh>=4AYDr0a zY=&b|ReIl1D_2Dn(gLQmuiRQwEGJD1Nr#IAPK!l{{WMGp7sb-3@Y4M(b`lFv8!O!I z@#Nz+6{?ayMO2M-^-ZKKlU0=!CX21*Mdt`c%(lFXwT4RT-{PWve;I&%td2^bDIxUH ztr3Oek_1fv$t$4`4n-9Dn7l5Wkg5-Ol4;Lt{e4k>3;gnll>AL zwz0y4JiusKUnsWU#n3xs)>|kKgw{rV(A#9 z{;+m*E(NA9Wh~gy#jCL#SQXMkAIy8HnwHDM0cm zqkux%86p+~taRvgCuv?9bVaP+jnoK4f(&jCL-zO%mQc_JRn7y7g{!d6RAwKW0*b_d z)i*t&F5lH{AZ0;ewy80y3 z_AxEQFRe(>fWjlsk_)JmsG}45SrIH-G+QTi5Z$RHEjg)Q;v=Vj$3I}x# zPk`$JJ^)Jx+)}Gf&#ZG+XL=i+7aj1ofAtu}7uk1IM}<{mRE0=2l8!~`TBWPg(!uD~q{`J}=?QX+-JnlXD&6}i7MUDN zBt52tQu*yQ{V%Uo0Ul)EUTaciGYHH&y>8HxdJ7gEqL0Sl16A~m(ekDr+n~qnVAQ7! zCC}*k6nbjrqh{T49p)w!Vo`9sF;Y`!Oq_6a5-kKH`%*QCP4rOA8AqgENuz$1AGFkc z`@6sYe}Ck@@p$`3fBLE0AA0>}`@jBs|FdV#|BtdT?Qh>*ef60)KHmQK*MIh9|Ld=& z`yc)7zdZIg3$LjA+dq|0w7=JVyz`mA_}BLLU;4slUi#2Kwcq^GH+P*M?zi9k>OUW@ zzP6WL|Lfm6^(Vjio%X*ye%T8jeqXwO^O+z2?2kWoOWpsw2_fx!{?a%8^5tK&@4f8& z;m();u6^%YU;E+ze(}}zy>EZxH@^0hH?r%0_S>)j<9~jxeebo;U0{o!&-UN^|M&~v I`HH&#msQ*e6aWAK literal 0 HcmV?d00001 diff --git a/test/tests/meshgenerators/file_mesh_generator/test_2edge.e b/test/tests/meshgenerators/file_mesh_generator/test_2edge.e new file mode 100644 index 0000000000000000000000000000000000000000..1509697bc9aa3e762f2a9720cc3120c845b0862d GIT binary patch literal 3172 zcmeHJJ#W-N5OuyTXWtP5p$J`m0V049-3gK)Q3WA_(yi>f&1Dte+VX8KAL2()@?-b` z6qL|IK}ADN#e1_}W2B^_JVu(E-8auOv*X?M?d(2Zh@yy^K>N^cuC21^M0+@jotbDg zb-IMJmQ{N?Q6n?qbO{=BZi|tgdcrbHmzd5qNsG|L&sCP62#Gmw5=)&^1}SR?@wJ)@ z80fa}8AJHUqe@*Nwm3S#ds$v&N?T{#1<8!D&KR>aKh;wNDKe68L;J*qKTE5c%XP0p zdueIZOjlD=SWKf8=)%+J@*C{xHTdRHUZl?0u__B=ojNP3k(!!I^e>Z3d!vftp7o2?Vr68+~kQvX8xyMS4?hSxW*TQIr;Rr z1sH3I-{Yq{EUeW-$6k}mJzjdy{`;Cc1xvnK;Mlj7@l%~~}2XV&hG2d_|M~>Jn&PUo#HBn=a z3Zgy+`CLC^vkqwT5%z|hi`Jn9Cb41L0*oaj+`}6<`~KwP7O1wv6(G+Esr4XpzC)>> zVd-mR@a2GVcko^8{@`ch`ak%G4si3&v-qoQemmf===oiLxJL#*UhY2}e0W>iUjgSG zd8W Date: Wed, 11 Oct 2023 16:41:30 -0500 Subject: [PATCH 07/17] Add a discontinuous IGA test with exact solution Of course *this* one works fine in parallel. But maybe it'll be useful for debugging any future regressions more quickly. --- .../exact_discontinuous_iga.i | 51 ++++++++++++++++++ .../gold/exact_discontinuous_iga_out.e | Bin 0 -> 34852 bytes .../file_mesh_generator/test_quadratic.e | Bin 0 -> 4032 bytes .../meshgenerators/file_mesh_generator/tests | 14 +++++ 4 files changed, 65 insertions(+) create mode 100644 test/tests/meshgenerators/file_mesh_generator/exact_discontinuous_iga.i create mode 100644 test/tests/meshgenerators/file_mesh_generator/gold/exact_discontinuous_iga_out.e create mode 100644 test/tests/meshgenerators/file_mesh_generator/test_quadratic.e diff --git a/test/tests/meshgenerators/file_mesh_generator/exact_discontinuous_iga.i b/test/tests/meshgenerators/file_mesh_generator/exact_discontinuous_iga.i new file mode 100644 index 000000000000..7a551f002ec0 --- /dev/null +++ b/test/tests/meshgenerators/file_mesh_generator/exact_discontinuous_iga.i @@ -0,0 +1,51 @@ +[Mesh] + [cyl2d_iga] + type = FileMeshGenerator + file = test_quadratic.e + discontinuous_spline_extraction = true + [] + allow_renumbering = false +[] + +[Variables] + [u] + order = SECOND # Must match mesh order + family = RATIONAL_BERNSTEIN + [] +[] + +[Kernels] + [diff] + type = Diffusion + variable = u + block = 0 # Avoid direct calculations on spline nodes + [] + [null] + type = NullKernel + variable = u + block = 1 # Keep kernel coverage check happy + [] +[] + +[BCs] + [left] + type = DirichletBC + variable = u + boundary = left + value = 1.0 + [] + [right] + type = DirichletBC + variable = u + boundary = right + value = 2.0 + [] +[] + +[Executioner] + type = Steady +[] + +[Outputs] + exodus = true +[] diff --git a/test/tests/meshgenerators/file_mesh_generator/gold/exact_discontinuous_iga_out.e b/test/tests/meshgenerators/file_mesh_generator/gold/exact_discontinuous_iga_out.e new file mode 100644 index 0000000000000000000000000000000000000000..b19cbda6a59a69737b3ec00245a6909ea408194c GIT binary patch literal 34852 zcmeHQJ&YVjb{_do5|^S*XX=6j3>XE=z}YaC%Uv!hlE6qj-t5* zF1xy>`s-Kkz509gs`~D|Pkyi0>+$a*{@%ymg-}LLTf3wwz;`n@SZ=F^ywY44c#MJn5T zT8UTYJ8=^{rt&Gi1ko_YEqSB(9R?Bj-TMjtF5&NO{JrB&1Ls60X%I)aRXlv-;EjL% zC4X;8|KTu<-7JVk&NL39%y}LsKFaeovA>Hy(tZJdTf-m(Z}ktg-}w9c6ZbTjP@d2M z2g)$Y#z-!o7l>a1_>Fxig`4FRo54o3_huGkAzz3O4pu`;{t&6@Xv&|5Fz~tzxZqXwMH_#4KN zf11QoBU|PrzxmvEwazP=zjJgCDY{(d`Cs9FD@qsMN19BizwhsncmetA?_=I|-}ev5 zyOl1!kGf$({eAz4#0$t@e;@O%`@Vld-nYCM8gr`P$rIKU6G;1y_05HU!Byvdp;DIQ zWm8b({d#3y;!{Qx1@(&hW8`VdanboC2*c0aFxSe%JbXU`;=EY8I>sPI9_{HDick~a#*h!>^*2L7~-$(jo{Fwlp z^sn^I5BM%(nwOinxte#@@71Fp|5{yoA5&Vj@bSOvoAUQx7x({JT<7xZ)bsz(IdIN_ za}K;b4uo=;y}Y--{5L1TX#83;;N^?poPu)>yi5*ojK}dFCpa7za{kBhCFfgap2Bpd z@i(R`dlRnpXX|l;Qy)L_9ZcWK@!#tlK4=JM;>QQ#ZX(@G`eN`xJb3iz>7n@GY=9f5 z@B`2oPbO~Ui~BIWiGyU6Ph^y(0otJoRG6S-}Kx0NgSqsUWZK4R7tq>o>0r{lN|Y)~?*Y;qDD?UU%JVH*VbA z-}~^oyuN#L|JDo|ZXO>VKlK*jgt&HP_tuqWBIYq>?w{U2eB1_j^GabGnR#Nx&D}gn zpmpLHCibelVC7$+@;ve~vAZjFZ`|I!d3*o5V&4k9-r=eA^2{BCQenbCJDDor<6LUM zs|39I;qBf1wE$PW@a6wbnq5KbNiU**z)WgXR6SrV%X zC`4Zz`2G_as!E=YgXse|az`?0qeGy=Mp+VvPUeo%I#Lnui=&61AKWj@;%zWhx)4D` z+D>Ex7?b+qPomiI<`$u)+hPO^vcL_UVG>UqwSK|H0DYB-}&nlE!%pNd|pJP315a;No8iGBb8v@da-MQPywb?|{;ttX{%rl7<7_3lY0pH88)K_$? zU3jRPe>ZtiD`*w%3pLcSMX_57tKs|NN9hX;ohE`~ zs2|O0al;*-nqn#|o9W|l&R!^ODyT^8B9v5O*3=gI!fwJxoII8ETT#mGI~Z(^Yv^47 zRhv-5n2C1bAx5l0YhCFJJJS}aX;EoSXkS1F0zY@d1|4ZgVy|~)XH#8X!wb%wr21vG z=&&Iw$*3;bo=Lz^UldiA5mh0Sai9;x!C#NtU zMwmXH!tQ)3{ld9E$Yi1n>_A<6ZZQuy4Llb<`kk;8&4t>y{)roTO=ds)7Mps=^+C@` zJr_F@>PiS*vI?LceI78(;;^jCw(BW{VI1Zqr-Cs}gQ_WV z`}*2y*aL;ybU)S0Cl%6JDNH0*k;;NrxDblmzIwBzqAezKK#aC9bEZ>mORJ!C%%av6 zH3f5dPL^kCajA<@oJ<^VEWKxIDIp+9p$05!@folR3YS{Gy0nM&2CG4Llu)eHW(F`E z#YvEjCzWXWLi1u#p3Nyj)B({Nc=TUSCgs|Xy_D;KY_o%Efs1t7Hx4s^*m|~6z=rB# zUsHk+ct(@6wg#mR80oMXHx0O5n=M4>G$>15%Bag$Y7?sNbB=Q}fL9Tkb)+Gc6gD+7 zC8Flx^@c~AlK@yX{Uu5*(ig>Ib@THJsWX7i6s|g0fUn^}52x7M=q)_hmLLR``E-bl zIr-FwdsQ8b@xIvKyS0Do#)o^i)+#puheKcO@|lOy_kHnv+D1)Nx)6Gh6vYvl=_bw; zgNJEc4j>>D3mLP4jkF^A!AYaIDdj)2(!&CfZUaq-_N`3;*rtq14{a%ID=x*NirnHs z)gl!5Ot4+$D6?!Z56v5TQ~Rk_EirR**#SgcV~zx8eD8=0sx!Z(bhD9MKJ)Pc`b)uGkB;e zwyth!9^ss%&X;aE2geWs9=rt)ZR>J5F@eJ~LE}N-M(F=WG8@aJL*oFz%uO=whc|H3 znTp58!M?~`Y^@u=Kvxy)^)Qi(fNBg{3{BYCU^c(jNJ;f6#?~-O;+!+mVyJCsw6A*U zvta7TQ<|6P!>Om9297+mm2w)%bTY)YTb8^@D7K0XPZoHIw3^g96d6tq8PBJvidm1@ z1JZq=jxZOu=OO$=q7<#+M49>045ul0GEL*c8PQ&`>wpVN)1{{9!M~mexS3;|U|~xe z(Vxe@$kz@N(`lHs@5Wl*^#_`z_vk);aNe63ioQ0{`7FM5@N zE`!e$p9*g6ke;#^qcMKUMNhp?eCwFvgzs1UEEyWsEEyAOst4kzZt*>eAi zc=ttl7=D|&QPNX#l=b1f8#T>>zV^jKut<+ON=7Y-@TuQ*H)=bRI`2m3Hx_moAAbA2 z`%2piHYl34{!+WOscK=Ha;s|d!4+=&=lyak_w*0s;WyyM|629SufU9+2M27blc(aW zm@U)HqP_d>ne8EulNhcg6Fqq*IkSqnqMPqx2r@fU8n;~$(02J_+5Co@hV4)r(wnFl;k)um@B-p`D;@ACkRLG~oQ&mLaf%l+Av^ z79#DZW+h-Vk@kQX5T{v^!}UvTjShnp3k}@$uTIV6EsDLqaa;(cB{1j~lm>Ek`O%)j zML;OyP#?Dd(b!9Q>Au!3H6>ME#BICYt@x+gL$zArV>i88tVkr5qj1nO7cSf(=lF>n z@*tgc1%vrKmp@9Bg*Ab7flMAqzxKZ66kZvrYbcYMHX*wb0^~WJf$`MXUNQ2q&N3G2Eb(wAOJ%7dwiK?Anf=ddjoy0^?!m8i=~85Q8vFR}kG9~OGuiMJA&FB7>4;}#t0ILw z0g%_m&|l*`R-VL~?o1?4TTRF1Fe%+pXA7d29#>;0b%Lw&Rqp!S@#&e*7B%QgkwmJN zcJ*DPtc_KT1@1C=ESA)tV*}U10g`x${d^gKV;nEYZ+|oq*`>%!jRai)*(+fV4nkDW9FUUp&%Vd^Rw4)x@kAzp_u^ zsq1m1U#v&31rd$K!Psj$mDsb;wBM19wtb6T17t^l_4S#)%5o&q!(fhE#?|K*XMoI& zJ%wHz&hJSvV7%IozG$7p7ZvLu^>`=1<%dPv?*1rdDfMCSEQMAUf*ypDp40*6_6(vYLf)g5#9v&SZ9o!ekM-OnW_2FaZ@Zr7IqZ)a}B}X)M zjG4H-l+`$RaT|$2wXZ3o$%ch8it=d~&sJbV9%#XO*rc7=#7+l?K1 zTO$3%p)2bA9hgQS5?F9YShB~5WR!w-C~_UpY=A6<+@TxQ6e`1jH+CaZSD)&3kg}q% zIxS(L)OLk7JgflcN^J+Hq^u}ZvJT_LO?&o5rh%ZU7z;7`s;73(&f6^H382zNF61)dilffMb910UAmJm7K4Dr#A4yhcm9UDY-M%s z@tL7J3BuZ9Bh$*z7f%k3j~+cdxbNIKeDV;3_M?Z3bX6&QaOnuqzO^BsK4cF?*QZF9 z`+&z82)8h z*%I5MJ}OJ0u$0Px_-^F`YG^mLpBj(&Vk&dQ)DsQNhUfkZD(S=3ZB@GVAjt+orll%R zq0+cj+RmcnNEQf|4QS8&XnT#HMR|n{A)DZKG*x0I^Fn1QhTaNz;emJd6|(sk{rV`$ zqEe+}G0jw3isLMgk&`pz-&>|V*artAl^{=M!kOdz+9Wzr6Rel#CUAh`@&M}Eegp3z zoZsEXK}O{fQ*C|Y)my>F69xn^pZ>%Bm(Mn~Zr{Os8+;So(sgirkX1m&C0k;?vCiHO zWFQCUfiNbHxv<38;KnORAJ}Q0P3`#*$9&Re+0>z4)8dmoXw154`x@^-QT!2y>#zhI z*VK4J39K7O*9Y`h#i8w6U>d^=r*eWfgJ*qPqMvm@-a#6iE+OsrGTi$bccxu1{$;pN z&@(qo0ndx);^q)T>S*5hh6nfOV%)le7gOJ!zW0hGr#>Bf=+Cay3hd=S~hgq7+Vaf9+7CH9tJ^-NAfjrT@M#OTQOP#M}=oy4Q_>Nlc=y=*+QTGG(Jxj!aySE^H)0eKq@%!*A ze|0ah&WtbZZ`<*A749Bzq0@`(H^E$xnfE23VhFt_?vK2p&yrfR-!}6xoaw$r-3^Xo znHG{-m5wrMjSuFBbw%HH|M&6U(Y?fJbkGso1-k?tm(bn8V-5Y%bW3#3eABi~ZAVZC zY>734FRJ$#l5-?V%Z#O||D}>YjQsz+UIr35Dz9hpZ=rKFKA7hO;3J!BGS~6^rojBn znBV>Pt9#>j$vg2_&UJI(UCVpc-u;yB{cFekeL$Q$r|;(rzbkOv4-cz Date: Sat, 28 Oct 2023 09:33:28 -0500 Subject: [PATCH 08/17] Add a discontinuous-IGA case with an L2 problem *This* one works in parallel! --- .../2d_discontinuous_iga_l2.i | 53 ++++++++++++++++++ .../gold/2d_discontinuous_iga_l2_out.e | Bin 0 -> 36132 bytes .../meshgenerators/file_mesh_generator/tests | 16 +++++- 3 files changed, 68 insertions(+), 1 deletion(-) create mode 100644 test/tests/meshgenerators/file_mesh_generator/2d_discontinuous_iga_l2.i create mode 100644 test/tests/meshgenerators/file_mesh_generator/gold/2d_discontinuous_iga_l2_out.e diff --git a/test/tests/meshgenerators/file_mesh_generator/2d_discontinuous_iga_l2.i b/test/tests/meshgenerators/file_mesh_generator/2d_discontinuous_iga_l2.i new file mode 100644 index 000000000000..fc23b6a6c968 --- /dev/null +++ b/test/tests/meshgenerators/file_mesh_generator/2d_discontinuous_iga_l2.i @@ -0,0 +1,53 @@ +[Mesh] + [cyl2d_iga] + type = FileMeshGenerator + file = PressurizedCyl_Patch6_4Elem.e + discontinuous_spline_extraction = true + [] + allow_renumbering = false +[] + +[Variables] + [u] + order = SECOND # Must match mesh order + family = RATIONAL_BERNSTEIN + [] +[] + +[Kernels] + [time] + type = TimeDerivative + variable = u + block = 0 # Avoid direct calculations on spline nodes + [] + [rxn] + type = Reaction + variable = u + rate = -0.1 + block = 0 # Avoid direct calculations on spline nodes + [] + [null] + type = NullKernel + variable = u + block = 1 # Keep kernel coverage check happy + [] +[] + +[ICs] + [u] + type = FunctionIC + variable = u + function = '1.0' + [] +[] + +[Executioner] + type = Transient + num_steps = 2 + solve_type = NEWTON + dtmin = 1 +[] + +[Outputs] + exodus = true +[] diff --git a/test/tests/meshgenerators/file_mesh_generator/gold/2d_discontinuous_iga_l2_out.e b/test/tests/meshgenerators/file_mesh_generator/gold/2d_discontinuous_iga_l2_out.e new file mode 100644 index 0000000000000000000000000000000000000000..7b1f934b1aa49f9510a0daed5ab492639cb5fac4 GIT binary patch literal 36132 zcmeHQdyFKGc*`|f;xKnys+i3{xDfV0ngTOLg|w|im5y?nQJurH?5)7>-EzMk%J zchBCg1u=;tVj>iS5=4vuS&B)|ufDIUzOTNju0FYP>bquUX5N5!0P!ZoeV*eRzS(nTX7=E8 zW`FGV93zUIKAnkhgFomQw%emK%DIds!)hzG>E(ef_yn<(GS z@;}$7zZd!c*WvdTiD&ek5E%Tz{=)rB{Qg^u_wx%x{EZ4{!wl(W`j6A)`t1OxKTYRU zUeEPW-kApw4ciG}0R?{=dEL1!nB}9itz(er#L` zLfeSkCkUI$#h;IKf7|uIQ*j}4jjq43XKuPZ;wNNafVfWF#otMOhU0?(KYAEKG2;QT zTw&vIW^e4q9z76&c>LBS8={J=zz1Q2U=SH@$27bpBN)VU^1d5!-_3{zxU8>q0|`L% zP*mZ19PAS(PrU=_ey5cZfpfZ-_@n~SCGX$+UOLA8{V?d0K-Im(Z#uqC)HzP`8`f6v zk0_VQZ{z%WD*#6~`8CQz^;QMYw&z)Kve!c@?%m@5jJWp(dOp!bk0)Ld&r}%{Z3m1~ zuID{udIP~Qx`!|+;FwM+?E78(93UDKG;vYgpow>Y@D5VDGN>Z=OlpGpvJj8|8WG_j z-MlZmv-eY*G~%tkls&)afd`jYPT(d4!uth&G;`BWAys(PsOVmj z2erjfcsE=;|6{8cR&kfWdnZ4dxta1DZa$6lAa4g?f*~3)(Lj9_ z1x3Rc`6yA{{Hw@w(}w9q)L9W8?^#5-@_fzIr>CetLz`;{M(9|mfvTQRyCE7>0BxI+ zcel{92va?LPxG7NHksbt{OSj{p84OGexvbAPrYq>d;43>fBh1LS48?(_cp)s;@^J7 z{=i)0W9F?-Zf}31`EnbfbE`~GHvi_MPhD7iZlUqfyI(|pqxsK=C>)pRBh4>=`SW)! z{LXgc=lCsE?ZT{@noB#3*!mnH>(__tldh_iMY=8X2 zO_eWCzsi@VTjl%vBXhTX@hzWe{Njn{o7>x8s>aWi^J|~}?Wdnyxl}0Mi}&39#P&}< zlF*gMSNZaMQu%WHx&IHqrS0wiU7i-NHeEa#| zll=Ry`#(vSx1Zs1S$N@dLQY5CB;|G6{lCete~Zgy{@{B=dG@_emg`92`NZea z;=c>!BR;?S#f>Wbrp%Y;Q&Z*R^*QB_YGk|VE(dlw@Vevx=~B}FWcQJsKz5q4Z>iG@c>5_7=oj5x*Ppdl27 zC;Z)M{vr7cK+Nko+*J~0hHc21oJ}|S{4RDfbj;YXS!;`(Id|^jDx-=+p%RDNK-=57 z>(Z7jAy5G3ZQUZizV4g#maoOyw6NW*k(%AsM zfY7xa{NV?)9lJH#GOf#TXj%>x=ON)iU~3pX$U< zaPjo&`6|HUbJW3UhLeLZM6F}%aDAoa1uIX5E<0z%Z0Q(VT)J=Z#C?lLIs3-puWmZl zAU0c`!!hBw>Ge5q5pZE5z;O@!&R?g>zvlM`F@x{KDXaS5X&Ns*I%NtNGX}o<_`oqt zi<&X_0;LU*u?AbV?F){_8+p`q`)5qw>^NZ+9SjZDkHf$-VzU!vk>qiMt(|>n`LuMG zRl%hHGv@k4+lEsDMt?Qf+x@_>owhmfVlpMh1Ojwp*Yu2b81xM8*CZz*`HT4@8tnA) z`s&5?F*wQPGC(+LQgZ7QSx{)lDN5!@AJ7Y~z%Shj22fKaK_Afc0BQJ62Nh$(F&@9- zxScMPo(v5}hT7x-)n)*d*Mn0%ur4Ps3a*ZnB_MoX&}S)-I?^vPc>&SI8XQr;ge>KQW!aez4qxOIIS-9+GH_9i2o8cpx&@3EtEt;aA!AUKmL ztT@y!^4N;K9-4mSqPb^zUwI6nw1N;_gKtEEhyJi_x?Yyd2GcVV0ccrHpCe*+4O(i` zvE>w<8#^Ia$7s+Kqviqou4Q7TuND?(gRF<9y9Gd;WRWo#fqCc{%7Va{%(VwFC0b592pwvbkrNv&E1F0)WGo3l zSzuJ4AztcNY_h&!F7fGMBXZ2p>Ka^EqxyJe4+9(nW<8Nd9+j~{m~(kRFAh}gt|1hm z=*?DSq%c{bEIPACCa%{riRI82)?11sk7`pGIxu3Bf~L3-iZi=#d@S$O4UE|?dh6Ke z_c4Ko3E31%9O>6X%#PE@Y^YmE>Y^WnJ%d;Mi7Fw;Bnkzf>=t%_#Y+vqa?4c}UB6=6 zG5SH3y_k|v?meVaZ2;){MIiN$K823YN6QKzpVk4i4H}u`k>t64k+~a)lvH`y^(oA% zrIH3s6za9*AwaJagl^pJrR8rhQ9)JH8kaJb1wnGYVO%I$-Iv<|&m3BO;8G2b3@{Q}H4c=T*6U&!^NfA)C=B(Dzv8hyBL;greXp&s^~ zII<>$BNad#+-zf^&7g0?<(>tj;|4ptbmH)dqsNv`OjK?Fw%h2xDo2nFcBMZtdAN}= zpG(;fd}5{<8hx0E{U9+DAruW63qZIv!+6l=iv^QTK^0xUbMQsYw%v`=`a{*`UM<^3 z0+Mg446!zaRkfQ}Y+VYgiYuxOioB@P^iX8pQPbe0AQGyFlY>$9LbVaSLwn%EEH4IL zZ8R!Q0$55O@Ec5nfr>)vsHZCBlAb>)JD0MS3LvG4A=wk?!vb%>J;k=$ZgXX7w6Ubhw#u+bZ2~oTc&bMsxK5iKsGia<1sVwR>2Iz%s3~tCLno|u#kX2gO+RhQ28Av z?mA(O-UbwsZ6oHlT4prVlBu`R=iwmY9ZsAGuj`Shm!2ov7DFTP9YPNm*Bo3ioxT)!txdSKqUK?yM@(xRkL79E#}$leC!GqZ&eq1mDW+ zVJ;M1?_2~&C*kald;4qYsOFcj}1jzA^AUO{gC%0l5!-_6H zGElwq{fweecO=CDVcGQfTvQa}#VE^(qCn2C>+6&(V4=URqojoILijPvn0rq-OIuEd zoW1%LoBGur7M8dsy!S-|LLU%o?^;;TCqz*E)djo;rqtI?MpKzeAE#?-UuW#n-27mQ zG?G|HuB+m!s@B9Zq?2kIUrcVBHXI)ocs8F|O-z6^vsbyUO24-mv#v_NEdy}WEWQ3H z03BGR^AP-Kuv75pkRn_=If>|rhkSJr9h8MSJmf2?%y6p=4;1Q-ml+#5iByM&RRD1z4H%7rMS@)kqedK1PPh)qgH6+!Sc-w|63Aht z7Ax23vB&5X{oaP0&PHTn6JnwF;GV#2uxJYd zA==qwqCmN(ua_(8G$*l+lp_akhE;_+rU;GpRBN#fyc5vPpw|JMvgH6%WT?N&P)COi zhRJbqbW9{;nhTh+DKg=6C%=vB<#5xVmRd%UcdnN~Pg(*3+||p(g)}tS1qTyB7@DPg z8g)cDzukIS*Y`zhBzN_4Dr4$xY7B3$dcA!0t5GceuXG%bC;xBtGFckaQq07#osO4>EFS{g)yu?%G&I;cM7DxiwGAxS*cfXjQqEnyobr#3 z>hy=Z|uh>+d!^0hZ!9#R>j zMU1OcA1_}6u6x=@Wm-y}7(7PJQwJ+vh+dIiOZ$8n;H{IM@RAcb!*Uv<<_S=Dshe^E zl{!r8sDP<&VKLHCMHK3Gu&!@;1MD7v?H8~J5v!d$-5^SrZRvNgVzqM)hPp`+^Z~h_ z-EgZMe)@V@0Ol7tD*|gMdNmnbL8N{?OaRM?)CYt?9L3=P@1pXpB0M+30#Vuph40u) zrzrO1o)N})Ex~FaNvZ1$FNM;lun36C*cRK*Ks5SN9=Rx?PD_dS#tW6(Z0PV$^oJ2G z3bnDj2szBKloQK{u^~1cn{ApC51lq`EH)`H#`6T2pcIUzNr6;1E7x%HmBdNw$R6lR zdE~+l5(;g=RP5>%MfsF5Y+^*~(iE?<)X}kik*R~_OIPzy;#?n)*7Zb#mW>?_=q1nm zibnm4Eq36+reA|z3nRL`Xr>cr;t_Z?+VaS$);gU-Yv6h| z?T;Y)L)s-NGw-9%Lkl~Ggn^5#I#eS@0~W9U6HoPdnD*F77cF}1qyq=&NwaW z$Sl0>Qiq2f6FW#wv*{`84fQa_4j}a1hAb!br7XVcKpwi%jUDs>>3t^8G1~>&1f3Ad zU(4|lb zXKE>*^tZOWC9H(*m{V{*P}#;$E|9(XZq=u!-iG;vlC*#XxI4K(Tu4J|v?okV8CG&c zIY-*C`V|{Ofl9$DV!U&P?Ka3_$~HfJ6~H$}B03eCT3DJ&-Yq|mhPsu z)BJh7s2UAI7%2A1EzdYl_R3M?@G3lh>RV4km4~Prsp{)USvado3rsz>mKL2M7(S(b z1)CQYHq^yM{qZsYyV!1(zUM>e`B00(QBHz7fTWc$21iBY^-;cmEJ@V|obt-F*P*9< zsM$>8HfE9U2^@Y-<@6A*P{UT3&N}YyGAzlirAW~n-k3|(VP&`4OhuEZqe`c2! ztzayWCs!YE6APmmv)MX-33b(dext_E)8{sZwL<82t)!8ZU`$nVdJV)!_)=jjsdG{`b1rJs9QtIoWgXc&_b#8XjD8*0du9+gX2=> z6mnU6L2*^1ev!#RkT#6Mk$%%lHhS!0Sx8XW;-ChFTfQZBWGPTbf9QKfun*F7mDEOj zxs){dr+$$)ZFAay**(_^*fl8ZRVh3Ut_ye#xLF=j+gML+a8}jz*8HCjz`y<5(I~!% zv7pVl5+Gfx7GLwxQNd4fA!sYd~b7z-N8}C}Za27`U z+Q>cVDBnm6QNOn#K=D;B6kU9HM^Q(qeh2#+Y?w3hi&71CarNZ6vnxu9<5DVBvfVA6 z#z^Wl z*u{vxLW3_}(ML;*LEpB)h}p%gPZ3J|qpMRW*VJdt+M@={btuFx<9K7FqRd!0;psS9 z2u6&hDiG_`L%wDlmU=0T`bEBmT=&U0d%r2VldJqZ`4(^`og=?-jBg|7zf~;CJt=Pf z8^`%?DHp$uJesfSxhk#B|L-r{@%RgOJkk{Dhh%EY^l_O!A=8#j)qSQsk7fD+nLaAh zr)B!EOgl1tN~S9^RnJ#ADpk+#%ljgkUdf%S`&7Qm^1Aw8 Date: Sat, 28 Oct 2023 09:35:42 -0500 Subject: [PATCH 09/17] Put test_2edge.e coords in *projective* space This is what I get for trying to write IGA files by hand. --- .../file_mesh_generator/test_2edge.e | Bin 3172 -> 3172 bytes 1 file changed, 0 insertions(+), 0 deletions(-) diff --git a/test/tests/meshgenerators/file_mesh_generator/test_2edge.e b/test/tests/meshgenerators/file_mesh_generator/test_2edge.e index 1509697bc9aa3e762f2a9720cc3120c845b0862d..c2ae2e9b5e03b8cf468ee45a7d88083c8dcb5b3f 100644 GIT binary patch delta 43 wcmaDN@kC;SA}jwT(=H{g+8g$d!1Vde2CU&MoV)E`2lyM?@7cY10_PPb0DHg@V*mgE delta 43 tcmaDN@kC;SA}jyZP!N!`e+s7kHyf~qvvAhfzYg#>w%_w>^90T-OaNi`5NrSd From d86958d48541b5bd59975cb1fff619379bac7d4c Mon Sep 17 00:00:00 2001 From: Roy Stogner Date: Sat, 28 Oct 2023 09:36:47 -0500 Subject: [PATCH 10/17] Better IGA tests for exposing bugs --- .../2d_discontinuous_iga.i | 1 - .../exact_discontinuous_iga.i | 33 ++++++++++--------- .../meshgenerators/file_mesh_generator/tests | 1 - 3 files changed, 17 insertions(+), 18 deletions(-) diff --git a/test/tests/meshgenerators/file_mesh_generator/2d_discontinuous_iga.i b/test/tests/meshgenerators/file_mesh_generator/2d_discontinuous_iga.i index de02a4ba24f8..60620c017cdd 100644 --- a/test/tests/meshgenerators/file_mesh_generator/2d_discontinuous_iga.i +++ b/test/tests/meshgenerators/file_mesh_generator/2d_discontinuous_iga.i @@ -5,7 +5,6 @@ discontinuous_spline_extraction = true [] allow_renumbering = false - parallel_type = replicated [] [Variables] diff --git a/test/tests/meshgenerators/file_mesh_generator/exact_discontinuous_iga.i b/test/tests/meshgenerators/file_mesh_generator/exact_discontinuous_iga.i index 7a551f002ec0..d7d67edd7c22 100644 --- a/test/tests/meshgenerators/file_mesh_generator/exact_discontinuous_iga.i +++ b/test/tests/meshgenerators/file_mesh_generator/exact_discontinuous_iga.i @@ -14,12 +14,25 @@ [] [] +[ICs] + [u] + type = FunctionIC + variable = u + function = 'x' + [] +[] + [Kernels] [diff] type = Diffusion variable = u block = 0 # Avoid direct calculations on spline nodes [] + [./time] + type = TimeDerivative + variable = u + block = 0 + [../] [null] type = NullKernel variable = u @@ -27,23 +40,11 @@ [] [] -[BCs] - [left] - type = DirichletBC - variable = u - boundary = left - value = 1.0 - [] - [right] - type = DirichletBC - variable = u - boundary = right - value = 2.0 - [] -[] - [Executioner] - type = Steady + type = Transient + num_steps = 1 + solve_type = NEWTON + dt = 1 [] [Outputs] diff --git a/test/tests/meshgenerators/file_mesh_generator/tests b/test/tests/meshgenerators/file_mesh_generator/tests index 5ca58127bd28..f296c20ee432 100644 --- a/test/tests/meshgenerators/file_mesh_generator/tests +++ b/test/tests/meshgenerators/file_mesh_generator/tests @@ -172,7 +172,6 @@ expect_out = 'Solve Converged' recover = false exodus_version = '>=8.0' - max_parallel = 1 exodiff = '2d_discontinuous_iga_out.e' requirement = 'The system shall have the ability to find continuous PDE solutions on discontinuous Bezier-Bernstein quads from IsoGeometric Analysis meshes.' From f9770d37e45dc370df5b2bfce994ec7848bc2d08 Mon Sep 17 00:00:00 2001 From: Roy Stogner Date: Sat, 28 Oct 2023 09:37:00 -0500 Subject: [PATCH 11/17] Updated gold outputs for Exodus IGA test changes --- .../gold/1d_discontinuous_iga_out.e | Bin 34352 -> 34352 bytes .../gold/2d_diffusion_iga_nosplines_out.e | Bin 32484 -> 32484 bytes .../gold/2d_diffusion_iga_out.e | Bin 35640 -> 35640 bytes .../gold/2d_discontinuous_iga_out.e | Bin 35480 -> 36044 bytes .../gold/exact_discontinuous_iga_out.e | Bin 34852 -> 36064 bytes 5 files changed, 0 insertions(+), 0 deletions(-) diff --git a/test/tests/meshgenerators/file_mesh_generator/gold/1d_discontinuous_iga_out.e b/test/tests/meshgenerators/file_mesh_generator/gold/1d_discontinuous_iga_out.e index 56bb10bb4718f35bcf4b3d6c47f3b17bfb34126c..f527d5937e21b8818cc98a44c2ecf57d3f9f8f1e 100644 GIT binary patch delta 521 zcmdnc!?Xd2=CV3o4F!Q1`=?+!$o>NZ3b4Nl7LT!i01<}MaP^ZFbLUPr$n7JD1^DimFf!IPG|)9NpUlW7I=O^jQX@4b+0ek$!ptlwG1)BD z)Z8QmC}3e^Y-*I2oMK>Xl$4TaXqjY?I(aMqX$iNYOa=es5(Oi31w$h%LrW_|bFi5} z(-bD}5Rk;8zyP8EGoAo`MekzmbtFUA)Og$|4>~Gf}NS@Pr*Zz)vK8w8eSNr>pXK!s- zz0Lmqy(IJf8^RFsF%fU>P35n!zq?~ze){XHaQZ%fTl>C0b6^@Apf0h$+rE7AVvf1? zcc+A$Qs88|yec_8m56ES9o|(Qx&X6?5lKHpuNG zhz0ogHn(#hWEQhDGD}QJHn%iX$j?(SGB7gMH8jvQnCz%8Ho1gfQp40J&C)2<)WSF| zCCMzs(lFJ;EYT>*!Zb0_*v!nr#K_z<)zr|!#B%ah{?ig6rKt-3$t4Pg1`39zRwl+k z1U7OqBcHb1hu!$QygcA98X@1)E2cWRdP z)+n8^zrXgBf9BH=`}_OORhX;^M##@qWN(nWUSWSv)FQ8#MF39U_kQuG)*}q2!2#+L z`@8uD$ELUZx4+wxnc%l&qy6102cCXhHO2n!)nhUT-fTn2cdBkU`FyJV-8EB{74#P& O=u_WImtEflQx5>liO?1R diff --git a/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_nosplines_out.e b/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_nosplines_out.e index c15f6d9f39d2f6eb71fd373f6e89e0b69c622951..ca488a11c42a58a3fc4a60832de811e6c5cad0f1 100644 GIT binary patch delta 751 zcmb7=eJqp#9LFQ6VTVecOfKQj*gem=H@CZ=idJ3oHYJNT_hPg|s2Lk3rBh+&@sc)# zy1dm+lH8a!*;K^L74Dv=M|UbE(-O(cb*n$?pZ@%Q_u1#WeU0lD#&rvIl`Ud`F(tF^ z-Wm$yBT>RpBLP#Vvt85WG2;Fn|SUH%<69Jb5zR1@6n;d(Gw6)+MNmtm&$T-n!4wJ_^wIm?%I z!G~ws&js>Ff(nJF$z=L2Q1be~WEz4un+xT}%tZkfZ&7M8{r6R1tVx~}O!WiZAdR35 zj4ld_UiX;?mIWmqOQNeif4>-<Aez}*P;3K8ADKdpYJj}*t+U<91xRj{ z z^uaClOAV=@KU6;NzP2w1093qC8qS zGM@$bgvLc3c>r+zwK{{367Xf&#-1+&fKQIC+)7sgzM@ThD9Zu-2UV;+lLCY&IW9Bw zF%YT*kA~?ebMBtb@-h2Lg1?upo&}<5$EMhcWS|qJ)hS;5B*3Sm>877Xz~43(^&dw7 SUt65p&o2UeRp;5VkoXI14Pg%e delta 750 zcmaFzm+{G8#tqZC4WHOwj@czD;+|!Hb#npx9h;Z-w<4-!i*gj~9|`t~FLjKyzn9j# zcbCI-`v*@Z7=QXRZSo0jdBF?zXRTlEn}3?a{vspm*Wb?1Cr{*2+bqY!#U!3$X=svS znw)HAqL81bU}RuqtZQJYYc$zWUTktGpQMJVQJSSus;Px>T1t``P;sh>S)x&rg=u1< zv6-2LiIKT!s;QxciRI*Je5WNsGD;QvlS>p7Oce~xtqhH zF%TI+6a-pNnY=+vZ*z};P)_}fyLsuSm)bvOX*S5{XtjSjHz@1*(gXG{1({x2r|z3 zul_jf#T`-KYTM3h`+*jGx9wB;=H0#rGZ0;Bwd$>Blq|{Ww7<8~jMZrJfBOe&D*6X6 zD%n41{k(OfsEqvsQ}@&^Gk5#jYD)G`H=MPM3;UD6wZwC(TS z>6-sobE5s-AKQ%H`Ao9Ed-KxISw&0j@BS96`@hBB{_f8dla#P#`+MHGvE_c-?C&j0 zQ(5#d1jyg{uu@aW{$8(58Lz+@`+F;d8Deha*sI^`>s38=O3wc7!*KQgzn<9Ny}z>a Yu$O`T-G9o9_i(PUzx!9f@aD�NBM^VE_OC diff --git a/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_out.e b/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_out.e index 30b2416848d15006039e0924a454337bba769e91..168243484bfedca60b18e79f67afa863adf660cc 100644 GIT binary patch delta 859 zcmbV|eJs>*9LF8P)^<8KVTrRMwaMM@x;qb@50Qs?Y-{Z*THKw7Iz*ags|yV+&hhgw z*J%FWCN#P{WSN%IgfS1}oZr_ECo`=i-GuzQ&3~iysQ_1@lpyt=DQyQ@sw40Mae z8Lym?dVH6DO5S{<{99ffEVv<;YLx<}%?f1t>`h^1kt;<}KVWk7qTiSv3GsMZDMc;8 zcp3Mp`q*D$a)JJSY!gsoNpFvT4J>2{%fL#A<(RAJ&W<@ziGi59oiwcWwY_Sk3 zjvK-YUr&O?uvML|^_HgD`W4MyD&UFvVu`@Zo0Xo%!q^zsgTwZ~1UQ3fi_&8&T8JwBIiWX z9vi1IUGXU!hd(xsw-1M{p)s{t@_jR)ADl?qbSQ;YhUu!Wv-Rt>>NG5<+h8twP3v3Z zh?9?0^2%<_pk4+j^s$OD;0dU;L)~l{0;pngIEdIXaChPzX#A3>b+se>y3BQ|< delta 858 zcmbV|ZAeoA6o%U}wE}A`9V5ysrHN#@yWP!gtz%65NG*`G^$FYD%rx!8`ccpXHI*KzpK z7^;_|IxJqO&(VmPXPro_Mm6!6S}hcd<3*?h(_w;mkraQ$>E~|VeF)j4Ek+Ow5lCbL zRE7y~7B>wX7 zawX{wNJy0-G4bs_pj0IfFV=?xwXkPHuJ0+38=KqZ`AndyMY&2sz^JOr|IT(5+*^NE zbuK-RQ&m<&eXH~HgRelQoUn49tO2s)hOI3q6-fLM@AlkD#@9L#4dY*Ocb25ufvh{c z(jnRjWHQ};+IJR+$xN(eR~Qh^!2*9sDx-hScQ=&*;qoYhYhC~`)e-AR>jQ$Yv$ikW z2Bd9}vQHyG*3bM~aB03hgw>mPX+fq0)c5asRv4Lv)VZvX%Q diff --git a/test/tests/meshgenerators/file_mesh_generator/gold/2d_discontinuous_iga_out.e b/test/tests/meshgenerators/file_mesh_generator/gold/2d_discontinuous_iga_out.e index c6d6edce714398cb36699fa70c0b0a265860412b..8905ebd6b21b1c1af994a4f58a3bc79cb6b3d3ff 100644 GIT binary patch delta 1516 zcmb7?eM}Q)9LH&SRTC*CQkbrmAc$mbukA`HSRdYGke0X1EleSmwny7pIJ}lY*g{n% zbsJ#ip>eRt#vi~Ouw-`NvWzfaHbYVaOpFwIw<0)jE*mf;n>xD!a!XXQ?f!V4g(;$6&z>o zezR%eL>rijAtTSWHiDU<8``Ca)8LkGonfYONoG40?eTc-gK587M0e_4FgAQ(-$>MR z9v{)iU5^mi8R;PKcpTuTrRK|n9u~OsbpE+gr)=U1nAYTP{cdClz>|zv#SdpzU?ByQ zB@F{#C%1r+oUGi=o&Z_nu4TF9n=vp^Q2ioq;BfdltQ6%NO_&zf@wIvrzZAoaD7OSh zjXL-{1KG^u@qD-_j$=6FhR9(}NJTg*5v$|X8cBkbgOzcFh)~28AY7pYPBFJZ4l^Pa z)u;uCI3Ye>RjiIj#gdmb-a&#;Bo>yaHHb*4(i97%Dg=ea%$Mj=9;6{WIM4Klmzivu z$eSV$&A|1XT(yZKlyC&XM1eF>Ac5ADY|k4 zMx7QdLt&I&4nIlsvsdZ4JK(1$-#4$>N2zu_9_>xj>Lxui}qBEon z7Sm39JkylgYnCB0^Z#{pUB}xY4XlooBH*6?CE)!2*mWJBD)?8885HOUYMbmWuYk#x zR&$y)Z>m)WNDFx#TwoJAe~R#!Ohl>e&?MlH0DKOgWh3Enm#tcZlL#iW8;bu?4_9nYKN9 ze9^XHfxFoGdfgodaPMd$9%!wAEHJ01|ELG9x$WI_M(R97#90$sWHoJlP#ICp9{tsOGn!_Bs@>_amsATvK`4f#G6 zkiC6f-GRw~Y#wu@6ARwC+R0Jv0pw?{FOFsP0J71@NK)&MXXHSD%OjsT?9 zPq(-y&>L`eoL}1q$WFFm?~n>#xVmUi};Qv$=)< delta 1393 zcmb7?Z%k8H7{+MbVCy~H?#O2m}f!6s)b1MUhZ2yBn%Z`&}#_LYxel^bf^Yy6XgN!H~wMx z^sRg?y-0IZeF+%m2K+M9Ja~x_!h^;AJ4#)e(LuPM`_R2I>mxulyr{I1-G9UTSePoT z@9#@)gVrDn4gM}P2S zK{2?+Vzx*m5+9@p&Op=rW#~X~h0sulDRI3jR-qJ|%fxb-ToEHv#K;lNUxmbcQBoXy#pTeT_MHzRCQ1Yt0wafpsW zr(cNFhFfiOlT})ZTTMp144?aqy~P>GkrG6lD~?%%T($M~`3SInPX%+QlJGv5>^(WC zJ)>v+Ued?^l9P>78ia8g8D@gb*Syt58vFkZ%0edY1=kCmDfvK8I24{|#}+cNNZ?* zmWTnf575!+BK-wFph|olUuAv3a&8@#cDryK#v(@@tB!=RoO|bw-JdvOVN5c6@!hxb zSkArkZdZ_K8th}*(5yS2mL&yf%$0la!8M>oC#ue+T?9&JPffjB1hlT|@|Vr4f#wao zZS)KQ^-f1mQgP@@0t@{YIJd4-yFulGFNsIIedQYDYT|;J*N|;)>-V$ zPS>wLOaZl}u08ll6i{cs>FN&H1QZ$)e!G0edtTrY>)!62I7CqD*qbV3{0aE0Z_XJ=#jGX}v5kZ7T0W)A84B;oUC@u>D_f j)1g$({-%zWu9qj9_ zohp$?ohGI}6pmDh2P85L?PU)GzfM{uCQaJY`Y;J0ZA+duAs&E9dz-Xt_Xl=g?DL)P zJNKOLp0C&0`aSl$TMV<}R+tpSt_{22vl{_+Z5#^0UUq5FtpxvQ?MNn)YO;7%#3C{c zB*cY8jE~27sKvwXTDXOYj9D7WD+KWg#EU5*o=PTt8r(dH;Fpd=|0Bdx2{FZs+k`v% z1E?H~x=$ipll%NkkbSRU$%d_AgTPkoWV?Lbw1{LHl5Lh2u~d|`5>)5Txlw+U3C?3> zL6V&%XpnKYDD=L|ssy^`jzE1<45&nxmuPetenzLwMsUTE#j2v7nflwN>g)j_d)6z z+htn1T+%+nUds;AJhfIVN(RNH>&vthG$v#4X?iATIJ&QQtLX?Az{Jb(&V;)$@+ZS7 z^H>M@a0Du+xsaxpOPZpQ7-8Em5Je0B6!2 zupkNn3hazrIW0E%U3+sLg8EA0eS&3m*&BsmNyP*!Ub}E-b};;QW?Muw9YfoNRG3cu zMnz28k8CIy&ZP%vU9LqcV@oqti5hA{?Sl~y4DEt`DI2QUO!E!|BP0y#SuyGvDeaX* zZ=W>cci*aocFC%{mZ8y|hO6hNTNSzw#tpi>3NbC% z>dV09hj|#*U%UT|_3wwuiK)9M1DtR0%}xVCpP&47sJ;_72HGE8xw^53|Ni5P+&2#w_VT;D`u*Eazm0u7w8Q@o DNQIE- delta 612 zcmaDblWEBWrVaIsjJcZ|7_AtY+FLg7VVcIu)ZQ~WghPv|xM6bx$9-OLOJfsbV}qo$ zBqN3VJOv{IBV%1d16{+(fm&jld&GDcCBjov6#SD*6bua&49%L*cqPHo^IHQa|M5Awf=J|djNN1T(^E2~dV z2 Date: Thu, 2 Nov 2023 12:30:35 -0500 Subject: [PATCH 12/17] Re-gold IGA tests affected by libMesh caching bug --- .../gold/1d_discontinuous_iga_out.e | Bin 34352 -> 34352 bytes .../gold/2d_diffusion_iga_nosplines_out.e | Bin 32484 -> 32484 bytes .../2d_diffusion_iga_nosplines_out_002_0.vtu | 6 +++--- .../gold/2d_diffusion_iga_out.e | Bin 35640 -> 35640 bytes .../gold/2d_diffusion_iga_out_002_0.vtu | 16 ++++++++-------- .../gold/2d_discontinuous_iga_out.e | Bin 36044 -> 36044 bytes 6 files changed, 11 insertions(+), 11 deletions(-) diff --git a/test/tests/meshgenerators/file_mesh_generator/gold/1d_discontinuous_iga_out.e b/test/tests/meshgenerators/file_mesh_generator/gold/1d_discontinuous_iga_out.e index f527d5937e21b8818cc98a44c2ecf57d3f9f8f1e..fe4f6078d25d3b0c737144693691377661bc2791 100644 GIT binary patch delta 326 zcmdnc!?Xd2-1#L;43iU$%#)Lh4HWY86pRdvjCBnSbd4+~GxCW}F5#EdNHVufF-)^C zH84&#N;5J@H8ZzNO*Tn1GO#o=NU=yxF;6u}F-bJFoV=C)v_wcose)gAnSz3mf`O%# ziK&&T1=viWX$q5f2uPwUFtIW;gD6nUojX~GxCW}F5#EdNKHvLG%&R= zGfPTLHcK@%H%S2sSQr_b8l@$t7#JHRr6d|!CK;qo-pYSk!mTJ%!9Te~!N^>}(8$Wr z(#p^rY$nh&g~>YvB(W$kfGAMRojX~=gAEb%ooeas I8QWkQ030}e%m4rY diff --git a/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_nosplines_out.e b/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_nosplines_out.e index ca488a11c42a58a3fc4a60832de811e6c5cad0f1..11ea3c6aec0682d01ca5820153ded39dd9f594a8 100644 GIT binary patch delta 568 zcmaFzm+{G8#tj;L5+;VpiALthNyY{W`FRRP21drZh6cJu7LzCPh)xdWlhjBuw@fii zvoJL!lX<_fEe4pZkZwzWzaNHG_%6B>Q_;Dz`dUowa{p+PpT)S;YQ9+f~VL2N>-i z*#1dSmbz?zhcT-$UgCrO?XUV5f4qEZf9JkXSCL4a{hd2M<0k)I;J^R@_V=Ye-Y;I2 zW`A%0+e^~Ar-5|H_o|;x_V>lT6hqp3?e8;f`h4JQ4^aHzk(E4E_4fDvpHJj4*R;Pc z`K-ow&ouk{jQ4B*wZ+)qi}ie3S|VeAFX1D1(rZ)udzGFYCI9^G@5O~HEj1{yzgIn> zhEXU0DE^^I{hOivy?ra*R!qHMe_v>K-qLUCK>mZ(6N;DG-}_^=UNnBH{e2O?HGvNz ufgxE|a;HGw{$6U^tN2BdGwtsszgg(q+hl*Qt~&bSek=QXwIw2V!}|e(4F?$j delta 568 zcmaFzm+{G8#tj;L5~=2<$!5tZ<`$L;`FRRP21drZh6cJu=94G#h)xdWlhjB}Nj5Yv zwJn41IN{_au$Q6z-9KO6C|h^h-~D}Cb=8Lq0CZ9R*#H0l diff --git a/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_nosplines_out_002_0.vtu b/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_nosplines_out_002_0.vtu index cd8c86ae4056..01c0055492df 100644 --- a/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_nosplines_out_002_0.vtu +++ b/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_nosplines_out_002_0.vtu @@ -10,9 +10,9 @@ 1.000000000000002 0.853553390594143 0.8535533905941495 0.8535533905941719 1.000000000000019 0.853553390594156 1.000000000000001 - - 0.4794383698264516 0.47785813166466456 0.47345648342332447 0.47937076789790656 0.477804004039832 0.4734428002629079 - 0.4792017393054785 0.47766875832572725 0.47340889522507873 0 0 0 + + 0.4860195810457867 0.48453973583237864 0.4804644157236139 0.48595583094624695 0.48448931984066745 0.48045391278373695 + 0.485796482609765 0.48436340785780146 0.48042797248844604 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 diff --git a/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_out.e b/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_out.e index 168243484bfedca60b18e79f67afa863adf660cc..9cfb506a413b8af0524afc9e3ad62961e2e6e89b 100644 GIT binary patch delta 808 zcmYk1YcN!C7>C;(MVToqYsHL}Tcc5DkF^&jZ=>00lxU<2sb!n*yGnmrG zU=mTbKImejb*siUN}*|Kec0oFj!@BUchX(y%rxlz_&)Fa-sh*UX6vijw+Tm{P@s_G zVnuQg#?Q*+;}{+kD8K@7iDrW1qN#Lr+n6kt5`q*-Fczf1Q*cZv5=)f|p&ZAgB8-qI z2(c0)gz{ji#_0Hr6PK3FkIu^B^Krfa55c7&IIdxFd^Inf-2P{nG(;#25r}jw&XsKS ztuT+q+RH%e;tR{|CxOgAr1#_f0&|~3zG-tB(5)BRe;Yc1zEilvCd>=ST!Two&qFXj z62+yB*8+W-Ccfl!16}l_le@Zw;R!!DVdG%_QD3E6X$N#ceZ=(8aiH@oS50eFOuqjC z>lPW1tggM&8hmwI%b$K7`p+P-f#EcFb^a$Z$&W+OX4JI&E;M90`t>@{Vk$yAgRtHGOjI{(afC1 z8Vhvs^c!Qfy&C9k3B#?c^neWa)E{Ki%f{lyn|B>}fVymAqnkDZIdO7IS26_Vxo>1s z!&4w9<5{}dy2U2jR{R^F+`TO)c@Kc%wl@zXMgdC2AD=Z_0VU;MYVF&2`T(I_%C-c~&i{<7F6x;Qp zTr$M$UxYtbSB3)Rn(jNZsUOgSg==NMVt}fC{sAA5s_30}dEt^d%JtR{k0k**T-V^o{vx42E+zQ_XR?~_>Zeu_BB59&k))>)617;mL&XL%JWz}uB~-#d zJYAJ0lBWyRY-Qp)K z*H|Ib{tF*pX@POCK?@oZ=p?N|F_pb;EVx%$JWb z%n4si!PneZKMAqfNN9xK2evHa+( zB0DXYEi}dN9hn0(X|J3sV+;DNVHfk-fo>tm#zZ$jpEMlDj$%O5sfMxMOrTpm;&u)f z0)6|u{8r@vpwAI@iatG{_Zee8v!m+*SNaXZGg&}KUuxhF1+WKu^>^%}fUbS}LV@Pl zoG@Y9bQW}&dx~QWK-ZM`e0LrKI%=i!x}_P={L-0)!8dE>=yxkxH+j8F4^r} W$S@$L?{Zmk2%u%cDaMXm{p2s&xqle| diff --git a/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_out_002_0.vtu b/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_out_002_0.vtu index cd3587801499..7d6e762e7216 100644 --- a/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_out_002_0.vtu +++ b/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_out_002_0.vtu @@ -12,14 +12,14 @@ 1.0000000000000115 1.000000000000019 0.8535533905941581 0.853553390594165 1.000000000000002 0.853553390594143 0.8535533905941495 0.8535533905941719 1.000000000000019 0.853553390594156 1.000000000000001 - - 0.47941867657053144 0.4793470160152781 0.46700861382638914 0.4664190830301079 0.4794061486286514 0.4793335948962325 - 0.4670174842023429 0.46643240648164813 0.4789510392927619 0.4788559857818215 0.46737074950821267 0.46688048895779016 - 0.4788890594965306 0.47879152897639765 0.4674236823510491 0.46694417032597063 0.47941867657053144 0.47788827474102535 - 0.4731778149208336 0.47935239194713514 0.4778339381041117 0.47316833690414034 0.47917859396070667 0.47769192507288627 - 0.47314445359715235 0.4683219129028348 0.4664190830301079 0.4683625381259725 0.4664850859282808 0.4684695154461442 - 0.46665644771971915 0.4789924330106902 0.47754011391283846 0.4731196986377274 0.4788890594965306 0.47745638288945796 - 0.4731076056637234 0.4685851008080247 0.4668403989903176 0.46865155015628907 0.46694417032597063 + + 0.48603151780879367 0.48576107973844235 0.47504872784287416 0.47445598718343496 0.48601958977655924 0.4857486294487372 + 0.4750570768596925 0.4744689533746331 0.4855891196012813 0.48531027773763513 0.4753867178447443 0.47489867992056944 + 0.4855307957818521 0.48525155364613637 0.47543616234385805 0.4749595435801812 0.48603151780879367 0.48455306103882423 + 0.4804049037906583 0.4859687630127081 0.4845026921292515 0.4803965713929476 0.48580435468892025 0.4843712066686249 + 0.4803756754727022 0.476166228893085 0.47445598718343496 0.4762048058286064 0.4745194276450756 0.4763062311269181 + 0.47468381664760123 0.48562834741833366 0.48423076218860617 0.4803541322625197 0.4855307957818521 0.4841535162208571 + 0.4803438579949972 0.47641578987894656 0.4748601800172304 0.47647876717561843 0.4749595435801812 diff --git a/test/tests/meshgenerators/file_mesh_generator/gold/2d_discontinuous_iga_out.e b/test/tests/meshgenerators/file_mesh_generator/gold/2d_discontinuous_iga_out.e index 8905ebd6b21b1c1af994a4f58a3bc79cb6b3d3ff..75f2ef90ec4adacb597b594da3601f4e224bfd20 100644 GIT binary patch delta 985 zcmZ{gdq`7p6vuOD@2obpnL@WTm7=`fZFifIM~wBLx0pS&Z0^pAN{KR3%Q%c-1coc4rDU{BYAg}WH&UX&e>G`z6lpFI{m9hc zj0vtND*@J=m|HaEEzsH99)yYez}`F~+qQTour-FeLzWs~Yjc--$M^xAW${mFtOfg1 zGGX^S1F*+f>iglFz~()#3z%f$^rS&i%uBEjno4!vL%`;kR(xpZ?i_JmwoN|3-4Ea5 zic|pYRUcoK{{-|b|3F&HujWJ~tbBIn9?B(Ud{ z#X6G&*bN<&rRx?T+vjEB^Bf+IM)sv>!&grLo}*Dsmbh97_7_`Qu99Y+(*w=arKtP| z!Hyhmr)$NS25iBHk6oq1z{VzZ)QrCY$Z}u*MP(aJD6hD2XNE7JdhdGz3V=DuRSqS$dG158ke^N{#DFadWjwI<5}?INBXxAy?r zo_!@?VjwWHU7uc`e$DB5do`^-oW99QV!51dp0NLf4j55y>Yrx@V5Z00y1rHd+SjzdocW z-IIU^QC7%tvj z;zAW}PeAlLWy1OG{eYU>YZfdKP8-vJ2T{(sD(P@6*HqAeEs|@d$vwgyt+1BlJmDMq E3nu{LlK=n! delta 985 zcmZ{gYfwya7{)7_Win!SP-=3Rsn&jV8izR&ynKhN`;t3Axs9ybT( z`e-F#xER++q%tgPHzpE@LPLat5RrtmbLWv5*O#x=;6g!|R4h(Y<6>=?WQB$uNJXJx zqI6s%2opYZs+L14~B+`oL11j-UV9JP5AD8572 z_#g$CO)sRXRt_*3XYtD1K312_8tIe+wMpUc75oYqRsPc}IJS7MZ5idK;LLHM4n|EU%%j&xcC;JcX>XY+b;sL zcI6*FE z`K)1+{u(D5=+J)CjP5A`>!&u*cQD3LU!aNU+Z;UWMJWE^``UOsw B*v9|> From 29b8e4ae6acfb7a30723b59ace04cd37c207dfb7 Mon Sep 17 00:00:00 2001 From: Roy Stogner Date: Wed, 8 Nov 2023 13:24:37 -0600 Subject: [PATCH 13/17] IGA on 32-bit ids is fine with the new libMesh --- modules/solid_mechanics/examples/cframe_iga/tests | 2 -- test/tests/meshgenerators/file_mesh_generator/tests | 6 ------ 2 files changed, 8 deletions(-) diff --git a/modules/solid_mechanics/examples/cframe_iga/tests b/modules/solid_mechanics/examples/cframe_iga/tests index d88382346c2a..fe450b76da32 100644 --- a/modules/solid_mechanics/examples/cframe_iga/tests +++ b/modules/solid_mechanics/examples/cframe_iga/tests @@ -11,8 +11,6 @@ method = '!DBG' # See #24187 petsc_version = '>=3.12.0' - # See #25819 - dof_id_bytes = 8 valgrind = none [] [] diff --git a/test/tests/meshgenerators/file_mesh_generator/tests b/test/tests/meshgenerators/file_mesh_generator/tests index f296c20ee432..d4e9bdb38369 100644 --- a/test/tests/meshgenerators/file_mesh_generator/tests +++ b/test/tests/meshgenerators/file_mesh_generator/tests @@ -31,8 +31,6 @@ exodus_version = '>=8.0' # See #24187 petsc_version = '>=3.12.0' - # See #25819 - dof_id_bytes = 8 # Set to >=9.0 once #21449 is resolved vtk_version = '>=9.1' # VTK diffs via XMLDiff are really fragile @@ -69,8 +67,6 @@ exodus_version = '>=8.0' # See #24187 petsc_version = '>=3.12.0' - # See #25819 - dof_id_bytes = 8 # Set to >=9.0 once #21449 is resolved vtk_version = '>=9.1' # VTK diffs via XMLDiff are really fragile @@ -107,8 +103,6 @@ exodus_version = '>=8.0' # See #24187 petsc_version = '>=3.12.0' - # See #25819 - dof_id_bytes = 8 # Set to >=9.0 once #21449 is resolved vtk_version = '>=9.1' # VTK diffs via XMLDiff are really fragile ... and we can't add an From 935b4683331f3c03cad85a88619f678b6e764566 Mon Sep 17 00:00:00 2001 From: Roy Stogner Date: Wed, 8 Nov 2023 13:25:21 -0600 Subject: [PATCH 14/17] IGA tests can all be enabled with the new libMesh --- test/tests/meshgenerators/file_mesh_generator/tests | 2 -- 1 file changed, 2 deletions(-) diff --git a/test/tests/meshgenerators/file_mesh_generator/tests b/test/tests/meshgenerators/file_mesh_generator/tests index d4e9bdb38369..f621307a40d8 100644 --- a/test/tests/meshgenerators/file_mesh_generator/tests +++ b/test/tests/meshgenerators/file_mesh_generator/tests @@ -23,7 +23,6 @@ [../] [./exodus_iga_test] - skip = 'libMesh issue, #25819' type = XMLDiff input = '2d_diffusion_iga.i' expect_out = 'Solve Converged' @@ -59,7 +58,6 @@ [] [./exodus_iga_nosplines_test] - skip = 'libMesh issue, #25819' type = XMLDiff input = '2d_diffusion_iga_nosplines.i' expect_out = 'Solve Converged' From dcd05f0a6372eed8d088beb0b4d31b9c18e2e31b Mon Sep 17 00:00:00 2001 From: Roy Stogner Date: Wed, 17 Jan 2024 12:33:11 -0600 Subject: [PATCH 15/17] Update VTK golds for new-libMesh compatibility We've added a libmesh_node_ids array so we can preserve those in future reads, but that's enough to confuse the XML differ. --- .../gold/2d_diffusion_iga_nosplines_out_002.pvtu | 1 + .../gold/2d_diffusion_iga_nosplines_out_002_0.vtu | 7 +++++++ .../gold/2d_diffusion_iga_out_002.pvtu | 1 + .../gold/2d_diffusion_iga_out_002_0.vtu | 9 +++++++++ 4 files changed, 18 insertions(+) diff --git a/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_nosplines_out_002.pvtu b/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_nosplines_out_002.pvtu index 5f1498f2fbbb..4d8974ad9f8f 100644 --- a/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_nosplines_out_002.pvtu +++ b/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_nosplines_out_002.pvtu @@ -2,6 +2,7 @@ + diff --git a/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_nosplines_out_002_0.vtu b/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_nosplines_out_002_0.vtu index 01c0055492df..716eff514797 100644 --- a/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_nosplines_out_002_0.vtu +++ b/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_nosplines_out_002_0.vtu @@ -3,6 +3,13 @@ + + 16 17 18 19 20 21 + 22 23 24 25 26 27 + 28 29 30 31 32 33 + 34 35 36 37 38 39 + 40 + 1.000000000000002 0.8535533905941428 0.8535533905941497 1.000000000000005 0.8535533905941461 0.8535533905941527 1.000000000000012 0.8535533905941521 0.8535533905941588 0.8535533905941567 1.000000000000002 0.8535533905941592 diff --git a/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_out_002.pvtu b/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_out_002.pvtu index e268c3ebf7b2..d541020fcda0 100644 --- a/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_out_002.pvtu +++ b/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_out_002.pvtu @@ -2,6 +2,7 @@ + diff --git a/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_out_002_0.vtu b/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_out_002_0.vtu index 7d6e762e7216..365c26219c24 100644 --- a/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_out_002_0.vtu +++ b/test/tests/meshgenerators/file_mesh_generator/gold/2d_diffusion_iga_out_002_0.vtu @@ -3,6 +3,15 @@ + + 0 1 2 3 4 5 + 6 7 8 9 10 11 + 12 13 14 15 16 17 + 18 19 20 21 22 23 + 24 25 26 27 28 29 + 30 31 32 33 34 35 + 36 37 38 39 40 + 1.000000000000002 0.8535533905941428 0.8535533905941567 1.000000000000002 1.000000000000005 0.8535533905941461 0.8535533905941592 1.000000000000004 1.000000000000019 0.8535533905941581 0.8535533905941719 1.000000000000019 From f16abd7a9ebf64ebf2af2e1a9e03fc6ce1af1a7f Mon Sep 17 00:00:00 2001 From: Roy Stogner Date: Mon, 11 Mar 2024 16:09:19 -0500 Subject: [PATCH 16/17] Just allow_renumbering on new tests If we switch these back to XMLDiff then we'll have to cope with it's limitations, but that can wait until then. --- .../meshgenerators/file_mesh_generator/1d_discontinuous_iga.i | 1 - .../meshgenerators/file_mesh_generator/2d_discontinuous_iga.i | 1 - .../meshgenerators/file_mesh_generator/2d_discontinuous_iga_l2.i | 1 - .../meshgenerators/file_mesh_generator/exact_discontinuous_iga.i | 1 - 4 files changed, 4 deletions(-) diff --git a/test/tests/meshgenerators/file_mesh_generator/1d_discontinuous_iga.i b/test/tests/meshgenerators/file_mesh_generator/1d_discontinuous_iga.i index b3fcf6a83548..fe0c77d608de 100644 --- a/test/tests/meshgenerators/file_mesh_generator/1d_discontinuous_iga.i +++ b/test/tests/meshgenerators/file_mesh_generator/1d_discontinuous_iga.i @@ -4,7 +4,6 @@ file = test_2edge.e discontinuous_spline_extraction = true [] - allow_renumbering = false [] [Variables] diff --git a/test/tests/meshgenerators/file_mesh_generator/2d_discontinuous_iga.i b/test/tests/meshgenerators/file_mesh_generator/2d_discontinuous_iga.i index 60620c017cdd..e90a1358f132 100644 --- a/test/tests/meshgenerators/file_mesh_generator/2d_discontinuous_iga.i +++ b/test/tests/meshgenerators/file_mesh_generator/2d_discontinuous_iga.i @@ -4,7 +4,6 @@ file = PressurizedCyl_Patch6_4Elem.e discontinuous_spline_extraction = true [] - allow_renumbering = false [] [Variables] diff --git a/test/tests/meshgenerators/file_mesh_generator/2d_discontinuous_iga_l2.i b/test/tests/meshgenerators/file_mesh_generator/2d_discontinuous_iga_l2.i index fc23b6a6c968..ac1e6b6b9221 100644 --- a/test/tests/meshgenerators/file_mesh_generator/2d_discontinuous_iga_l2.i +++ b/test/tests/meshgenerators/file_mesh_generator/2d_discontinuous_iga_l2.i @@ -4,7 +4,6 @@ file = PressurizedCyl_Patch6_4Elem.e discontinuous_spline_extraction = true [] - allow_renumbering = false [] [Variables] diff --git a/test/tests/meshgenerators/file_mesh_generator/exact_discontinuous_iga.i b/test/tests/meshgenerators/file_mesh_generator/exact_discontinuous_iga.i index d7d67edd7c22..29ff09443b49 100644 --- a/test/tests/meshgenerators/file_mesh_generator/exact_discontinuous_iga.i +++ b/test/tests/meshgenerators/file_mesh_generator/exact_discontinuous_iga.i @@ -4,7 +4,6 @@ file = test_quadratic.e discontinuous_spline_extraction = true [] - allow_renumbering = false [] [Variables] From b922a8224860ff37d7f11558cf94916c169c34a1 Mon Sep 17 00:00:00 2001 From: Roy Stogner Date: Mon, 11 Mar 2024 16:26:48 -0500 Subject: [PATCH 17/17] Make SQA requirements unique --- test/tests/meshgenerators/file_mesh_generator/tests | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/test/tests/meshgenerators/file_mesh_generator/tests b/test/tests/meshgenerators/file_mesh_generator/tests index f621307a40d8..22de96a6c9d1 100644 --- a/test/tests/meshgenerators/file_mesh_generator/tests +++ b/test/tests/meshgenerators/file_mesh_generator/tests @@ -37,7 +37,7 @@ vtk = true xmldiff = '2d_diffusion_iga_out_002.pvtu 2d_diffusion_iga_out_002_0.vtu' - requirement = 'The system shall have the ability to read IsoGeometric Analysis meshes with rational bases in ExodusII format.' + requirement = 'The system shall have the ability to read IsoGeometric Analysis meshes with rational bases and output spline results in VTK format.' design = 'meshgenerators/FileMeshGenerator.md' issues = '#18768' [] @@ -52,7 +52,7 @@ # See #24187 petsc_version = '>=3.12.0' exodiff = '2d_diffusion_iga_out.e' - requirement = 'The system shall have the ability to read IsoGeometric Analysis meshes with rational bases in ExodusII format.' + requirement = 'The system shall have the ability to read IsoGeometric Analysis meshes with rational bases and output spline results in ExodusII format.' design = 'meshgenerators/FileMeshGenerator.md' issues = '#18768' [] @@ -72,7 +72,7 @@ vtk = true xmldiff = '2d_diffusion_iga_nosplines_out_002.pvtu 2d_diffusion_iga_nosplines_out_002_0.vtu' - requirement = 'The system shall have the ability to remove spline nodes from IsoGeometric Analysis meshes.' + requirement = 'The system shall have the ability to remove spline nodes from IsoGeometric Analysis meshes and output Bernstein-Bezier results in VTK format.' design = 'meshgenerators/FileMeshGenerator.md' issues = '#18768' [] @@ -88,7 +88,7 @@ petsc_version = '>=3.12.0' exodiff = '2d_diffusion_iga_nosplines_out.e' - requirement = 'The system shall have the ability to remove spline nodes from IsoGeometric Analysis meshes.' + requirement = 'The system shall have the ability to remove spline nodes from IsoGeometric Analysis meshes and output Bernstein-Bezier results in ExodusII format.' design = 'meshgenerators/FileMeshGenerator.md' issues = '#18768' []