| 56 |
|
#define G_USEINIT_XS_STATE 35 |
| 57 |
|
#define G_ORTHOBOXTOLERANCE 36 |
| 58 |
|
#define G_LJRCUT 37 |
| 59 |
+ |
#define G_MINIMIZER 38 |
| 60 |
+ |
#define G_MIN_MAXITER 39 |
| 61 |
+ |
#define G_MIN_WRITEFRQ 40 |
| 62 |
+ |
#define G_MIN_RESETFRQ 41 |
| 63 |
+ |
#define G_MIN_FTOL 42 |
| 64 |
+ |
#define G_MIN_GTOL 43 |
| 65 |
+ |
#define G_MIN_LSTOL 44 |
| 66 |
+ |
#define G_MIN_LSMAXITER 45 |
| 67 |
|
|
| 68 |
|
|
| 69 |
|
Globals::Globals(){ |
| 117 |
|
addHash( "useInitialExtendedSystemState", G_USEINIT_XS_STATE); |
| 118 |
|
addHash( "orthoBoxTolerance", G_ORTHOBOXTOLERANCE); |
| 119 |
|
addHash( "LJrcut", G_LJRCUT); |
| 120 |
< |
|
| 121 |
< |
|
| 122 |
< |
// define some default values |
| 120 |
> |
addHash( "minimizer", G_MINIMIZER); |
| 121 |
> |
addHash( "minMaxIter", G_MIN_MAXITER); |
| 122 |
> |
addHash( "minWriteFrq", G_MIN_WRITEFRQ); |
| 123 |
> |
addHash( "minResetFrq", G_MIN_RESETFRQ); |
| 124 |
> |
addHash( "minFTol", G_MIN_FTOL); |
| 125 |
> |
addHash( "minGTol", G_MIN_GTOL); |
| 126 |
> |
addHash( "minLSTol", G_MIN_LSTOL); |
| 127 |
> |
addHash( "minLSMaxIter", G_MIN_LSMAXITER); |
| 128 |
> |
|
| 129 |
|
|
| 130 |
|
strcpy( mixingRule,"standard"); //default mixing rules to standard. |
| 131 |
|
usePBC = 1; //default periodic boundry conditions to on |
| 1231 |
|
strdup( "Global error. zconsForcePolicy unrecognized.\n" ); |
| 1232 |
|
return 0; |
| 1233 |
|
break; |
| 1234 |
+ |
} |
| 1235 |
+ |
break; |
| 1236 |
+ |
|
| 1237 |
+ |
case G_MINIMIZER: |
| 1238 |
+ |
switch( the_type ){ |
| 1239 |
+ |
|
| 1240 |
+ |
case STRING: |
| 1241 |
+ |
strcpy(minimizer_name, the_event->evt.asmt.rhs.sval); |
| 1242 |
+ |
|
| 1243 |
+ |
for(int i = 0; zconsForcePolicy[i] != '\0'; i++){ |
| 1244 |
+ |
zconsForcePolicy[i] = toupper(zconsForcePolicy[i]); |
| 1245 |
+ |
} |
| 1246 |
+ |
have_minimizer= 1; |
| 1247 |
+ |
return 1; |
| 1248 |
+ |
break; |
| 1249 |
+ |
|
| 1250 |
+ |
case DOUBLE: |
| 1251 |
+ |
the_event->err_msg = |
| 1252 |
+ |
strdup( "Global error. minimizer_name is not a double or int.\n" ); |
| 1253 |
+ |
return 0; |
| 1254 |
+ |
break; |
| 1255 |
+ |
|
| 1256 |
+ |
case INT: |
| 1257 |
+ |
the_event->err_msg = |
| 1258 |
+ |
strdup( "Global error. minimizer_name is not a double or int.\n" ); |
| 1259 |
+ |
return 0; |
| 1260 |
+ |
break; |
| 1261 |
+ |
|
| 1262 |
+ |
default: |
| 1263 |
+ |
the_event->err_msg = |
| 1264 |
+ |
strdup( "Global error. minimizer_name unrecognized.\n" ); |
| 1265 |
+ |
return 0; |
| 1266 |
+ |
break; |
| 1267 |
+ |
} |
| 1268 |
+ |
break; |
| 1269 |
+ |
|
| 1270 |
+ |
case G_MIN_MAXITER: |
| 1271 |
+ |
switch( the_type ){ |
| 1272 |
+ |
|
| 1273 |
+ |
case STRING: |
| 1274 |
+ |
the_event->err_msg = |
| 1275 |
+ |
strdup( "Global error. minimizer_maxiteration is not a double or int.\n" ); |
| 1276 |
+ |
return 1; |
| 1277 |
+ |
break; |
| 1278 |
+ |
|
| 1279 |
+ |
case DOUBLE: |
| 1280 |
+ |
minimizer_maxiteration = the_event->evt.asmt.rhs.dval; |
| 1281 |
+ |
have_minimizer_maxiteration = 1; |
| 1282 |
+ |
return 1; |
| 1283 |
+ |
break; |
| 1284 |
+ |
|
| 1285 |
+ |
case INT: |
| 1286 |
+ |
minimizer_maxiteration = the_event->evt.asmt.rhs.ival; |
| 1287 |
+ |
have_minimizer_maxiteration = 1; |
| 1288 |
+ |
return 1; |
| 1289 |
+ |
break; |
| 1290 |
+ |
|
| 1291 |
+ |
default: |
| 1292 |
+ |
the_event->err_msg = |
| 1293 |
+ |
strdup( "Global error. minimizer_maxiteration unrecognized.\n" ); |
| 1294 |
+ |
return 0; |
| 1295 |
+ |
break; |
| 1296 |
|
} |
| 1297 |
|
break; |
| 1298 |
+ |
|
| 1299 |
+ |
case G_MIN_WRITEFRQ: |
| 1300 |
+ |
switch( the_type ){ |
| 1301 |
+ |
|
| 1302 |
+ |
case STRING: |
| 1303 |
+ |
the_event->err_msg = |
| 1304 |
+ |
strdup( "Global error. minimizer_writefrq is not a double or int.\n" ); |
| 1305 |
+ |
return 1; |
| 1306 |
+ |
break; |
| 1307 |
+ |
|
| 1308 |
+ |
case DOUBLE: |
| 1309 |
+ |
minimizer_writefrq= the_event->evt.asmt.rhs.dval; |
| 1310 |
+ |
have_minimizer_writefrq = 1; |
| 1311 |
+ |
return 1; |
| 1312 |
+ |
break; |
| 1313 |
+ |
|
| 1314 |
+ |
case INT: |
| 1315 |
+ |
minimizer_writefrq= the_event->evt.asmt.rhs.ival; |
| 1316 |
+ |
have_minimizer_writefrq = 1; |
| 1317 |
+ |
return 1; |
| 1318 |
+ |
break; |
| 1319 |
+ |
|
| 1320 |
+ |
default: |
| 1321 |
+ |
the_event->err_msg = |
| 1322 |
+ |
strdup( "Global error. minimizer_writefrq unrecognized.\n" ); |
| 1323 |
+ |
return 0; |
| 1324 |
+ |
break; |
| 1325 |
+ |
} |
| 1326 |
+ |
break; |
| 1327 |
+ |
|
| 1328 |
+ |
case G_MIN_RESETFRQ: |
| 1329 |
+ |
switch( the_type ){ |
| 1330 |
+ |
|
| 1331 |
+ |
case STRING: |
| 1332 |
+ |
the_event->err_msg = |
| 1333 |
+ |
strdup( "Global error. minimizer_resetfrq is not a double or int.\n" ); |
| 1334 |
+ |
return 1; |
| 1335 |
+ |
break; |
| 1336 |
+ |
|
| 1337 |
+ |
case DOUBLE: |
| 1338 |
+ |
minimizer_resetfrq= the_event->evt.asmt.rhs.dval; |
| 1339 |
+ |
have_minimizer_resetfrq = 1; |
| 1340 |
+ |
return 1; |
| 1341 |
+ |
break; |
| 1342 |
+ |
|
| 1343 |
+ |
case INT: |
| 1344 |
+ |
minimizer_resetfrq= the_event->evt.asmt.rhs.ival; |
| 1345 |
+ |
have_minimizer_resetfrq = 1; |
| 1346 |
+ |
return 1; |
| 1347 |
+ |
break; |
| 1348 |
+ |
|
| 1349 |
+ |
default: |
| 1350 |
+ |
the_event->err_msg = |
| 1351 |
+ |
strdup( "Global error. minimizer_resetfrq unrecognized.\n" ); |
| 1352 |
+ |
return 0; |
| 1353 |
+ |
break; |
| 1354 |
+ |
} |
| 1355 |
+ |
break; |
| 1356 |
+ |
|
| 1357 |
+ |
case G_MIN_FTOL: |
| 1358 |
+ |
switch( the_type ){ |
| 1359 |
+ |
|
| 1360 |
+ |
case STRING: |
| 1361 |
+ |
the_event->err_msg = |
| 1362 |
+ |
strdup( "Global error. minimizer_ftol is not a double or int.\n" ); |
| 1363 |
+ |
return 1; |
| 1364 |
+ |
break; |
| 1365 |
+ |
|
| 1366 |
+ |
case DOUBLE: |
| 1367 |
+ |
minimizer_ftol= the_event->evt.asmt.rhs.dval; |
| 1368 |
+ |
have_minimizer_ftol = 1; |
| 1369 |
+ |
return 1; |
| 1370 |
+ |
break; |
| 1371 |
+ |
|
| 1372 |
+ |
case INT: |
| 1373 |
+ |
minimizer_ftol= the_event->evt.asmt.rhs.ival; |
| 1374 |
+ |
have_minimizer_ftol = 1; |
| 1375 |
+ |
return 1; |
| 1376 |
+ |
break; |
| 1377 |
+ |
|
| 1378 |
+ |
default: |
| 1379 |
+ |
the_event->err_msg = |
| 1380 |
+ |
strdup( "Global error. minimizer_ftol unrecognized.\n" ); |
| 1381 |
+ |
return 0; |
| 1382 |
+ |
break; |
| 1383 |
+ |
} |
| 1384 |
+ |
break; |
| 1385 |
+ |
|
| 1386 |
+ |
case G_MIN_GTOL: |
| 1387 |
+ |
switch( the_type ){ |
| 1388 |
+ |
|
| 1389 |
+ |
case STRING: |
| 1390 |
+ |
the_event->err_msg = |
| 1391 |
+ |
strdup( "Global error. minimizer_gtol is not a double or int.\n" ); |
| 1392 |
+ |
return 1; |
| 1393 |
+ |
break; |
| 1394 |
+ |
|
| 1395 |
+ |
case DOUBLE: |
| 1396 |
+ |
minimizer_gtol= the_event->evt.asmt.rhs.dval; |
| 1397 |
+ |
have_minimizer_gtol = 1; |
| 1398 |
+ |
return 1; |
| 1399 |
+ |
break; |
| 1400 |
+ |
|
| 1401 |
+ |
case INT: |
| 1402 |
+ |
minimizer_gtol= the_event->evt.asmt.rhs.ival; |
| 1403 |
+ |
have_minimizer_gtol = 1; |
| 1404 |
+ |
return 1; |
| 1405 |
+ |
break; |
| 1406 |
+ |
|
| 1407 |
+ |
default: |
| 1408 |
+ |
the_event->err_msg = |
| 1409 |
+ |
strdup( "Global error. minimizer_gtol unrecognized.\n" ); |
| 1410 |
+ |
return 0; |
| 1411 |
+ |
break; |
| 1412 |
+ |
} |
| 1413 |
+ |
break; |
| 1414 |
+ |
|
| 1415 |
+ |
case G_MIN_LSMAXITER: |
| 1416 |
+ |
switch( the_type ){ |
| 1417 |
+ |
|
| 1418 |
+ |
case STRING: |
| 1419 |
+ |
the_event->err_msg = |
| 1420 |
+ |
strdup( "Global error. minimizer_ls_maxiteration is not a double or int.\n" ); |
| 1421 |
+ |
return 1; |
| 1422 |
+ |
break; |
| 1423 |
+ |
|
| 1424 |
+ |
case DOUBLE: |
| 1425 |
+ |
minimizer_ls_maxiteration = the_event->evt.asmt.rhs.dval; |
| 1426 |
+ |
have_minimizer_ls_maxiteration = 1; |
| 1427 |
+ |
return 1; |
| 1428 |
+ |
break; |
| 1429 |
+ |
|
| 1430 |
+ |
case INT: |
| 1431 |
+ |
minimizer_ls_maxiteration = the_event->evt.asmt.rhs.ival; |
| 1432 |
+ |
have_minimizer_ls_maxiteration = 1; |
| 1433 |
+ |
return 1; |
| 1434 |
+ |
break; |
| 1435 |
+ |
|
| 1436 |
+ |
default: |
| 1437 |
+ |
the_event->err_msg = |
| 1438 |
+ |
strdup( "Global error. minimizer_ls_maxiteration unrecognized.\n" ); |
| 1439 |
+ |
return 0; |
| 1440 |
+ |
break; |
| 1441 |
+ |
} |
| 1442 |
+ |
break; |
| 1443 |
+ |
|
| 1444 |
+ |
case G_MIN_LSTOL: |
| 1445 |
+ |
switch( the_type ){ |
| 1446 |
+ |
|
| 1447 |
+ |
case STRING: |
| 1448 |
+ |
the_event->err_msg = |
| 1449 |
+ |
strdup( "Global error. minimizer_ls_tol is not a double or int.\n" ); |
| 1450 |
+ |
return 1; |
| 1451 |
+ |
break; |
| 1452 |
+ |
|
| 1453 |
+ |
case DOUBLE: |
| 1454 |
+ |
minimizer_ls_tol= the_event->evt.asmt.rhs.dval; |
| 1455 |
+ |
have_minimizer_ls_tol = 1; |
| 1456 |
+ |
return 1; |
| 1457 |
+ |
break; |
| 1458 |
+ |
|
| 1459 |
+ |
case INT: |
| 1460 |
+ |
minimizer_ls_tol= the_event->evt.asmt.rhs.ival; |
| 1461 |
+ |
have_minimizer_ls_tol = 1; |
| 1462 |
+ |
return 1; |
| 1463 |
+ |
break; |
| 1464 |
+ |
|
| 1465 |
+ |
default: |
| 1466 |
+ |
the_event->err_msg = |
| 1467 |
+ |
strdup( "Global error. minimizer_ls_tol unrecognized.\n" ); |
| 1468 |
+ |
return 0; |
| 1469 |
+ |
break; |
| 1470 |
+ |
} |
| 1471 |
+ |
break; |
| 1472 |
+ |
|
| 1473 |
|
// add more token cases here. |
| 1474 |
|
case G_SEED: |
| 1475 |
|
switch( the_type ){ |