Fixed difussion approximation in all formal solvers by removing the term with the source function derivative. This term can be negative and lead to negative intensities at the boundary
Some cleaning up and remove debug code