From 035a3a4b74baeae717085b1938868333ea498ba3 Mon Sep 17 00:00:00 2001 From: ilhamv Date: Mon, 20 Jan 2025 03:17:53 -0800 Subject: [PATCH] add time census regression test --- test/regression/azurv1_census/answer.h5 | Bin 0 -> 196072 bytes test/regression/azurv1_census/input.py | 50 ++++++++++++++++++++++++ 2 files changed, 50 insertions(+) create mode 100644 test/regression/azurv1_census/answer.h5 create mode 100644 test/regression/azurv1_census/input.py diff --git a/test/regression/azurv1_census/answer.h5 b/test/regression/azurv1_census/answer.h5 new file mode 100644 index 0000000000000000000000000000000000000000..2e7bc4d504801fbb68f7dbeb9f6c7426a0f8da5c GIT binary patch literal 196072 zcmeEP2V4|67hj5s9Xpm&#D;ybVtr%BjvXr?psdoQ2x7yIdLoDg3%WD~6loSjkP)Sv z9tzm7qn?T#&xZY**(4yl>~Qd*?&^#`c`0v_mw(=S$t0QCwFUzRSFO~bl7PK5Gz97b zO-Yda`HHWc&#ICn5-0g4IG}>5JEp;iU7#ZQClF}i^cuK)X}<0!JWU^ER~y5vSVR7Fq@@!}g)1V{tP{c4P6v*G$GDiY0Rtq0{HzH$sd8&_3D z$r{)a)JQ>6_9hVc$dN=blE;l2Jy7KvFX2)S$wlxG06Vy?JFBiu(pQChb z+&s%oh=cH%q@s}98AK7P0)e)gV&G9)DdhP@Y5b61RF&%&wHqP{EU_SVbrDWfX$9l& zQz``H>5)6Ozu*Sk%_+}hYht5fHD}0x(IX@pgr@;Wn*;)ukG+H)^z=IDG5!oi^Z!bl zBOb2v^RZhHpoxY+U~e(Y%EHW1G~WWTkn$_wJM64XZ7kZm2uNJseZ(ISK5}10Tn@>% zHg&YH7nxd-a!7t#2cdRL?e-2#lO$JVOo{xntejlfq$)U_)HBY)(NSbGTTm6J6ZvhO z9GQVfAgG4Z38_v3m|}PZ1UEA~CrM5XoKNH%Z((L>Bbw)gVG`8D`Gh*yS&1Z)w70i) zvSNh7jI}Tyk;~H5%8?ls1+{TH>5o|=7YlPHtqx8nayi&K*_&Btix|{Q9oYdDUpb_~ zB!Qvb|WNu+B#C&A^ ziYg^QXuZ3xiR%nD62lA`k0R%UYO4(8pXw~i27hYH>vUr8=!LRu)P^# zJ|o;{ZviKMk%`|fyR1DQb55zs;UI? zf_+rCrlM9F>)y9Yz~B;}9{c2i-NM`_B;qRgpJn8;e!@lLiqEc1wsf zfZY~h^CgUgK@OFziQOC`C6SF6IX818%a;3%#H?V6-^3;5eA@z6hO%jav_Y{|k(xi_ z^Es00nLpCUb=Z3e66{tM`=7*FX{EOhnWQ2wQeTmoWu3TWWru!jnaeR12uN2HP%Ka` z9FM`^F%?Oa0Mq`c;@kvGfwe&bHQCKXdF)@( z5&?zMG%(e~v^J(%m^Q*x8&e%jn`7D%Q(a8kVA>AT4wwot)x)$CrkyeE!cysW6^egr zTfWkHv(_`r_*=Mf({IJqm7zR2o_u0%c+^mvsoE5U;#ka+3wk)A4&l3hsPr_kc? zAnle%OOJ=RDQd8^pZQ(!2v!n23i}Ufr7gWa$g>Yc_a8`qiG7gvm;79O1nz$_(;>+E zkbpu`GLJ)ER3OM}&n>19XKQri=F`X7TI2>f>=aJd)sw@ZB+0UVxUiiO{UVg;S;@y~ z8ta*~-&HI-BkSEpxc?9u0s)Qpb4a(6^UY-Ai{L{Mh7wYF>*K*9RFEbSNUB`f{*(3# zx{}9mEoKFXm63hAM6zlK2KO5=j@_OaKXwAzrC~DU&$CjSqN>r|1!5>x|H&oS2oqi!@mWjCGl+>OwS%nga6=d})lq?hR#3?zwWwMki;m7%%mr9{m zr1S7GEWyx$C}kMVFf+Atbh5V)m|2R9*$o<&vSTt!9T>{I!EU;Mr6LDMd(liM@-E&i zk%QztcT7!ZT3Ik$%(lOk>1_5s8&l>bca}Oyie>369G49FOq^vCIGC9-?+n@t9PBJ4 z?*Tg6GH=Q==}z<5TrxKo2*i?@;{Ew?`u1Ka^zFytp`wGdLD|NAB(DlqWi>%rJ&#tz zc#x8bt&kCb#D6YD<6OD?uT13i;g&~?vkfV3zF?e9_vhw|z}e{oxcOpmb~g5x+}uAt zx484_x%m8`v~~P=N%zd%Gur?TFCzwoMXV#EJs|yG@^d@=aeK+gPP*NP z1Qe3W>c{u;A7ka{gJSs)@?2tsdkc9F=J}H(;B*Hk`&p)D1@BptI5~bq^hx+1dtKgh zOzK=zK#q0M1(A6tsfQk0G#^RT$AR>5B&BFd;pFV`Mdy=L102YmKvI;G7H0?2ts;J+ zD3&i(5sT;QPOocdpS6GDx|8E z91rw;`jU^+u5xfu^l|!ojZ4dw0uQ3J^3I!#GR9v>!i?XZKVf}n{<#x*=S|4@l48!4 zknb=J9Vj_w!kj-r-(@nluyGJMvcJTDuK#-eL{Y2}jdL0TG8@yMjvFF}xZ!k9V>t{s znOK$BgqzO?XGfTF^F`ondL}nt49<3uoDP!DOXDo}+46bVoK&3cLk4a+>_5(rty2oS z{c(6~Qc^rLB$`B?v&!z=&JK(RDVZoNu}BjCxwN?RCMNg`{b=<2xw7)L|9Aw7`)*zg zJGe8Kw&X9fp0^MMC>)P*Lje;U7q3A9Vgv*r6yZAUo_N|OBj^eR9p z_3HcKv7g6-v5A6nC(%3}h447e<538YQ#>Ap@HorkQ3#I zTqGYqHst?WE{)H>BQ9p$+QQ0;{q_lYUN9gTWlYJuQx}&jbmFdn5lM1U^SSx-akil| zHy?$wO}9tejc9mBi;;SlODP z(1e zpp`xaK?MwbKRo*Icrg84uq zHy@3&wf^MhOU2oC9^8CcINNO_H{XBFm9;Z1$*p8QDP2zK=S#}^o)Pk!Wb`2IFZuhJ zBe9(+*Yj6#=a^`3q#@c^QS^7C0EOdp01BAkcyK5R5Hl2FgCly4*?6KP2*~=6qV1*R z?>kIa;&Dkm|9n-%Chm?T-aumuHy`rK0)f_6Za!_Ct-FnzPtkL~pI%ugrCxnMJgj*< z7-J}y4?FRA6vAUMk4GUqR`GZg!eaxEMunfXFSE+d?q+s?8nWg=()w)S!vI~Bw_s`D=!2I zW%oN}wO@J(<3naD1Rt_uOX5G5%Ib%fL?E(nL8K>?tlO7-9rCgqf28PmpqZfb_tDgF zw-w-m5=xKHEsPJszd$Qo3W7f|)DPFqpMxpeP9Wo)B`!9?pSyfJoUIYW&F6-*wf1uJ z5$n_4$IYkcxw3vtLn3m4eU_D9f2*W^T^f(lK5x8}0MTJOq?*5{33(_J;FpLC>rZG5K^I*~WB$dVi?QoJ(G)=|H zE|Ep&lhiC6aE~G>{~wk0!&;I8$$V0}^wQ6le)n@Q)v=u^*Yj6#-+RhBiZn#pdr0dI zAOAHGK>DF7_xHtwv8dg;xcR}SQNaDzgoAID6ar|RF8G}<4g7sLGS8R(u9MX5d=w)f zJ5D2yOK>dt?|e0pGw%~C1?J>AoQyl#$GF8Nqpr(wZay^oVZIL*$IVClnNK`7pQ7iM zy`&FORu|CHuInwuIFS*Yn2S3JC=?w!5Q#(H zCpNaXm@TrkF&=9;N|G^4#C%7HVUket5-^j4*a&4g&}DDyXzIumfTIHPJ6PmRVwv}f z?M-bQtO{-+^GIp*WbM0f|Eb~%2`%|}xDEH4wBJ+A@sO-L3URxGPY}e5rW8(&Nhms> zq#EErDos*~rbak9>txaSB-IiJd`^*+qNxi`&P^&hpQPqnQ3^@<|7dY~C0+XOx9(Wj z&I{ALqS=_p#OB%vKLt0Vz_-}!Lyzhh?G)5zt(HR5-Ij>CP0V)^|bh`af zKwde)vrxXq7;e7o?YO{Qn94!FIDJk(k8YG61?3dSC;wsy5IrU1yEbm0HuiI*xn#m8 z`OYMn+%cI%`}7Ipf27juvqCiNI??MNUc$A2vTw_%f=&!CZs^h?R-#R+(vmLB26 zN`WcTtAgKKpwDrOP5ho#1~=dDmMg1Qh3ADtzX&b4o~6jqGsS%VLgT1F?tbI%rXPHJ z3nL&sjuP)_d9g4bu@}R_$)%3Tn@8s0dNr<~qKlgcQr8e4k-O9a@i!_Z0@C`6r(*(L ztUsht4z^DA-~8--{(Oxh3?-!U^fwsSn|c)$i1b%ouK)G9f%vp>{hAL^fWmQ~rzl{F z(t#INVx;^xyP z+1VtZoV0koDy#QBu3&v58cXnUCjo_^vUn7Jt|R+tgqD82KMUhQG?O&{cf})z$AeLU zg7-i2loStLi6WEtJc*siqoubWjTXpav=Pbw>yjV(8G+*b2l0M{mcISExcz7zP|#LI z_a8{!g7%ly_#(u3kdg@^iiHQ6hZ|w>$b7>N&##Hp))o$yl26);EF8Z5-71POl#t4s z-)GfAc;(fR#tl^u4J*}{9= zd=$>6{^90R^julJLY6?@Z<8J`O7ACoU_Dc|c~)83*-)jlGg@Nz z(-v70Lyg%GOWB{SW2q}cEo?08XEO=2?QNY{+Y~w2I@;UXiOiTVd)e9B&Sa^h$l8LX zPS`M=*h1JiTL^PGTZ$Z61_@I)4KgNk@CefArCuoxJQ|XVxBZ|j^>FI*?r^H_Z58_1$CHsaY|9#ZV zUjmOh*qR9?c$ECR^LH^GWJQa(3FY!1v!4nWMjQ9^!Rf&cP8EIOxs!Zz0xGSGySJ(wX}z>o$26RHxt7LccxEw3yX2C z*_qBxUvv;|b)rAb84*|^v=gmWjw#8h6U`R0tzReFu*dOL##KAfRu|SSY;{$ScDJ4~ z??s><9ay2w?4x#iH0}3(+S1;7^o?GBjUA$)N9SGgF{pk~NY`d7-Xo-il^BWU3TdYF zHnY16X;z$$YC^hVDsN zHuT*|W;%Tg+t3Xp0JNsZN(63AuX{}~HMgb{by);-=?~ME##I@tOLy{OD`-W}VQOvX z)Qa|Jnip546@9Hei%?6t=FFWFZ}e(O58K+aztx2n^m10rIW6eeFD^_ATG0P?Vi9ak zGw9^**K1Bo*#Bxqv)yynq#1odZAGG0}|AU^tGHY+0ca7)=ESxPH z(Gw@IS3^3PZN|z5bhnn7OwIM_y&js^58c(G>%JQmb+LvPy_!MoHC2y(%Scz-t}Z=G z(qVPzET$XBpQ}w9sIUOnqGQ-@XjhBYWs8|wllEbozpe(Y5#To9)iF&vm#sFlI-UOT z!#xXOH98{3c~g)_721<2{r=NRv>l^1&4*W{4Vhjy52-*4B??!b)=Xn>P^WX5)Mlqu z>2yYI&2&^~6Sg~4z9wihtxbymn2?9exG~bGZMl#jkaS2ap@7o?^8vmOa5;d>1y58Dq@#L46UhOFNG?c4a#A@+ zZb}2SgVIIqq6|8+SAF0}?U#URU&r}!G@6>6e2h@0^7t}vUPbf#EHj#gi1wzh4MsxhB88ahtfuVh{{2JiHcy|mw5+@ zwR*;HQLHy%{1_!MTi=cp_Hz{WdldGA6!wc0_LCI$n-un=6!xnW_OlfByGT8&cA~If zrm&x;u-~SzAE&Tir?8)=B+Ym2OyO~W!s7yk#|a9L8x$T#sB%m|Q$rFH)Mmrq#MwCp znRjh?c3Au_IifG(%mqXR0`~rbi=@k=1)N2Ow2j58iA)Yl`(yh5#F-*;@%Nx|L96}1 z@gJL3q%ozK%+g%AX~l`k6wO(+^!@);lppszwXP+SNwJ?V32tN`ZItm$kqwh*>MD|C z&BhGYPAK4DW@~TZ!0eOZ^y2or$T}xx|Gx=ZcPcN?#U;6^5eqB2XdzDaROjZ?$9&>) z+fC&Wm@l;^Hy>szGxFo?V24uJ0ouRezt>P^3igW?MD}lpKf^=5 zK&x(HyCpGV*Z(gJ$lq09)}_!zS3s21F2bjo zisHdizvGNj<^SdFTRt{?jzA#KzVqNVCkvK|xdTQ|pfLo+4HU2GyhQ&20oms>4ar@o zR|NvHPgC-JyihsMV??w5gK$5Ah0KHTG*OV(mmpEJForl=gW~2R^HCQ)ZayKdCbtDQ zA8|Z(J-GR3T%NEiH=h{uSqi!NFk7*tc)cpC_nTs|o=MB?UaX8ilHteUar_s-qon7N z;&8i>|0?KPaY0alL&fVyN%1(%`bWDpK~$<4II2&QTXFNb5dQYue5sg^>c-7Ss#R24yk3=+ zUwEbze&NUA@%$ITqqO7b3*3I``KMz0o#OSQr1Q>qa&QnUs>>)OSJ~t7mB)iYxL_Xs zq!j{O+D^oU`aw}6^4_15-dB{zM@>TGt0A6uk#S!OH#of+_xMX@HM9;lpF1wkpaVCb zCFaYPng=Q>{_^B1d%LZmyll2m^qBCyb33G={l%=m5lnX?`AqpS^{d zlf8qj1O~xOQ=7TC3C1&>X3Zk?O44!V-=;g7+RwINa6mqa@mnN7Zv4$@@&3NF>#?uo z*tN2cf2Exdd{jz*kDf!y3AC}@WcB7YAaZ;urw=zD@%EOaU^#Xu#>M0f;x5ks-><35 z%_pat@1Bg#fmV=v4%BTplAFfj3y-5f43;eTFp5)nUs)p=#l<+^(;nyJc-{mYM?|2z z00rD}K2H=2f&c1Fi-zhI3>lHqSup-OpHaDD|dxV=W4`;i_bMu)H zzG!YfPm&!&0!0f*>{>gPo3H4*OP^Z2UX|7JTtlp9((O|$kMra3XvE{esBOXX{vW>x z9;KcCHOB2%E45uBgH(m7vOmw1mH+Rir1otf=|>g(-3X#prKP3cKj*D9J-~?_yZJeFw#LZ`it4STp%|}L&)Dhf#6wcOi;pQv4b)`>5 z&$}vepLgx1AbGS(Uj=`s5+fk%caSEKeOYW)0zymgSBGFfL@b7Y=1u|%L6P67AP{0{ z@{GCdRUh+dm~!*s#tU+0aPw*7^0KCL^AW#jGMSr?pdB-nn~&Ik>NIXXEMBqX$LZB1 zrO>M%hsSIl4@Nr+o_C#jJPP5lhR3529@}_43gNL&N%7ENT2ZjSTz2PUj$%Bd$3=HZ zeTwpjtPh!BF|)B?o*Pht?D%6K;rRK_+a}eJgsA@t@H8q2eut2}bu9tjNYLI?xW8rb zL-G=?lK&n}j4o=Hbi2jhp|xdSWZC0ED>WX)J!fZ$kfxA%cG>p`U!o!-AZss36UzQM zI~CgzdGaHuxRZcFQ0eDv6IzMrz@<|E!Y)smZ! zpshWZn~&JP(2AQ+PV;1vf1F-rDTQABI6QKAJQ!^(c+Sq_@hF6cMr&lO#4`}nbteIZ zpg+!Dw0Jz29xt#LT_wdsL!!52^eXvtk3McUDyP5*#05bG4v~FHLkxhOD?y@Yp|MjG zWR2wJGs4+!_&f*UrlihUnBZz+q|SXv*KULp@b#60#)TkD8s$y`3PFpmkNABAYQJR!ax9+Q<44>T^R3xrsWHJqjOxiFG$IGQz%n~!J#E%o~m(%)|(){--x zyS&0uC|<01y(%leW2zK>hx8D6uhN-$x6#tk!ra)xTIArsyqJg?MAmjzA~TVrv6-oz zqmw;;%aM71($>z2`4&C<2BVp+jibG-6}rLB44p`|Gc%sWya|clDijD%`Nqy78*^J{ zhSiz%n$;NxP)h{$M(B!P|F6pk;D%3mOr?<$q= z`fWb(E zxVSPS-%qbJ@qW2->DBkcL;DxSqom*Ou$S}PA}%wX{+iS zKJV+x!G(Byf|Fv=;`O7ncu>D69woIG135VSD)v%RJPef@4>XS}CkV!7Y#`@7O=9I) zI6&cZZ)CQVyBCEt@pF(C=8IWQI7kVFsU?nTkPjFZ&d!e!?`-16U0(hTKR6)kcZ?(w zsFXwU5~<6Imj3-3BRTqN@LdpdFl82eESNPIM@w6Fp@aR-7F*FQkqfgpgugE3U~9$v zD50eOx%LaJZ08#@DXTeg=^& zBRSbP7Q`LUs>QeUvK5lwO+>3e%dj*A@7{G^tdLE&mqNmW$#D*srz%gjj#PDAwZso$Z9Tyd9fO-aRk;xi~f__Fc3ya3P7V=!N?2jRd8Qp5zu^Hy{7N!_M! zQCWY7kS_c0K8BfY91(8uDrL~-}82Ifl-<>nJ&KFdgMK7x*J7&jlW z2mOQGeA>8PgK%y>;=hd|xcP8##gZSVSD{LwS3eGqSRM~XI}7HWXLvjc;c<<}qYxgC zc{~c?@j*%PKqE03y~=Jqyi#>kCMlT#OL8gJyhV0=k;bE>=S?f)cB67gs|5maK~RB1 z($5Dz7SjD!Som}tb&+~L2%d%R)4<=i(8R^$IC7Uqevd8nBsX7SDHJbO_Ig!Isr9OO zJW6}+zo8sld=yV9zLDtr;n9r8gTc9A9@|za@j!mKydW66LCs^_{xAz?yTozxY2qJ% z@HxTFmxXoJZ5}sYF6PT};O2|LeC~L?2}i#>`F?iRSt<4E`{B{+7saEjz7sM4<0EY+ z$~GS^D?A1&DIRDAk*xDdJFwEf_wyBB(j_Xoe(-omEE=1&ih#$1$D_c8Wyjm!TR+IW zR*2`JF|_n}Sb!ZHGQoT~r?~kHF<(|9H=i5kGfd&;v%`FOr@8sWSX`}{+U}`={b6HAQ!A?l#*(dPV%1bBfhkTe z?J`geyk3?PU?smzo?hx>e?^)%+g>F34k0-|`~T`C8i!Swbux5m<8wTuTq0l$asExp@Gi7|mpxB5*6Tu#KWB4i(!I{Kon_xEO~O0Vf2vkXd+6Gk4&HS$ zF??`m`gFIj80VUu>Fo4H2jNyH`qP{dffYhK(OTu0lAJoxY%$yVb)pS>9A9NzwG(Z1 zVco)3SM_Lj>nZbI1nSX&723=`YNtoje($F(?X5@O==ImwAsTvg-X$M{>KBD{ZMNb) zLRwgfk!Y@vW=d}}yStEP#p$Rfq$@^_Z2R1|Bi--*)4EH#cBB{7XC>@FOOnky(2tj$ zx!V3yds^p+C@J2%Jw4^I29wpE-qe_FK|A`)*ddH4?P%3C4;((fY)kiK#N4p9Ep7jz z`S3rQwxvUxu`O&v-<@Qp)5ow4-9Q3BYkI6i;MVlI*A!E8YdTSvML?JSFl}jEmBG4n zCoi^wR`eXE)^<*XeX#~}OH=v@+uY~cbhAEe4`|a<=dlgbrdu+tJ0faAZ)HlkZbhYj3(z7HTR)@}Fx^euu+O&ZR3vewuhV6!SwP;2F%ctzTf>2>pv3bat7aOG*u zH1-B{I+sapc3PEAXVlhAM};1BBaPaY3kd>Ahr|*JI2|w_;QIiV1Grr9MD;*Asuwho z9AJp#f>b0Ym4oD_G*CMzUDPhh0JW2{MD3>B5e`%^!i7pjI8iwWH%cAhNNq;AQq2+0 z)FFgBH4ya!brJOoWsdrZdX4&xT8a9Rs*U=U3Pk-(bwT}3okn^y`@-fSlCL3V$J=@Mjc8b zex446^`63ZKg$}ODYT59VJF~NLeDgq=d*$sXS!2 z6piedaz}Pe86rEUw2<9X*~kx2!N@OAcF0dq`p9ok0^~=iRODAEPvmDPBjk4|ZRCfj z9ORd%2-ba>cc56SXZ#k$dK1QvQ4+KD?MPuiM`6E5VLwP=zer&}NnyW9VLwV?ze-_0 zOJTo@)U#?Q3j1XW`)LaMZ3_Ew3j1{m`*}*zeAmtt9tS8qE>L)!pzydsNd_9XE{Nuj zI=O&@I9&r%O-yTJs)cDIOtmr9!L&K1Eiu)_v<;^1FztY;5K}!&J7L-xQ)xK`zwK4r zyb>=z3kWU!ymAoUXC*To0=Z(%D}Njw6MhLi>fml9wDkI6avjZRgNb5b`n&3f#V?^B z#BUHp{-Ov- zuV0dNkSx4zNSbYJVJi8ZBL{Pf`D_|lC7}pI390OSKvzR^Nay|GGNlEG`({1H=*`ICo4lZ5BpF2bD~~{XlOjLLiNBp zVsK7aSFTd!t*}Eqylazv)5to*?*ky;QXT;w0Um+h41p>`)k0sS?t<=3#b3L3ItY%Q zi8s%WJOIP5o^a2y2!c`T1D97D91q^@?R#9io&uB4ecs`)@g%&^8`NTq)^!+n(^O}M z;UySMx$c{~Js!@hPMp$BBOD&2RHlyA^?`}n&vFO;vl-@BU;M`b$8bfg_K@hT92 zmwUUncKkT=o9P?hTpj@)fnN^+2ZM~Bi5Cuo-o5hmbvhq{K*ySIZ&p48x3%)LM+#_Y zy}{$~+6%{^{aTZqi`t)q#qhSCLr4Z#L>{wGG)ja7XMxTKo0G6$YufG+p-Hf)N5tI+ zdrrc#MLUiKj*bGI3SKiWh8=^$;mectoFie@y(%`F+Id6F7~Q%{9=XA)HRCQle6e1J zhMEXB`>a~Sk2AlX%JD7Y5#SL}as;lQ@8{ca+)fDjLp12wxIh>1OD%WZ8?UCK1qmmvg_2h*(IAyW%uybS6xDJLuy<`!ty0_watqgH&if zZospkxo6?szKO<*Rz||mSEY3_Jp@oEv5YJlhYe{ zP7y;|Rr-}v&s}iv`MDZPZ8zjoVFism>wU|>Z)wV(nF~!IaKr@53PFJ&(u>*}+%AA0 zSCm};c=34zcm$LWfulB)pB`v9AME#~J=q|NhU#~g?sE+HgI%*nz0n(+08iFE@q4{1 z7P?*-(ndEV4AyCi?$3F50Gb^1T6kI)b=Tw6^d*RHP zJ|^95!(egG!SAiBo`m=7BM#nMoea)x^-UeeoP?fJpN|^-_!x}y5TDyx^Du0V)tvLf zCjt)kKXY*F-4J;Fan;_>ZmZzfomK|6A3XEvkcfzs0rtV*x~KB8)bR&lZhZ4&q7z%8 zv8{U8&ePFgJtnQHN0k-)c%po&$V<#4z$2g}2qfEO-klo|2xA&YTEu5ZgYzbh7E|{A z4VrhJKYe6*6duN(6s>#`18{rhI1ucCRd;Gyt)1fpPSZyRPF}VJ1j1e?w4eIJ#_^Fh zK6fMGsP)x#D#k})@l5wc!ox98-ZW-_(2tDI_taIH6quplgYX3YEeAB#BdiL{%%X0_nEpZBh%YwCQ7M&Nv z-fK;-UEV&4A77M2PkAwU1b75~J_58^Ww-9>T|j-lrok+Yt*}x1^531#C%}d|(>~r! zIRrtkPsWC?V)iFqCLa6SEDlEOkG<9FnIANqbNbeuun4%@xOo+E;u%o&-}yCbO$Z!0 z_+e0`#4QkdYfyYl)E3y(bZ_%?w@vx9l6piIsPBcj^XE1mlNJxQ15RA(+V?n&o*62r zx_&Qw`1jx>M>T(#yr%j&+nE6{#;)B<`cnXC)><^MZY?o%9J9gtXgx7Je7LXtlek2v zyY7TR*T3Uo--)C|kBULCub$J_2hR%jL!90wb(yn&3B>oVF|W71mkj;7;hM|R<^|6s zeJs+yJj!+C$C00}n!E%&0z3i=hCs%>HmfXM*Ffah5aZpKyz*)9_{lHZjPZxn$^TYc zdo>ZRozOZfzIO_yIUip(!E`^=Q8$VR?e7Jt9`*EViz6UWZ|9*_*(YK3+(z4*HQxu} z>ju0_%uIx-{(~B9d>fqaf0M21mTw~Pgm+O5DpZRupw&ldK3LI940IQ>8EvX1Lc;-m zNw0Jw;FYLF@B0u77ly4pR%z}DI6k4}sS_^);KPIp+Rxv`LFl!)YIWZ-^4!_iB=&VM zJZxOCUUJJlFwt_&qv3}8U}L8p!)uN?4Iy@$vNE)i;hE2ZNm?FPWN!Q4t2 zy>hPXg(|scY95;w0rOj2zL@Up2i0S}dj|A846(H=%`T3Kg4Qs9ldsni$Qt6fH^V9d zu58)9(*Ir*_-XfEbMI>+2+Zg9JvJa72Kw|(JS;v9qWYQFBcFx9i!<&it8N^DNe;7a zsoe>K7d?Gd`mdh#pZ)K&UBV$rQ&+>?ZhjZ)UOoV$2Dt67IW`eSG#)iIbaNc2+byqm zSuF~dk6XDqE-4ybR1V+XuudWjALQ1)&-!p^;pC9fVUZUEK(#Fo9-MADD!lvc$#Ybhy>e*g$>BZ%~>bVW)$;dnN$f6o!8?AsgR}Prh5pR~Eb+q+6bl=OI zKjtqU0UiM!0VPAATV|!)oYRg_IqYUq&*MA6Z^oLMN5Z^iXiwLg2m3dkp6};g41QAm z#p#XkR5+=s2^@xIbk&;3)XJev9o-2OfdM6}{^;c@PbcH2ZhS%!r1gRS&sE zUkrpr^{I0^;hxf74B@~kps zIV_qnWAVl}8)0R;ffv?S2R@IXD`oy^8?D?Bsg`JP(6;-0lZCXO2Uq z#_i5tn6wYB?t3}oz^y&dyUKL)%h6%balk#hum(q67yzI67;C>H_-{5Os0VwAocpL%7?L+=#~sFs?%HCYg}|kQOjC_H^=gnRd#< z??vRXc1>FK24S{p-G9ezg!g+a-|ki0ZdjMTw^zbHF>iX;Sr-J` z^HitXEQOgz3((8W+XjE~&8K}wz1<3yzB2rzpQEma z)SowGGeqSau<4N=1=RdcL%v2{gRLE&_Dvjf35Hq>{ZP*+2G-6?7qEgLX?=%(SFJw2}Nk)flWHUBtCe}`PZr057$TODioufPp@ z&1pQ%%HbLedbPA!|H@~YIW-m5F328d(lHsn{M&Q*otSuNJiNxM zr?(Q}!GwLsHU-7Nj8$*%F7=Lvk>i)0Ib{+8=I*C$&x+DOZ3eJ(F~cM5I| z-ST4drxRdrk=v!Kg%}RKJG*hRFaqo}46FKvhd{r>GvYkAg~H>5T^F>{2!z&MUKyB9 ziO#3}Z+$Uc{4x@D_zqV&u*O4a^`ki4_|7kNYN!KdQC?AR%np8pcl z_|5&vzlpU2bK^wU89y7KZpaH3f0~We0>N;c+4B_KW<;KxGs3K zewXz=KX|hH+MmlShJ$If7AGD%L_wP>@uNouEiCqac{+Y~e;73~_@qXQ zMWyybKfVRrHICa9F>Ye6EBGGW=XEJD?W~bn`?`Vf=9y{UQ$~z~|Een83M~^jH zo*Mx-b)M@iw>t`nX>%ObWFCMjk)tNJT78Up59Gnz$^BB`h8pXE%N!mm%izIX8lTG} zz$5S@5m=lTn)Dzs7FHRqs{)jLpB_C0qQS53FQ8)c%O!nv^2F=3_YXEg=ay=BZ07EQ7hzdn zt7@fycf0H1o*fRu#wo%R#y&Byzj0*fKbw+Zk3geQmEnKG;?T|?FYEb3|Ihn3`_wxP zE4yu)Q15XPr2M_KTJJh>(7eIT{SNeL__u=T6c1GzF6L~h?%Qa&O#S9v>#nrBCj#|b z*2&t27jLKt+}G~>lWV&5#dQgBe$u{`{|vt!ZdZl;ln_q z=O*ggVQ8{?ccTpX`v-xV&drBY%i%Bm^8M%a3(ThnawZO$+%ObEJ3pJ}XUyyySLicg zd)WAVKR154Ws85hhCodiyFK$$2!vnRuyEk*C2*(q$XSD?9)O{t2Aew8OMxXBfiKO3 z=U`RBg~Txzk3oaCsweH1#lg(&lPs2~hrxZPFMIQzorDaR8}<7M_QUz425r|f`|$*) z;kUZi*yS7ogVt^hoH;Qa&fcls%YJM!%-neRa{ts6*tNo`hRLSCVf30ke&;e`!D4va z_ztG=aB^DZh_Q`N!pkjJU$yXyh9N8b7j76A1F_8}Zl&npLyA$D^Hm4ns$NZfRs=_AMfE!e zFP!T3h(4PIJLmlEcf>gpoNeasRC$;H_iGNX@p|4_2=_hw;PvOj&~upA?32DRFgWby z(r!)TpyHw>Cj0W@;f%%nE6vVDz^qxTUk2nIhHVvgr?>gr50=}wY%VQ_;+!9(Ak1b?3#(P-wqd27?Pk564Y0#t*wPg{;S4b48? z`X@l=08|_CaOgJX`ABv0=o>BvePrl^wWGTS4B87F+l{Lw*4YkoIvH;XDi;AJ2Mp9F zRF44PM{Ro>99G_Q6Ut`+yu>^LJObYffdBop-B_ zT>q*^&?y)f-n{?F3x}X}NY=Sco^kN7M%+fK74zO9rU~L36?G1q=06`4CwJq%3pZYz zy&uB2x0yEF{0JLs`|B|h+#K_+w>QC(wp<7{~U zA@KFNKd!>6FV&9epS=g`hKJbB?S2t@Z0xBuYF#)a9&g?1<{dxyvSEF*J}G-(=E9r% zy@dh!H08SSONT|#5Rx1n*Y(p;=sGa1YeZ5Mq+3;}8EhL11GeqFx?;_7F!S;`9+vSJ z44m|_Qun24(A?cawR~DMWcFNiexSSj@2SXZ_nxQsD);N+3XxmxM~Ex$hmCU#+CB_B z0vo^-zMhDK4%2M6UcMs#ys*4-`E(uu9swSKAB}*hdC!&e8id07l{vvP??geZ>Wl2J z_Ya4*&KFKpZ>1?S4vc8kaX{m7?qEBn?^2^-N#JocY`(emaT&hi+A-tPCY674qhWJK zLl2uru`t*E$`)&tFlg;LWlo!BM?rg=)4NA|PeSLJyDx{E-hm0Oo99`cxd}R7{Tl43 zbQfxkTbs3PT^77q+iKym?tekOYh4DsYk3B2|9EY9J|Pv#C9gP`9C;aP)tb||QCKWQ zI(J*!`cVw5u1|Hy{dg7}n1{)BqX zo*afr7ltJbZ5RtR4z9m2@!kIaTuyN~^ZCCA0xE)X0?9AiolCwjU=M)Um>C$Bf+ zpTEkwgZT1z1b75^1QZN`F_wA>y7f1MPTfi^(=PhKS2KfGBaiHaZ9_5}CT!RO?=JLM zS-tiFn4f*|>VO>)U_ZK=`h<{RP&?PA+sBl6$kKjM@4D78$XmI2WvyF5(6*al#cI>{ zK;|9`=L!bv;Y{Y5*`q6O0F$#GX}9b{LH($IKj+2FIce~6Z}--YA7>OyA9!(i1b75^ z1WJm)$B(e0-K}tV)oYjkV~t>VH6r5Kkknx4xcQ~!!|Bm5IjWwem(CGbS}*zSv-|?p&{SjMpuB{25Hm)%?vh7tuzJn7OAlYHm!Y91!p%Ob3Vx5Qr0wGQ zc?5U_cm#L^z7GP|&-e3f$h>!d`u>GhzJ9(ivbS~L_kyFa_tdrnW7mtJE;TL5_4A*Q zG~lJ}ESEhH_iE_9ZR-2rz~4SS8ZSKtm4tsz9e(i;Y>&E{cu@5am~K#6e&AdmHI5{+FP?bP9 zw&QMWH5Ki`9c*(%UZs`h^F<+9!|P z>n{!f_Zm~w)_qt3$L_Q;u>IhfPa{?bFF##>H{4j)Pv^-HKlsaM%u1&(UXYv-mGWd? zI9QKKtLjnZzkfGVE--usj{uJVj{uKAkq9K)W!{||5D4>(T34;{ApnkdpK>*XqT%wY zUMeH)0^pTV*oD^*Lm>WMRz($KFIau2rq$XxPT(|sbl~J=%zlV>#nscPZSZb?wV9Ti z{lOxjS?>jhgTc$9ed|YanC}Rz7iH+jF~5hSK4a^PoKT20y>Z{i*dN-t)!5eQ;a&&~ zt=Yz9q!`X#8UCmG&LB|nny@NVBLeQ$uI?~QJrt^HR(N${;Q@HvJNa0Tu^XYMMM{m> zS{r2OZSN)-BS*M_&@;ZrI@gUb+HJLBOB&ZPu&;r%N@t`(GX# znC1r_{U1$uHZ}mBwpeF>JSqw*1k}7($BX&y<{yKfG`|o4bFzKT)SiA4`r27eeYD6I z`i1R0a%(a>HleWUD9-oDOP2UPu=Ut93 zW6pc^nyhl8?T0W(xVPD1`O=-Rub$J_2hR%jL*~Ct>N02l5{Tx#=!N>onGzR8US_FjUqz(dx8CQ_c>W!fzW%6$+JwgBanAVV_(fG zyI@kk8Xeqn4ngYJ6q@oZcyB-1s&4rv0#C?`YEYqCbOHTZafH>xdfrg$`O%C0t3<*1 zIj_%YuVu!c9TqVZ6$pNT-Rf=241>HT<1eQhZiR{SQX)GYKL#F|LFa<5heG$4i39FG z_5*v3kd5QAncumOO0zxfzY|t%$cRlDco`2 zUs!FqUVqB0#V|DUY}g_LPZ=tHFW&9zx=ilMuHlwo$Wr$kTh7}L zs@4-ye;)Az@$oBne5@m(wdL|2A4Y|Ok8bmTc{7i|a{Jry@yYvO+RPZ=gSMeCXK=%` zK`Vkl=x5gB#wnBk?0={25)Mh4x*D2zwQN@1%MWbZ_=ffWI~o?+hB@0ighNv`ovy0w z1EGHH-K%!r3xtE8hcQb)_euO_tf_e<%v*-`bgg-?f8*)-{_(}& zC)Hn^-UyxIP7HOv5eY7HYOHD9ZWrXfhWRHyZ-=keS2|Xy?g!5z#=pH+Z7)==6uEl3 zQ6S7;Keh72%WL7EwRd~EGV4#Ke!7}(GJ;`rgDEL?M}i=?M|%0RErFmuKtE)4`aU>+ z$R(uvyInAF^4SJY5ATEruMIMWdu)N8aq*pez1QcO9lV9X)t zYqD_X;rKWRsMLDx{JuM3+2!;$hMz=0e^{qxXmA)>xMrBP&^rhXJH(qU8Ib@to@{t( ztm6+hpKgTZW*vq#W5%U#ix$J}P{-`dSF+!OBlTagJ@meIzjboQnTTg9S7tqy{T_E^ zkET0cZI?gJSQ~f;Zqen(8G$?|^4DwRNjkXoyMRj1qv1@|9GyO0`j?3w*5|G60!p6$Icyjciz9B|JrticiZ(x9tLzgF9!Rrh*LQ<{cA z%HArUmvo7SyYtVs{WK*4vh-Sx^o$RJ2DiP-P3Rp0jh0>3x$kuZf@=?d)qDFv7|_b% zbIq-B5ZZ3ZfgrmO2z9vDU2Aw0Ob;|nqIVyGJweY#n|+9ZrR{bH`P=RU(Y?Gt2dhmo z?UaYxi^yZ`nzZT-!fe&L|Bl-T^}XL%KdRyh@AUoq-oD@m9TP4uZFtv3ro3JBhKwDq zE5KlUX2hxf@o>Dp%J~Wh{b79Gpp~t5t(M`JS3jT5BfulTBOr~ya=lLXH{bDwcO#v3 z>WR0*(yp2R&PZo|$4tHU)g|RaL09eWfT~AA;OMzY8COC}KcDS0H>^wF+biKJRk`n` z7GmgE-fsB1Dks3mdFO|YU7{d$;N_H&TbKOb{_Jlz_k(fg{lF~tM)|GGd+S2AUh(1D z5umfi_mY14A-M0+<(bvUaOiZ$J>^QqK8PK)>CEnNevo~<>)nUFd_Zm2fNRm06XAXD zxr002i-zG{+nLs>d=&N!Nn7-^c{2F5Z2iF_@ECLs*rLC5br@)#I=Z@cvjgxh-r@P7 zal653m35s;=5g@i?fr}KW`6L#$#RDm%R->f_xkGi2-fL?CiHu#foKJ7c|?N+dq|Gnj? z>ml{$4cQD)IR|Wdq(?!^R%tL^`w9$u)7sK?Knk3*)5^OXwI6yv+H?CKZ9h1Ebhl7- z{LcToAEfn+kMjud2>dz-9LtEeY{2-Tn}&Dx^tiT1hK_pH{Np5j<^7$JYO7-n{}s4_ z|4oO9JCZKL+#Qx*D?f~dKclPnnIVh?&BHemo83!>^=1pYn6-$7f!Df}znvHXuTPI! zUTIDg^mDUcv19B(XnEt5M@Hg4m|;49j>fdTpeG)@=J;nC+W8GOS^oYoETMG%^ocwJ z6)m@(xa}4LniKl10*xRz;IzM8g~~yYy8o-eh!pwjSG6sN>h%rtgAQ8Wt0LM5!=*Fc zLsb6R57FtzwBFwi1j{xfsn`VneEMK-Tf)RQ#76JLS^zUt$`+ zxiA0sxsnDc6aF2Pz{m1Pq*1s-z zw0@WMK0oL+_R^rS7em0bT8k5p9ipI3mH5#kgBBKhzr6lAwjk5^to3^MBkYjHk`+<# zNbNuymx+mRJ#47c>sH6%?$m)F1~ii0SCLnqbh_@1HMb3uUBUUP+pXUI0ia#+(U*}A z17MovghnUs1i+tTCtn)n%sj{P7momsfYKmvnT~LIxg!!jQ%~FX8R!qAMh2hMXtC&* zvE%=>cP`LSZ+je{@#y80Nyp>R!$nCdTtuQDl^e<QEUp%uPiiA*w-P6vn$@ zGG;tO)F@**BFCeIO2@l8Qo5x=cQB^4a^>8!R_<}ne`c-u&VK&){``OIKl}gR`?tq3 zf^M<1?-;}tTnTTe;`oVh)uVWW_xhn&2#7vI3DjjlI@wR+?b=*8E>MOQ#9s$N9S!tS zt8y?T?RlDwVgg>kvUUFAY{;{AzjOsf0%->wQ~SS1L*1`3F^Y>QFbD02FTR@rwfoFS zK~+hh{jfK4WoafH)KIv8GNcgRlwQ;zX6FHUWmhZBkq+(aUEeDfXTv@cE9*6>3{K+* zvA~%RIrujrB~*MeoM1D(ZqIU@a2^KeVL>l~$dE1Te1H3i0=OXWwsqSM%7lDV8Q=q) z-TB#^eO4mPI#iYUNYJkul08@PHUtH{Rn_&OLWTWFDuJ_4i~>pAMbd;0`>vWwig4CH z^bZ#dm8--9-^z2wM?NOOFgrcGq~S8iG3&MDx-P*yMkUe2dgAl6soFE1F5R@^stG$7 zP936aE#8pG)4q)o;snk*w()$3A_Ncue?9^}0~zI=WknEZc3RZu)yx?W3gwCkZ#hi_ z&7HjpNA<}Ql-kggHt%3G@A@<3D;v_l@ds><_%$Q3H4Ro|eh`tqRtD9gdnJGNDuO4^ zPi=GhF$aqDMfVt7PlRP@`1%XC62LH)U}2>jLKUdX|5*0B29jVKqwPTR2#(|e&hI&>=usMADa#X!8NS&n+l+sU&Q`&<4wpU z7&_e)ivzVGr~5xXDFg|>WkwzuIS?!BI}ph6d#7c$c}f}Pf^9NmSlcuY>boO_;}-i* zIQ~9$u>kkB10e9YHT)c3?4+vMzpCCNKb%+Id_|yOGWjNmt`YAQRQcocSBDaOJ+#Ot zfVyrV=943trys7x;}kPhc*h&(g(yM*A@Ki4z~JmXZOo}Z{@mhQ5m9EE*__-R1w_@x zOpHtxeAr^s>mU#V1}p&$n_niMKS;qGk=v&04|gV7%qIb zc|EQK{8Fxvb=|K*?}gqFrL8%jrD?GQtB?YVss@Mlv=l>!HcQW3fdYe?t4)u&e4fwt z!dj)ao)`_3V9SmYz!2@N zf?ZM}*kEA((#wb+{Fj`NFhT$!Fmnh5U8=p@FKY!^;lJGKWK&@t6V`4COM;*zOhg!e z1jsIMIXh1vdQxr3>A0wm%YlmPUV;W19KSxCiubbj?%*7H{8>sf1B|=1R)%u+Rc!uw zj^KGyW7ufE$~LkogQv|YFXp1MC!hZbFOb)h;rKrq2CG=!!IVK+U;Phdhi^m8(dyvR z#v5=M-%tKQi~-t8>Q%Ms6|n3(jQx{t8aU_X2eOzusSKNR0UlNyS5Nkc) zM46Qc>V_AeRbQh+DRHY$5;hSUdg8Kt?qA!ntf{GxLbsM_3_%HOh{-zKO-8`z^38IhON5!(6#V7 zQCB{Wr=@eWpF(Hh8REb;OtXREY%d6x8%Q{w;;c06Vg*V2JID0 zuZNP!Fo^ZaNb#XVwMc#Ta#TL7wR4NQ|KRie)l~eDbc6sx03qwNG140A3 zDUMIGVUe&mzIjs?C?2UR6_&?L_#N0Tue5ooMKHMS*yw9xR1W7KX8v;@u6RPesrs{L zFarI$F9-HUT;F78TL4o0(P5%Fsh|*MdvKL(0hHfvqna8OLl1hFgD|=pif))5s_VW5 z4t^nhVo7)5Y~=jf`b$jk_IsJP+@}QIH5GO1Jt+k`CPJh3w<3rxUnQ|!oCR2tw*NwZ zI=l}%6!Hp1gL1lp=8zu?tUM~!;>FUTMrOTPQcWDt8nVgNhcuSV|7#FHqek8+lmxfC3sIQSZ+*}CMitXLrk3XV`Y2A$ zPhY&a#ZaV|WBKNcME?-b;RM|A``xj4tmgqY z*U==812}@q=OD_B;OS2AM4=q9u7~V#SXZpme%yXn&rw?)DKKWGBNiz*7B^zVoEMsF z$!PRg1UH&jj-Qi061PwqP-C{Li&7IpseR%fEzYTJ&S*LxAKz%vIo$gA(d_>JY~jIU z1^*A?A5JZN17`-9jV!hkbeB;QqjtRkRCmu?S5zce v>YivlEcba|+i%qlq+WyoLI5Fv5I_hZ1g4L`^o{FUXT-Q}`gnBD2p<0isWI-p literal 0 HcmV?d00001 diff --git a/test/regression/azurv1_census/input.py b/test/regression/azurv1_census/input.py new file mode 100644 index 00000000..6a09f1bf --- /dev/null +++ b/test/regression/azurv1_census/input.py @@ -0,0 +1,50 @@ +import numpy as np +from pprint import pprint +import mcdc + +# ============================================================================= +# Set model +# ============================================================================= +# Infinite medium with isotropic plane surface at the center +# Based on Ganapol LA-UR-01-1854 (AZURV1 benchmark) +# Effective scattering ratio c = 1.1 + +# Set materials +m = mcdc.material( + capture=np.array([1.0 / 3.0]), + scatter=np.array([[1.0 / 3.0]]), + fission=np.array([1.0 / 3.0]), + nu_p=np.array([2.3]), +) + +# Set surfaces +s1 = mcdc.surface("plane-x", x=-1e10, bc="reflective") +s2 = mcdc.surface("plane-x", x=1e10, bc="reflective") + +# Set cells +mcdc.cell(+s1 & -s2, m) + +# ============================================================================= +# Set source +# ============================================================================= +# Isotropic pulse at x=t=0 + +mcdc.source(point=[0.0, 0.0, 0.0], isotropic=True, time=[1e-10, 1e-10]) + +# ============================================================================= +# Set tally, setting, and run mcdc +# ============================================================================= + +mcdc.tally.mesh_tally( + scores=["flux"], + x=np.linspace(-20.5, 20.5, 202), + t=np.linspace(0.0, 20.0, 21), +) + +# Setting +mcdc.setting(N_particle=30) +mcdc.time_census(np.linspace(0.0, 20.0, 21)[1:-1]) +mcdc.population_control() + +# Run +mcdc.run()