This work aimed at generating homogenized group constants using the Serpent code and then using the CriMR diffusion code to model the large SFR OECD 3600 MWth MOX core. The results were compared with a full core reference Monte Carlo solution by Serpent. Reactivity feedback parameters were also considered. Generating the group constants from separate fuel assemblies allows for simultaneously carrying out calculations and then using the results as input in diffusion codes rather than waiting so long for a 3D full core Monte Carlo calculation to be completed. From the results of the integral parameters we see a close agreement in the calculation codes. The differences can be attributed to the errors that could arise from generating the constants from individual sub-assemblies. The differences in the underlying physics and approximations used in development of the codes could also be a factor. Another way the errors could be reduced is by checking to see that the sub-assembly configurations used in the non-multiplying zones are as close as possible to the real layout in a full 3D core.