Finding optimum balances between conflicting interests in multipurpose reservoirs often represents an important challenge for decision makers. This study assesses the use of different computational tools to obtain optimal reservoir operations applied to the Hatillo dam in the Dominican Republic. A multiobjective optimization approach is used, in which non-dominated sorting genetic algorithm II (NSGAII) and multi-objective evolutionary algorithm based on decomposition (MOEA/D) optimizers are applied to models that simulate reservoir operations. Three different Machine Learning (ML) models, namely, the multilayer perceptron (MLP), the radial basis network (RBN) and the linear function (LF), are employed to learn the current operation of the system. Subsequently, a general model is proposed to simulate daily reservoir operations (2009-2019), integrating water balances, physical constraints of the dam components and the ML models, the latter defining daily controlled discharges. In the optimization process, the ML parameters are the decision variables, while the objectives evaluated are irrigation, hydropower generation and flood control. The results are compared with the actual operation of the reservoir. Three dimensional Pareto fronts are obtained, from which, the wide variety of operations can be evidenced. The flood control objective was found to have a wide room for improvement over the current operation of the reservoir, and several of the solutions found improve the current operation for the three proposed objectives. The MLP models tend to generate the best results for this case study and the NSGAII optimizer generates the best optimization results.