added phi0_guess_j in the input argument.  This is needed because we need to use the scalar flux in the outer iteration from the past iteration to initialize the phi0_old_j in the inner iteration.  This means that the initial scattering source was constructed with the constructed with the same scalar flux as that used for the fission source
added phi0_guess_j in the input argument.  This is needed because we need to use the scalar flux in the outer iteration from the past iteration to initialize the phi0_old_j in the inner iteration.  This means that the initial scattering source was constructed with the constructed with the same scalar flux as that used for the fission source

jipuwang/FS_LS_MoC1D