Commit b9992edc authored by COULAUD Olivier's avatar COULAUD Olivier

Add new data to check the periodicity

parent e0bf4435
Periodic Data
To check the periodic Algorithm, we consider a crystal of NaCL. The computation are done by
DL_POLY molecular dynamic code with the Ewald method. We only have the energy and the forces.
To compare with our approach we removed the Electric moment and the new data are in the file
called forceNacl_2000.bfma. The potential is wrong
The new energy is Energy 7.4014682651e+01
Tests/Release/removeMoment -Ewald2FMM -fin ../Data/forceNacl_2000_dlpolyPer.bin -dlpoly -fout ../Data/forceNacl_2000.fma
Energy EWALD = -3.6662545891e+05
Energy New = 7.4014682651e+01
\ No newline at end of file
......@@ -232,7 +232,7 @@ public:
file >> index;
std::getline(file, line); // needed to skip the end of the line in non periodic case
std::cout << "line: " << line << std::endl;
// std::cout << "line: " << line << std::endl;
if ( levcfg == 0) {
file >> x >> y >> z;
}else if ( levcfg == 1) {
Markdown is supported
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment