######################################################################## # # Optimization Control Parameter # # model : select model (1:2-series tank, 2:storage function) # algorithm : select algorithm (1:SCE-UA, 2:PSO) # objfunc : objective function (1:RMSE, 2:Nash-siutcliffe efficiency) # &ctrl_opti model = 1, algorithm = 1, objfunc = 1, / ######################################################################## # # Optimizer Parameter File (SCE-UA) # # nc : number of complex # nm : member in a complex (0:set default 2*n+1) # ne : evolving member in a complex (0:set default n+1) # alpha : number of alpha iteration (0:set default 1) # beta : number of beta iteration (0:set default 2*n+1) # it : maximum iteration count (0:unlimitted) # ie : maximum evaluation count (0:unlimitted) # ir : interval of message printing (0:no message) # flag_quiet : restrain messages # flag_initp : use initial parameters # flag_widemutate : wide range mutation (recommended) # flag_qmcinit : initialize by Quasi-Monte Carlo # &opti_sceua nc = 5, nm = 0, ne = 0, alpha = 0, beta = 0, it = 0, ie = 100000, ir = 10, flag_quiet = 0, flag_initp = 1, flag_widemutate = 1, flag_qmcinit = 0, / ######################################################################## # # Optimizer Parameter File (PSO) # # np : population size (number of particle) # w : inertia parameter (0.9 is recommended) # c1 : self intention parameter (2.0 is recommended) # c2 : swarm intention parameter (2.0 is recommended) # c3 : random search parameter (very small value is recommended) # vmax : limit of velocity length (0.0-0.1 or smaller is recommended) # vrand : velocity perturbation parameter (0.0-0.2 is recommended) # it : maximum iteration count (0:unlimitted) # ie : maximum evaluation count (0:unlimitted) # ir : iteration of message printing (0:no message) # flag_quiet : restrain messages # flag_initp : use initial parameters # flag_lrand : local random search algorithm (use c3) # flag_vlimit : velocity limit algorithm (use vmax, strongly recommended) # flag_vrand : velocity perturbation algorithm (use vrand) # flag_gcenter : evaluate center of gravity (recommended) # flag_gjump : move to center of gravity (unrecommended) # flag_qmcinit : initialize by Quasi-Monte Carlo # &opti_pso np = 20, w = 0.9, c1 = 2.0, c2 = 2.0, c3 = 1.E-7, vmax = 0.5, vrand = 0.15, it = 0, ie = 100000, ir = 100, flag_quiet = 0, flag_initp = 1, flag_vlimit = 1, flag_lrand = 1, flag_vrand = 1, flag_gcenter = 1, flag_gjump = 0, flag_qmcinit = 0, /