| 48 |
|
#define G_TAUBAROSTAT 27 |
| 49 |
|
#define G_ZCONSTIME 28 |
| 50 |
|
#define G_NZCONSTRAINTS 29 |
| 51 |
< |
#define G_ZCONSTOL 30 |
| 51 |
> |
#define G_ZCONSTOL 30 |
| 52 |
> |
#define G_ZCONSFORCEPOLICY 31 |
| 53 |
> |
#define G_SEED 32 |
| 54 |
> |
#define G_RESETTIME 33 |
| 55 |
> |
#define G_USEINITTIME 34 |
| 56 |
|
|
| 57 |
|
Globals::Globals(){ |
| 58 |
|
|
| 59 |
|
int i; |
| 60 |
+ |
|
| 61 |
+ |
hash_size = 23; |
| 62 |
+ |
hash_shift = 4; |
| 63 |
|
|
| 64 |
|
components = NULL; |
| 65 |
|
|
| 83 |
|
addHash( "boxY", G_BOXY ); |
| 84 |
|
addHash( "boxZ", G_BOXZ ); |
| 85 |
|
addHash( "sampleTime", G_SAMPLETIME ); |
| 86 |
+ |
addHash( "resetTime", G_RESETTIME ); |
| 87 |
|
addHash( "statusTime", G_STATUSTIME ); |
| 88 |
|
addHash( "electrostaticCutoffRadius", G_ECR ); |
| 89 |
|
addHash( "dielectric", G_DIELECTRIC ); |
| 98 |
|
addHash( "tauBarostat", G_TAUBAROSTAT); |
| 99 |
|
addHash( "zconsTime", G_ZCONSTIME); |
| 100 |
|
addHash( "nZconstraints", G_NZCONSTRAINTS); |
| 101 |
< |
addHash( "zconsTol", G_ZCONSTOL); |
| 101 |
> |
addHash( "zconsTol", G_ZCONSTOL); |
| 102 |
> |
addHash( "zconsForcePolicy", G_ZCONSFORCEPOLICY); |
| 103 |
> |
addHash( "seed", G_SEED); |
| 104 |
> |
addHash( "useInitTime", G_USEINITTIME), |
| 105 |
|
|
| 106 |
+ |
|
| 107 |
|
// define some default values |
| 108 |
|
|
| 109 |
|
strcpy( mixingRule,"standard"); //default mixing rules to standard. |
| 110 |
|
usePBC = 1; //default periodic boundry conditions to on |
| 111 |
|
useRF = 0; |
| 112 |
+ |
useInitTime = 1; // default to pull init time from the init file |
| 113 |
|
|
| 114 |
|
have_force_field = 0; |
| 115 |
|
have_n_components = 0; |
| 129 |
|
have_box_z = 0; |
| 130 |
|
have_sample_time = 0; |
| 131 |
|
have_status_time = 0; |
| 132 |
+ |
have_reset_time = 0; |
| 133 |
|
have_thermal_time = 0; |
| 134 |
|
have_ecr = 0; |
| 135 |
|
have_dielectric = 0; |
| 143 |
|
have_n_zConstraints = 0; |
| 144 |
|
have_zConstraints = 0; |
| 145 |
|
have_zcons_tol = 0; |
| 146 |
< |
} |
| 146 |
> |
have_seed = 0; |
| 147 |
|
|
| 148 |
+ |
} |
| 149 |
+ |
|
| 150 |
|
Globals::~Globals(){ |
| 151 |
|
int i; |
| 152 |
|
|
| 264 |
|
|
| 265 |
|
case STRING: |
| 266 |
|
return current_zConstraint->assignString( the_event->evt.asmt.lhs, |
| 267 |
< |
the_event->evt.asmt.rhs.sval, |
| 268 |
< |
&(the_event->err_msg)); |
| 267 |
> |
the_event->evt.asmt.rhs.sval, |
| 268 |
> |
&(the_event->err_msg)); |
| 269 |
|
break; |
| 270 |
|
|
| 271 |
|
case DOUBLE: |
| 272 |
|
return current_zConstraint->assignDouble( the_event->evt.asmt.lhs, |
| 273 |
< |
the_event->evt.asmt.rhs.dval, |
| 274 |
< |
&(the_event->err_msg)); |
| 273 |
> |
the_event->evt.asmt.rhs.dval, |
| 274 |
> |
&(the_event->err_msg)); |
| 275 |
|
break; |
| 276 |
|
|
| 277 |
|
case INT: |
| 278 |
|
return current_zConstraint->assignInt( the_event->evt.asmt.lhs, |
| 279 |
< |
the_event->evt.asmt.rhs.ival, |
| 280 |
< |
&(the_event->err_msg)); |
| 279 |
> |
the_event->evt.asmt.rhs.ival, |
| 280 |
> |
&(the_event->err_msg)); |
| 281 |
|
break; |
| 282 |
|
|
| 283 |
|
default: |
| 735 |
|
} |
| 736 |
|
break; |
| 737 |
|
|
| 738 |
< |
case G_THERMALTIME: |
| 738 |
> |
case G_RESETTIME: |
| 739 |
> |
switch( the_type ){ |
| 740 |
> |
|
| 741 |
> |
case STRING: |
| 742 |
> |
the_event->err_msg = |
| 743 |
> |
strdup( "Global error. resetTime is not a double or int.\n" ); |
| 744 |
> |
return 0; |
| 745 |
> |
break; |
| 746 |
> |
|
| 747 |
> |
case DOUBLE: |
| 748 |
> |
resetTime = the_event->evt.asmt.rhs.dval; |
| 749 |
> |
have_reset_time = 1; |
| 750 |
> |
return 1; |
| 751 |
> |
break; |
| 752 |
> |
|
| 753 |
> |
case INT: |
| 754 |
> |
resetTime = (double)the_event->evt.asmt.rhs.ival; |
| 755 |
> |
have_reset_time = 1; |
| 756 |
> |
return 1; |
| 757 |
> |
break; |
| 758 |
> |
|
| 759 |
> |
default: |
| 760 |
> |
the_event->err_msg = |
| 761 |
> |
strdup( "Global error. resetTime unrecognized.\n" ); |
| 762 |
> |
return 0; |
| 763 |
> |
break; |
| 764 |
> |
} |
| 765 |
> |
break; |
| 766 |
> |
|
| 767 |
> |
case G_THERMALTIME: |
| 768 |
|
switch( the_type ){ |
| 769 |
|
|
| 770 |
|
case STRING: |
| 898 |
|
strdup( "Global error. tempSet was not \"true\" or \"false\".\n" ); |
| 899 |
|
return 0; |
| 900 |
|
break; |
| 901 |
+ |
|
| 902 |
+ |
case G_USEINITTIME: |
| 903 |
+ |
if( the_type == STRING ){ |
| 904 |
+ |
|
| 905 |
+ |
if( !strcasecmp( "true", the_event->evt.asmt.rhs.sval )) useInitTime = 1; |
| 906 |
+ |
else if( !strcasecmp( "false", the_event->evt.asmt.rhs.sval )) useInitTime = 0; |
| 907 |
+ |
else{ |
| 908 |
+ |
the_event->err_msg = |
| 909 |
+ |
strdup( "Global error. useInitTime was not \"true\" or \"false\".\n" ); |
| 910 |
+ |
return 0; |
| 911 |
+ |
} |
| 912 |
+ |
return 1; |
| 913 |
+ |
} |
| 914 |
|
|
| 915 |
+ |
the_event->err_msg = |
| 916 |
+ |
strdup( "Global error. useInitTime was not \"true\" or \"false\".\n" ); |
| 917 |
+ |
return 0; |
| 918 |
+ |
break; |
| 919 |
+ |
|
| 920 |
|
case G_USEPBC: |
| 921 |
|
if( the_type == STRING ){ |
| 922 |
|
|
| 1040 |
|
} |
| 1041 |
|
break; |
| 1042 |
|
|
| 1043 |
< |
case G_ZCONSTIME: |
| 1043 |
> |
case G_ZCONSTIME: |
| 1044 |
|
switch( the_type ){ |
| 1045 |
|
|
| 1046 |
|
case STRING: |
| 1069 |
|
} |
| 1070 |
|
break; |
| 1071 |
|
|
| 1072 |
< |
case G_ZCONSTOL: |
| 1072 |
> |
case G_ZCONSTOL: |
| 1073 |
|
switch( the_type ){ |
| 1074 |
|
|
| 1075 |
|
case STRING: |
| 1098 |
|
} |
| 1099 |
|
break; |
| 1100 |
|
|
| 1101 |
+ |
case G_ZCONSFORCEPOLICY: |
| 1102 |
+ |
switch( the_type ){ |
| 1103 |
+ |
|
| 1104 |
+ |
case STRING: |
| 1105 |
+ |
strcpy(zconsForcePolicy, the_event->evt.asmt.rhs.sval); |
| 1106 |
|
|
| 1107 |
+ |
for(int i = 0; zconsForcePolicy[i] != '\0'; i++) |
| 1108 |
+ |
{ |
| 1109 |
+ |
zconsForcePolicy[i] = toupper(zconsForcePolicy[i]); |
| 1110 |
+ |
} |
| 1111 |
+ |
have_zcons_force_policy = 1; |
| 1112 |
+ |
return 1; |
| 1113 |
+ |
break; |
| 1114 |
+ |
|
| 1115 |
+ |
case DOUBLE: |
| 1116 |
+ |
the_event->err_msg = |
| 1117 |
+ |
strdup( "Global error. zconsForcePolicy is not a double or int.\n" ); |
| 1118 |
+ |
return 0; |
| 1119 |
+ |
break; |
| 1120 |
+ |
|
| 1121 |
+ |
case INT: |
| 1122 |
+ |
the_event->err_msg = |
| 1123 |
+ |
strdup( "Global error. zconsForcePolicy is not a double or int.\n" ); |
| 1124 |
+ |
return 0; |
| 1125 |
+ |
break; |
| 1126 |
+ |
|
| 1127 |
+ |
default: |
| 1128 |
+ |
the_event->err_msg = |
| 1129 |
+ |
strdup( "Global error. zconsForcePolicy unrecognized.\n" ); |
| 1130 |
+ |
return 0; |
| 1131 |
+ |
break; |
| 1132 |
+ |
} |
| 1133 |
+ |
break; |
| 1134 |
|
// add more token cases here. |
| 1135 |
+ |
case G_SEED: |
| 1136 |
+ |
switch( the_type ){ |
| 1137 |
+ |
|
| 1138 |
+ |
case STRING: |
| 1139 |
+ |
the_event->err_msg = |
| 1140 |
+ |
strdup( "Global error. seed is not a string.\n" ); |
| 1141 |
+ |
return 0; |
| 1142 |
+ |
return 0; |
| 1143 |
+ |
break; |
| 1144 |
+ |
|
| 1145 |
+ |
case DOUBLE: |
| 1146 |
+ |
have_seed = 1; |
| 1147 |
+ |
seed = (int)the_event->evt.asmt.rhs.dval; |
| 1148 |
+ |
return 1; |
| 1149 |
+ |
break; |
| 1150 |
+ |
|
| 1151 |
+ |
case INT: |
| 1152 |
+ |
have_seed = 1; |
| 1153 |
+ |
seed = the_event->evt.asmt.rhs.ival ; |
| 1154 |
+ |
return 1; |
| 1155 |
+ |
break; |
| 1156 |
+ |
|
| 1157 |
+ |
default: |
| 1158 |
+ |
the_event->err_msg = |
| 1159 |
+ |
strdup( "Global error. seed unrecognized.\n" ); |
| 1160 |
+ |
return 0; |
| 1161 |
+ |
break; |
| 1162 |
+ |
} |
| 1163 |
+ |
break; |
| 1164 |
+ |
// add more token cases here. |
| 1165 |
|
|
| 1166 |
|
} |
| 1167 |
|
} |