 (\epsilon_{\alpha'})^2 \sin(\vartheta_{\beta'}) 
\sin(\varphi_{\gamma})(\Delta \epsilon)_{\alpha'}  (\Delta\vartheta)_{\beta'}.\end{aligned}$$ The bin center values $\epsilon_{\alpha'}$ and $\vartheta_{\beta'}$ are given just as $r_{i'}$ and $\theta_{j'}$ in (\[rcenter\]) and (\[thetacenter\]). Also, the momentum bin widths are given by the difference of the bounding edge values, just as in the case of the spatial zones. The definitions of $(1/r)_{i'}$ and $(\cot\theta)_{j'}$ will be given below. The values of ${\cal N}$ on momentum bin edges are given by “upwind” or “donor-cell” differencing (see [@liebendoerfer02]).

Finally, we show the remaining definitions needed to ensure that the second and fourth terms in (\[example\]), as well as the third and fifth terms, cancel for ${\cal N}$ spatially homogeneous. Given the differencings defined above, this is accomplished with the following definitions: $$\begin{aligned}
(\sin\vartheta)_{\beta'} &=& \sin(\vartheta_{\beta'}), \\
(\cos\vartheta)_{\beta'} &=&  {\sin(\vartheta_{\beta+1})^2 - \sin(\vartheta_{\beta})^2
\over 2 \sin(\vartheta_{\beta'}) (\Delta\vartheta)_{\beta'}}, \\
(\cos\varphi)_{\gamma'} &=&  {\sin(\varphi_{\gamma+1}) - \sin(\varphi_{\gamma})
\over (\Delta\varphi)_{\gamma'}}, \\
\left(1\over r\right)_{i'} &=& {(r_{i+1})^2 - (r_i)^2 \over 2 (r_{i'})^2 (\Delta r)_{i'}}, \\
(\cot\theta)_{j'} &=& {\sin(\theta_{j+1}) - \sin(\theta_{j})
\over \sin(\theta_{j'}) (\Delta\theta)_{j'}}.\end{aligned}$$

Radiative transfer algorithm and distributed-memory implementation {#sec:5}
==================================================================

In this section classes of operators involved in radiation transport are identified, and a strategy for implementing them on distributed-memory computer architectures is described.

The equations of neutrino radiative transfer are the transport equation for whatever radiation particle distribution function is used, together
(\epsilon_{\alpha'})^2 \sin(\vartheta_{\beta' }) 
 \sin(\varphi_{\gamma})(\Delta \epsilon)_{\alpha' }   (\Delta\vartheta)_{\beta'}.\end{aligned}$$ The bin center values $ \epsilon_{\alpha'}$ and $ \vartheta_{\beta'}$ are given just as $ r_{i'}$ and $ \theta_{j'}$ in (\[rcenter\ ]) and (\[thetacenter\ ]). besides, the momentum bank identification number width are given by the difference of the bounding edge value, just as in the case of the spatial zone. The definition of $ (1 / r)_{i'}$ and $ (\cot\theta)_{j'}$ will be given downstairs. The values of $ { \cal N}$ on momentum bin edge are given by “ upwind ” or “ donor - cell ” differencing (visualize [ @liebendoerfer02 ]). 

 Finally, we show the remaining definitions needed to ensure that the second and fourth terminus in (\[example\ ]), as well as the third and fifth terms, cancel for $ { \cal N}$ spatially homogeneous. give the differencings defined above, this is accomplished with the following definitions: $ $ \begin{aligned } 
 (\sin\vartheta)_{\beta' } & = & \sin(\vartheta_{\beta' }), \\ 
 (\cos\vartheta)_{\beta' } & = &   { \sin(\vartheta_{\beta+1})^2 - \sin(\vartheta_{\beta})^2 
 \over 2 \sin(\vartheta_{\beta' }) (\Delta\vartheta)_{\beta' } }, \\ 
 (\cos\varphi)_{\gamma' } & = &   { \sin(\varphi_{\gamma+1 }) - \sin(\varphi_{\gamma }) 
 \over (\Delta\varphi)_{\gamma' } }, \\ 
 \left(1\over r\right)_{i' } & = & { (r_{i+1})^2 - (r_i)^2 \over 2 (r_{i'})^2 (\Delta r)_{i' } }, \\ 
 (\cot\theta)_{j' } & = & { \sin(\theta_{j+1 }) - \sin(\theta_{j }) 
 \over \sin(\theta_{j' }) (\Delta\theta)_{j'}}.\end{aligned}$$ 

 Radiative transportation algorithm and distributed - memory execution { # sec:5 } 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

 In this section classes of operators involved in radiotherapy transport are identified, and a strategy for implementing them on distributed - memory calculator architecture is described. 

 The equations of neutrino radiative transfer are the transport equation for whatever radiation particle distribution function is used, in concert
 (\epdilon_{\alpha'})^2 \sin(\vartheta_{\beua'}) 
\sin(\varphi_{\gamma})(\Bwlta \e'silon)_{\ampha'}  (\Deuta\vartheta)_{\beta'}.\end{aligned}$$ Thx bib cenuvr values $\epsilon_{\alphx'}$ and $\varnheta_{\beta'}$ are tiven just as $r_{i'}$ akb $\thefw_{j'}$ iu (\[ccenter\]) and (\[thejacenter\]). Alsm, the momentum bkn widths are given by the difference of the blunding edge vwluex, jusf as in the case of the spatial zohes. The definitions pf $(1/r)_{i'}$ and $(\cot\theta)_{j'}$ will bf gigen below. The valufs of ${\cal N}$ on iimentum bin ddges are given by “upwjnd” or “donor-cell” differencing (sde [@lizbendoerfer02]).

Dibalpi, we show thx remapning definitlpns neaded to ensure that tme serond and fourth terms in (\[xxample\]), as well as tre third dnb fifth terms, cancel doe ${\cal N}$ s[atixoly hojojensous. Glvei the diffedencings dedined above, this is asbpmplished wifh the fjllowing definitions: $$\begin{aligned}
(\sin\varuheta)_{\geta'} &=& \sin(\vartheta_{\beta'}), \\
(\cis\vartheta)_{\beta'} &=&  {\sin(\vwrtheta_{\beea+1})^2 - \sin(\vartheta_{\beta})^2
\over 2 \sin(\vartheta_{\beta'}) (\Delta\vardheta)_{\uega'}}, \\
(\eis\varoyi)_{\hamma'} &=&  {\sin(\varphi_{\gamma+1}) - \sin(\varphi_{\gamma})
\over (\Dejfa\fagphi)_{\gamma'}}, \\
\left(1\ovev r\right)_{i'} &=& {(r_{i+1})^2 - (r_i)^2 \pvfr 2 (r_{i'})^2 (\Delta r)_{i'}}, \\
(\zot\thecz)_{j'} &=& {\sin(\theta_{j+1}) - \sin(\thfta_{j})
\ovet \sin(\tyeta_{j'}) (\Delua\theya)_{j'}}.\end{aligned}$$

Radiative tranwfer algoritkm qnd distributed-memlry implemeutatiom {#sec:5}
==================================================================

On this section classes of kperators ijvolved ih radiation transoorn ara identified, and a strategr for impoemeuting thdm om distwibuted-memlry computer architecturfs is dascribed.

Thf equations of neutrino radiatitx transfer arg tve nransport equabion for whatevqr radiation pcrticle bistricution funbtion is nsed, togethew
(\epsilon_{\alpha'})^2 \sin(\vartheta_{\beta'}) \sin(\varphi_{\gamma})(\Delta \epsilon)_{\alpha'} (\Delta\vartheta)_{\beta'}.\end{aligned}$$ The bin $\epsilon_{\alpha'}$ $\vartheta_{\beta'}$ are just as $r_{i'}$ (\[thetacenter\]). the momentum bin are given by difference of the bounding edge values, as in the case of the spatial zones. The definitions of $(1/r)_{i'}$ and will be given below. The values of ${\cal N}$ on momentum bin edges given “upwind” “donor-cell” (see [@liebendoerfer02]). Finally, we show the remaining definitions needed to ensure that the second and fourth in (\[example\]), as well as the third and terms, cancel for ${\cal spatially homogeneous. Given the differencings above, is accomplished the definitions: (\sin\vartheta)_{\beta'} &=& \sin(\vartheta_{\beta'}), (\cos\vartheta)_{\beta'} &=& {\sin(\vartheta_{\beta+1})^2 - \sin(\vartheta_{\beta})^2 \over 2 \sin(\vartheta_{\beta'}) (\Delta\vartheta)_{\beta'}}, \\ (\cos\varphi)_{\gamma'} &=& {\sin(\varphi_{\gamma+1}) - \sin(\varphi_{\gamma}) \over (\Delta\varphi)_{\gamma'}}, \left(1\over r\right)_{i'} - (r_i)^2 2 (\Delta \\ (\cot\theta)_{j'} &=& \sin(\theta_{j}) \over \sin(\theta_{j'}) (\Delta\theta)_{j'}}.\end{aligned}$$ Radiative transfer implementation {#sec:5} ================================================================== In this section classes of involved in transport are identified, and a strategy implementing them on distributed-memory computer architectures is described. equations of neutrino radiative transfer are the transport equation for whatever radiation particle distribution function together
 (\epsilon_{\alpha'})^2 \sin(\vartheta_{\beTa'}) 
\sin(\varphI_{\gammA})(\DeLta \EpSiloN)_{\alpHa'}  (\Delta\varthetA)_{\Beta'}.\End{aligned}$$ The bin center ValueS $\ePSiloN_{\AlPha'}$ anD $\vartheTA_{\bETA'}$ arE gIvEn jUsT As $R_{i'}$ and $\TheTa_{j'}$ in (\[rcEnter\]) and (\[thEtaCeNter\]). Also, the mOMeNtum bin widThs Are given by thE diFferenCe Of tHE bounDinG edge Values, JUst as iN the case oF tHE spatiAL zones. THE DeFiniTions of $(1/r)_{i'}$ and $(\cot\thETa)_{J'}$ Will be given belOw. The vAlUEs OF ${\Cal n}$ on Momentum biN eDges aRE given bY “UpWIND” or “DOnor-cell” diffeRencing (see [@lIEbeNdoerfEr02]).

finALly, we sHow thE rEMaiNing definitIons Needed to eNsure tHAt the seCOnd and fOurth tErmS in (\[ExamPLe\]), As WelL aS The THiRd aND fiFth terms, CaNcEl for ${\Cal N}$ SPATIallY hoMogeNeous. given the diffeRenCingS DefIned aBove, tHis iS aCcompLished With tHe Following definiTionS: $$\begin{aliGneD}
(\sIn\vArTheta)_{\BEta'} &=& \sin(\VarTheTa_{\beta'}), \\
(\cOs\varthETa)_{\bEtA'} &=&  {\SIN(\vArtheta_{\beta+1})^2 - \sin(\vartHeTA_{\BeTa})^2
\over 2 \siN(\varthETa_{\BeTA'}) (\Delta\vaRtHetA)_{\betA'}}, \\
(\COs\varPhi)_{\gAMmA'} &=&  {\sin(\varpHi_{\gammA+1}) - \SiN(\vArphi_{\gaMmA})
\over (\DElTa\vArpHi)_{\gamMA'}}, \\
\lefT(1\over r\Right)_{i'} &=& {(r_{i+1})^2 - (R_i)^2 \oveR 2 (R_{i'})^2 (\Delta r)_{i'}}, \\
(\cot\thETa)_{j'} &=& {\sin(\theta_{j+1}) - \sIN(\tHETa_{J})
\Over \Sin(\Theta_{j'}) (\Delta\ThetA)_{J'}}.\end{AligNEd}$$

radIAtive TransFeR AlGOrithm and distributeD-mEmory iMplemEntation {#sec:5}
==================================================================

In This sectioN CLAsses of oPeraTOrS Involved in radiAtion Transport aRE identifIed, anD a strateGy for implEMEnting thEm oN diStrIbuTED-mEmory computer ARChitEcTures is DesCribed.

THe eQuaTioNs oF nEutrino raDiative tRaNsFeR aRe tHe traNSport equAtIon FoR whAteveR RadiatIon paRticLe DiSTriBution fUNcTIOn is UsEd, TogeTheR
 (\epsilon_{\alpha'})^2 \s in(\varthe ta_{\ bet a'} )
\si n(\v arphi_{\gamma} ) (\De lta \epsilon)_{\alpha' }  (\ De l ta\v a rt heta) _{\beta ' }. \ e nd{ al ig ned }$ $  T he bi n c enter v alues $\ep sil on _{\alpha'}$a nd  $\varthet a_{ \beta'}$ are  gi ven ju st  as $r_{i '}$  and$\thet a _{j'}$  in (\[rc en t er\])a nd (\[t h e ta cent er\]). Also, them om e ntum bin width s aregi v en b y t hedifference  o f the boundin g  e d g e  va l ues, just asin the case ofthe sp at ial zones.  Thede f ini tions of $( 1/r) _{i'}$ an d $(\c o t\theta ) _{j'}$will b e g ive n be l ow .The  v a lue s  o f $ { \ca l N}$ on  m om entum  bin e d g es a regive n by“upwind” or “ don or-c e ll”  diff erenc ing(s ee [@ lieben doerf er 02]).

Finally,  weshow therem ai nin gdefin i tionsnee ded  to ens ure tha t  th es e c on d and fourth terms  i n (\ [example \]), a s  w el l  as theth ird  and f ifthterm s ,cancel f or ${\ c al  N }$ spat ia lly ho mo gen eou s. Gi v en t he dif ferencin gs de f ined above, th i s is accompli s he d wi t h th e f ollowing de fini t ions : $$ \ be gin { align ed}
( \s i n\ v artheta)_{\beta'} & =&  \sin( \vart heta_{\beta'} ), \\
(\co s \ v artheta) _{\b e ta ' } &=&  {\sin(\ varth eta_{\beta + 1})^2 -\sin( \varthet a_{\beta} ) ^ 2
\over2 \ sin (\v art h e ta _{\beta'}) (\ D e lta\ va rtheta) _{\ beta'}} , \ \
( \co s\v ar phi)_{\ga mma'} &= & { \s in (\v arphi _ {\gamma+ 1} ) -  \ sin (\var p hi_{\g amma} )
\o ve r( \De lta\var p hi ) _ {\ga mm a' }},\\\l eft(1 \ove r  r\ right)_ {i'} &=&{(r _ {i+1 }) ^2  - (r_i )^2 \over 2 ( r_ {i'})^2 (\ De lta  r)_{i ' } }, \\
(\ cot\theta)_{j'} &=& {\s i n(\thet a_{ j+1})  - \ sin(\thet a_{ j})
\o ver \sin(\ theta_ {j'})  ( \De l t a\the t a )_ {j' }} .\end{alig n e d}$ $

Ra di ativ e trans fer algorithm andd ist ributed-memor y i mple m e nt ati o n{ #se c: 5 }
= = = =============== ========== == = == ========== = === == ======= ======= ===
I n thissection c lasses of  o pera t o rsinvolved i n radiat ion trans p ort a r eident ifi ed, an da s trate gy for imp lemen ting t he m on d istri bu ted-memo ry computer architectur es isdescr ibe d.

The e qua t ion s of neut rino  radiative  tr ans fer a ret he tr ansp o rt  eq u ation  for whateverr ad iat i o nparticle di s t r ibu tionfun c tion i s us ed, together
 (\epsilon_{\alpha'})^2_\sin(\vartheta_{\beta'}) 
\sin(\varphi_{\gamma})(\Delta_\epsilon)_{\alpha'}  (\Delta\vartheta)_{\beta'}.\end{aligned}$$ The_bin center_values_$\epsilon_{\alpha'}$ and_$\vartheta_{\beta'}$_are given just_as $r_{i'}$ and_$\theta_{j'}$ in (\[rcenter\]) and_(\[thetacenter\]). Also, the_momentum_bin widths are given by the difference of the bounding edge values, just as_in_the case_of_the_spatial zones. The definitions of_$(1/r)_{i'}$ and $(\cot\theta)_{j'}$ will be_given below._The values of ${\cal N}$ on momentum bin_edges_are given by_“upwind” or “donor-cell” differencing (see [@liebendoerfer02]).

Finally, we show the_remaining definitions needed to ensure that_the second and_fourth_terms_in (\[example\]), as well_as the third and fifth terms,_cancel for ${\cal N}$ spatially homogeneous._Given the differencings defined above, this is_accomplished with the following definitions: $$\begin{aligned}
(\sin\vartheta)_{\beta'}_&=& \sin(\vartheta_{\beta'}), \\
(\cos\vartheta)_{\beta'} &=& _{\sin(\vartheta_{\beta+1})^2 -_\sin(\vartheta_{\beta})^2
\over 2 \sin(\vartheta_{\beta'}) (\Delta\vartheta)_{\beta'}}, \\
(\cos\varphi)_{\gamma'}_&=&  {\sin(\varphi_{\gamma+1})_- \sin(\varphi_{\gamma})
\over_(\Delta\varphi)_{\gamma'}}, \\
\left(1\over r\right)_{i'}_&=& {(r_{i+1})^2 - (r_i)^2 \over 2_(r_{i'})^2 (\Delta r)_{i'}},_\\
(\cot\theta)_{j'} &=& {\sin(\theta_{j+1}) - \sin(\theta_{j})
\over \sin(\theta_{j'})_(\Delta\theta)_{j'}}.\end{aligned}$$

Radiative_transfer algorithm and_distributed-memory_implementation_{#sec:5}
==================================================================

In this_section classes of_operators_involved in_radiation_transport are identified, and a strategy_for_implementing them on distributed-memory computer architectures is_described.

The equations of neutrino_radiative_transfer are the transport_equation for whatever radiation particle_distribution function is used, together