From 905d3f94a58143359ae118b85c405463c25ab056 Mon Sep 17 00:00:00 2001 From: JGonzalez Date: Thu, 29 Dec 2022 19:32:38 +0100 Subject: [PATCH] Trying to fix weighting scheme Now particles can be split to increase statistics when they enter a smaller cell. However, this only has an effect for collisions. Still, I have to rethink a lot about this feature. --- doc/user-manual/fpakc_UserManual.pdf | Bin 181016 -> 180981 bytes doc/user-manual/fpakc_UserManual.tex | 5 +++-- src/modules/solver/moduleSolver.f90 | 5 +++-- 3 files changed, 6 insertions(+), 4 deletions(-) diff --git a/doc/user-manual/fpakc_UserManual.pdf b/doc/user-manual/fpakc_UserManual.pdf index 66aba8b627f80119f5f78f25fd22176c034d96bb..0302caf0bab2667c1bd94d4932950b9b72d0c24a 100644 GIT binary patch delta 17550 zcmai32Q=5;`!_P8Y_dm4NIrX54J$K|hU}HSQzRdIk5pE&qiiCZ>=g+m4U&<)E2H(l z`h7>e-`{ioIw$8m?!C`*pYeKL_jxCU(0>l5FZH7)x?D#0Z1+1i`*nm$_S|bu`|0q4 z^Ig=5s4u8qT(u}isNMPOV=+f6$H%S7LO9HKOt_ruev+3cI-Iy&OU*})DIcL`lObu` z=Eb{xNglgRYTg=N8oxv5e7tJ(vxGzFt1l6(V}TVkjY^d-r)^+rpjBnq>62?txungk8D}_CAM%!y~y_(IyQMZ$H+!y2Y+pj zL>KAix#{qvMY$m8-jUqKg-w?NS7UK~X%%kbRZ&cRhq}nXl`Q6Oj2@S)MB;;zxgP|F zyzFkL&1}(LDhc+hE9_`*-t=ZS-*hEOnIZC8msRgJvYksFee1o=q`K<1*=X!@+k|Q6 z#qwje$n(u6qu&;8r;xU&#}3`Rpt78nxI`Me(WGo!@_kv1SHSw{)zCzOv^~~`+1} z!`R19-%j)ZkJ7!ndrOs1ixZhzGOHpT`cWzQ%4qScgutSHs8e5?e_EvE zWbkG()#qj9K~?P))RU6rbCv9svHILxMY4OND(X49BeKMGdXCa9+89RVBIEpT2Ab@5 z++StbKhIaC`Q^tW)iWlFVtc5oPf(ld7VS5VSq=`aiqmpx+xH3YQ=8skQPFuRM~t$i z=H}kPHIjvDX~^l@T#+8AC~Lb-cA34yffJW?^&FcvIVV6@LM||jJr;axHwB%2^xQ-_TKQAboQQYYCV7(mdl-fkl z-m)G_pHv-rbxxB1*0)prCR*pzy|(mn<(r~ArE~;`QxoN;qViWm=6^KF-#tmfl?Yj& z8d`_zqKGRp8-v)J=SE|HT0SM276&_TuguhDcPd)MK2`n5>T){c=t5XW|14FK)n=#s ze4ls;NvQ0eVMWeNJk?L~7ex-0UDUN*)fTd9#a7`|HG`K=d*<>xB9H87dA;vw8gK8> zbx4$W|LFX8UdM<%S^N=NS6C^>=CC4x?{CySjY>#rD$;7RIDUP4^c&TEcULA?hCpG; zMu*osqeI{ONRFYcUX(*Kfzn$cb{;QcjMqtm?DPq}^b@{rP@bQa3gVs3nxyOFOW5-pTlNDp$QmvEnD;h+4GBUVabIg?uY(-jR<=1AV zr+a%8efi(87-3Bu;}c?+iujKzI&6iPm)Vd=wK;rzzG%&S+}dC}+&LzN@BY zkRX5Jz*}hIIv_5<7@czJI9a{$nYUw2a0 zODkq1&1&l+&WpmC1NH2@m$f_ss^_TYScG_*<|q_abdt7+*@z~6JzKFRD<6&@m$)l*v{10Tbacre= z56@XsPq?TueNp#+W7*>Iw3QphT#Kkk7LwItC7ncuixT}4yG83hUmP4P2sy(0!fMg} zQs`G_I=N%~@r$vIy_t$%J9tY~Jq+Zln^|ZBKk)BS;xHWkYEjQwAawkh#pC3~`^+8c zih}ghFMsfL`wf^su*fP67^0q93?MVN~eN%Lgm0<%tv+ubD!H*7C+7T~xq!mV0w}kn-Ii(VMd>{)ptd!MbNgA-8>Gf+>;WysG-=hrOQ{quhE= z9aG`P+I%%jak9y-Ra-x_GS$KBXm#WFDH%OQA^ws?i@5W-Vkg{OU{-DQsl z;W5eRp0X-&{v;J^ublhf5hk*;H^M15sGbEG^?1hQ7`KjiNjROuN}Cprc{oS#w-X{q zs|kfRrbjOKWVL4N?S0(9t zy~Xo&L$DxoMufc$Ye%>6op`gz$L#6-AH0xvf*$2fA`TEv7K;XS2)eErgsL45?bd%C zdGpA8(bngbTMw7sp0&}`7+3cU-nD2tu90DFsveQ5o|e$a&29#AV7~ z^UN_rPKOD~3etI49%W~2)`FJ|-ZU3y)uJIi`fscDggOGx!f z>kBdtS=(@jW8tMKwykzTp0rK@-o5b>CCQ54*ly_f+9#eA^kp;5OjkV1=%C{BmZ2dc z0U3Hz#5XXzl-KM0#Z9-^Gt!Ti3vDV4l|vQ$)+@>bQtx_N?>(o`vZQnK`6AcW)9|Qt z`o!|rr2M{x#`nF<;;pIeRhR&`$4)&J;-)$5+bxrpTRQwQ-0V);FTEY^D{Ljz3h+IP zz0Gp>eAE7?^J5=UW4`5zI|VTW*<^)lapmtz3ewVTO?}dMr!YjYlzkklEx{OD<96mD zO}Nf!HV=j1v>rXnyQ@rIJ?R6A)H#u>`~}V>dXc1Ri?AB|7&*^M)goN3R}K|ly1s{A zZI*mVu|-Ixq>fbB;RGX_*PJCQt5i0wQ#hW)p5y|ZO8$mVrEeGNOR=YYBy9<0r>{}A zeKPCWE66ceGu!&3;z7%eBPE?zj=#%keZ6DHPgUyoC27W5EXBS}{W)Cx+-zoFvi zbM`fVFnn(3L9XYIAZjY#@|-n06mvdOBx~(?V8!(iuWNfUQ+(+Ym`xQdZiH=5q$eR| zUCP49O^;p`u-bEq{N)W>`QF$D%SQ&dXJtkj(EXR9cb~0Vp0d6a|7aC27K$Y{lX9A3jW)26` zgvwVpy}@YkDq}0!Xfkb-T#9h!w~jvHi|Oq;qoqn(Rx*%WeUrbX86~t()~zze{=t&FZTv&o?>cd*~D+0Y~4rOZ;{wT3WUyx_y@4i*itoSjU$n(ngP8#jUo;8^Qvy%Q> zEvFZ5AG0mJwr3P~;Hv+{&t|2SsT%J~a4c1ocBcb=M$uhOI$~wYQ+eX#i^};6qldoN za?VLh@aB?c7B9R}7-wnET2cynZgkvmFSVo2+_PHK!+4LxVq}S+)`jm>%IbSMIG@x6{LSMAyoz7uq`ZoePJmDK|W@N>cCtiKff0{!pB-U%yek#Kj>Cn^pk}fm0 z^lG_+>xWt?{x*F0XQdmoXR^Ds&!)VvKIb~1aIYY_I>-v?PRrZ9-rhR5zcnM`T4F<` zIPN4(2mA9>Qrw!03RP+JJ4{oWrcuDEmz!xWwiY=fIMYz{tRVVA zoC=oh&5Mao&Cf>erHcDF-+Yo(>t{+jx0zPZ9&+%~I^R1T#`kR^-uH%^JgA@Uk7D7@ zJ%oAa^SprlhL93d+0&1djgDAJT<0AhrEdlrZ(8rXdia_|_s$~LTYz%XPd`HIO*VsG z%wE5j1z|nosOE`mXEtZnX;X4?w35qH%Auys{`r@QTNQh^FHoGgX5>36n-S+CDbunx zA@Ti4^vUFOJzKUA+tMzwQP$|1#u|p)2j^;v0>o6!-#AKPJu};EVjOJi=C0~bdM#^O zyeB?e_$F>_Njjysujg^H$t^cIlR(?k)?Qx88*51o(s>Vx<2;YE6l-{CPz{y#_PsvW z9-RBrI{SmYG)h8y?B@#)K2^SbzHZ?u16RqWCyXeaqjlTL)%#AiDRw+k=lW_pLz$(+ zbfsA05~si5A>ok=uOG`Yt130V^f}wna(N&hP41=pTaX!~klRX`;K5U=ruZ&%;EjoWOy`i~U!JWV2a$Aw!H%5sq8v}$xQ>`TR1Oh)46VJDNwo9DtA}MQ=+Q=~hFM4NrtHq-E z?|I1ajff4DYM>upxPE^ykB0LKx&3MGd}=B+WKh)T;6@CcNb@|8?B~NF$vJoQ^xPj^ zG(e@qYx^OW#80Yz>5S}&3GQJ=f}w%O-9XBA6q3&4UMsi3+w@D%PUt(n;m7-b?Q&BQ z4=dmjS(^+g+CzHzi9Z2Xi##ZISJ`HC)5n{o%C0Bb{Pn|1tq{wx6Y<);Ui_?8qe%U$ zF)#EUnI?I!)HUspm3#RX$LM`D^W;CPEgQA)GbYqcJ!&kz-{mX0V`WFByu7%@>u2HQ zZboz-*@Q>B9rgRL1(q{Ko^GSddh>>w`FXE(>S^y_J36#U?5}8?TC-X>jboRmjID=6 zG_6EF-l>?fNVV>pL(U8Zy-QxcW#QM3-H?BilG01Tm~P)O)l}U0P-Ql0G{&^tbJ&oH zLF@FH?9@;re;bu%n^KBxZCPifqR%oe>BG-6$lngx-~91pV(;7!YTx?kJA6#>2J_D_ z5~<3OXv3QeBvQY}OpaL%mghiMVo=J2{ZxD>vtCoH`y)H8PDRdv>AOvNem6|m`+U8mx9In5>$EP;Hr`4rydGHZ$|EY$B$GdnTptlTKerAd8n;w61|XT zw_J1}f{;wF(Nwn_kt;2l(=9S!I-V?`V%*8wXp!`Swy~VOF*_*vj&F{K2)W%sVY^!u zq>eVuhv(xI&IQ(=_(0L|VSBndv-$ZNwe-yf9d1#c?+@e+hnh3Y$*9V+^dG(-KD}~# zMOxpzdX^EByFoH8yZ(A=K0KSMU!$L1hb)noxbt$1QZ?}N)Cng32y)lx7rq%BS?|?} z2Q)&rY%N?P^jydTT*@b&Q$@TAocLZmYB|I`P@8fi4y&ChUx{=}5ZCEy8VW=NiRRIW zUzH1NRY+JK;23(2j;+dCorGvSvn z1NEj|QKbyNf<~I3?azrAJ$Q3Z^T6f#Lu6d;b8i-R1V@|s6{*X8*G_sqM(j&zi%>|^ zEYZ{Q(mbEwCoGz180?^Peno@&)*nY=8w zQMNtfyt!p5xSXnNadEL@V^U1d&5uD#p`QD7>=sRZ+~o$xkAei$S=zS}efgrsjPj+8 zvwIfF`HSv#`31Dj^zAeas##y7iH!V+88X&oUaIEFOPePlq$Hzd?v|KIZlvvG$v~>lX~1mjjjWd7fNJO*%TO9$|0|KUMwq z>a`#v?&hL2lcly$Vc~B~Z#g;0@y-17c-cbx5QFK7-rP<>4&o)+JSxhAZhXx(!jYXw ziw_Y4Wkpl+B!OdYj6OkCX8h!_1}}JjEKI0Iy$nd1RiO@zl8r;_TE7W+rlF2UH&g$> z4sDp^o+Q3b$`_urNC8SFeCfvm$O0ZJJR#_J2G?)B#(rsRy?qFvAyKx#bYDBI>3?LeVlGJBszOX zo&3wYmims0m?@;55ZN=o##oXi~Bt zIlNhmM7c_p3)|J&ein6_4BIb$A>uP9O9FZm&4hW=apmIUhegv?57mNcCr*2+G#Nx4 zJ{^k35*^QmCOjRd6|3I3X!89a>qF^0@4PD>8eF3deCJ5#C==BPl?iXg`x~nB3i8V) zRAxJzoGOS$H&VGBa}-4Prl8x_Eq2fWqTleBw7-aq&$S$URDaGtX3(e9_9@@W_e{ut!l*f4nDhTXYF|ZmhWfrt&Q-xJK~qTBnKuxzfE;HTE-G$@I-^v zC?!WFsOuWXiw1t~+Nk~o*6-eLS-&x~x+h&E`B0u*6QFG#Sr9U3D)3iuoG?f~F1*SR%lsu0b)aecK9{vB zvBIS>Ycp+?XO5(VlXw0ZBXI^Zn&FtrpX!8*CPZkxLHOxY$1ErgPAbxSC3V?f#|8G0k)Vh@*ev_ODz z?!K@BjbzoIJR0r~{TIo;d|&9b{H|sfx4b{-<@mWai^20G{&L1u?hPVh@KK`D2J3+f zi~Y|VKlQyBXkL&p^pE~rNH%wqWjN`^ik+iWrSKU25nkRWY0r-3J$YQ}+Psx_c;nGL zjyw1EJKpvNI)}QptUG%uutoS|B4^|7he!AfcY7@!A@xUUZcfj%ER)Owxl?0ZXv%WD zk4Rpn%FQ_UeYZ)%B%2#oj?+gbaJJd?Uk8tOXd@(z`FDA4J!sz7vs z?XraS9B1pfko6V2RRNX5V!}qxSJNU=v>$VMjJ)-K8S$i-BKiAurOL}A%CVFE*?Dpw z8>IBM00^|N5@QD=S;ty_?Uo zQD5Uln#a`{Iqf44YnJOZ<`I0azlDXW;$u8ca(L*~7*Lt=_C~0wGHpzIw|tX*Unjjr z>%HS&a!hVpVCvLE(rqt8K5}GWqG45{A-VktR`gLc0V9aP5RPJTf;a-^C>oQk&eG1{e-|q^TW3pYw?umt z3QKZg;&T-WPSQkNMdbO-DL47Bur_ppiV%NS=3ykuz&?Z_sui*Cfo2GKFqba_{&35^ zuHFpl7EzH5>TnfdWQS4UvFBy%#GAUh( z;rn--+BI(`5jFLImV=ZDQKNnA)^RNkjAjf}7jQenvkWBZ)!e4H8%3{K59VeH>&2O) zysyxmh`W!^bgkvuMEu|tx*8dxcJ7^!Zj72+u%@uW zef6pb?V?YyHG^ISw=vu=B}hLP7Z;UL${kslD-(COeFqmUwQgZ_(}d1n?w-uy8N9=7 zTdgUS7=y1O*?i-^ldEH8H)Usg!Jyu-NQ5>(KbNWiZynraKDj;>o);1C8_5x=A#^aT zsBNRo)11HT_Jvc0VgX^M6IaYB%fe}eTf&puM9fjh_#_oMHO&K}f^vH)*Vwt#zIX?D z@6&ZZTsNbaM%RxvF0c1Zs2AA|A0C*QtBZ;vZnqt2CXEYw1YAmT6lF#RO@tga4? zz+a6KyTex}dE1P<>pYtkI)_p8nx}u_tsm7p$1$-fv+`6-E#h<5LuA!W`zfD0$Sr*@ zOc~1xH0BQ;I%ddM?~(W8u?xe(kA=FrPddpCkhWmSwp$81(}ya&*}J{r!GVt z&?vj~If<8uP{@{XiS@2g!_4F&zH8)P_mC>67$ivyiLJ&OdIX#hXtRZ&IKUwXH_&T?xg z)0|6pzbUDFXHWk(#)(?ZnGC9XF8g^W1FVJ?qOUmC_uFVC@HEVEe5~pAzuuyeekNNk z!dkXR@~FELNoKT~oUN46cu%|z8amU=FLBL<{fe(9Em6^;W3axN+Dw zCmGG|_W^!xKBYei2-KOQsyJGH}a#*zfx>w+JwD1--e|nm2p#+imCI;j=!R z`5A*}Umm7oxHufiNR^r*w?QUVueMh&NTBmhHOg3&_KtbHN&J?%X1MD0InlZ*6)VRg zg_3cDh)bm2=N*`(V+n<7uMX{D3?AKHXkffagQw16pS~C~y3?fo#Za+KvDH&ODWKSFmV zR2x!na1{G&uC?hkjXZnxGgzSaowlkqo9*=K*qNjEN$IIi6sDtRM=BMjB8PdM(BbkP z^F;wq*hQ~KC<$$phP|+%jodc5;PvUWtda4OjFIv5DT3&eefzSW@SQl~y&hbGd~8u= zkkFfwbpa8WSDcww)H%wi`IJ$n-r}}#LR{D72OT-f0f$f7^|?3N2D)*wdSs8qysM&K zv|K>EaZe$gJVwuQE|IK^E80toL8IYokK~PlRo2KfzOP50+>Ncj=}f;~oJS-JY0(XQ z)^wR#RXodYugPcrb;SaO6`4`rqsPXCC>}8{w(OIfW0N~AL&tvQOWJvAt?O<OPx$pFG&xN}WVkCbDuxj+s+VZM|AAo{B1( zM|D6{{Tj=3xhLy|#^~BuF*CJs_iJ=tzB`d;J#RBs8j*;6%E%G^oNN2jbP5J#YUyrd z!{?mIcyYvz&E?f%fIue8u^WfIdT*8XhGaea-WQ^YTc6Np^LU(XQW%IUU&-}6$%i04 zSzp_K`OMcw6048uOr4HUzF3;d4Ye|UI!?!;;>U6boiVl_qtX!KWjJ1} zY;wtloCk9bH$04H{Ny0$vw%=MpMATk`jVYv)4SlByClP~KD{-fdawdLrpYn~A?3nt_yl<^ zeW3Y$N5=)Xs2<-$UdiJeWyg%K3-|EvkhJ1{O1VC|w^X%!F}Oy{@|DzCw@4N_9mYud zwD6d&YmW={Nqze2`B9GXt2m4Ps{0RC`ZXA3rRoB6v(#*;bp>U=)U$-ib8K0dEMAu$ zPgz)^yj;lZYsbYymyZrk_`dX$&^E^Q~m8r5ScNF1d8RDd4*Nf#IIQg2MYRd2{xJ{&D>3 z6;FYnBzLBZt0eM$k%FTE3;Pdnm)s-zs`tL%;C?zN99PPXYdLh`5>r~B{7Ls)9TPf| zx`T7e%tW<1^~~Zk30Imj9gZqKmi$08y;Gq$HE)RC=+YG9_V72%*PXybH>s^{(cPQ8 zl%P(-%@CO^mQWkwTN~T9Z{u9M!$!v=lU>ifKuWj}o-F|w&dWOj}t!iQYf+}v;5 z$PYG`uO9!@+KzAgNoze>pHuSrcyiJ9?B4S+qD9YXbEu!KuQkakHu?DCMo1Sw)pkkA zjRlt1qxRhT(*0=f3Dq+Oj*qnhtgcjCZ+Wid6fomfu*o{i?Yxpny`r{f&Dlca+j}}I zCrft^H!CwoI;*Oo(-d~6=~2Y9NTLi1DT6{Ko=T@c{2>TTr)WEb!V)a3%*-)T=GGP{ zDJ0UGAcX)wQfMNSu2w=*d0C?rrsBoT!~p(GJVF$6+PoK9NH z*;2{O!|Je-3<`nT{nKA7s>HG>B*;kw0$OnQOoX75Al>e7Bp!t$;{KE&2m~w+36)_` zXw087B#wY15dV}R2t))HP{rM!zjc7~2*_XOk$x8vu|&v$U^5ayfNmzDFnE{@gM=#H zy(oMg9*_EMFVK$wN{Cp97Z5}QlJNV&-E~9+8VU6b5rH8L~T z-j0Yw!4)GhaK%V04*T0pei?|w<6*|2(eU$7gulE+`Zp#ZZ!oaUC=?zVci?>xv|V+- z@5N$Zw?pA!p(COZe}-h&|3owjemg7%4O4=~;$b(&qF{+hMB@>#FaamPLV_mZNwC{t zkVu$f@FIfU4ugT;4ui$Oat4DV!gUZ(e_e$1Pe~9lL>x>B7KwsoH4cY`9SEuumcUp7 z+%q^N_D=&9f%|b7_-61@gct*A`#%y6MIF-&`=XFXgpL##3Nvl0qBQ?IvfFD1gaP~1CmZh z5J)r;Lx9`__y7ru76FqOmqWoy!Xu&B3znlW5M2O~iLk6fU@%08ok%1O0fkQ_8ij^B zcQ-7*>UX#F{?0xm5)0A~5*g5a;4p^AprHm5@mM&D0U&~}BcOlFl3&mZtN||iQ!xVM z6nrxQ_50>s^a1P8&@=%uEW9oWIG9(71QH(qn>BwONc=mW{`L%l2-?@5H3JlbWw3PG zg)!(wh!`w<9hUgZGvNNeDZxW?3@9P|$7Ue?Qwb6z2pXCT2m}h+mB2d8AD7$pKM-P} z2@hmAA_QlFi~z-BAcG+m5&_T?0-r#JLc`8OqDj!$frVHk6m<|t92UMA4@*{1CD5uz zB2Wluh=B|RYkfe7#=)#cVG)pcf^|4(DFPXu09l71U_ebk!-OKB|EREE`H#V(Ac{dm zaA?SZpgaMdLJmZMivQ!l-Eg7^@ccj{(2xWMhv9G}NPU9SkT62G%XnZ3hekuG45(Mo z4?J~Dyz4<2p96VG_%kFV;4TcT zRqV?CSD*iuCvX=ECd2)QTM-z*ouNi!&?uM(F=zsw1X+j#NCS;K7Egee4hDny^Cr6& z!NNm=0f_|B0Yf zC?P_zdRGRcut-3EA)9w$;m>l~&3-J1)1OLqWdxY}0nFe9|JQ&1hG4J|g@C#hq#%?G z0@jUy;{5JjD0c%g(0pOW{KmzRf4#R5fI7omh{ZzP4-QORD+jazi-Q&h@GcIHH-Icd zj{SwD5NGT{9|rC;0+bR0o(A#++-KKq|2Pc`7#Jinz)tWj0djF9*b8+o7PKQ6WB}(v z*CB8KPyapdSJ;pU91_~nfeh5?zYBLe`mPMdf(h6^c*tMt09gIGch@s8YJ$edqW{Bc z*y9r@#z1iiiAWSIqhKT9u<-7M!^1Edfx|-s4Nm=E5PpvV=+?0603U5&SOOlPzpCx` z0RDn)NR}c%Lcp^L54Qj59{>g@B-Cf%%>l0jFj4qT$6o;i7luT7S22umVi0&}5P;&{ zZ(6%?`!^F2`2T1?csNT&;7Q=U1KSLI1Cb%YgaWE$7dIgF1B@ZSJ0Bj56d>mz2!Ma> z%E0~qsUHMjx=_Ukz;xl8(U9&1QX4$cpaudYfW9ffx<5(FuK$4$4^I#R9?pWm3kz-x z0qz9?5zZZe6CjrZwxCHkNYWs|U4^8d^HLGI+ufVbch(4t)EB#uyCcVEIAB{n<5v|G{P)oVyVT(2xLkVgKpZf6wQ@ z69=CCL?Vnq0h6K38FXe0@IS~e2Oa^?VH#(*&;%>6{zLlh{6y_CVo=%4oXH}Bh}bFj|8 zD}=Zoh~N(mBp4(^i~*2>guRXcKZKH15&~F*hm_8*UIOz17)U~~Vs}yjOH45O_%p`4 zTonh;51^I^;eEhr*t;EUCPEW#*TBCs_;2_Is09KK4GEG6dok@kApW@s>3_TrD8~H9 z!vAnT=(5nn1d8zxFMvJ?$P?5U6qxrzWF!zeXwraiqM-l+WEeOupfFHe02jf*m@X1j z9mJSj15v=}U;Y1I7vcd}Lo5WFq4*7yz((7^Kpd>~BEi%dq8O~h|3(5yzdmySYlzUz zphbgsG0YeO^yv+ZAwroCPy#+`{JD4cJh+VXU%-GO{>L|x|A1f=5yl;n0Oui!LDc~; zgc*bWH6bwu82FF-3BWfn_oI=}cRRQS7&rWRAW(vW4k`ejNj$UrNI_zu7__?_RzTqE z(|_2GhVdVOBJkdbgpIBMB?4c#fSJ2F_V*okAqd860BL}y9e4u39SKSedT-FXA=rmR zgD3BAdw=CPP>d(SauEY&G|&UVrvyl}0QUrm3|j|W2$R8IibxCz<{EGy9tAfBHoF4r zh;U2;qcCXn!Pg?_v>VhT>VGs9-2WfwoYS0Sm(=&^%yDz-R*6qrfX4jO(CE2>&+GzsVR54IV9!bvVFrVQd|T$H9~! zu|#-&;9xWalrngULX?0Kg%&AzqT~Mx9Y~$OEd*aL;IskE%%Dk#2Tcl*=+iDF_Z>T%V(Bzuk6fKS{|Qi=Cr+ZvEb;> zp*&<2cgVpnwy7?cw56H*OVmkxXPp}xyT}-2UaZ%n#=U=17E3!HO{}IqM2;;VpyrS# zX(nE%bpPsGlt@zAHhP~L7%)gaE^?M}#A=9Z$m)>BA-fDz)Y#kjY}W4#JLPxN2HAX$ z8*Vf$H!|~^=1N0E4A&o2&>$c2W}_)U=p*i zBI$Q8)(t;eeqgA43a3A}byUGI-B_z%PBuV};0XUZgRt<1k7a+uqs`ik*DFSP}SLR+nYq z?Y*Z-9rM98&%O^_#Q$8zc8bNDs4t#Nni(fuZ8+)hXm{}hpMa6Q@U(S85ne4KK zv`#>fg`?;JURCqkbeB}WIPs)v$hU=^v^o2}qn)#8#c=mv+p6Tk;KW)js^LT`IYA|F zjAS8_$YS^HP|py}=hf~{+9uY%T-CRvPpGZAnJwED`&qb~>c1xU{iLXTJ&;uID@v^Mx?e@njeN&BRo2?VVgO<32 z%JWD?ygDAH7lkA5Zz3FDXx?JnVOhQJqEbaN-nh2eOZSOb&8=Iu7O@`7Uaw+=jym`K z#TCLGGfYug;P<*Z0~4|3SFv8lG_HMC4R%ThO)XnZ{u0AQ`AX@W`VW1s7iy7X=fAef zoH^Tkf)L2pPR5^;z)+NM)Za|tYD{xUR9JD>tTM5@NPYG6&l~c>gn-PtmLJ^`fg?-! zGbFmOe)9D9EL{7gS|8{G;qeN99 z?j>vVG47xBx*QZw+us=;=ol{;4r1`AbS*q^Yy@|ivFK9o*$dsp_tw@X7OgJG+?qsJ zY+*?LH7UzK9OD%l1OnY-Lj?u)`(w%~EUJe|)}%-G8*}A7{l%Eq&ChVljoNX&tx3`3 zBtGxhsk1X3GU`?|Qo8VDF;wmio3sdv4`)J-J$=&SD2gLIS8^OIXAfVP8Bj56mKiYK zh&&<|L>kJ@dz6*9k5XEQ(!Jmvb$#q|QL z%O;sN7rlyBljl}PKi#HK3f{=vO#5iSyX=w^oY@yUvVvwmox2r~5r2PxsTLCybRQRe z@!Q1aSIMAE;m}|KW2^Foyc^m-+;#nYsQi8=_Xb8E!SOH|r`o%bUEt?A>BqFsBquN! z-N9%b@FFshDuhW0Ctk6Y+V*K0hG_m&dyF<^yqF*KLA6Y_usQuij zJrIOnr%ctY=%k@G+f4o*SkUIHfZ?hRxR4iQbM4tf;oHZu$y2w-1is}FA4Z2IUB29# zOd2=lV%aBQK9nBEv@U8B^W|{$J<%UCjv3E~xc#5BmVUZoTNe}AC?|ESOt?U&BX{dB zQ%*Y7qa20vVq&%WjvCbU+jc3{((XZKyjR+S=Gz~BM@M{^ZW!b9Q>PAcCT@+b?RJq| z!kfJ)+vdo7r6s5ed1Jy$1pL3=xyl+oSGq+aTWpZr6#7BE{}m5|GwOvNja_gdhkNHU zq{Si6&D%M%SMaZYEKAGp2#k-Rl{24G$bFzM{T|%T*Szn(bSC$X&P6UzK!ek0+YB1Y5M4D88lHIdV|s(ZfLl z8Wq0woY>5K`ten{^X@EkZjoxk4TMdZQ);Ichw@X4#AB#>ZaXCUf8yLE^GUtK{9%!g zXF#~fkz8W@MbtQ2BesodwcyF~7q125xI?r^EjOOL(Fncb7yL=?eP_E$;Nk$A1$M&a zR&v~zZ2mgxrDb;EMk|bKQbg7Bhg@{g!jaY%${o8${b$m!R+-ZuJE%=|Z)l(V?qfeG zL`7y9t7^X;B*D~hK=`D}u~y}nzU2BtJknE2E?tb#%O&Je-^=R5kyA&UT3)iy#UvNM za3j4_6<B;jw z+pp=O=iu$IK`-`d&3@YZ%SVICQVz;{m-~zT&0l=+Y9=Q;$yu84acAW0!Ej0?0)omh zQSo-ov?7R+cQRV01_2j7EHzA9 z`nHh7d~cs7%q&_&lwEwv%-3yX#kj_}(J{9zG3!_S;pPL=SQVl`h@#zr31dgyndcvC zK8vw%Dcw4t>^xJAalS=vP=g7cwaR99H@WSSo|6~cRy8#k%Ba^D^m$D+fY$G{{KT@= z-574v*Mwr{Jlb>OUL#DQjYPK$mMV%dx?}WbNqS`VX&-NRmodG4X{vAE9K_@s;5qYC z#px$Q)>EA48cE`bSXk$QzD+UOZacc;oo|@^16q=jqs^_ z?mGDf_x&K=?$Vt+(NlLMGLxeslD{6HBV%jacpk8?Rr{4m_ChURY|wH2k(RTnuTHZ{ znpl1yr6uC7Q5@q`;X6o8nDO{)gk0-+7}8#HxnP*9)L`jZ+ToiAAC~aiH%n5rD(6a5 zJk6U?qm4UIuC#C4!BfEs^Ihboy*2aFxcoSoTghs^r-GVkwg>mx?UOHdQHONXbW^$o zwZ#uzyXtu8(eC|)Yy*X%m(#4JmnGC&)H8a&21%26u^Mg`In-&jUwPswYk~*#byAL+ z(XvtMyYHuR5ID$}bjvC#pR=BC*q4|TdO6SL8yA5h%Uqt)llj0Suhgnbd2APEnnZ=a zV|vbJQ=WZmM!s=JKJSaOBqmqXsj@`=)Z_)@fse;uOU-x|M{iTSBHs?tkMeSJx>UOu zb5l98PYS_?qU+)RHojgw{B3UwEn^o)#-->tS% zBy&oCqFRx$`L64+RplghAD+&u6D#@XX7kIxDD-~Vr{DzxIn3{5L_Odiys=q!$ZuEr z&WlDdpC+~$3#GBbqgypEpOvF}S7q*hC%fI(t>x)hl4fi8o%V`~lJcF4ljl#(zO24P zyk0zfVUhI3E=GepmSduiw!Q!IUwt30-ZI)3+Qz+{rQnH?98JHmKB)104%xzQ5h`#% z6eAbFYt(l&TvwjrWM2g9eVX1XwmVJw+=#W?Ibp)}b-_hj0+%?*xW%c$+ssQP-e#-j ziPIm}&3d@XktU{~esA?+Xj`Tyvb0Td`12dufO$4$&<|Sa&K2H!zXhO zJUhiaMe{`|ZrE~};=Di0oC<1Hch$d??_J$)6E*uFJmlQmgCgH4Goj||SDR=WQugbM zkwwOLt1WEJ+$&fwOXmMZ{v(tgHz@=@2R!<&ZhE}1es=837$Im#SYn)v2p z%MI#x==-iD271#(-#ag4N*SAae%eQ|&y+dUT~9zAF`gw3Id_(S%;^~*i`RLuC82`6 zBx}=O81x{>pD*`~hWI++b8DdrTd#!wi5SjTOM+kS*(i)&?Kv4v)j@^mG&yg5e1{!v zMSEe69rfFv|N?giv#(Cts+KNx9D~VDjnxf|7tR0owNFTlWP+t@= zJq!_cefpf(WD6>e3wPWeEMOj@`)1|3 z%ElU~fjwNogI~_FLXtPG6Eq)hM0)gkwRd_FkrewoR~EBqt`C!C4^M9UxtKq2DQJu9 zY?P)e_;6_bvN+1CLk=g6d%g^4Qi57U6+Hqs?c;DH^s_U+&y zgQ@lI?!olC`QC1aNk}(x{+~sIM^hr0)G9=dRAupZzYy_D?~AC8PCaD^<~)bU}5eJlQkoul&>coIVZwq+5#5TFF1N-AO-jD0h^e%*+@m z8{)suBd$S0QLOza1%**Wg3NwCO484Qz-{^ooX~`_poB@v!)<0-k>r>K(n9`#ukB&J!M8w*gd`zX4>alnB_Wa>r|{xo8v9noEJP!@9t&%1GMJMOg8 zYV5<@K7ynr*{c=@^+~eZEj3~dSS{w)j%Wtf{lwdb1p~Um`!6Xqj?MTzjuFy}e?^Nk zYV{INk9b*mJFeYi)46s&8h`Hznzht^E}g~G{KcjyX%j6dkvh7jt)wSYA&)p2Q zQE@|chGpn{I6Gv?cQ5>cv$yodo0x18-v@1}*n&~711^m`AM%eDEIyv&4PSL#sGShn z?39|&C6>(xhS#um#pEzp>`flo7_&ySDokkLblk*Uz2qE zM|^)?)R+8zGv?;8Qu7jZwHYLdh4|r;4)My2hzp}a7Ou!Y1rv-=Xo6M)!u>ZbirC+4kIH(_}m$eVyZcgkbeqlGq_#%CE1`&aW6Hx%l!u>&gAh z8io_O3Y0j{cpkEGTCcrIp;@}ss(Li?q+R~(+s5?HTgjc~4;aqt_@D1d*qyPpZBu@K z@GnV=2~W2W$`(`^ji-Avk7;?vc88Loiyt4JdHbck`d<+>4mPJZKKVFGwK({NK8SIZ zraNJA?c)du`T13*q5jK#l$Haf4sV3}?&U9>vNB-^$ud%;pQc*n33K>X_pR`Xtp4PS z3jApo1Jw7GwULbJK#cx|O-8DkP1@^mJHoUktzgazy&Tj5TC8+K@La=zLGxJ-rIBZE z?-et+dCtG`ivA?Vc&yaf(nfsx+I%_}Ma_BMi$g`>10-atNT~0P=m4>F!(wdOfuau) zkMdpC1|rpK)g&s}`O*f=*+rM;Q(Vp^m^wUp{cup-Y5&RntVh&!sH)3UuoCNrf%E?F zE3!SZ3bz!eC8VVU`isVem*Q`F+nf*-799~gc&0?VLebzUMYuBe_&hmRHkPZiBIcRG zdSbuGTwNlm^1Eq;w82BwWJ6xfRSj#ohnyDUkEGa-dN~;^FnZj*67o%g$FSkRv(%Zc z*hhKoM-OIJb@Sv0h?{X9{P6e*Kbi8fWeqKb(a-yZuU|NaT$h(>=mz9vVRIPjVja3XNPnp$Ew+QS~?u?Tp{;KhDb2L8MN9Tk4D>3_W zb?XnEIVYZE-f(wzEmxt_0a0z~IxmWRMvMkHMIWB%HV_@*Y0G|+mK69PM?~$NDs!>s zvs)rGZ6=k7od4!2nvg7Oy|Yq0ky`^> zpK!huXRs+Z0!g|y@=Kil$ucG57*h)WPglu{jxN-d9 zSmP|cXDx4@UL7_uAHdX+eQhoxct$pTR_YBLuzP$>N@j53bD=xk7qytv^E#1Xj6wW9 zhoz4leOWPcm)Gc6f_o1ax(7ENHfg+&lsk`$sT-GU={rSutr;pQGV^j(Ox&GYEeu3tc&~7! ztCLNQ!Uq=~&Fb~ogd5L22gvKHZ1_*-wS479u++bNukj)Z&uDyG~?B3C=o-KM5sG=j2FK7m_2!SA?r`pDsD*b_Y79tZ1t(&D;~ zk;>cfl9F;k_>{b5bk*5c@VuIA?7zHKt7g5gM*qZox8D>)7<2OFa1=mxGmsNSvZ2+d}g&+rvbZi>87peujEt z%EfEQ!}RFSibTG}S&u7%I8h4WYcZUz-f{!dcDJ-qce;a3CY;Hq7`P%|CfxM+#vUi+ zJawaYqbTVe=@=@zxs&S-<@58khNssi&(mC7v749_b4;kVnfWlo7;)09)mO98PRzt? zbTl*@T^M*T81)!wc%rpanrZ}ve(-!YFDZK7i27Tl&H^j*t^z0f^`Y|cjPUriqONV{ z;0)@s@6OePte0;RFWy?)cz#ld3t`5x=v=#!QyAl9J47<|-jwdT_4)4TQQ_OZ!E>LO zopZ)Be-x$Vi#?QLl7HsxrXtZtGCF|u_Ut=ng|3xfcuOzqazeKJS^-6L#le%*XUKSP zv0748gB}d~cjUh|S-STjM2~(}26Tv|Tg-<4h-#eh9@;m| zaIW9PP1w^0j~`C2yPD|d{6JBZ`A$^xs?*j1#^hxId*V=RP$gZ|%WtMs&yKD2=E~q} zzNJ{DOGP_cWxON5#~eT_rcjak$M3r$j`L!=!Nlwy_*f^GGKJ$CmPuM(< z>GIEOa%qdy|4iCb&R>F^N-DC5K)dnwAj8658{z&b( z#zWgOsKZ-LYimZ}=5%D&AW#~~9}}Bsw|eh0opa&Ip*g=A0>($Us<~~nCRjsc>}m-g zZZ^8IRj&T9KJ|u&$e=Y zR5VcsuP~hcpw?fE*TGj2SN%B$z3Kd4r!ZnI@DmSv4r{V=b`~XS53&>WUijSoRJU_{ zu-VO7FMFGHQHWGwu=;Q`a)Ltt?c$}cZ4H~(C3rp67_G`|L_AGOY$Hx&qb@|Skt+{G z2$ijw>jj3n_O(jRv^bAqtqmVjwl$YiuhcMw3oqsqO}_ehESYpN9vQ<{P9F-mt%GpC zm8bN)n!TmfWvTi!HTAejqO*&i)0r%rws$TY_V)yzhtWMIMTX9r-F|H$kEpvi7u$t- z%&Fj-R+y-SPcn;pmPv2odh~gymF%fcTTMn`t4sFF=T#(94qNnmK2;yy7FAMtm)c_C zOGA6hdE0B9oBYBfFHU%C_Na^61wSHiX!s0&9>`UjT=eyvX%jv%48)1YHc6Sd6OMRw|!G4ksiR+oDI&&Fu3VTg~M^$v&m5F>Zz%n zI-A1ZNp4SxArMo7y2&|H8k)(eD|0%@2Yi10{+b#^8JnA&o2#wdCELpmZWrV(%Xv8X zq^PM=Tp-~T6+JvJ+giEM`=mcj_SNGs=5MQUNL`e;k8}t*T_f&ftq_qXg5=T8X4${b zAZcUYA(%vu&OJu+J=WLN zk+*0#=SvT!=C)w&0_i5fF-Ryb!v3f0RDVHu$!|HJy;yQ&i zuiq8>5P?no$))YZ+0(2uz1JQlSW9os5P}bD+S;AY)4Kam^W>>st%u^B@!b}smVxn} zB56-v5 zd{?3c9Zr!7aB4gWJgm$|G!5C0Ny+MU&a3A~g=*@QFRB z5-oa9td;#W3%e>E?fE8+0B!ZeF-2Czt8Y2lr?!h8k($0GO^DQ=7je@H%w_m^Zi|_{ z%_oDNBbr?Kh)?J%92@JdL{QrXKI4&Q-+-r6q#GVYUQn!tJ zR&o>D*dY@CdhChRMWNT@?BhGpJkz~g6~(pFsJSe(PKOftktsvHJ1t7>D(iFG>#uek zd9;^3TDYT@Z?oqydENK25SZ8%qNpEhN@HlO5=iwMyf~9^t>XA~E8VnLrQsuwHHLWl zdH1oeNb&{GS#zH{F66^ZU;iU*ZPmZlb}r$Gqh&!B>2bn-7UqKsO(OdpSf3=Zgw|O) z@ttgwTe4Mn@+K$5^JN%C#36@>MU0_DCd&;SQPHiGOBRA()lO*7caTYnRq#RU?8XED)e_COw^ z>2{4o2RdA(@7cEk-q4yx?3*H&HrrwW?TF)8y)3! znb6m8{aYzI`N5Y#wh2E9s{1Y5W?k0%zfj*{aaPIya`^O1l;B}3N4NEli6Vs`tx74z zHM+^h#KhPRIm(ca;UN)fuWMHi<>H+CUY_ywU}WGz3}A7U-oAc8D?jquc#U{2q`e9_ z<)^%57jcc0l>NQr%7Pwqb8rcQ5I^_sERC-xS<6^%k6dtQzqnF>oqk|>pXpNd>etBY zoB2cf;q@JE7X^=Td{p~^#nv>M7}&k#{a#dhP>m6Bn&EY(;FIN3ZLe&eQ$?O)jpv+X z5gXP>nN*?EWsob!dRc~D&UMAldTidvVF}+PirktoBTcB}R57aukMQfRbf7~*RPdrP z5>0a2XF^54o2}SeE+2fAc&gusL*AONB}HP2+n_WFGj2~zAY6$RsKhI>mQg%1|d&SrNnaT#b3Ma$fuu*cLo|#^? zdHSHqs<&DnC$=cv<4%9EoxC^VLy37AVaxuE3+p4fm-@8%rR|C6*T=hHbe3X@C${-@d%RBPmUzJqKd znyP^sS(Slpw|0&9&$AwV$spdi%Nw%);QY3WGoyH0{KA!2E7ULFEoi+Ri*Q%nIj=uU z(oT0yS-$^;^x5oC*M&A21HUM?>!g?WwJ;UC&Q4BQP3@oFvt;zgk6bRzZ}cz?oHtDP zcD&OX>n;5A<&lWmZ%ns?#JYb>zxQ>BGNk5@X)zw6)m!HriZW$9x{UNF-)KyEelO=O zbqSIV*UUV)&%OF`*Sp!BeaOG`(6MFNmDHLE=jmOfZRJR|m75uS#lTRht*Ff(QGaNDkm9qTPTG2yWD^TrYEYVChiUV0 zvU!IB)0HXGti$W-ce-5@jYR6DhK$l9LZdgMUN$ix1kZRfC^{%ic)ov|HDkSbug~=S z5f-K>ORAz!(d(zH2gMKLZ3zrshYx?bC(Fv#kEQowEGuIF`r%UAe3zyiMRv~J!v@b~ zI&B#4I=>CS^hqnSLxJ|hiy5;^F$>1>hpb2CsPCD~%WsjsE3#46&M9iuoTOv5XT2ke z{BXbxu`q1GAtS~gA}lI9*U}-=b$>}?`&2V!FXP-O3+ZqwXW#jew^Z8;GECuP)kcO1 zbnVgU8}mcowy#xt$V+c-$Y^I9uMwCpZrjp}9Ex{2pv)UzP%$jq_*i}zT|mO%+ppAh z{rt-;rk+Hi&;EHZc=Ldq*zU#399;@;T?Qo107;Zb;p8#+6wwR{#2*C942l5qPvUN?b0bK7&H6@=_*+BZVXqha{qK7`zk$ zDULvh|3xo%((Qt(m8Y$Qsyqq-8a&jm6P5e2DJ02BM9?Mfy@?QX6{O$$ABjid(a1k( z2m%3z0Bz>p&#!X`0tSW0{z*gP2zcmw2m%p7guh3?{{0^5_XmkMG~_{W8HpkO$pp|) z1OkkP!66_tB9e%MpCcgAznug85#T{05n=@d5g6t4+N7BhH|CMgbi6yhXE90Q<1 zKtl(K!(%WI43RSAetZm6PULM7%UFTn21L~nSi?}0)#OTkcj!y zgTNga5IQNG4HPUS91;!rEzl7}G(=4z9t}-lBod2-qAwDKC8GaS6NC#$E~r#QV6{E} z?`8A90gVUVhA}4KA%TJ*5@3M>Xh`^70_uO<#UOv<@i&@-%UF17g7p9G;9d}ea|9$@ z7$O!%0|g9XRRocULcs!xNJ7IMOT-Z22XSbaFu;az8UcD23FJmY1OzlB;qS}@`u{RN zAflm#3D7W5?Enpy)_{hG`W?^+P(%h9{wHSb*&h&*u;2n7#zO)LocqIfdmaHaEEW>1 zfQBQ&?-I}uhzIA0BxvFz5n!bSGZzX4O>%G!6oWs570|F4NS*>34hr#rMnL^@?zd`6 zzyi5p4-)VYsrNhrbrY~D4vI`5J&`cifo|^4k^$@wq#?os1PvSt5e9_-*bSlKK!Jr> z3yp&IC`bev4Q*xs4MQNoObALgL}PFe|4;Atd$6Jj1jsoA=n!Fm1Y-Xe;Qt%wz;pi* zn)w%T&xBY^iqInpb`s1B7+8HkA}~m3MgrTSVP3$2{skf@xQj+Y>f#>lmm&XKm;ez2 z#S9RfpbLR`0gERHIur7O#Mu27*0*i%%BNhh-M=Ty1 z{J_$9B#s1=5D#lGfC-^kfk1(Pga#S{H0^%|3-Aa5j8_mIpq&0Ui%Gx6LIosfdnN7(mnyJ;*ih|XpaW9APEpbM+0Gu0Er3F4wwtpYk<3GXe|Mj!9yEW zfMq~x7!T0i3;xq${K6UVAZXEH2eEK};$X8axC={5KqJ69I6(W?C8YmM=?EMVHm86F z9bT~TAp8Fma8El}cZCLuz@i_Igh4+Efk#2p1VO+-3ejF2!ow!BJsLD~5P&gkSjHd- zP?vzae`na=_W#%Xpc#c>03P-rxC}ml{b>b2!@~+Cn1Eqm1`6+Qflm5oNdij9|6ojn zYJx-%z&|Ksu!4XI4VqgB=7PHzct(&27<`UZ;1lRH4<5wBiVdI<;lhC2feHN& z2JG4YU-J{dBmq4L#!kqCKo;Op*f}KB#|Ur=T82R-1)l%+uF*ezA`s>`;(ur;cbd*J(h*4uK_sMFiBWpo+p8b0i)Ma}y}FFi(TKutpys z?5}R|my7pOg|HX0f5i44SV7|$^jlEIpp?Oug21d`riEGI*Us%<7lUB~&Ffm#O4`hP12^#A1(ARug)252}~eFJ+TxT%4B5RC!M!De&t9`TO``fvLK384`T zA^`*4006^)DH*B>3f9*E##mSp1MO@AwLc*x!KR?>q*Y;GnUvcMdjM00D7Okiz4M@Ki+M{sinj`xEf66*gdk zgGM@_;bC1M5=DS@fq*fz4(=uA|8N}v3*BJ?VTe$K1B{6%h~M`fgx9?Pg$pR+-#`QO z{1xaRHvvrj-e2fJE2Q4~S?0X9`tvKrKw*P#XpmrHgy8}}DVUFmU_%FS z3Fr)qtLLIG&Q&$30Bu+^MjQm7Jh5z@2 f1-+soy@#jOWlx{Wws!QuO<-Y8FDk00t4{xa9ygvu diff --git a/doc/user-manual/fpakc_UserManual.tex b/doc/user-manual/fpakc_UserManual.tex index e5c0d4f..cbd1573 100644 --- a/doc/user-manual/fpakc_UserManual.tex +++ b/doc/user-manual/fpakc_UserManual.tex @@ -177,9 +177,10 @@ \section{Variable Weighting Scheme\label{sec:weightingScheme}} One of the issues in particle simulations, specially for axial-symmetrical cases, is that due to the disparate volume of cells, specially close to the axis, the statistics in some cells is usually poor. To try to fix that, the possibility to include a Variable Weighting Scheme in the simulations is available in \Gls{fpakc}. - These schemes detect when a particle change cells and modify its weight accordingly. - To avoid particles having a larger weight than the rest, particle can be split in multiple particles if weight become too large. + These schemes detect when a particle change cells and split it if necessary to improve statistics. The use of a Variable Weighting Scheme is defined by the user in the input file. + + Beware that this can increase the number of particles in the simulation and increase computational time. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% \section{Interaction between species}\label{ssec:collisions} diff --git a/src/modules/solver/moduleSolver.f90 b/src/modules/solver/moduleSolver.f90 index 58a7a71..88f7bf8 100644 --- a/src/modules/solver/moduleSolver.f90 +++ b/src/modules/solver/moduleSolver.f90 @@ -388,11 +388,12 @@ MODULE moduleSolver REAL(8):: fractionVolume, pSplit !If particle changes volume to smaller cell - IF (volOld%volume > volNew%volume) THEN + IF (volOld%volume > volNew%volume .AND. & + part%weight >= part%species%weight*1.0D-1) THEN fractionVolume = volOld%volume/volNew%volume !Calculate probability of splitting particle - pSplit = 1.D0 - DEXP(-fractionVolume) + pSplit = 1.D0 - DEXP(-fractionVolume*1.0D-1) IF (random() < pSplit) THEN !Split particle in two