  When multiple fission modes coexist in a given nucleus, distinct fragment
yield distributions appear. Multimodal fission has been observed in a number of
fissioning nuclei spanning the nuclear chart, and this phenomenon is expected
to affect the nuclear abundances synthesized during the rapid neutron-capture
process (r-process). In this study, we generalize the previously proposed
hybrid model for fission-fragment yield distributions to predict competing
fission modes and estimate the resulting yield distributions. Our framework
allows for a comprehensive large-scale calculation of fission fragment yields
suited for r-process nuclear network studies. Nuclear density functional theory
is employed to obtain the potential energy and collective inertia tensor on a
multidimensional collective space defined by mass multipole moments. Fission
pathways and their relative probabilities are determined using the nudged
elastic band method. Based on this information, mass and charge fission yields
are predicted using the recently developed hybrid model. Fission properties of
fermium isotopes are calculated in the axial quadrupole-octupole collective
space for three energy density functionals (EDFs). Disagreement between the
EDFs appears when multiple fission modes are present. Within our framework, the
UNEDF1$_{\textrm{HFB}}$ EDF agrees best with experimental data. Calculations in
the axial quadrupole-octupole-hexadecapole collective space improve the
agreement with the experiment for SkM$^{*}$. We also discuss the sensitivity of
fission predictions on the choice of EDF for several superheavy nuclei. Fission
fragment yield predictions for nuclei with multiple fission modes are sensitive
to the underlying EDF. For large-scale calculations in which a minimal number
of collective coordinates is considered, UNEDF1$_{\textrm{HFB}}$ provides the
best description of experimental data.
