Equation in plant.utils.R (Rs.Comp) inconsistent with Ruffault et al. (2022) Eq. 45

The code in question reads:

A <- max((-1 * (psi + PiFT - Esymp) - sqrt((psi + PiFT - Esymp)^2 + 4 * (psi * Esymp))) / (2 * Esymp), 1 - PiFT / psi)

The term 4 * (psi * Esymp) is inconsistent with Eq. 45 in Ruffault (2022) that is 4 * (PiFT * Esymp). Can anyone comment which one is correct?