Finally, I found this discussion on UDN, which mentions a fix introduced in UE5 (CL 17679918).
The solution is to modify UnrealMathSSE.h and UnrealMathNeon.h, specifically the VectorSinConstant values:
static const float a = 7.58946609f; // 16 * sqrtf(p)
static const float b = 1.63384342f; // (1 - p) / sqrtf(p)
After applying this change, my Cook process is now working fine.
Hope this helps!