  A general method for classifying the possible quark models of a
multi-Higgs-doublet model, in the presence of Abelian symmetries, is presented.
All the possible sets of textures that can be present in a given sector are
shown, thus turning the determination of the flavor models into a combinatorial
problem. Several symmetry implementations are studied for two and three Higgs
doublet models. Some models implementations are explored in great detail, with
a particular emphasis on models known as Branco-Grimus-Lavoura and
nearest-neighbour-nnteraction. Several considerations on the flavor changing
neutral currents of multi-Higgs models are also given.
