Parallel Shooting Method for Boundary-value Problems: Application to the Neutron Transport Equation
by J. Canosa, H. R. Penafiel
A direct method is given for the solution of the spherical harmonics approximation to the Boltzmann equation for neutron transport in slab geometry. The roundoff instability of the problem is eliminated by performing linear transformations of the matrices involved, which ensure that the matrix columns are linearly independent. The novelty of the method lies in that only a minimum number of matrix transformations is performed, the precise number being determined dynamically, efficiently, and in a new way by the program itself in the course of the computation.