Discussions and Closures Closure to Application of the Firefly Algorithm to Optimal Operation of Reservoirs with the Purpose of Irrigation Supply and Hydropower Productionby Irene Garousi-Nejad, Omid Bozorg-Haddad, Hugo A. Loáiciga, and Miguel A. Mari ˜ no DOI: 10.1061/(ASCE)IR.1943-4774.0001064 Irene Garousi-Nejad, S.M.ASCE 1 ; Omid Bozorg-Haddad 2 ; and Hugo A. Loáiciga, Ph.D., P.E., F.ASCE 3 1 Ph.D. Student, Dept. of Civil and Environmental Engineering, Utah State Univ., Logan, UT 84322. E-mail: I.Garousi@aggiemail.usu.edu 2 Professor, Faculty of Agricultural Engineering and Technology, Dept. of Irrigation and Reclamation Engineering, College of Agriculture and Natural Resources, Univ. of Tehran, Karaj, 31587-77871 Tehran, Iran (corresponding author). E-mail: OBHaddad@ut.ac.ir 3 Professor, Dept. of Geography, Univ. of California, Santa Barbara, CA 93016-4060. E-mail: Hugo.Loaiciga@ucsb.edu The writers thank the discusser for their interest in the original pa- per and for giving the writers the opportunity to clarify a few issues raised therein. The firefly algorithm (FA) was proposed by Yang (2008) to solve continuous optimization problems. Continuous algorithms can be used to solve discrete problems; however, it is recommended to use for that task algorithms that are inherently designed for dis- crete problems, such as the ant-colony optimization algorithm. The discussed paper did not apply the storage carryover con- straint in the simulation modeling because the simulation periods of the irrigation supply and hydropower production problems were long enough so that there was no need to use the carryover con- straint. The carryover constraint is appropriately applied when the simulation period is shorter than 10 years (see, e.g., Aboutalebi et al. 2015a). The discussed papers irrigation supply problem (Aydoghmoush Reservoir) and hydropower production problem (Karun-4 Reservoir) had simulation periods equal to 10 years (120 months) and 42 years (504 months), respectively. The discussed paper set the value of the tailwater to be constant in the hydropower generation problem because fluctuations of the tailwater are negligible compared with changes in hydraulic head (Aboutalebi et al. 2015b). The discussed paper evaluated the performance of the FA in solving reservoir operation problems with two different purposes. The standard form of the FA was used for this purpose. The choice of the standard FA is thoroughly addressed in the paper by Garousi-Nejad et al. (2016). The role of penalty functions is to modify the objective function so that the decision variables take values within the feasible range, that is, constraints are not violated. Therefore, the penalty functions adds a positive (or negative) value to the objective function in minimization (or maximization) problems. This was done in the discussed paper. The reason to use evolutionary and metaheuristic algorithms (FA in the discussed paper) in problems that can be solved with classic methods (such as linear programming, nonlinear program- ming, or dynamic programming) is to evaluate the performance of the algorithms and to compare the results with those obtained via the classic methods. This approach has been discussed at length in other papers (see Aboutalebi and Garousi-Nejad 2015). Independent runs of an evolutionary or metaheuristic algorithm are necessary because these algorithms start with a randomly gen- erated population that changes from one run to another. Multiple runs allow the characterization of near-optimal solutions in terms of averages and standard deviations. The discussed paper applied five runs to generate a range of results. Considering a larger number of runs (say 10, 25, or 100) depends on the problem complexity (Aboutalebi et al. 2015a). The discussed paper established that five runs provided a suitable characterization of the variability of inter- run results. The number of functional evaluations has been used in several studies (Aboutalebi et al 2016a, b) because it is more informative than computational time or number of iterations of an optimization algorithm. In fact, the number of functional evaluations is indepen- dent of the computing equipment and thus provides a fairer comparison. The number of functional evaluations is calculated by multiplying the population size by number of algorithmic iterations. References Aboutalebi, M., Bozorg-Haddad, O., and Loáiciga, H. A. (2015a). Appli- cation of the SVR-NSGAII to hydrograph routing in open channels. J. Irrig. Drain. Eng., 10.1061/(ASCE)IR.1943-4774.0000969, 04015061. Aboutalebi, M., Bozorg-Haddad, O., and Loáiciga, H. A. (2015b). Optimal monthly reservoir operation rules for hydropower generation derived with SVR-NSGAII.J. Water Resour. Plann. Manage., 10.1061/(ASCE)WR.1943-5452.0000553, 04015029. Aboutalebi, M., Bozorg-Haddad, O., and Loáiciga, H. A. (2016a). Multi- objective design of water-quality monitoring networks in river-reservoir systems.J. Environ. Eng., 10.1061/(ASCE)EE.1943-7870.0001155, 04016070. Aboutalebi, M., Bozorg-Haddad, O., and Loáiciga, H. A. (2016b). Simulation of methyl tertiary butyl ether concentrations in river- reservoir systems using support vector regression.J. Irrig. Drain. Eng., 10.1061/(ASCE)IR.1943-4774.0001007, 04016015. Aboutalebi, M., and Garousi-Nejad, I. (2015). Discussion of application of the water cycle algorithm to the optimal operation of reservoir sys- tems.J. Irrig. Drain. Eng., 10.1061/(ASCE)IR.1943-4774.0000947, 07015029. Garousi-Nejad, I., Bozorg-Haddad, O., and Loáiciga, H. A. (2016). Modi- fied firefly algorithm for solving multireservoir operation in continuous and discrete domains.J. Water Resour. Plann. Manage., 10.1061 /(ASCE)WR.1943-5452.0000644, 04016029. Yang, X. S. (2008). Firefly algorithm, nature-inspired meta-heuristic algorithms, Wiley, New York, 7990. © ASCE 07017020-1 J. Irrig. Drain. Eng. J. Irrig. Drain Eng., 2018, 144(1): 07017020 Downloaded from ascelibrary.org by 54.161.69.107 on 07/04/20. Copyright ASCE. For personal use only; all rights reserved.