 smoothness indicators in are calculated as: $$\label{eq:weno_betal_k2}
    \begin{aligned}
        \beta_0 &= (U_{i} - U_{i+1})^2 \\
        \beta_1 &= (U_{i-1} - U_{i})^2.
    \end{aligned}$$ For higher-order reconstructions, the derivatives in can be computed by using a general expression of the reconstruction polynomial. For $k=3$, after rearranging terms, we find (see [@Jiang1996; @Shu1998]): $$\label{eq:weno_betal_k3}
    \begin{aligned}
        \beta_0 &= \frac{1}{4}(3U_{i} - 4U_{i+1} + U_{i+2})^2 + \frac{13}{12}(U_{i} - 2U_{i+1} + U_{i+2})^2 \\
        \beta_1 &= \frac{1}{4}(U_{i-1} - U_{i+1})^2 + \frac{13}{12}(U_{i-1} - 2U_{i} + U_{i+1})^2 \\
        \beta_2 &= \frac{1}{4}(U_{i-2} - 4U_{i-1} + 3U_{i})^2 + \frac{13}{12}(U_{i-2} - 2U_{i-1} + U_{i})^2,
    \end{aligned}$$ see also the work of [@Henrick2005]. For orders $k=4,5,6$, the corresponding expressions were obtained in [@Balsara2000] and extended to even higher orders $k=7,8,9$ in [@Gerolymos2009]. From the expressions in and, we see that $\beta_\ell = O({\Delta x}^2)$ as soon as the solution is smooth over the stencil, while we have $\beta_\ell = O(1)$ for stencils containing a discontinuity. Using equation, this results in weights $\omega_\ell \approx d_\ell$ for smooth solutions and $\omega_\ell = O({\Delta x}^4)$ for discontinuous solutions, as desired.

Evaluating the Boltzmann collision operator with a fast spectral scheme {#sub:FastSpectral}
================================================================
smoothness indicators in are calculated as: $ $ \label{eq: weno_betal_k2 } 
     \begin{aligned } 
         \beta_0 & = (U_{i } - U_{i+1})^2 \\ 
         \beta_1 & = (U_{i-1 } - U_{i})^2. 
     \end{aligned}$$ For higher - order reconstruction, the derivative instrument in can be computed by using a general formulation of the reconstruction polynomial. For $ k=3 $, after rearranging terms, we discover (see [ @Jiang1996; @Shu1998 ] ): $ $ \label{eq: weno_betal_k3 } 
     \begin{aligned } 
         \beta_0 & = \frac{1}{4}(3U_{i } - 4U_{i+1 } + U_{i+2})^2 + \frac{13}{12}(U_{i } - 2U_{i+1 } + U_{i+2})^2 \\ 
         \beta_1 & = \frac{1}{4}(U_{i-1 } - U_{i+1})^2 + \frac{13}{12}(U_{i-1 } - 2U_{i } + U_{i+1})^2 \\ 
         \beta_2 & = \frac{1}{4}(U_{i-2 } - 4U_{i-1 } + 3U_{i})^2 + \frac{13}{12}(U_{i-2 } - 2U_{i-1 } + U_{i})^2, 
     \end{aligned}$$ see besides the workplace of [ @Henrick2005 ]. For orders $ k=4,5,6 $, the represent expressions were obtained in [ @Balsara2000 ] and gallop to even higher orders $ k=7,8,9 $ in [ @Gerolymos2009 ]. From the expression in and, we see that $ \beta_\ell = O({\Delta x}^2)$ as soon as the solution is politic over the stencil, while we have $ \beta_\ell = O(1)$ for stencils containing a discontinuity. use equation, this results in weights $ \omega_\ell \approx d_\ell$ for smooth solution and $ \omega_\ell = O({\Delta x}^4)$ for discontinuous solutions, as desired. 

 Evaluating the Boltzmann collision operator with a fast apparitional scheme { # sub: FastSpectral } 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = =
 smlothness indicators in ave calculated as: $$\label{ex:weno_bstal_k2}
    \begin{aligned}
        \beta_0 &= (U_{m} - U_{u+1})^2 \\
        \beta_1 &= (U_{i-1} - U_{i})^2.
    \evd{aligned}$$ For higyer-ocder reconstructmkns, the derivzbives mn can be compujed by using a general expsersnon of the reconstruction polynomial. For $k=3$, sfher rearrangind tegmf, we find (see [@Jiang1996; @Shu1998]): $$\label{eq:weno_betzl_k3}
    \uegin{aligned}
        \beta_0 &= \frac{1}{4}(3U_{i} - 4U_{i+1} + U_{i+2})^2 + \fraf{13}{12}(U_{i} - 2U_{i+1} + U_{i+2})^2 \\
        \beta_1 &= \frax{1}{4}(U_{i-1} - Y_{i+1})^2 + \frac{13}{12}(U_{i-1} - 2J_{i} + U_{i+1})^2 \\
        \beta_2 &= \ftac{1}{4}(U_{i-2} - 4U_{i-1} + 3U_{i})^2 + \frac{13}{12}(U_{i-2} - 2U_{i-1} + U_{i})^2,
    \end{akigned}$$ see aoso jhe work of [@Ienrich2005]. For orders $k=4,5,6$, the cmrrespomding expressipns weee obtained in [@Balsare2000] and extended to evgn higher mrbers $k=7,8,9$ in [@Gerolymos2009]. Feon the expsesskins in aid, se see thet $\beta_\ell = O({\Delta x}^2)$ aw soon as the solutoog is smooth ovsr the seencil, while we have $\beta_\ell = O(1)$ for stetcims containing a discontunuity. Using equation, this restlts in weights $\omega_\ell \approx d_\ell$ for smooth smlutikvs ckq $\imfga_\ell = O({\Delta x}^4)$ for discontinuous solutions, wa cexired.

Evaluatinn the Boltzmann cokllsojn operator wkth a yzsf spectral scheme {#dub:FastFpecteal}
================================================================
smoothness indicators in are calculated as: $$\label{eq:weno_betal_k2} &= - U_{i+1})^2 \beta_1 &= (U_{i-1} reconstructions, derivatives in can computed by using general expression of the reconstruction polynomial. $k=3$, after rearranging terms, we find (see [@Jiang1996; @Shu1998]): $$\label{eq:weno_betal_k3} \begin{aligned} \beta_0 &= - 4U_{i+1} + U_{i+2})^2 + \frac{13}{12}(U_{i} - 2U_{i+1} + U_{i+2})^2 \\ \beta_1 &= - + - + U_{i+1})^2 \\ \beta_2 &= \frac{1}{4}(U_{i-2} - 4U_{i-1} + 3U_{i})^2 + \frac{13}{12}(U_{i-2} - 2U_{i-1} + U_{i})^2, see also the work of [@Henrick2005]. For orders the corresponding expressions were in [@Balsara2000] and extended to higher $k=7,8,9$ in From expressions and, we see $\beta_\ell = O({\Delta x}^2)$ as soon as the solution is smooth over the stencil, while we have = O(1)$ containing a Using this in weights $\omega_\ell for smooth solutions and $\omega_\ell = discontinuous solutions, as desired. Evaluating the Boltzmann collision with a spectral scheme {#sub:FastSpectral} ================================================================
 smoothness indicators in are Calculated As: $$\labEl{eQ:weNo_BetaL_k2}
    \beGin{aligned}
        \beta_0 &= (u_{I} - U_{i+1})^2 \\
        \bEta_1 &= (U_{i-1} - U_{i})^2.
    \end{aligned}$$ For hiGher-oRdER recONsTructIons, the DErIVAtiVeS iN caN bE CoMputeD by Using a gEneral exprEssIoN of the reconsTRuCtion polynOmiAl. For $k=3$, after rEarRanginG tErmS, We finD (seE [@JianG1996; @Shu1998]): $$\laBEl{eq:weNo_betal_k3}
    \bEgIN{alignED}
        \beta_0 &= \frAC{1}{4}(3u_{i} - 4u_{i+1} + U_{i+2})^2 + \Frac{13}{12}(U_{i} - 2U_{i+1} + U_{i+2})^2 \\
        \beta_1 &= \fraC{1}{4}(u_{i-1} - u_{I+1})^2 + \frac{13}{12}(U_{i-1} - 2U_{i} + U_{i+1})^2 \\
        \betA_2 &= \frac{1}{4}(U_{I-2} - 4U_{I-1} + 3u_{i})^2 + \FRAc{13}{12}(U_{I-2} - 2U_{i-1} + u_{i})^2,
    \end{alignEd}$$ See alSO the worK Of [@hENRicK2005]. for orders $k=4,5,6$, the CorrespondiNG exPressiOnS weRE obtaiNed in [@baLSarA2000] and extendeD to eVen higher Orders $K=7,8,9$ In [@GerolYMos2009]. From The expResSioNs in ANd, We See ThAT $\beTA_\eLl = O({\dEltA x}^2)$ as soon As ThE soluTion IS SMOoth OveR the StencIl, while we have $\BetA_\ell = o(1)$ For StencIls coNtaiNiNg a diScontiNuity. usIng equation, this ResuLts in weigHts $\OmEga_\ElL \apprOX d_\ell$ fOr sMooTh solutIons and $\OMegA_\eLL = o({\deLta x}^4)$ for discontinuoUs SOLuTions, as dEsired.

eVaLuATing the BOlTzmAnn cOLLisioN opeRAtOr with a fAst speCTrAl Scheme {#sUb:fastSpEcTraL}
================================================================
 smoothness indicators inare calcul atedas:  $$ \l abel {eq: weno_betal_k2}    \begin{aligned}
       \be ta _ 0 &= (U _{i}- U_{i+ 1 }) ^ 2  \\ 
    \ be ta_1&=(U_{i-1 } - U_{i}) ^2. 
   \end{alig n ed }$$ For hi ghe r-order reco nst ructio ns , t h e der iva tives  in ca n  be co mputed by  u s ing ag enerale x pr essi on of the reconst r uc t ion polynomial . For$k = 3$ , aft errearrangin gterms ,  we fin d  ( s e e  [@ J iang1996; @Sh u1998]): $$ \ lab el{eq: we no_ b etal_k 3}
  \ beg in{aligned} 
       \bet a_0 &= \frac{1 } {4}(3U_ {i} -4U_ {i+ 1} + U_ {i +2} )^ 2  +\ fr ac{ 1 3}{ 12}(U_{i }-2U_{i +1}+ U _ {i+2 })^ 2 \\ 
       \beta_1 & = \ frac { 1}{ 4}(U_ {i-1}  - U _{ i+1}) ^2 + \ frac{ 13 }{12}(U_{i-1} -  2U_ {i} + U_{ i+1 }) ^2\\ 
      \be ta_ 2 & = \frac {1}{4}( U _{i -2 } -  4 U_{i-1} + 3U_{i})^ 2+ \f rac{13}{ 12}(U_ { i- 2} - 2U_{i- 1}  +U_{i } ) ^2,
   \ e nd {aligned }$$ se e  a ls o the w or k of [ @H enr ick 2005] .  For  order s $k=4,5 ,6$,t he correspondi n g expressions we r e  o b tain edin [@Balsar a200 0 ] an d ex t en ded to ev en hi gh e ro rders $k=7,8,9$ in[@ Geroly mos20 09]. From the  expressio n s in and,we s e et hat $\beta_\el l = O ({\Delta x } ^2)$ assoonas the s olution i s smooth o ver  th e s ten c i l,  while we hav e $\be ta _\ell =  O( 1)$ for  st enc ils  co nt aining adisconti nu it y.  U sin g equ a tion, th is  re su lts  in w e ights$\ome ga_\ el l\ app rox d_\ e ll $ forsm oo th s olu ti ons a nd $ \ ome ga_\ell  = O({\De lta x}^4 )$  f or disc ontinuous sol ut ions, as d es ire d.

Ev a l uating t he Boltzmann collisiono perator  wi th afast  spectral  sc heme { #su b :FastS pectra l}
== == === = = ===== = = == === == ========== = = === ===== == ==== ======= ========
 smoothness_indicators in_are calculated as: $$\label{eq:weno_betal_k2}
_ __\begin{aligned}
 __  _  \beta_0_&= (U_{i} - U_{i+1})^2_\\
  __    \beta_1 &= (U_{i-1} - U_{i})^2.
    \end{aligned}$$ For_higher-order_reconstructions, the_derivatives_in_can be computed by using_a general expression of the_reconstruction polynomial._For $k=3$, after rearranging terms, we find (see_[@Jiang1996;_@Shu1998]): $$\label{eq:weno_betal_k3}
 _  \begin{aligned}
      _ \beta_0 &= \frac{1}{4}(3U_{i} - 4U_{i+1}_+ U_{i+2})^2 +_\frac{13}{12}(U_{i}_-_2U_{i+1} + U_{i+2})^2 \\
_     _ \beta_1 &= \frac{1}{4}(U_{i-1} - U_{i+1})^2_+ \frac{13}{12}(U_{i-1} - 2U_{i} + U_{i+1})^2 \\
_     _ \beta_2 &= \frac{1}{4}(U_{i-2} -_4U_{i-1} +_3U_{i})^2 + \frac{13}{12}(U_{i-2} - 2U_{i-1}_+ U_{i})^2,
 _ _\end{aligned}$$ see also_the work of [@Henrick2005]. For orders_$k=4,5,6$, the corresponding_expressions were obtained in [@Balsara2000] and_extended_to even higher_orders_$k=7,8,9$_in [@Gerolymos2009]._From the expressions_in_and, we_see_that $\beta_\ell = O({\Delta x}^2)$ as_soon_as the solution is smooth over the_stencil, while we have_$\beta_\ell_= O(1)$ for stencils_containing a discontinuity. Using equation,_this results in weights $\omega_\ell \approx_d_\ell$ for_smooth solutions_and $\omega_\ell = O({\Delta x}^4)$ for discontinuous solutions, as desired.

Evaluating the_Boltzmann collision operator with a fast_spectral scheme {#sub:FastSpectral}
================================================================