Abstract When solving the first-order form of linear Boltzmann equation, a common misconception is that matrix-free computational method “sweeping mesh”, used in conjunction with Discrete Ordinates method, too complex or does not scale well enough to be implemented modern high performance computing codes. This has led considerable efforts development matrix-based methods are computationally exp...