#endif
#endif
- /* Mask for the cut-off to avoid overflow in gmx_simd_exp_r */
+ /* Mask for the cut-off to avoid overflow of cr2^2 */
cr2_S0 = gmx_simd_mul_r(lje_c2_S, gmx_simd_blendzero_r(rsq_S0, wco_vdw_S0));
#ifndef HALF_LJ
cr2_S2 = gmx_simd_mul_r(lje_c2_S, gmx_simd_blendzero_r(rsq_S2, wco_vdw_S2));