1565 lines
63 KiB
C++
1565 lines
63 KiB
C++
/* ----------------------------------------------------------------------
|
|
LAMMPS - Large-scale Atomic/Molecular Massively Parallel Simulator
|
|
https://www.lammps.org/, Sandia National Laboratories
|
|
Steve Plimpton, sjplimp@sandia.gov
|
|
|
|
Copyright (2003) Sandia Corporation. Under the terms of Contract
|
|
DE-AC04-94AL85000 with Sandia Corporation, the U.S. Government retains
|
|
certain rights in this software. This software is distributed under
|
|
the GNU General Public License.
|
|
|
|
See the README file in the top-level LAMMPS directory.
|
|
------------------------------------------------------------------------- */
|
|
|
|
// unit tests for pair styles intended for molecular systems
|
|
|
|
#include "error_stats.h"
|
|
#include "test_config.h"
|
|
#include "test_config_reader.h"
|
|
#include "test_main.h"
|
|
#include "yaml_reader.h"
|
|
#include "yaml_writer.h"
|
|
|
|
#include "gmock/gmock.h"
|
|
#include "gtest/gtest.h"
|
|
|
|
#include "atom.h"
|
|
#include "compute.h"
|
|
#include "force.h"
|
|
#include "info.h"
|
|
#include "input.h"
|
|
#include "kspace.h"
|
|
#include "lammps.h"
|
|
#include "modify.h"
|
|
#include "pair.h"
|
|
#include "universe.h"
|
|
#include "utils.h"
|
|
|
|
#include <cctype>
|
|
#include <cstdio>
|
|
#include <cstdlib>
|
|
#include <cstring>
|
|
#include <mpi.h>
|
|
|
|
#include <map>
|
|
#include <string>
|
|
#include <utility>
|
|
#include <vector>
|
|
|
|
using ::testing::HasSubstr;
|
|
using ::testing::StartsWith;
|
|
|
|
using namespace LAMMPS_NS;
|
|
|
|
static void delete_file(const std::string &filename)
|
|
{
|
|
remove(filename.c_str());
|
|
};
|
|
|
|
void cleanup_lammps(LAMMPS *lmp, const TestConfig &cfg)
|
|
{
|
|
delete_file(cfg.basename + ".restart");
|
|
delete_file(cfg.basename + ".data");
|
|
delete_file(cfg.basename + "-coeffs.in");
|
|
delete lmp;
|
|
lmp = nullptr;
|
|
}
|
|
|
|
LAMMPS *init_lammps(int argc, char **argv, const TestConfig &cfg, const bool newton = true)
|
|
{
|
|
LAMMPS *lmp;
|
|
|
|
lmp = new LAMMPS(argc, argv, MPI_COMM_WORLD);
|
|
|
|
// check if prerequisite styles are available
|
|
Info *info = new Info(lmp);
|
|
int nfail = 0;
|
|
for (auto &prerequisite : cfg.prerequisites) {
|
|
std::string style = prerequisite.second;
|
|
|
|
// this is a test for pair styles, so if the suffixed
|
|
// version is not available, there is no reason to test.
|
|
if (prerequisite.first == "pair") {
|
|
if (lmp->suffix_enable) {
|
|
style += "/";
|
|
style += lmp->suffix;
|
|
}
|
|
}
|
|
|
|
if (!info->has_style(prerequisite.first, style)) ++nfail;
|
|
}
|
|
delete info;
|
|
if (nfail > 0) {
|
|
cleanup_lammps(lmp, cfg);
|
|
return nullptr;
|
|
}
|
|
|
|
// utility lambdas to improve readability
|
|
auto command = [&](const std::string &line) {
|
|
lmp->input->one(line.c_str());
|
|
};
|
|
auto parse_input_script = [&](const std::string &filename) {
|
|
lmp->input->file(filename.c_str());
|
|
};
|
|
|
|
if (newton) {
|
|
command("variable newton_pair index on");
|
|
} else {
|
|
command("variable newton_pair index off");
|
|
}
|
|
|
|
command("variable input_dir index " + INPUT_FOLDER);
|
|
|
|
for (auto &pre_command : cfg.pre_commands) {
|
|
command(pre_command);
|
|
}
|
|
|
|
std::string input_file = INPUT_FOLDER + PATH_SEP + cfg.input_file;
|
|
parse_input_script(input_file);
|
|
|
|
command("pair_style " + cfg.pair_style);
|
|
|
|
for (auto &pair_coeff : cfg.pair_coeff) {
|
|
command("pair_coeff " + pair_coeff);
|
|
}
|
|
|
|
// set this here explicitly to a setting different
|
|
// from the default, so we can spot YAML files for
|
|
// long-range interactions that do not include these
|
|
// settings. they will fail after restart or read data.
|
|
command("pair_modify table 0");
|
|
command("pair_modify table/disp 0");
|
|
|
|
for (auto &post_command : cfg.post_commands) {
|
|
command(post_command);
|
|
}
|
|
|
|
command("run 0 post no");
|
|
command("variable write_data_pair index ii");
|
|
command("write_restart " + cfg.basename + ".restart");
|
|
command("write_data " + cfg.basename + ".data pair ${write_data_pair}");
|
|
command("write_coeff " + cfg.basename + "-coeffs.in");
|
|
|
|
return lmp;
|
|
}
|
|
|
|
void run_lammps(LAMMPS *lmp)
|
|
{
|
|
// utility lambda to improve readability
|
|
auto command = [&](const std::string &line) {
|
|
lmp->input->one(line.c_str());
|
|
};
|
|
|
|
command("fix 1 all nve");
|
|
command("compute pe all pe/atom pair");
|
|
command("compute sum all reduce sum c_pe");
|
|
command("thermo_style custom step temp pe press c_sum");
|
|
command("thermo 2");
|
|
command("run 4 post no");
|
|
}
|
|
|
|
void restart_lammps(LAMMPS *lmp, const TestConfig &cfg, bool nofdotr = false, bool newton = true)
|
|
{
|
|
// utility lambda to improve readability
|
|
auto command = [&](const std::string &line) {
|
|
lmp->input->one(line.c_str());
|
|
};
|
|
|
|
command("clear");
|
|
if (newton)
|
|
command("newton on");
|
|
else
|
|
command("newton off");
|
|
command("read_restart " + cfg.basename + ".restart");
|
|
|
|
if (!lmp->force->pair) {
|
|
command("pair_style " + cfg.pair_style);
|
|
}
|
|
if (!lmp->force->pair->restartinfo || !lmp->force->pair->writedata) {
|
|
for (auto &pair_coeff : cfg.pair_coeff) {
|
|
command("pair_coeff " + pair_coeff);
|
|
}
|
|
}
|
|
|
|
for (auto &post_command : cfg.post_commands) {
|
|
command(post_command);
|
|
}
|
|
if (nofdotr) command("pair_modify nofdotr");
|
|
|
|
command("run 0 post no");
|
|
}
|
|
|
|
void data_lammps(LAMMPS *lmp, const TestConfig &cfg)
|
|
{
|
|
// utility lambdas to improve readability
|
|
auto command = [&](const std::string &line) {
|
|
lmp->input->one(line.c_str());
|
|
};
|
|
auto parse_input_script = [&](const std::string &filename) {
|
|
lmp->input->file(filename.c_str());
|
|
};
|
|
|
|
command("clear");
|
|
command("variable pair_style delete");
|
|
command("variable data_file delete");
|
|
command("variable newton_pair delete");
|
|
command("variable newton_pair index on");
|
|
|
|
for (auto &pre_command : cfg.pre_commands) {
|
|
command(pre_command);
|
|
}
|
|
|
|
command("variable pair_style index '" + cfg.pair_style + "'");
|
|
command("variable data_file index " + cfg.basename + ".data");
|
|
|
|
std::string input_file = INPUT_FOLDER + PATH_SEP + cfg.input_file;
|
|
parse_input_script(input_file);
|
|
|
|
for (auto &pair_coeff : cfg.pair_coeff) {
|
|
command("pair_coeff " + pair_coeff);
|
|
}
|
|
for (auto &post_command : cfg.post_commands) {
|
|
command(post_command);
|
|
}
|
|
command("run 0 post no");
|
|
}
|
|
|
|
// re-generate yaml file with current settings.
|
|
|
|
void generate_yaml_file(const char *outfile, const TestConfig &config)
|
|
{
|
|
// initialize system geometry
|
|
const char *args[] = {"PairStyle", "-log", "none", "-echo", "screen", "-nocite"};
|
|
|
|
char **argv = (char **)args;
|
|
int argc = sizeof(args) / sizeof(char *);
|
|
LAMMPS *lmp = init_lammps(argc, argv, config);
|
|
if (!lmp) {
|
|
std::cerr << "One or more prerequisite styles are not available "
|
|
"in this LAMMPS configuration:\n";
|
|
for (auto prerequisite : config.prerequisites) {
|
|
std::cerr << prerequisite.first << "_style " << prerequisite.second << "\n";
|
|
}
|
|
return;
|
|
}
|
|
|
|
const int natoms = lmp->atom->natoms;
|
|
std::string block("");
|
|
|
|
YamlWriter writer(outfile);
|
|
|
|
// write yaml header
|
|
write_yaml_header(&writer, &test_config, lmp->version);
|
|
|
|
// pair_style
|
|
writer.emit("pair_style", config.pair_style);
|
|
|
|
// pair_coeff
|
|
block.clear();
|
|
for (auto pair_coeff : config.pair_coeff) {
|
|
block += pair_coeff + "\n";
|
|
}
|
|
writer.emit_block("pair_coeff", block);
|
|
|
|
// extract
|
|
block.clear();
|
|
for (auto data : config.extract)
|
|
block += fmt::format("{} {}\n", data.first, data.second);
|
|
writer.emit_block("extract", block);
|
|
|
|
// natoms
|
|
writer.emit("natoms", natoms);
|
|
|
|
// init_vdwl
|
|
writer.emit("init_vdwl", lmp->force->pair->eng_vdwl);
|
|
|
|
// init_coul
|
|
writer.emit("init_coul", lmp->force->pair->eng_coul);
|
|
|
|
// init_stress
|
|
auto stress = lmp->force->pair->virial;
|
|
// avoid false positives on tiny stresses. force to zero instead.
|
|
for (int i = 0; i < 6; ++i)
|
|
if (fabs(stress[i]) < 1.0e-13) stress[i] = 0.0;
|
|
block = fmt::format("{:23.16e} {:23.16e} {:23.16e} {:23.16e} {:23.16e} {:23.16e}", stress[0],
|
|
stress[1], stress[2], stress[3], stress[4], stress[5]);
|
|
writer.emit_block("init_stress", block);
|
|
|
|
// init_forces
|
|
block.clear();
|
|
auto f = lmp->atom->f;
|
|
for (int i = 1; i <= natoms; ++i) {
|
|
const int j = lmp->atom->map(i);
|
|
block += fmt::format("{:3} {:23.16e} {:23.16e} {:23.16e}\n", i, f[j][0], f[j][1], f[j][2]);
|
|
}
|
|
writer.emit_block("init_forces", block);
|
|
|
|
// do a few steps of MD
|
|
run_lammps(lmp);
|
|
|
|
// run_vdwl
|
|
writer.emit("run_vdwl", lmp->force->pair->eng_vdwl);
|
|
|
|
// run_coul
|
|
writer.emit("run_coul", lmp->force->pair->eng_coul);
|
|
|
|
// run_stress
|
|
stress = lmp->force->pair->virial;
|
|
// avoid false positives on tiny stresses. force to zero instead.
|
|
for (int i = 0; i < 6; ++i)
|
|
if (fabs(stress[i]) < 1.0e-13) stress[i] = 0.0;
|
|
block = fmt::format("{:23.16e} {:23.16e} {:23.16e} {:23.16e} {:23.16e} {:23.16e}", stress[0],
|
|
stress[1], stress[2], stress[3], stress[4], stress[5]);
|
|
writer.emit_block("run_stress", block);
|
|
|
|
block.clear();
|
|
f = lmp->atom->f;
|
|
for (int i = 1; i <= natoms; ++i) {
|
|
const int j = lmp->atom->map(i);
|
|
block += fmt::format("{:3} {:23.16e} {:23.16e} {:23.16e}\n", i, f[j][0], f[j][1], f[j][2]);
|
|
}
|
|
writer.emit_block("run_forces", block);
|
|
|
|
cleanup_lammps(lmp, config);
|
|
return;
|
|
}
|
|
|
|
TEST(PairStyle, plain)
|
|
{
|
|
if (test_config.skip_tests.count(test_info_->name())) GTEST_SKIP();
|
|
|
|
const char *args[] = {"PairStyle", "-log", "none", "-echo", "screen", "-nocite"};
|
|
|
|
char **argv = (char **)args;
|
|
int argc = sizeof(args) / sizeof(char *);
|
|
|
|
::testing::internal::CaptureStdout();
|
|
LAMMPS *lmp = init_lammps(argc, argv, test_config, true);
|
|
|
|
std::string output = ::testing::internal::GetCapturedStdout();
|
|
if (verbose) std::cout << output;
|
|
|
|
if (!lmp) {
|
|
std::cerr << "One or more prerequisite styles are not available "
|
|
"in this LAMMPS configuration:\n";
|
|
for (auto &prerequisite : test_config.prerequisites) {
|
|
std::cerr << prerequisite.first << "_style " << prerequisite.second << "\n";
|
|
}
|
|
GTEST_SKIP();
|
|
}
|
|
|
|
EXPECT_THAT(output, StartsWith("LAMMPS ("));
|
|
EXPECT_THAT(output, HasSubstr("Loop time"));
|
|
|
|
// abort if running in parallel and not all atoms are local
|
|
const int nlocal = lmp->atom->nlocal;
|
|
ASSERT_EQ(lmp->atom->natoms, nlocal);
|
|
|
|
double epsilon = test_config.epsilon;
|
|
// relax test precision when using pppm and single precision FFTs
|
|
#if defined(FFT_SINGLE)
|
|
if (lmp->force->kspace && lmp->force->kspace->compute_flag)
|
|
if (utils::strmatch(lmp->force->kspace_style, "^pppm")) epsilon *= 2.0e8;
|
|
#endif
|
|
auto f = lmp->atom->f;
|
|
auto tag = lmp->atom->tag;
|
|
ErrorStats stats;
|
|
stats.reset();
|
|
const std::vector<coord_t> &f_ref = test_config.init_forces;
|
|
ASSERT_EQ(nlocal + 1, f_ref.size());
|
|
for (int i = 0; i < nlocal; ++i) {
|
|
EXPECT_FP_LE_WITH_EPS(f[i][0], f_ref[tag[i]].x, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][1], f_ref[tag[i]].y, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][2], f_ref[tag[i]].z, epsilon);
|
|
}
|
|
if (print_stats) std::cerr << "init_forces stats, newton on: " << stats << std::endl;
|
|
|
|
auto pair = lmp->force->pair;
|
|
auto stress = pair->virial;
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(stress[0], test_config.init_stress.xx, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[1], test_config.init_stress.yy, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[2], test_config.init_stress.zz, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[3], test_config.init_stress.xy, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[4], test_config.init_stress.xz, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[5], test_config.init_stress.yz, epsilon);
|
|
if (print_stats) std::cerr << "init_stress stats, newton on: " << stats << std::endl;
|
|
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_vdwl, test_config.init_vdwl, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_coul, test_config.init_coul, epsilon);
|
|
if (print_stats) std::cerr << "init_energy stats, newton on: " << stats << std::endl;
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
run_lammps(lmp);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
f = lmp->atom->f;
|
|
tag = lmp->atom->tag;
|
|
stress = pair->virial;
|
|
|
|
const std::vector<coord_t> &f_run = test_config.run_forces;
|
|
ASSERT_EQ(nlocal + 1, f_run.size());
|
|
stats.reset();
|
|
for (int i = 0; i < nlocal; ++i) {
|
|
EXPECT_FP_LE_WITH_EPS(f[i][0], f_run[tag[i]].x, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][1], f_run[tag[i]].y, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][2], f_run[tag[i]].z, 5 * epsilon);
|
|
}
|
|
if (print_stats) std::cerr << "run_forces stats, newton on: " << stats << std::endl;
|
|
|
|
stress = pair->virial;
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(stress[0], test_config.run_stress.xx, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[1], test_config.run_stress.yy, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[2], test_config.run_stress.zz, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[3], test_config.run_stress.xy, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[4], test_config.run_stress.xz, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[5], test_config.run_stress.yz, epsilon);
|
|
if (print_stats) std::cerr << "run_stress stats, newton on: " << stats << std::endl;
|
|
|
|
stats.reset();
|
|
int id = lmp->modify->find_compute("sum");
|
|
double energy = lmp->modify->compute[id]->compute_scalar();
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_vdwl, test_config.run_vdwl, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_coul, test_config.run_coul, epsilon);
|
|
// skip comparing per-atom energy with total energy for "kim"
|
|
if (std::string("kim") != lmp->force->pair_style)
|
|
EXPECT_FP_LE_WITH_EPS((pair->eng_vdwl + pair->eng_coul), energy, epsilon);
|
|
if (print_stats) std::cerr << "run_energy stats, newton on: " << stats << std::endl;
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
cleanup_lammps(lmp, test_config);
|
|
lmp = init_lammps(argc, argv, test_config, false);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
// skip over these tests if newton pair is forced to be on
|
|
if (lmp->force->newton_pair == 0) {
|
|
|
|
f = lmp->atom->f;
|
|
tag = lmp->atom->tag;
|
|
stats.reset();
|
|
for (int i = 0; i < nlocal; ++i) {
|
|
EXPECT_FP_LE_WITH_EPS(f[i][0], f_ref[tag[i]].x, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][1], f_ref[tag[i]].y, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][2], f_ref[tag[i]].z, epsilon);
|
|
}
|
|
if (print_stats) std::cerr << "init_forces stats, newton off:" << stats << std::endl;
|
|
|
|
pair = lmp->force->pair;
|
|
stress = pair->virial;
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(stress[0], test_config.init_stress.xx, 3 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[1], test_config.init_stress.yy, 3 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[2], test_config.init_stress.zz, 3 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[3], test_config.init_stress.xy, 3 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[4], test_config.init_stress.xz, 3 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[5], test_config.init_stress.yz, 3 * epsilon);
|
|
if (print_stats) std::cerr << "init_stress stats, newton off:" << stats << std::endl;
|
|
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_vdwl, test_config.init_vdwl, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_coul, test_config.init_coul, epsilon);
|
|
if (print_stats) std::cerr << "init_energy stats, newton off:" << stats << std::endl;
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
run_lammps(lmp);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
f = lmp->atom->f;
|
|
tag = lmp->atom->tag;
|
|
stress = pair->virial;
|
|
stats.reset();
|
|
for (int i = 0; i < nlocal; ++i) {
|
|
EXPECT_FP_LE_WITH_EPS(f[i][0], f_run[tag[i]].x, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][1], f_run[tag[i]].y, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][2], f_run[tag[i]].z, 5 * epsilon);
|
|
}
|
|
if (print_stats) std::cerr << "run_forces stats, newton off:" << stats << std::endl;
|
|
|
|
stress = pair->virial;
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(stress[0], test_config.run_stress.xx, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[1], test_config.run_stress.yy, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[2], test_config.run_stress.zz, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[3], test_config.run_stress.xy, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[4], test_config.run_stress.xz, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[5], test_config.run_stress.yz, epsilon);
|
|
if (print_stats) std::cerr << "run_stress stats, newton off:" << stats << std::endl;
|
|
|
|
stats.reset();
|
|
id = lmp->modify->find_compute("sum");
|
|
energy = lmp->modify->compute[id]->compute_scalar();
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_vdwl, test_config.run_vdwl, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_coul, test_config.run_coul, epsilon);
|
|
// skip comparing per-atom energy with total energy for "kim"
|
|
if (std::string("kim") != lmp->force->pair_style)
|
|
EXPECT_FP_LE_WITH_EPS((pair->eng_vdwl + pair->eng_coul), energy, epsilon);
|
|
if (print_stats) std::cerr << "run_energy stats, newton off:" << stats << std::endl;
|
|
}
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
restart_lammps(lmp, test_config);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
f = lmp->atom->f;
|
|
tag = lmp->atom->tag;
|
|
stats.reset();
|
|
ASSERT_EQ(nlocal + 1, f_ref.size());
|
|
for (int i = 0; i < nlocal; ++i) {
|
|
EXPECT_FP_LE_WITH_EPS(f[i][0], f_ref[tag[i]].x, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][1], f_ref[tag[i]].y, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][2], f_ref[tag[i]].z, epsilon);
|
|
}
|
|
if (print_stats) std::cerr << "restart_forces stats:" << stats << std::endl;
|
|
|
|
pair = lmp->force->pair;
|
|
stress = pair->virial;
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(stress[0], test_config.init_stress.xx, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[1], test_config.init_stress.yy, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[2], test_config.init_stress.zz, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[3], test_config.init_stress.xy, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[4], test_config.init_stress.xz, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[5], test_config.init_stress.yz, epsilon);
|
|
if (print_stats) std::cerr << "restart_stress stats:" << stats << std::endl;
|
|
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_vdwl, test_config.init_vdwl, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_coul, test_config.init_coul, epsilon);
|
|
if (print_stats) std::cerr << "restart_energy stats:" << stats << std::endl;
|
|
|
|
// pair style rann does not support pair_modify nofdotr
|
|
if (test_config.pair_style != "rann") {
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
restart_lammps(lmp, test_config, true);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
f = lmp->atom->f;
|
|
tag = lmp->atom->tag;
|
|
stats.reset();
|
|
ASSERT_EQ(nlocal + 1, f_ref.size());
|
|
for (int i = 0; i < nlocal; ++i) {
|
|
EXPECT_FP_LE_WITH_EPS(f[i][0], f_ref[tag[i]].x, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][1], f_ref[tag[i]].y, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][2], f_ref[tag[i]].z, epsilon);
|
|
}
|
|
if (print_stats) std::cerr << "nofdotr_forces stats:" << stats << std::endl;
|
|
|
|
pair = lmp->force->pair;
|
|
stress = pair->virial;
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(stress[0], test_config.init_stress.xx, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[1], test_config.init_stress.yy, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[2], test_config.init_stress.zz, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[3], test_config.init_stress.xy, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[4], test_config.init_stress.xz, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[5], test_config.init_stress.yz, epsilon);
|
|
if (print_stats) std::cerr << "nofdotr_stress stats:" << stats << std::endl;
|
|
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_vdwl, test_config.init_vdwl, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_coul, test_config.init_coul, epsilon);
|
|
if (print_stats) std::cerr << "nofdotr_energy stats:" << stats << std::endl;
|
|
}
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
data_lammps(lmp, test_config);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
f = lmp->atom->f;
|
|
tag = lmp->atom->tag;
|
|
stats.reset();
|
|
ASSERT_EQ(nlocal + 1, f_ref.size());
|
|
for (int i = 0; i < nlocal; ++i) {
|
|
EXPECT_FP_LE_WITH_EPS(f[i][0], f_ref[tag[i]].x, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][1], f_ref[tag[i]].y, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][2], f_ref[tag[i]].z, epsilon);
|
|
}
|
|
if (print_stats) std::cerr << "data_forces stats:" << stats << std::endl;
|
|
|
|
pair = lmp->force->pair;
|
|
stress = pair->virial;
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(stress[0], test_config.init_stress.xx, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[1], test_config.init_stress.yy, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[2], test_config.init_stress.zz, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[3], test_config.init_stress.xy, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[4], test_config.init_stress.xz, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[5], test_config.init_stress.yz, epsilon);
|
|
if (print_stats) std::cerr << "data_stress stats:" << stats << std::endl;
|
|
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_vdwl, test_config.init_vdwl, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_coul, test_config.init_coul, epsilon);
|
|
if (print_stats) std::cerr << "data_energy stats:" << stats << std::endl;
|
|
|
|
if (pair->respa_enable) {
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
cleanup_lammps(lmp, test_config);
|
|
lmp = init_lammps(argc, argv, test_config, false);
|
|
lmp->input->one("run_style respa 2 1 inner 1 4.8 5.5 outer 2");
|
|
run_lammps(lmp);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
// need to relax error by a large amount with tabulation, since
|
|
// coul/long styles do not use tabulation in compute_inner()
|
|
// and compute_middle() so we get a significant deviation.
|
|
pair = lmp->force->pair;
|
|
if (pair->ncoultablebits) epsilon *= 5.0e6;
|
|
|
|
f = lmp->atom->f;
|
|
tag = lmp->atom->tag;
|
|
stats.reset();
|
|
for (int i = 0; i < nlocal; ++i) {
|
|
EXPECT_FP_LE_WITH_EPS(f[i][0], f_run[tag[i]].x, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][1], f_run[tag[i]].y, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][2], f_run[tag[i]].z, 5 * epsilon);
|
|
}
|
|
if (print_stats) std::cerr << "run_forces stats, r-RESPA:" << stats << std::endl;
|
|
|
|
stress = pair->virial;
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(stress[0], test_config.run_stress.xx, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[1], test_config.run_stress.yy, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[2], test_config.run_stress.zz, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[3], test_config.run_stress.xy, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[4], test_config.run_stress.xz, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[5], test_config.run_stress.yz, epsilon);
|
|
if (print_stats) std::cerr << "run_stress stats, r-RESPA:" << stats << std::endl;
|
|
|
|
stats.reset();
|
|
id = lmp->modify->find_compute("sum");
|
|
energy = lmp->modify->compute[id]->compute_scalar();
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_vdwl, test_config.run_vdwl, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_coul, test_config.run_coul, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS((pair->eng_vdwl + pair->eng_coul), energy, epsilon);
|
|
if (print_stats) std::cerr << "run_energy stats, r-RESPA:" << stats << std::endl;
|
|
}
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
cleanup_lammps(lmp, test_config);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
};
|
|
|
|
TEST(PairStyle, omp)
|
|
{
|
|
if (!LAMMPS::is_installed_pkg("OPENMP")) GTEST_SKIP();
|
|
if (test_config.skip_tests.count(test_info_->name())) GTEST_SKIP();
|
|
|
|
const char *args[] = {"PairStyle", "-log", "none", "-echo", "screen", "-nocite",
|
|
"-pk", "omp", "4", "-sf", "omp"};
|
|
|
|
// cannot run dpd styles with more than 1 thread due to using multiple pRNGs
|
|
if (utils::strmatch(test_config.pair_style, "^dpd")) args[8] = "1";
|
|
|
|
char **argv = (char **)args;
|
|
int argc = sizeof(args) / sizeof(char *);
|
|
|
|
::testing::internal::CaptureStdout();
|
|
LAMMPS *lmp = init_lammps(argc, argv, test_config, true);
|
|
|
|
std::string output = ::testing::internal::GetCapturedStdout();
|
|
if (verbose) std::cout << output;
|
|
|
|
if (!lmp) {
|
|
std::cerr << "One or more prerequisite styles with /omp suffix\n"
|
|
"are not available in this LAMMPS configuration:\n";
|
|
for (auto &prerequisite : test_config.prerequisites) {
|
|
std::cerr << prerequisite.first << "_style " << prerequisite.second << "\n";
|
|
}
|
|
GTEST_SKIP();
|
|
}
|
|
|
|
EXPECT_THAT(output, StartsWith("LAMMPS ("));
|
|
EXPECT_THAT(output, HasSubstr("Loop time"));
|
|
|
|
// abort if running in parallel and not all atoms are local
|
|
const int nlocal = lmp->atom->nlocal;
|
|
ASSERT_EQ(lmp->atom->natoms, nlocal);
|
|
|
|
// relax error a bit for OPENMP package
|
|
double epsilon = 5.0 * test_config.epsilon;
|
|
// relax test precision when using pppm and single precision FFTs
|
|
#if defined(FFT_SINGLE)
|
|
if (lmp->force->kspace && lmp->force->kspace->compute_flag)
|
|
if (utils::strmatch(lmp->force->kspace_style, "^pppm")) epsilon *= 2.0e8;
|
|
#endif
|
|
auto f = lmp->atom->f;
|
|
auto tag = lmp->atom->tag;
|
|
|
|
const std::vector<coord_t> &f_ref = test_config.init_forces;
|
|
ErrorStats stats;
|
|
stats.reset();
|
|
for (int i = 0; i < nlocal; ++i) {
|
|
EXPECT_FP_LE_WITH_EPS(f[i][0], f_ref[tag[i]].x, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][1], f_ref[tag[i]].y, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][2], f_ref[tag[i]].z, epsilon);
|
|
}
|
|
if (print_stats) std::cerr << "init_forces stats, newton on: " << stats << std::endl;
|
|
|
|
auto pair = lmp->force->pair;
|
|
auto stress = pair->virial;
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(stress[0], test_config.init_stress.xx, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[1], test_config.init_stress.yy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[2], test_config.init_stress.zz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[3], test_config.init_stress.xy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[4], test_config.init_stress.xz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[5], test_config.init_stress.yz, 10 * epsilon);
|
|
if (print_stats) std::cerr << "init_stress stats, newton on: " << stats << std::endl;
|
|
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_vdwl, test_config.init_vdwl, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_coul, test_config.init_coul, epsilon);
|
|
if (print_stats) std::cerr << "init_energy stats, newton on: " << stats << std::endl;
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
run_lammps(lmp);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
f = lmp->atom->f;
|
|
stress = pair->virial;
|
|
tag = lmp->atom->tag;
|
|
|
|
const std::vector<coord_t> &f_run = test_config.run_forces;
|
|
ASSERT_EQ(nlocal + 1, f_run.size());
|
|
stats.reset();
|
|
for (int i = 0; i < nlocal; ++i) {
|
|
EXPECT_FP_LE_WITH_EPS(f[i][0], f_run[tag[i]].x, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][1], f_run[tag[i]].y, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][2], f_run[tag[i]].z, 5 * epsilon);
|
|
}
|
|
if (print_stats) std::cerr << "run_forces stats, newton on: " << stats << std::endl;
|
|
|
|
stress = pair->virial;
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(stress[0], test_config.run_stress.xx, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[1], test_config.run_stress.yy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[2], test_config.run_stress.zz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[3], test_config.run_stress.xy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[4], test_config.run_stress.xz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[5], test_config.run_stress.yz, 10 * epsilon);
|
|
if (print_stats) std::cerr << "run_stress stats, newton on: " << stats << std::endl;
|
|
|
|
stats.reset();
|
|
int id = lmp->modify->find_compute("sum");
|
|
double energy = lmp->modify->compute[id]->compute_scalar();
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_vdwl, test_config.run_vdwl, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_coul, test_config.run_coul, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS((pair->eng_vdwl + pair->eng_coul), energy, epsilon);
|
|
if (print_stats) std::cerr << "run_energy stats, newton on: " << stats << std::endl;
|
|
|
|
// skip over these tests if newton pair is forced to be on
|
|
if (lmp->force->newton_pair == 0) {
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
cleanup_lammps(lmp, test_config);
|
|
lmp = init_lammps(argc, argv, test_config, false);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
f = lmp->atom->f;
|
|
tag = lmp->atom->tag;
|
|
stats.reset();
|
|
for (int i = 0; i < nlocal; ++i) {
|
|
EXPECT_FP_LE_WITH_EPS(f[i][0], f_ref[tag[i]].x, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][1], f_ref[tag[i]].y, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][2], f_ref[tag[i]].z, epsilon);
|
|
}
|
|
if (print_stats) std::cerr << "init_forces stats, newton off:" << stats << std::endl;
|
|
|
|
pair = lmp->force->pair;
|
|
stress = pair->virial;
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(stress[0], test_config.init_stress.xx, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[1], test_config.init_stress.yy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[2], test_config.init_stress.zz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[3], test_config.init_stress.xy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[4], test_config.init_stress.xz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[5], test_config.init_stress.yz, 10 * epsilon);
|
|
if (print_stats) std::cerr << "init_stress stats, newton off:" << stats << std::endl;
|
|
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_vdwl, test_config.init_vdwl, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_coul, test_config.init_coul, epsilon);
|
|
if (print_stats) std::cerr << "init_energy stats, newton off:" << stats << std::endl;
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
run_lammps(lmp);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
f = lmp->atom->f;
|
|
tag = lmp->atom->tag;
|
|
stats.reset();
|
|
for (int i = 0; i < nlocal; ++i) {
|
|
EXPECT_FP_LE_WITH_EPS(f[i][0], f_run[tag[i]].x, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][1], f_run[tag[i]].y, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][2], f_run[tag[i]].z, 5 * epsilon);
|
|
}
|
|
if (print_stats) std::cerr << "run_forces stats, newton off:" << stats << std::endl;
|
|
|
|
stress = pair->virial;
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(stress[0], test_config.run_stress.xx, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[1], test_config.run_stress.yy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[2], test_config.run_stress.zz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[3], test_config.run_stress.xy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[4], test_config.run_stress.xz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[5], test_config.run_stress.yz, 10 * epsilon);
|
|
if (print_stats) std::cerr << "run_stress stats, newton off:" << stats << std::endl;
|
|
|
|
stats.reset();
|
|
id = lmp->modify->find_compute("sum");
|
|
energy = lmp->modify->compute[id]->compute_scalar();
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_vdwl, test_config.run_vdwl, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_coul, test_config.run_coul, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS((pair->eng_vdwl + pair->eng_coul), energy, epsilon);
|
|
if (print_stats) std::cerr << "run_energy stats, newton off:" << stats << std::endl;
|
|
}
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
restart_lammps(lmp, test_config, true);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
f = lmp->atom->f;
|
|
tag = lmp->atom->tag;
|
|
stats.reset();
|
|
ASSERT_EQ(nlocal + 1, f_ref.size());
|
|
for (int i = 0; i < nlocal; ++i) {
|
|
EXPECT_FP_LE_WITH_EPS(f[i][0], f_ref[tag[i]].x, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][1], f_ref[tag[i]].y, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][2], f_ref[tag[i]].z, 5 * epsilon);
|
|
}
|
|
if (print_stats) std::cerr << "nofdotr_forces stats:" << stats << std::endl;
|
|
|
|
pair = lmp->force->pair;
|
|
stress = pair->virial;
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(stress[0], test_config.init_stress.xx, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[1], test_config.init_stress.yy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[2], test_config.init_stress.zz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[3], test_config.init_stress.xy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[4], test_config.init_stress.xz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[5], test_config.init_stress.yz, 10 * epsilon);
|
|
if (print_stats) std::cerr << "nofdotr_stress stats:" << stats << std::endl;
|
|
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_vdwl, test_config.init_vdwl, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_coul, test_config.init_coul, 5 * epsilon);
|
|
if (print_stats) std::cerr << "nofdotr_energy stats:" << stats << std::endl;
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
cleanup_lammps(lmp, test_config);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
};
|
|
|
|
TEST(PairStyle, gpu)
|
|
{
|
|
if (!LAMMPS::is_installed_pkg("GPU")) GTEST_SKIP();
|
|
if (!Info::has_gpu_device()) GTEST_SKIP();
|
|
if (test_config.skip_tests.count(test_info_->name())) GTEST_SKIP();
|
|
|
|
const char *args_neigh[] = {"PairStyle", "-log", "none", "-echo",
|
|
"screen", "-nocite", "-sf", "gpu"};
|
|
const char *args_noneigh[] = {"PairStyle", "-log", "none", "-echo", "screen", "-nocite", "-sf",
|
|
"gpu", "-pk", "gpu", "0", "neigh", "no"};
|
|
|
|
char **argv = (char **)args_neigh;
|
|
int argc = sizeof(args_neigh) / sizeof(char *);
|
|
|
|
// cannot use GPU neighbor list with hybrid pair style (yet)
|
|
if (test_config.pair_style.substr(0, 6) == "hybrid") {
|
|
argv = (char **)args_noneigh;
|
|
argc = sizeof(args_noneigh) / sizeof(char *);
|
|
}
|
|
|
|
::testing::internal::CaptureStdout();
|
|
LAMMPS *lmp = init_lammps(argc, argv, test_config, false);
|
|
|
|
std::string output = ::testing::internal::GetCapturedStdout();
|
|
if (verbose) std::cout << output;
|
|
|
|
if (!lmp) {
|
|
std::cerr << "One or more prerequisite styles with /gpu suffix\n"
|
|
"are not available in this LAMMPS configuration:\n";
|
|
for (auto &prerequisite : test_config.prerequisites) {
|
|
std::cerr << prerequisite.first << "_style " << prerequisite.second << "\n";
|
|
}
|
|
GTEST_SKIP();
|
|
}
|
|
|
|
EXPECT_THAT(output, StartsWith("LAMMPS ("));
|
|
EXPECT_THAT(output, HasSubstr("Loop time"));
|
|
|
|
// abort if running in parallel and not all atoms are local
|
|
const int nlocal = lmp->atom->nlocal;
|
|
ASSERT_EQ(lmp->atom->natoms, nlocal);
|
|
|
|
// relax error for GPU package depending on precision setting
|
|
double epsilon = test_config.epsilon;
|
|
if (Info::has_accelerator_feature("GPU", "precision", "double"))
|
|
epsilon *= 7.5;
|
|
else if (Info::has_accelerator_feature("GPU", "precision", "mixed"))
|
|
epsilon *= 5.0e8;
|
|
else
|
|
epsilon *= 1.0e10;
|
|
// relax test precision when using pppm and single precision FFTs, but only when also
|
|
// running with double precision
|
|
#if defined(FFT_SINGLE)
|
|
if (lmp->force->kspace && lmp->force->kspace->compute_flag &&
|
|
Info::has_accelerator_feature("GPU", "precision", "double"))
|
|
if (utils::strmatch(lmp->force->kspace_style, "^pppm")) epsilon *= 2.0e8;
|
|
#endif
|
|
const std::vector<coord_t> &f_ref = test_config.init_forces;
|
|
const std::vector<coord_t> &f_run = test_config.run_forces;
|
|
ErrorStats stats;
|
|
|
|
auto f = lmp->atom->f;
|
|
auto tag = lmp->atom->tag;
|
|
stats.reset();
|
|
for (int i = 0; i < nlocal; ++i) {
|
|
EXPECT_FP_LE_WITH_EPS(f[i][0], f_ref[tag[i]].x, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][1], f_ref[tag[i]].y, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][2], f_ref[tag[i]].z, epsilon);
|
|
}
|
|
if (print_stats) std::cerr << "init_forces stats, newton off:" << stats << std::endl;
|
|
|
|
auto pair = lmp->force->pair;
|
|
auto stress = pair->virial;
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(stress[0], test_config.init_stress.xx, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[1], test_config.init_stress.yy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[2], test_config.init_stress.zz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[3], test_config.init_stress.xy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[4], test_config.init_stress.xz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[5], test_config.init_stress.yz, 10 * epsilon);
|
|
if (print_stats) std::cerr << "init_stress stats, newton off:" << stats << std::endl;
|
|
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_vdwl, test_config.init_vdwl, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_coul, test_config.init_coul, epsilon);
|
|
if (print_stats) std::cerr << "init_energy stats, newton off:" << stats << std::endl;
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
run_lammps(lmp);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
f = lmp->atom->f;
|
|
tag = lmp->atom->tag;
|
|
stats.reset();
|
|
for (int i = 0; i < nlocal; ++i) {
|
|
EXPECT_FP_LE_WITH_EPS(f[i][0], f_run[tag[i]].x, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][1], f_run[tag[i]].y, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][2], f_run[tag[i]].z, 5 * epsilon);
|
|
}
|
|
if (print_stats) std::cerr << "run_forces stats, newton off:" << stats << std::endl;
|
|
|
|
stress = pair->virial;
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(stress[0], test_config.run_stress.xx, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[1], test_config.run_stress.yy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[2], test_config.run_stress.zz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[3], test_config.run_stress.xy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[4], test_config.run_stress.xz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[5], test_config.run_stress.yz, 10 * epsilon);
|
|
if (print_stats) std::cerr << "run_stress stats, newton off:" << stats << std::endl;
|
|
|
|
stats.reset();
|
|
auto id = lmp->modify->find_compute("sum");
|
|
auto energy = lmp->modify->compute[id]->compute_scalar();
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_vdwl, test_config.run_vdwl, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_coul, test_config.run_coul, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS((pair->eng_vdwl + pair->eng_coul), energy, epsilon);
|
|
if (print_stats) std::cerr << "run_energy stats, newton off:" << stats << std::endl;
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
cleanup_lammps(lmp, test_config);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
};
|
|
|
|
TEST(PairStyle, intel)
|
|
{
|
|
if (!LAMMPS::is_installed_pkg("INTEL")) GTEST_SKIP();
|
|
if (test_config.skip_tests.count(test_info_->name())) GTEST_SKIP();
|
|
|
|
const char *args[] = {"PairStyle", "-log", "none", "-echo", "screen", "-nocite",
|
|
"-pk", "intel", "0", "mode", "double", "omp",
|
|
"4", "lrt", "no", "-sf", "intel"};
|
|
|
|
// cannot use more than 1 thread for dpd styles due to pRNG
|
|
if (utils::strmatch(test_config.pair_style, "^dpd")) args[12] = "1";
|
|
|
|
char **argv = (char **)args;
|
|
int argc = sizeof(args) / sizeof(char *);
|
|
|
|
::testing::internal::CaptureStdout();
|
|
LAMMPS *lmp = init_lammps(argc, argv, test_config, true);
|
|
|
|
std::string output = ::testing::internal::GetCapturedStdout();
|
|
if (verbose) std::cout << output;
|
|
|
|
if (!lmp) {
|
|
std::cerr << "One or more prerequisite styles with /intel suffix\n"
|
|
"are not available in this LAMMPS configuration:\n";
|
|
for (auto prerequisite : test_config.prerequisites) {
|
|
std::cerr << prerequisite.first << "_style " << prerequisite.second << "\n";
|
|
}
|
|
GTEST_SKIP();
|
|
}
|
|
|
|
// relax error a bit for INTEL package
|
|
double epsilon = 7.5 * test_config.epsilon;
|
|
// relax test precision when using pppm and single precision FFTs
|
|
#if defined(FFT_SINGLE)
|
|
if (lmp->force->kspace && lmp->force->kspace->compute_flag)
|
|
if (utils::strmatch(lmp->force->kspace_style, "^pppm")) epsilon *= 2.0e8;
|
|
#endif
|
|
|
|
// we need to relax the epsilon a LOT for tests using long-range
|
|
// coulomb with tabulation. seems more like mixed precision or a bug
|
|
for (auto post_cmd : test_config.post_commands) {
|
|
if (post_cmd.find("pair_modify table") != std::string::npos) {
|
|
if (post_cmd.find("pair_modify table 0") == std::string::npos) epsilon *= 1000000.0;
|
|
}
|
|
}
|
|
|
|
EXPECT_THAT(output, StartsWith("LAMMPS ("));
|
|
EXPECT_THAT(output, HasSubstr("Loop time"));
|
|
|
|
// abort if running in parallel and not all atoms are local
|
|
const int nlocal = lmp->atom->nlocal;
|
|
ASSERT_EQ(lmp->atom->natoms, nlocal);
|
|
|
|
auto f = lmp->atom->f;
|
|
auto tag = lmp->atom->tag;
|
|
|
|
const std::vector<coord_t> &f_ref = test_config.init_forces;
|
|
ErrorStats stats;
|
|
stats.reset();
|
|
for (int i = 0; i < nlocal; ++i) {
|
|
EXPECT_FP_LE_WITH_EPS(f[i][0], f_ref[tag[i]].x, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][1], f_ref[tag[i]].y, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][2], f_ref[tag[i]].z, epsilon);
|
|
}
|
|
if (print_stats) std::cerr << "init_forces stats:" << stats << std::endl;
|
|
|
|
Pair *pair = lmp->force->pair;
|
|
double *stress = pair->virial;
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(stress[0], test_config.init_stress.xx, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[1], test_config.init_stress.yy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[2], test_config.init_stress.zz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[3], test_config.init_stress.xy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[4], test_config.init_stress.xz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[5], test_config.init_stress.yz, 10 * epsilon);
|
|
if (print_stats) std::cerr << "init_stress stats:" << stats << std::endl;
|
|
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_vdwl, test_config.init_vdwl, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_coul, test_config.init_coul, epsilon);
|
|
if (print_stats) std::cerr << "init_energy stats:" << stats << std::endl;
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
run_lammps(lmp);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
f = lmp->atom->f;
|
|
tag = lmp->atom->tag;
|
|
stress = pair->virial;
|
|
|
|
const std::vector<coord_t> &f_run = test_config.run_forces;
|
|
ASSERT_EQ(nlocal + 1, f_run.size());
|
|
stats.reset();
|
|
for (int i = 0; i < nlocal; ++i) {
|
|
EXPECT_FP_LE_WITH_EPS(f[i][0], f_run[tag[i]].x, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][1], f_run[tag[i]].y, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][2], f_run[tag[i]].z, 5 * epsilon);
|
|
}
|
|
if (print_stats) std::cerr << "run_forces stats:" << stats << std::endl;
|
|
|
|
stress = pair->virial;
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(stress[0], test_config.run_stress.xx, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[1], test_config.run_stress.yy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[2], test_config.run_stress.zz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[3], test_config.run_stress.xy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[4], test_config.run_stress.xz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[5], test_config.run_stress.yz, 10 * epsilon);
|
|
if (print_stats) std::cerr << "run_stress stats:" << stats << std::endl;
|
|
|
|
stats.reset();
|
|
int id = lmp->modify->find_compute("sum");
|
|
double energy = lmp->modify->compute[id]->compute_scalar();
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_vdwl, test_config.run_vdwl, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_coul, test_config.run_coul, epsilon);
|
|
|
|
// rebo family of pair styles will have a large error in per-atom energy for INTEL
|
|
if (test_config.pair_style.find("rebo") != std::string::npos) epsilon *= 100000.0;
|
|
|
|
EXPECT_FP_LE_WITH_EPS((pair->eng_vdwl + pair->eng_coul), energy, epsilon);
|
|
if (print_stats) std::cerr << "run_energy stats:" << stats << std::endl;
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
cleanup_lammps(lmp, test_config);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
};
|
|
|
|
TEST(PairStyle, opt)
|
|
{
|
|
if (!LAMMPS::is_installed_pkg("OPT")) GTEST_SKIP();
|
|
if (test_config.skip_tests.count(test_info_->name())) GTEST_SKIP();
|
|
|
|
const char *args[] = {"PairStyle", "-log", "none", "-echo", "screen", "-nocite", "-sf", "opt"};
|
|
|
|
char **argv = (char **)args;
|
|
int argc = sizeof(args) / sizeof(char *);
|
|
|
|
::testing::internal::CaptureStdout();
|
|
LAMMPS *lmp = init_lammps(argc, argv, test_config);
|
|
|
|
std::string output = ::testing::internal::GetCapturedStdout();
|
|
if (verbose) std::cout << output;
|
|
|
|
if (!lmp) {
|
|
std::cerr << "One or more prerequisite styles with /opt suffix\n"
|
|
"are not available in this LAMMPS configuration:\n";
|
|
for (auto prerequisite : test_config.prerequisites) {
|
|
std::cerr << prerequisite.first << "_style " << prerequisite.second << "\n";
|
|
}
|
|
GTEST_SKIP();
|
|
}
|
|
|
|
EXPECT_THAT(output, StartsWith("LAMMPS ("));
|
|
EXPECT_THAT(output, HasSubstr("Loop time"));
|
|
|
|
// abort if running in parallel and not all atoms are local
|
|
const int nlocal = lmp->atom->nlocal;
|
|
ASSERT_EQ(lmp->atom->natoms, nlocal);
|
|
|
|
// relax error a bit for OPT package
|
|
double epsilon = 2.0 * test_config.epsilon;
|
|
// relax test precision when using pppm and single precision FFTs
|
|
#if defined(FFT_SINGLE)
|
|
if (lmp->force->kspace && lmp->force->kspace->compute_flag)
|
|
if (utils::strmatch(lmp->force->kspace_style, "^pppm")) epsilon *= 2.0e8;
|
|
#endif
|
|
auto f = lmp->atom->f;
|
|
auto tag = lmp->atom->tag;
|
|
|
|
const std::vector<coord_t> &f_ref = test_config.init_forces;
|
|
ErrorStats stats;
|
|
stats.reset();
|
|
for (int i = 0; i < nlocal; ++i) {
|
|
EXPECT_FP_LE_WITH_EPS(f[i][0], f_ref[tag[i]].x, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][1], f_ref[tag[i]].y, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][2], f_ref[tag[i]].z, epsilon);
|
|
}
|
|
if (print_stats) std::cerr << "init_forces stats:" << stats << std::endl;
|
|
|
|
Pair *pair = lmp->force->pair;
|
|
double *stress = pair->virial;
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(stress[0], test_config.init_stress.xx, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[1], test_config.init_stress.yy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[2], test_config.init_stress.zz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[3], test_config.init_stress.xy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[4], test_config.init_stress.xz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[5], test_config.init_stress.yz, 10 * epsilon);
|
|
if (print_stats) std::cerr << "init_stress stats:" << stats << std::endl;
|
|
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_vdwl, test_config.init_vdwl, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_coul, test_config.init_coul, epsilon);
|
|
if (print_stats) std::cerr << "init_energy stats:" << stats << std::endl;
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
run_lammps(lmp);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
f = lmp->atom->f;
|
|
tag = lmp->atom->tag;
|
|
stress = pair->virial;
|
|
const std::vector<coord_t> &f_run = test_config.run_forces;
|
|
ASSERT_EQ(nlocal + 1, f_run.size());
|
|
stats.reset();
|
|
for (int i = 0; i < nlocal; ++i) {
|
|
EXPECT_FP_LE_WITH_EPS(f[i][0], f_run[tag[i]].x, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][1], f_run[tag[i]].y, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][2], f_run[tag[i]].z, 5 * epsilon);
|
|
}
|
|
if (print_stats) std::cerr << "run_forces stats:" << stats << std::endl;
|
|
|
|
stress = pair->virial;
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(stress[0], test_config.run_stress.xx, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[1], test_config.run_stress.yy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[2], test_config.run_stress.zz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[3], test_config.run_stress.xy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[4], test_config.run_stress.xz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[5], test_config.run_stress.yz, 10 * epsilon);
|
|
if (print_stats) std::cerr << "run_stress stats:" << stats << std::endl;
|
|
|
|
stats.reset();
|
|
int id = lmp->modify->find_compute("sum");
|
|
double energy = lmp->modify->compute[id]->compute_scalar();
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_vdwl, test_config.run_vdwl, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_coul, test_config.run_coul, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS((pair->eng_vdwl + pair->eng_coul), energy, epsilon);
|
|
if (print_stats) std::cerr << "run_energy stats:" << stats << std::endl;
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
restart_lammps(lmp, test_config, true);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
f = lmp->atom->f;
|
|
tag = lmp->atom->tag;
|
|
stats.reset();
|
|
ASSERT_EQ(nlocal + 1, f_ref.size());
|
|
for (int i = 0; i < nlocal; ++i) {
|
|
EXPECT_FP_LE_WITH_EPS(f[i][0], f_ref[tag[i]].x, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][1], f_ref[tag[i]].y, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[i][2], f_ref[tag[i]].z, 5 * epsilon);
|
|
}
|
|
if (print_stats) std::cerr << "nofdotr_forces stats:" << stats << std::endl;
|
|
|
|
pair = lmp->force->pair;
|
|
stress = pair->virial;
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(stress[0], test_config.init_stress.xx, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[1], test_config.init_stress.yy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[2], test_config.init_stress.zz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[3], test_config.init_stress.xy, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[4], test_config.init_stress.xz, 10 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(stress[5], test_config.init_stress.yz, 10 * epsilon);
|
|
if (print_stats) std::cerr << "nofdotr_stress stats:" << stats << std::endl;
|
|
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_vdwl, test_config.init_vdwl, 5 * epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(pair->eng_coul, test_config.init_coul, 5 * epsilon);
|
|
if (print_stats) std::cerr << "nofdotr_energy stats:" << stats << std::endl;
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
cleanup_lammps(lmp, test_config);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
};
|
|
|
|
TEST(PairStyle, single)
|
|
{
|
|
if (test_config.skip_tests.count(test_info_->name())) GTEST_SKIP();
|
|
|
|
const char *args[] = {"PairStyle", "-log", "none", "-echo", "screen", "-nocite"};
|
|
|
|
char **argv = (char **)args;
|
|
int argc = sizeof(args) / sizeof(char *);
|
|
|
|
// need to add this dependency
|
|
test_config.prerequisites.push_back(std::make_pair("atom", "full"));
|
|
|
|
// create a LAMMPS instance with standard settings to detect the number of atom types
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
LAMMPS *lmp = init_lammps(argc, argv, test_config);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
if (!lmp) {
|
|
std::cerr << "One or more prerequisite styles are not available "
|
|
"in this LAMMPS configuration:\n";
|
|
for (auto &prerequisite : test_config.prerequisites) {
|
|
std::cerr << prerequisite.first << "_style " << prerequisite.second << "\n";
|
|
}
|
|
test_config.prerequisites.pop_back();
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
cleanup_lammps(lmp, test_config);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
GTEST_SKIP();
|
|
}
|
|
test_config.prerequisites.pop_back();
|
|
|
|
// gather some information and skip if unsupported
|
|
int ntypes = lmp->atom->ntypes;
|
|
int molecular = lmp->atom->molecular;
|
|
if (molecular > Atom::MOLECULAR) {
|
|
std::cerr << "Only atomic and simple molecular atom styles are supported\n";
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
cleanup_lammps(lmp, test_config);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
GTEST_SKIP();
|
|
}
|
|
|
|
Pair *pair = lmp->force->pair;
|
|
if (!pair->single_enable) {
|
|
std::cerr << "Single method not available for pair style " << test_config.pair_style
|
|
<< std::endl;
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
cleanup_lammps(lmp, test_config);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
GTEST_SKIP();
|
|
}
|
|
|
|
if (!pair->compute_flag) {
|
|
std::cerr << "Pair style disabled" << std::endl;
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
cleanup_lammps(lmp, test_config);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
GTEST_SKIP();
|
|
}
|
|
|
|
// now start over
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
|
|
// utility lambda to improve readability
|
|
auto command = [&](const std::string &line) {
|
|
lmp->input->one(line.c_str());
|
|
};
|
|
|
|
command("clear");
|
|
command("variable newton_pair delete");
|
|
command("variable newton_pair index on");
|
|
|
|
command("variable input_dir index " + INPUT_FOLDER);
|
|
|
|
for (auto &pre_command : test_config.pre_commands) {
|
|
command(pre_command);
|
|
}
|
|
|
|
command("atom_style full");
|
|
command("units ${units}");
|
|
command("boundary p p p");
|
|
command("newton ${newton_pair} ${newton_bond}");
|
|
|
|
if (molecular == Atom::MOLECULAR) {
|
|
command("special_bonds lj/coul "
|
|
"${bond_factor} ${angle_factor} ${dihedral_factor}");
|
|
}
|
|
|
|
command("atom_modify map array");
|
|
command("region box block -10.0 10.0 -10.0 10.0 -10.0 10.0 units box");
|
|
|
|
auto cmd = fmt::format("create_box {} box", ntypes);
|
|
if (molecular == Atom::MOLECULAR) {
|
|
cmd += " bond/types 1"
|
|
" extra/bond/per/atom 1"
|
|
" extra/special/per/atom 1";
|
|
}
|
|
command(cmd);
|
|
|
|
command("pair_style " + test_config.pair_style);
|
|
|
|
pair = lmp->force->pair;
|
|
|
|
for (auto &pair_coeff : test_config.pair_coeff) {
|
|
command("pair_coeff " + pair_coeff);
|
|
}
|
|
|
|
// create (only) two atoms
|
|
|
|
command("mass * 1.0");
|
|
command("create_atoms 1 single 0.0 -0.75 0.4 units box");
|
|
command("create_atoms 2 single 1.5 0.25 -0.1 units box");
|
|
command("set atom 1 charge -0.5");
|
|
command("set atom 2 charge 0.5");
|
|
command("set atom 1 mol 1");
|
|
command("set atom 2 mol 2");
|
|
command("special_bonds lj/coul 1.0 1.0 1.0");
|
|
|
|
if (molecular == Atom::MOLECULAR) {
|
|
command("create_bonds single/bond 1 1 2");
|
|
command("bond_style zero");
|
|
command("bond_coeff 1 2.0");
|
|
}
|
|
|
|
for (auto &post_command : test_config.post_commands) {
|
|
command(post_command);
|
|
}
|
|
|
|
command("run 0 post no");
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
int idx1 = lmp->atom->map(1);
|
|
int idx2 = lmp->atom->map(2);
|
|
double epsilon = test_config.epsilon;
|
|
double **f = lmp->atom->f;
|
|
double **x = lmp->atom->x;
|
|
double delx = x[idx2][0] - x[idx1][0];
|
|
double dely = x[idx2][1] - x[idx1][1];
|
|
double delz = x[idx2][2] - x[idx1][2];
|
|
double rsq = delx * delx + dely * dely + delz * delz;
|
|
double fsingle = 0.0;
|
|
double epair[4], esngl[4];
|
|
double splj = lmp->force->special_lj[1];
|
|
double spcl = lmp->force->special_coul[1];
|
|
ErrorStats stats;
|
|
|
|
epair[0] = pair->eng_vdwl + pair->eng_coul;
|
|
esngl[0] = pair->single(idx1, idx2, 1, 2, rsq, spcl, splj, fsingle);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx1][0], -fsingle * delx, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx1][1], -fsingle * dely, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx1][2], -fsingle * delz, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx2][0], fsingle * delx, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx2][1], fsingle * dely, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx2][2], fsingle * delz, epsilon);
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
command("displace_atoms all random 0.5 0.5 0.5 723456");
|
|
command("run 0 post no");
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
f = lmp->atom->f;
|
|
x = lmp->atom->x;
|
|
idx1 = lmp->atom->map(1);
|
|
idx2 = lmp->atom->map(2);
|
|
delx = x[idx2][0] - x[idx1][0];
|
|
dely = x[idx2][1] - x[idx1][1];
|
|
delz = x[idx2][2] - x[idx1][2];
|
|
rsq = delx * delx + dely * dely + delz * delz;
|
|
fsingle = 0.0;
|
|
|
|
epair[1] = pair->eng_vdwl + pair->eng_coul;
|
|
esngl[1] = pair->single(idx1, idx2, 1, 2, rsq, spcl, splj, fsingle);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx1][0], -fsingle * delx, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx1][1], -fsingle * dely, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx1][2], -fsingle * delz, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx2][0], fsingle * delx, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx2][1], fsingle * dely, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx2][2], fsingle * delz, epsilon);
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
command("displace_atoms all random 0.5 0.5 0.5 3456963");
|
|
command("run 0 post no");
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
f = lmp->atom->f;
|
|
x = lmp->atom->x;
|
|
idx1 = lmp->atom->map(1);
|
|
idx2 = lmp->atom->map(2);
|
|
delx = x[idx2][0] - x[idx1][0];
|
|
dely = x[idx2][1] - x[idx1][1];
|
|
delz = x[idx2][2] - x[idx1][2];
|
|
rsq = delx * delx + dely * dely + delz * delz;
|
|
fsingle = 0.0;
|
|
|
|
epair[2] = pair->eng_vdwl + pair->eng_coul;
|
|
esngl[2] = pair->single(idx1, idx2, 1, 2, rsq, spcl, splj, fsingle);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx1][0], -fsingle * delx, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx1][1], -fsingle * dely, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx1][2], -fsingle * delz, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx2][0], fsingle * delx, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx2][1], fsingle * dely, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx2][2], fsingle * delz, epsilon);
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
command("displace_atoms all random 0.5 0.5 0.5 9726532");
|
|
command("run 0 post no");
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
f = lmp->atom->f;
|
|
x = lmp->atom->x;
|
|
idx1 = lmp->atom->map(1);
|
|
idx2 = lmp->atom->map(2);
|
|
delx = x[idx2][0] - x[idx1][0];
|
|
dely = x[idx2][1] - x[idx1][1];
|
|
delz = x[idx2][2] - x[idx1][2];
|
|
rsq = delx * delx + dely * dely + delz * delz;
|
|
fsingle = 0.0;
|
|
|
|
epair[3] = pair->eng_vdwl + pair->eng_coul;
|
|
esngl[3] = pair->single(idx1, idx2, 1, 2, rsq, spcl, splj, fsingle);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx1][0], -fsingle * delx, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx1][1], -fsingle * dely, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx1][2], -fsingle * delz, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx2][0], fsingle * delx, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx2][1], fsingle * dely, epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(f[idx2][2], fsingle * delz, epsilon);
|
|
if (print_stats) std::cerr << "single_force stats:" << stats << std::endl;
|
|
|
|
if ((test_config.pair_style.find("coul/dsf") != std::string::npos) &&
|
|
(test_config.pair_style.find("coul/wolf") != std::string::npos)) {
|
|
stats.reset();
|
|
EXPECT_FP_LE_WITH_EPS(epair[0], esngl[0], epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(epair[1], esngl[1], epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(epair[2], esngl[2], epsilon);
|
|
EXPECT_FP_LE_WITH_EPS(epair[3], esngl[3], epsilon);
|
|
if (print_stats) std::cerr << "single_energy stats:" << stats << std::endl;
|
|
} else if (print_stats)
|
|
std::cerr << "skipping single_energy test due to self energy\n";
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
cleanup_lammps(lmp, test_config);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
}
|
|
|
|
TEST(PairStyle, extract)
|
|
{
|
|
if (test_config.skip_tests.count(test_info_->name())) GTEST_SKIP();
|
|
|
|
const char *args[] = {"PairStyle", "-log", "none", "-echo", "screen", "-nocite"};
|
|
|
|
char **argv = (char **)args;
|
|
int argc = sizeof(args) / sizeof(char *);
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
LAMMPS *lmp = init_lammps(argc, argv, test_config, true);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
if (!lmp) {
|
|
std::cerr << "One or more prerequisite styles are not available "
|
|
"in this LAMMPS configuration:\n";
|
|
for (const auto &prerequisite : test_config.prerequisites) {
|
|
std::cerr << prerequisite.first << "_style " << prerequisite.second << "\n";
|
|
}
|
|
GTEST_SKIP();
|
|
}
|
|
|
|
auto pair = lmp->force->pair;
|
|
if (!pair->compute_flag) {
|
|
std::cerr << "Pair style disabled" << std::endl;
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
cleanup_lammps(lmp, test_config);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
GTEST_SKIP();
|
|
}
|
|
|
|
void *ptr = nullptr;
|
|
int dim = 0;
|
|
for (auto &extract : test_config.extract) {
|
|
ptr = pair->extract(extract.first.c_str(), dim);
|
|
EXPECT_NE(ptr, nullptr);
|
|
EXPECT_EQ(dim, extract.second);
|
|
}
|
|
ptr = pair->extract("does_not_exist", dim);
|
|
EXPECT_EQ(ptr, nullptr);
|
|
|
|
// replace pair style with the same.
|
|
// should just update setting, but not create new style.
|
|
|
|
int ntypes = lmp->atom->ntypes;
|
|
for (int i = 1; i <= ntypes; ++i) {
|
|
for (int j = 1; j <= ntypes; ++j) {
|
|
pair->cutsq[i][j] = -1.0;
|
|
}
|
|
}
|
|
|
|
// utility lambda to improve readability
|
|
auto command = [&](const std::string &line) {
|
|
lmp->input->one(line.c_str());
|
|
};
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
command("pair_style " + test_config.pair_style);
|
|
EXPECT_EQ(pair, lmp->force->pair);
|
|
|
|
for (auto &pair_coeff : test_config.pair_coeff) {
|
|
command("pair_coeff " + pair_coeff);
|
|
}
|
|
pair->init();
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
|
|
for (int i = 1; i <= ntypes; ++i) {
|
|
for (int j = 1; j <= ntypes; ++j) {
|
|
EXPECT_GE(pair->cutsq[i][j], 0.0);
|
|
}
|
|
}
|
|
|
|
if (!verbose) ::testing::internal::CaptureStdout();
|
|
cleanup_lammps(lmp, test_config);
|
|
if (!verbose) ::testing::internal::GetCapturedStdout();
|
|
}
|