1 |
|
/********************************************************************** |
2 |
|
|
3 |
|
This basic Periodic Table class was originally taken from the data.cpp |
4 |
< |
file in OpenBabel. The code has been modified to match the OOPSE coding style. |
4 |
> |
file in OpenBabel. The code has been modified to match the OpenMD coding style. |
5 |
|
|
6 |
|
We have retained the OpenBabel copyright and GPL license on this class: |
7 |
|
|
25 |
|
* @file ElementsTable.cpp |
26 |
|
* @author gezelter |
27 |
|
* @date 12/21/2007 |
28 |
– |
* @time 11:30am |
28 |
|
* @version 1.0 |
29 |
|
*/ |
30 |
|
|
31 |
|
#include "config.h" |
32 |
+ |
|
33 |
+ |
#include <iostream> |
34 |
|
#include <cstdlib> |
35 |
|
#include <string> |
36 |
|
#include <fstream> |
37 |
|
#include <cstdlib> |
38 |
|
#include "utils/ElementsTable.hpp" |
39 |
|
#include "utils/simError.h" |
40 |
< |
#include "io/basic_ifstrstream.hpp" |
40 |
> |
#include "io/ifstrstream.hpp" |
41 |
|
|
42 |
< |
#if !HAVE_STRNCASECMP |
43 |
< |
extern "C" int strncasecmp(const char *s1, const char *s2, size_t n); |
42 |
> |
#ifdef _MSC_VER |
43 |
> |
#define strncasecmp _strnicmp |
44 |
> |
#define strcasecmp _stricmp |
45 |
|
#endif |
46 |
|
|
47 |
+ |
|
48 |
|
#ifdef WIN32 |
49 |
|
#define FILE_SEP_CHAR "\\" |
50 |
|
#else |
55 |
|
#define BUFF_SIZE 32768 |
56 |
|
#endif |
57 |
|
|
58 |
< |
namespace oopse { |
58 |
> |
namespace OpenMD { |
59 |
|
|
60 |
|
ElementsTable etab; |
61 |
|
|
62 |
|
ElementsTable::ElementsTable() { |
63 |
|
init_ = false; |
64 |
< |
STR_DEFINE(dir_, FRC_PATH ); |
64 |
> |
dir_ = std::string("TO_STRING(FRC_PATH)"); |
65 |
|
envvar_ = "FORCE_PARAM_PATH"; |
66 |
|
filename_ = "element.txt"; |
67 |
|
} |
80 |
|
RealType red, green, blue; |
81 |
|
|
82 |
|
// skip comment line (at the top) |
80 |
– |
|
83 |
|
if (line[0] != '#') { |
84 |
|
sscanf(line,"%d %5s %lf %*f %lf %d %lf %lf %lf %lf %lf %lf %lf %255s", |
85 |
|
&num, |
100 |
|
elNeg, ionize, elAffin, red, green, blue, |
101 |
|
name); |
102 |
|
elements_.push_back(ele); |
103 |
+ |
|
104 |
|
} |
105 |
|
} |
106 |
|
|
111 |
|
return elements_.size(); |
112 |
|
} |
113 |
|
|
114 |
< |
char *ElementsTable::GetSymbol(int atomicnum) { |
114 |
> |
const char *ElementsTable::GetSymbol(int atomicnum) { |
115 |
|
if (!init_) |
116 |
|
Init(); |
117 |
|
|
296 |
|
buffer = getenv(envvar_.c_str()); |
297 |
|
buffer += FILE_SEP_CHAR; |
298 |
|
|
299 |
+ |
|
300 |
+ |
|
301 |
+ |
|
302 |
|
if (!subdir_.empty()) { |
303 |
|
subbuffer = buffer; |
304 |
|
subbuffer += subdir_; |
305 |
|
subbuffer += FILE_SEP_CHAR; |
306 |
|
} |
307 |
|
|
308 |
+ |
|
309 |
+ |
|
310 |
|
buffer += filename_; |
311 |
|
subbuffer += filename_; |
312 |
+ |
|
313 |
|
|
314 |
|
ifs1.open(subbuffer.c_str()); |
315 |
|
ifsP= &ifs1; |
316 |
< |
if (!(*ifsP)) { |
316 |
> |
if (!(ifsP->is_open())) { |
317 |
|
ifs2.open(buffer.c_str()); |
318 |
|
ifsP = &ifs2; |
319 |
|
} |
320 |
+ |
|
321 |
|
} else { |
322 |
|
sprintf( painCave.errMsg, |
323 |
|
"ElementsTable error.\n" |
330 |
|
if ((*ifsP)) { |
331 |
|
while(ifsP->getline(charBuffer,BUFF_SIZE)) |
332 |
|
ParseLine(charBuffer); |
333 |
< |
|
334 |
< |
if (ifs1) |
335 |
< |
ifs1.close(); |
336 |
< |
if (ifs2) |
337 |
< |
ifs2.close(); |
338 |
< |
if (ifs3) |
339 |
< |
ifs3.close(); |
340 |
< |
if (ifs4) |
341 |
< |
ifs4.close(); |
333 |
> |
|
334 |
> |
if (ifs1) |
335 |
> |
ifs1.close(); |
336 |
> |
if (ifs2) |
337 |
> |
ifs2.close(); |
338 |
> |
if (ifs3) |
339 |
> |
ifs3.close(); |
340 |
> |
if (ifs4) |
341 |
> |
ifs4.close(); |
342 |
> |
|
343 |
> |
if (GetSize() == 0) { |
344 |
> |
sprintf( painCave.errMsg, |
345 |
> |
"ElementsTable error.\n" |
346 |
> |
"\tCannot initialize database %s \n", filename_.c_str()); |
347 |
> |
painCave.isFatal = 0; |
348 |
> |
simError(); |
349 |
> |
} |
350 |
|
|
333 |
– |
if (GetSize() == 0) { |
334 |
– |
sprintf( painCave.errMsg, |
335 |
– |
"ElementsTable error.\n" |
336 |
– |
"\tCannot initialize database %s \n", filename_.c_str()); |
337 |
– |
painCave.isFatal = 0; |
338 |
– |
simError(); |
351 |
|
} |
340 |
– |
|
341 |
– |
} |
352 |
|
|
353 |
|
} |
354 |
|
} |