Positive P_N closures

We introduce a modification to the standard spherical harmonic closure used with linear kinetic equations of particle transport. While the standard closure is known to produce negative particle concentrations, the modification corrects this defect by requiring that the ansatz used to close the equations itself be a nonnegative function. We impose this requirement via explicit constraints in a quadratic optimization problem.