  We use an elliptic system of equations with complex coefficients for a set of
complex-valued tensor fields as a tool to construct infinite-dimensional
families of non-singular stationary black holes, real-valued Lorentzian
solutions of the Einstein-Maxwell-dilaton-scalar
fields-Yang-Mills-Higgs-Chern-Simons-$f(R)$ equations with a negative
cosmological constant. The families include an infinite-dimensional family of
solutions with the usual AdS conformal structure at conformal infinity.
