Internal conversion (IC) between the two lowest singlet excited states, $1^1B{}_u^+$ and $2^1A{}_g^-$, of s-$\textit{trans}$-butadiene and s-$\textit{trans}$-hexatriene is investigated using a series of single- and multi- reference wave function and density functional theory (DFT) methodologies. Three independent types of the equation-of-motion coupled-cluster (EOMCC) theory capable of providing an accurate and balanced description of one- as well as two-electron transitions, abbreviated as δ-CR-EOMCC(2,3), DIP-EOMCC(4h2p){No}, and DEA-EOMCC(4p2h){Nu} or DEA-EOMCC(3p1h,4p2h){Nu}, consistently predict that the $1^1B{}_u^+$/$2^1A{}_g^-$ crossing in both molecules occurs along the bond length alternation coordinate. However, the analogous $1^1B{}_u^+$ and $2^1A{}_g^-$ potentials obtained with some multireference approaches, such as CASSCF and MRCIS(D), as well as with the linear-response formulation of time-dependent DFT (TDDFT), do not cross. Hence, caution needs to be exercised when studying the low-lying singlet excited states of polyenes with conventional multiconfigurational methods and TDDFT. The multistate many-body perturbation theory methods, such as XMCQDPT2, do correctly reproduce the curve crossing. Among the simplest and least expensive computational methodologies, the DFT approaches that incorporate the contributions of doubly excited configurations, abbreviated as MRSF (mixed reference spin-flip) TDDFT and SSR(4,4), accurately reproduce our best EOMCC results. This is highly promising for nonadiabatic molecular dynamics simulations in larger systems."