A new mathematical adjoint for the modified SAAF-SN equations

Numerical results verify the correctness of the implementation and show its utility both for fixed source and eigenvalue problems.