ADAPTIVE SIMULATED ANNEALING /* $Id: asa.c,v 15.10 1997/06/20 21:36:22 ingber Exp ingber $ */ OPTIONS_FILE = 0 OPTIONS_FILE_DATA = 0 RECUR_OPTIONS_FILE = 0 RECUR_OPTIONS_FILE_DATA = 0 COST_FILE = 1 ASA_LIB = 1 HAVE_ANSI = 1 IO_PROTOTYPES = 0 TIME_CALC = 1 TIME_STD = 0 INT_LONG = 0 INT_ALLOC = 0 SMALL_FLOAT = 1e-18 MIN_DOUBLE = 1e-18 MAX_DOUBLE = 1e+18 EPS_DOUBLE = 1e-18 CHECK_EXPONENT = 0 NO_PARAM_TEMP_TEST = 0 NO_COST_TEMP_TEST = 0 SELF_OPTIMIZE = 0 ASA_TEST = 1 ASA_TEST_POINT = 0 ASA_TEMPLATE = 0 MY_TEMPLATE = 1 ASA_TEMPLATE_LIB = 1 ASA_TEMPLATE_ASA_OUT_PID = 0 ASA_TEMPLATE_MULTIPLE = 0 ASA_TEMPLATE_SELFOPT = 0 ASA_TEMPLATE_SAMPLE = 0 ASA_TEMPLATE_QUEUE = 0 ASA_TEMPLATE_PARALLEL = 0 ASA_TEMPLATE_SAVE = 0 USER_INITIAL_COST_TEMP = 0 RATIO_TEMPERATURE_SCALES = 0 USER_INITIAL_PARAMETERS_TEMPS = 0 DELTA_PARAMETERS = 0 QUENCH_PARAMETERS = 0 QUENCH_COST = 0 QUENCH_PARAMETERS_SCALE = 1 QUENCH_COST_SCALE = 1 OPTIONAL_DATA = 0 OPTIONAL_DATA_INT = 0 USER_COST_SCHEDULE = 0 USER_ACCEPT_ASYMP_EXP = 0 USER_ACCEPTANCE_TEST = 0 USER_GENERATING_FUNCTION = 0 USER_REANNEAL_COST = 0 USER_REANNEAL_PARAMETERS = 0 MAXIMUM_REANNEAL_INDEX = 50000 REANNEAL_SCALE = 10 ASA_SAMPLE = 0 ASA_QUEUE = 0 ASA_RESOLUTION = 0 ASA_PARALLEL = 0 ASA_SAVE = 0 SYSTEM_CALL = 1 FDLIBM_POW = 0 FDLIBM_LOG = 0 FDLIBM_EXP = 0 ASA_PRINT = 1 ASA_OUT = asa.out USER_ASA_OUT = 0 ASA_PRINT_INTERMED = 1 ASA_PRINT_MORE = 0 G_FIELD = 12 G_PRECISION = 7 OPTIONS->Limit_Acceptances = 1000 OPTIONS->Limit_Generated = 99999 OPTIONS->Limit_Invalid_Generated_States = 1000 OPTIONS->Accepted_To_Generated_Ratio = 0.0001 OPTIONS->Cost_Precision = 1e-18 OPTIONS->Maximum_Cost_Repeat = 5 OPTIONS->Number_Cost_Samples = 5 OPTIONS->Temperature_Ratio_Scale = 4e-07 OPTIONS->Cost_Parameter_Scale_Ratio = 1 OPTIONS->Temperature_Anneal_Scale = 100 OPTIONS->Include_Integer_Parameters = 0 OPTIONS->User_Initial_Parameters = 0 OPTIONS->Sequential_Parameters = -1 OPTIONS->Initial_Parameter_Temperature = 1 OPTIONS->Acceptance_Frequency_Modulus = 100 OPTIONS->Generated_Frequency_Modulus = 10000 OPTIONS->Reanneal_Cost = 1 OPTIONS->Reanneal_Parameters = 1 OPTIONS->Delta_X = 0.001 OPTIONS->User_Tangents = 0 OPTIONS->Curvature_0 = 0 start_asa:time: 0h 0m 5.75938s; incr: 0h 0m 0.001277s Initial Random Seed = 377 *number_parameters = 52 index_v parameter_minimum parameter_maximum parameter_value parameter_type 0 -10 60 5.927299 -1 1 0 2.8 0.120949 -1 2 -1.02 3.59 0.398564 -1 3 -0.5 4.5 3.259212 -1 4 -10 90 10 -1 5 0 2.8 0.337503 -1 6 -1.02 3.59 0.59 -1 7 -0.5 4 2.42153 -1 8 -52.9 -0.1 -10 -1 9 0.05 2.4 0.41717 -1 10 0.01 2.96 0.4778 -1 11 -1.5 4.7 1.5416 -1 12 2.6 32 7.01224 -1 13 -0.8 1.5 0.587316 -1 14 0.1 2.5 0.673489 -1 15 -1 1 0.248997 -1 16 0 55 7.166667 -1 17 -0.5 2.5 0.59288 -1 18 -0.51 2.65 0.176994 -1 19 0.1 3.5 2.235 -1 20 -22 -1 -4.382042 -1 21 -1.01 1.6 0.259147 -1 22 0.01 0.99 0.185202 -1 23 0.1 4 2.5 -1 24 0 200 14.59202 -1 25 -1.1 5 -0.755816 -1 26 -1 2.5 0.325556 -1 27 -1.1 5 4.304003 -1 28 0 200 15.0193 -1 29 0.1 5 3.974043 -1 30 -1 3.5 0.472427 -1 31 -1.1 5 0.587854 -1 32 0.6 32 2.2202 -1 33 -0.1 1.5 0.270855 -1 34 0 3.5 0.040061 -1 35 0 4 0.929117 -1 36 0 1.5 0.718176 -1 37 0 1 0.169737 -1 38 0 1 0.347722 -1 39 0 1 0.10192 -1 40 0 1.5 1.000203 -1 41 -1 1 -0.044163 -1 42 -0.5 1 0.6535 -1 43 -0.5 1 -0.018 -1 44 -0.5 1 -0.037159 -1 45 -0.5 0.5 -0.100041 -1 46 -0.5 0.5 0.03133 -1 47 -0.5 0.5 -0.099307 -1 48 -2.05969 -0.05969 -1.05969 -1 49 0 1 1 -1 50 0.2 5 0.75 -1 51 -1 1 0.1 -1 temperature_scale = 13.48324 temperature_scale_parameters[0] = 13.48324 *temperature_scale_cost = 13.48324 :time: 0h 0m 21.7529s; incr: 0h 0m 15.9935s *index_cost_acceptances = 0, *current_cost_temperature = 1.067118e+17 *accepted_to_generated_ratio = 1, *number_invalid... = 0 *number_generated = 0, *number_accepted = 0 best...->cost = 6.807714, last...->cost = 6.807714 index_v best...->parameter current_parameter_temp tangent 0 15.39165 1 0 1 0.531282 1 0 2 3.575981 1 0 3 2.71608 1 0 4 36.23843 1 0 5 0.036893 1 0 6 1.010751 1 0 7 1.514914 1 0 8 -46.9979 1 0 9 1.99214 1 0 10 1.51804 1 0 11 3.393068 1 0 12 3.185141 1 0 13 -0.4789843 1 0 14 0.969089 1 0 15 0.408756 1 0 16 42.76942 1 0 17 0.8063101 1 0 18 2.37605 1 0 19 2.150466 1 0 20 -8.474534 1 0 21 0.5871924 1 0 22 0.3481399 1 0 23 2.004986 1 0 24 140.7697 1 0 25 0.3596286 1 0 26 0.2157005 1 0 27 4.401157 1 0 28 96.34045 1 0 29 4.820492 1 0 30 1.646191 1 0 31 4.244539 1 0 32 3.315556 1 0 33 0.1150464 1 0 34 2.077153 1 0 35 1.32347 1 0 36 0.3357023 1 0 37 0.7748877 1 0 38 0.326149 1 0 39 0.8857665 1 0 40 0.5867356 1 0 41 0.47815 1 0 42 0.7689394 1 0 43 0.8993201 1 0 44 0.9795074 1 0 45 0.06169163 1 0 46 -0.3472331 1 0 47 -0.3201971 1 0 48 -1.885904 1 0 49 0.4467561 1 0 50 3.135538 1 0 51 0.9655843 1 0 best...->cost=5.293438 *number_accepted=2 *number_generated=2 best...->cost=5.187541 *number_accepted=3 *number_generated=3 best...->cost=5.133138 *number_accepted=4 *number_generated=4 best...->cost=4.14622 *number_accepted=5 *number_generated=5 best...->cost=3.656432 *number_accepted=12 *number_generated=14 :time: 0h 7m 31.7255s; incr: 0h 7m 9.97262s *index_cost_acceptances = 0, *current_cost_temperature = 4.280641e+10 *accepted_to_generated_ratio = 0.9673913, *number_invalid... = 0 *number_generated = 105, *number_accepted = 100 best...->cost = 3.656432, last...->cost = 6.405846 index_v best...->parameter current_parameter_temp tangent 0 37.14344 3.945065e-07 2.641054e-06 1 0.6221269 3.945065e-07 -5.420856e-05 2 3.504279 3.945065e-07 -0.0003304126 3 2.201857 3.945065e-07 -0.0006429994 4 87.69806 3.945065e-07 0.006218714 5 0.2634984 3.945065e-07 -0.2512975 6 1.815763 3.945065e-07 -0.9312221 7 1.21851 3.945065e-07 -0.6809777 8 -38.84437 3.945065e-07 -1.65543e-11 9 2.166033 3.945065e-07 -3.758095e-10 10 1.623778 3.945065e-07 -8.39126e-09 11 4.610327 3.945065e-07 -7.428478e-09 12 14.34061 3.945065e-07 0.01660983 13 -0.3651819 3.945065e-07 -0.4978449 14 1.439734 3.945065e-07 -0.3995085 15 0.2750655 3.945065e-07 -1.01591 16 48.40915 3.945065e-07 -4.6641e-06 17 0.1333343 3.945065e-07 0.0001197692 18 2.438705 3.945065e-07 0.0006981029 19 2.150276 3.945065e-07 0.0008896974 20 -16.68175 3.945065e-07 0.09613881 21 0.902882 3.945065e-07 0.6805685 22 0.2426506 3.945065e-07 3.503563 23 2.484928 3.945065e-07 0.1022041 24 145.5784 3.945065e-07 0.003303394 25 0.2157853 3.945065e-07 -0.2472618 26 1.124556 3.945065e-07 -1.376171 27 2.117616 3.945065e-07 -0.7371101 28 109.7975 3.945065e-07 0.003929876 29 4.702343 3.945065e-07 -0.1981155 30 0.7368416 3.945065e-07 -0.6688134 31 0.9957017 3.945065e-07 -0.2064132 32 11.78367 3.945065e-07 0.000281792 33 0.2743805 3.945065e-07 -0.009292187 34 2.078797 3.945065e-07 -0.009039811 35 0.4496174 3.945065e-07 -0.03349971 36 0.8393064 3.945065e-07 0.1249977 37 0.4892337 3.945065e-07 -0.03843324 38 0.3935065 3.945065e-07 -0.02900256 39 0.7902602 3.945065e-07 0.0003061408 40 1.014078 3.945065e-07 -0.2346984 41 0.5879797 3.945065e-07 0.1174917 42 0.232972 3.945065e-07 0.07664256 43 0.882173 3.945065e-07 3.912724e-06 44 -0.3374997 3.945065e-07 -1.725312 45 0.1003566 3.945065e-07 0.2524135 46 0.1754344 3.945065e-07 -0.2725218 47 0.2442378 3.945065e-07 0.1520637 48 -1.475948 3.945065e-07 -0.009384554 49 0.5055021 3.945065e-07 -0.2570595 50 3.868493 3.945065e-07 -3.577941e-08 51 -0.4785027 3.945065e-07 -2.498858e-08 best...->cost=2.338397 *number_accepted=107 *number_generated=118 best...->cost=2.119854 *number_accepted=108 *number_generated=169 best...->cost=2.079151 *number_accepted=109 *number_generated=232 best...->cost=2.004001 *number_accepted=110 *number_generated=305 best...->cost=1.954364 *number_accepted=111 *number_generated=352 best...->cost=1.904111 *number_accepted=112 *number_generated=425 best...->cost=1.870824 *number_accepted=113 *number_generated=455 best...->cost=1.855202 *number_accepted=114 *number_generated=495 best...->cost=1.6444 *number_accepted=115 *number_generated=542 best...->cost=1.580338 *number_accepted=116 *number_generated=645 best...->cost=1.389681 *number_accepted=117 *number_generated=704 best...->cost=1.389423 *number_accepted=118 *number_generated=722 best...->cost=1.367948 *number_accepted=119 *number_generated=770 best...->cost=0.9306805 *number_accepted=120 *number_generated=785 best...->cost=0.9261777 *number_accepted=121 *number_generated=1137 best...->cost=0.8874596 *number_accepted=122 *number_generated=1204 best...->cost=0.8737539 *number_accepted=123 *number_generated=1607 best...->cost=0.8650884 *number_accepted=124 *number_generated=1727 best...->cost=0.8576754 *number_accepted=125 *number_generated=1749 best...->cost=0.8398228 *number_accepted=126 *number_generated=1994 best...->cost=0.8366002 *number_accepted=127 *number_generated=2149 best...->cost=0.8225202 *number_accepted=128 *number_generated=2329 best...->cost=0.6877029 *number_accepted=129 *number_generated=2344 best...->cost=0.6761485 *number_accepted=130 *number_generated=2561 best...->cost=0.65366 *number_accepted=131 *number_generated=5639 :time: 7h 39m 54.3795s; incr: 0h 5m 6.94263s *index_cost_acceptances = 0, *current_cost_temperature = 3.557337e-06 *accepted_to_generated_ratio = 0.0002292526, *number_invalid... = 0 *number_generated = 10000, *number_accepted = 131 best...->cost = 0.65366, last...->cost = 0.65366 index_v best...->parameter current_parameter_temp tangent 0 41.70471 1.025741e-07 0.003604065 1 1.535058 1.025741e-07 -0.1605917 2 2.139116 1.025741e-07 -0.1996208 3 0.2666296 1.025741e-07 -0.4564135 4 32.25038 1.025741e-07 -0.001709917 5 1.184166 1.025741e-07 0.01814815 6 1.681267 1.025741e-07 0.08597013 7 1.510732 1.025741e-07 0.03661235 8 -37.87002 1.025709e-07 0.03731762 9 2.071735 1.025709e-07 1.283447 10 1.300241 1.025709e-07 1.394613 11 -0.02812291 1.025709e-07 1.620657 12 13.57843 1.025741e-07 -0.007827798 13 0.1646008 1.025741e-07 0.05288933 14 0.7701157 1.025741e-07 0.138246 15 0.625305 1.025709e-07 0.01750213 16 35.48538 1.025741e-07 0.00458072 17 0.9773088 1.025741e-07 -0.04153031 18 0.7778884 1.025741e-07 -0.2726003 19 1.237569 1.025741e-07 -0.1424105 20 -3.601587 1.025741e-07 0.01520563 21 1.465942 1.025741e-07 0.0397662 22 0.5116924 1.022427e-07 0.1504907 23 1.381081 1.025741e-07 0.05645354 24 83.21117 1.025741e-07 0.0003991088 25 3.07243 1.025741e-07 -0.03215623 26 1.595934 1.025645e-07 -0.05895398 27 0.6367598 1.025741e-07 -0.07956389 28 150.8911 1.025741e-07 0.0005790899 29 4.331202 1.025741e-07 -0.1353762 30 0.2746439 1.025741e-07 -0.5864849 31 1.641983 1.025741e-07 -0.1441606 32 11.5499 1.025741e-07 -0.002986789 33 0.2399144 1.025741e-07 0.03771472 34 0.6264627 1.025741e-07 0.2126323 35 2.004352 1.025741e-07 0.105643 36 1.436069 1.025741e-07 -0.06803593 37 0.6795059 1.025741e-07 0.6701851 38 0.7362047 1.025741e-07 0.04336908 39 0.1795409 1.025741e-07 0.2244795 40 0.8988927 1.025741e-07 -0.03370332 41 0.2960007 1.025741e-07 -0.5124151 42 0.6325087 1.025741e-07 -0.05151384 43 0.06631532 1.025741e-07 1.080041 44 0.2473832 1.025484e-07 0.09879578 45 -0.1488179 1.025741e-07 0.7040745 46 -0.02406017 1.025741e-07 -0.07207785 47 -0.2167138 1.025741e-07 -1.728853 48 -1.647116 1.025741e-07 0.003696711 49 0.09356013 1.025741e-07 0.3131578 50 1.080582 1.025709e-07 0.002114105 51 -0.4132841 1.025709e-07 0.002889524