This paper addresses reliable and efficient calculation of the mode of a multivariate sample, which is a classical fusion function. In particular, we focus on the inputs given on the unit simplex, when aggregating elements of Atanassov intuitionistic fuzzy sets, interval‐valued fuzzy sets and their extensions, as well as compositional data. We outline the use of a specially designed 2‐additive fuzzy measures and the Choquet integral for the purposes of reducing computational complexity in higher dimensions. We present computational analysis and benchmark four different methods of density‐based mode estimation.