The diffuse term was erroneously multiplied by (1-Fc)^2 after a previous refactoring. The documentation was correct.