1 |
< |
#include <cstdlib> |
2 |
< |
#include <cstdio> |
3 |
< |
#include <cstring> |
1 |
> |
#include <stdlib.h> |
2 |
> |
#include <stdio.h> |
3 |
> |
#include <string.h> |
4 |
|
|
5 |
|
#include "Globals.hpp" |
6 |
|
#include "simError.h" |
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 |
> |
#define G_USEINIT_XS_STATE 35 |
57 |
> |
#define G_ORTHOBOXTOLERANCE 36 |
58 |
|
|
59 |
+ |
|
60 |
|
Globals::Globals(){ |
61 |
|
|
62 |
|
int i; |
63 |
+ |
|
64 |
+ |
hash_size = 23; |
65 |
+ |
hash_shift = 4; |
66 |
|
|
67 |
|
components = NULL; |
68 |
|
|
86 |
|
addHash( "boxY", G_BOXY ); |
87 |
|
addHash( "boxZ", G_BOXZ ); |
88 |
|
addHash( "sampleTime", G_SAMPLETIME ); |
89 |
+ |
addHash( "resetTime", G_RESETTIME ); |
90 |
|
addHash( "statusTime", G_STATUSTIME ); |
91 |
|
addHash( "electrostaticCutoffRadius", G_ECR ); |
92 |
|
addHash( "dielectric", G_DIELECTRIC ); |
101 |
|
addHash( "tauBarostat", G_TAUBAROSTAT); |
102 |
|
addHash( "zconsTime", G_ZCONSTIME); |
103 |
|
addHash( "nZconstraints", G_NZCONSTRAINTS); |
104 |
< |
addHash( "zconsTol", G_ZCONSTOL); |
104 |
> |
addHash( "zconsTol", G_ZCONSTOL); |
105 |
> |
addHash( "zconsForcePolicy", G_ZCONSFORCEPOLICY); |
106 |
> |
addHash( "seed", G_SEED); |
107 |
> |
addHash( "useInitialTime", G_USEINITTIME); |
108 |
> |
addHash( "useInitialExtendedSystemState", G_USEINIT_XS_STATE); |
109 |
> |
addHash( "orthoBoxTolerance", G_ORTHOBOXTOLERANCE); |
110 |
|
|
111 |
+ |
|
112 |
|
// define some default values |
113 |
|
|
114 |
|
strcpy( mixingRule,"standard"); //default mixing rules to standard. |
115 |
|
usePBC = 1; //default periodic boundry conditions to on |
116 |
|
useRF = 0; |
117 |
< |
|
117 |
> |
useInitTime = 1; // default to pull init time from the init file |
118 |
> |
useInitXSstate = 1; // default to pull the extended state from the init file |
119 |
> |
orthoBoxTolerance = 1E-6; |
120 |
> |
|
121 |
|
have_force_field = 0; |
122 |
|
have_n_components = 0; |
123 |
|
have_target_temp = 0; |
136 |
|
have_box_z = 0; |
137 |
|
have_sample_time = 0; |
138 |
|
have_status_time = 0; |
139 |
+ |
have_reset_time = 0; |
140 |
|
have_thermal_time = 0; |
141 |
|
have_ecr = 0; |
142 |
|
have_dielectric = 0; |
150 |
|
have_n_zConstraints = 0; |
151 |
|
have_zConstraints = 0; |
152 |
|
have_zcons_tol = 0; |
153 |
+ |
have_seed = 0; |
154 |
|
|
133 |
– |
|
134 |
– |
strcpy( initial_config, "./fooBar" ); |
135 |
– |
|
155 |
|
} |
156 |
|
|
157 |
|
Globals::~Globals(){ |
404 |
|
break; |
405 |
|
} |
406 |
|
break; |
407 |
+ |
|
408 |
+ |
case G_ORTHOBOXTOLERANCE: |
409 |
+ |
switch( the_type ){ |
410 |
+ |
|
411 |
+ |
case STRING: |
412 |
+ |
the_event->err_msg = |
413 |
+ |
strdup( "Global error. orthoBoxTolerance is not a double or int.\n" ); |
414 |
+ |
return 0; |
415 |
+ |
break; |
416 |
+ |
|
417 |
+ |
case DOUBLE: |
418 |
+ |
orthoBoxTolerance = the_event->evt.asmt.rhs.dval; |
419 |
+ |
have_target_temp = 1; |
420 |
+ |
return 1; |
421 |
+ |
break; |
422 |
+ |
|
423 |
+ |
case INT: |
424 |
+ |
orthoBoxTolerance = (double)the_event->evt.asmt.rhs.ival; |
425 |
+ |
have_target_temp = 1; |
426 |
+ |
return 1; |
427 |
+ |
break; |
428 |
+ |
|
429 |
+ |
default: |
430 |
+ |
the_event->err_msg = |
431 |
+ |
strdup( "Global error.orthoBoxTolerance unrecognized.\n" ); |
432 |
+ |
return 0; |
433 |
+ |
break; |
434 |
+ |
} |
435 |
+ |
break; |
436 |
|
|
437 |
|
case G_ENSEMBLE: |
438 |
|
if( the_type == STRING ){ |
766 |
|
default: |
767 |
|
the_event->err_msg = |
768 |
|
strdup( "Global error. statusTime unrecognized.\n" ); |
769 |
+ |
return 0; |
770 |
+ |
break; |
771 |
+ |
} |
772 |
+ |
break; |
773 |
+ |
|
774 |
+ |
case G_RESETTIME: |
775 |
+ |
switch( the_type ){ |
776 |
+ |
|
777 |
+ |
case STRING: |
778 |
+ |
the_event->err_msg = |
779 |
+ |
strdup( "Global error. resetTime is not a double or int.\n" ); |
780 |
|
return 0; |
781 |
|
break; |
782 |
+ |
|
783 |
+ |
case DOUBLE: |
784 |
+ |
resetTime = the_event->evt.asmt.rhs.dval; |
785 |
+ |
have_reset_time = 1; |
786 |
+ |
return 1; |
787 |
+ |
break; |
788 |
+ |
|
789 |
+ |
case INT: |
790 |
+ |
resetTime = (double)the_event->evt.asmt.rhs.ival; |
791 |
+ |
have_reset_time = 1; |
792 |
+ |
return 1; |
793 |
+ |
break; |
794 |
+ |
|
795 |
+ |
default: |
796 |
+ |
the_event->err_msg = |
797 |
+ |
strdup( "Global error. resetTime unrecognized.\n" ); |
798 |
+ |
return 0; |
799 |
+ |
break; |
800 |
|
} |
801 |
|
break; |
802 |
|
|
934 |
|
strdup( "Global error. tempSet was not \"true\" or \"false\".\n" ); |
935 |
|
return 0; |
936 |
|
break; |
937 |
+ |
|
938 |
+ |
case G_USEINITTIME: |
939 |
+ |
if( the_type == STRING ){ |
940 |
+ |
|
941 |
+ |
if( !strcasecmp( "true", the_event->evt.asmt.rhs.sval )) useInitTime = 1; |
942 |
+ |
else if( !strcasecmp( "false", the_event->evt.asmt.rhs.sval )) useInitTime = 0; |
943 |
+ |
else{ |
944 |
+ |
the_event->err_msg = |
945 |
+ |
strdup( "Global error. useInitTime was not \"true\" or \"false\".\n" ); |
946 |
+ |
return 0; |
947 |
+ |
} |
948 |
+ |
return 1; |
949 |
+ |
} |
950 |
|
|
951 |
+ |
the_event->err_msg = |
952 |
+ |
strdup( "Global error. useInitTime was not \"true\" or \"false\".\n" ); |
953 |
+ |
return 0; |
954 |
+ |
break; |
955 |
+ |
|
956 |
+ |
case G_USEINIT_XS_STATE: |
957 |
+ |
if( the_type == STRING ){ |
958 |
+ |
|
959 |
+ |
if( !strcasecmp( "true", the_event->evt.asmt.rhs.sval )) |
960 |
+ |
useInitXSstate = 1; |
961 |
+ |
else if( !strcasecmp( "false", the_event->evt.asmt.rhs.sval )) |
962 |
+ |
useInitXSstate = 0; |
963 |
+ |
else{ |
964 |
+ |
the_event->err_msg = |
965 |
+ |
strdup( "Global error. useInitExtendedSystemState was not \"true\" or \"false\".\n" ); |
966 |
+ |
return 0; |
967 |
+ |
} |
968 |
+ |
return 1; |
969 |
+ |
} |
970 |
+ |
|
971 |
+ |
the_event->err_msg = |
972 |
+ |
strdup( "Global error. useInitExtendedSystemState was not \"true\" or \"false\".\n" ); |
973 |
+ |
return 0; |
974 |
+ |
break; |
975 |
+ |
|
976 |
|
case G_USEPBC: |
977 |
|
if( the_type == STRING ){ |
978 |
|
|
1154 |
|
} |
1155 |
|
break; |
1156 |
|
|
1157 |
+ |
case G_ZCONSFORCEPOLICY: |
1158 |
+ |
switch( the_type ){ |
1159 |
+ |
|
1160 |
+ |
case STRING: |
1161 |
+ |
strcpy(zconsForcePolicy, the_event->evt.asmt.rhs.sval); |
1162 |
|
|
1163 |
+ |
for(int i = 0; zconsForcePolicy[i] != '\0'; i++) |
1164 |
+ |
{ |
1165 |
+ |
zconsForcePolicy[i] = toupper(zconsForcePolicy[i]); |
1166 |
+ |
} |
1167 |
+ |
have_zcons_force_policy = 1; |
1168 |
+ |
return 1; |
1169 |
+ |
break; |
1170 |
+ |
|
1171 |
+ |
case DOUBLE: |
1172 |
+ |
the_event->err_msg = |
1173 |
+ |
strdup( "Global error. zconsForcePolicy is not a double or int.\n" ); |
1174 |
+ |
return 0; |
1175 |
+ |
break; |
1176 |
+ |
|
1177 |
+ |
case INT: |
1178 |
+ |
the_event->err_msg = |
1179 |
+ |
strdup( "Global error. zconsForcePolicy is not a double or int.\n" ); |
1180 |
+ |
return 0; |
1181 |
+ |
break; |
1182 |
+ |
|
1183 |
+ |
default: |
1184 |
+ |
the_event->err_msg = |
1185 |
+ |
strdup( "Global error. zconsForcePolicy unrecognized.\n" ); |
1186 |
+ |
return 0; |
1187 |
+ |
break; |
1188 |
+ |
} |
1189 |
+ |
break; |
1190 |
|
// add more token cases here. |
1191 |
+ |
case G_SEED: |
1192 |
+ |
switch( the_type ){ |
1193 |
+ |
|
1194 |
+ |
case STRING: |
1195 |
+ |
the_event->err_msg = |
1196 |
+ |
strdup( "Global error. seed is not a string.\n" ); |
1197 |
+ |
return 0; |
1198 |
+ |
return 0; |
1199 |
+ |
break; |
1200 |
+ |
|
1201 |
+ |
case DOUBLE: |
1202 |
+ |
have_seed = 1; |
1203 |
+ |
seed = (int)the_event->evt.asmt.rhs.dval; |
1204 |
+ |
return 1; |
1205 |
+ |
break; |
1206 |
+ |
|
1207 |
+ |
case INT: |
1208 |
+ |
have_seed = 1; |
1209 |
+ |
seed = the_event->evt.asmt.rhs.ival ; |
1210 |
+ |
return 1; |
1211 |
+ |
break; |
1212 |
+ |
|
1213 |
+ |
default: |
1214 |
+ |
the_event->err_msg = |
1215 |
+ |
strdup( "Global error. seed unrecognized.\n" ); |
1216 |
+ |
return 0; |
1217 |
+ |
break; |
1218 |
+ |
} |
1219 |
+ |
break; |
1220 |
+ |
// add more token cases here. |
1221 |
|
|
1222 |
|
} |
1223 |
|
} |