lammps/unittest/python/test_python_package.cpp

/* ----------------------------------------------------------------------
   LAMMPS - Large-scale Atomic/Molecular Massively Parallel Simulator
   https://lammps.sandia.gov/, 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.
------------------------------------------------------------------------- */

#include "atom.h"
#include "info.h"
#include "input.h"
#include "gmock/gmock.h"
#include "gtest/gtest.h"

#include <cmath>
#include <cstring>
#include <vector>

// location of '*.py' files required by tests
#define STRINGIFY(val) XSTR(val)
#define XSTR(val) #val
std::string INPUT_FOLDER = STRINGIFY(TEST_INPUT_FOLDER);

// whether to print verbose output (i.e. not capturing LAMMPS screen output).
bool verbose = false;

using LAMMPS_NS::utils::split_words;

namespace LAMMPS_NS {
using ::testing::MatchesRegex;
using ::testing::StrEq;
using ::testing::HasSubstr;

class PythonPackageTest : public ::testing::Test {
protected:
    LAMMPS *lmp;
    Info *info;

    void command(const std::string &line) { lmp->input->one(line.c_str()); }

    void SetUp() override
    {
        const char *args[] = {"PythonPackageTest", "-log", "none", "-echo", "screen", "-nocite"};
        char **argv        = (char **)args;
        int argc           = sizeof(args) / sizeof(char *);
        if (!verbose) ::testing::internal::CaptureStdout();
        lmp = new LAMMPS(argc, argv, MPI_COMM_WORLD);
        if (!verbose) ::testing::internal::GetCapturedStdout();
        ASSERT_NE(lmp, nullptr);
        info = new Info(lmp);
        if (!verbose) ::testing::internal::CaptureStdout();
        command("units real");
        command("dimension 3");
        command("region box block -4 4 -4 4 -4 4");
        command("create_box 1 box");
        command("create_atoms 1 single  0.0  0.0 0.0    units box");
        command("create_atoms 1 single  1.9 -1.9 1.9999 units box");
        command("pair_style zero 2.0");
        command("pair_coeff * *");
        command("mass * 1.0");
        command("variable input_dir index " + INPUT_FOLDER);
        if (!verbose) ::testing::internal::GetCapturedStdout();
    }

    void TearDown() override
    {
        if (!verbose) ::testing::internal::CaptureStdout();
        delete info;
        delete lmp;
        if (!verbose) ::testing::internal::GetCapturedStdout();
    }
};

TEST_F(PythonPackageTest, python_invoke)
{
    if (!info->has_package("PYTHON")) GTEST_SKIP();
    // execute python function from file
    if (!verbose) ::testing::internal::CaptureStdout();
    command("python printnum file ${input_dir}/func.py");
    if (!verbose) ::testing::internal::GetCapturedStdout();
    ::testing::internal::CaptureStdout();
    command("python printnum invoke");
    std::string output = ::testing::internal::GetCapturedStdout();
    if (verbose) std::cout << output;
    ASSERT_THAT(output, HasSubstr("2.25\n"));

    // execute another python function from same file
    if (!verbose) ::testing::internal::CaptureStdout();
    command("python printtxt exists");
    if (!verbose) ::testing::internal::GetCapturedStdout();
    ::testing::internal::CaptureStdout();
    command("python printtxt invoke");
    output = ::testing::internal::GetCapturedStdout();
    if (verbose) std::cout << output;
    ASSERT_THAT(output, HasSubstr("sometext\n"));

    // execute python function that uses the LAMMPS python module
    if (!verbose) ::testing::internal::CaptureStdout();
    command("variable idx equal 2.25");
    command("python getidxvar input 1 SELF format p exists");
    if (!verbose) ::testing::internal::GetCapturedStdout();
    ::testing::internal::CaptureStdout();
    command("python getidxvar invoke");
    output = ::testing::internal::GetCapturedStdout();
    if (verbose) std::cout << output;
    ASSERT_THAT(output, HasSubstr("2.25\n"));
}

TEST_F(PythonPackageTest, python_variable)
{
    if (!info->has_package("PYTHON")) GTEST_SKIP();
    if (!verbose) ::testing::internal::CaptureStdout();
    command("variable sq python square");
    command("variable val index 1.5");
    command("python square input 1 v_val return v_sq format ff file ${input_dir}/func.py");
    if (!verbose) ::testing::internal::GetCapturedStdout();
    ::testing::internal::CaptureStdout();
    command("print \"${sq}\"");
    std::string output = ::testing::internal::GetCapturedStdout();
    if (verbose) std::cout << output;
    ASSERT_THAT(output, MatchesRegex("print.*2.25.*"));
}

} // namespace LAMMPS_NS

int main(int argc, char **argv)
{
    MPI_Init(&argc, &argv);
    ::testing::InitGoogleMock(&argc, argv);

    // handle arguments passed via environment variable
    if (const char *var = getenv("TEST_ARGS")) {
        std::vector<std::string> env = split_words(var);
        for (auto arg : env) {
            if (arg == "-v") {
                verbose = true;
            }
        }
    }
    if ((argc > 1) && (strcmp(argv[1], "-v") == 0)) verbose = true;

    int rv = RUN_ALL_TESTS();
    MPI_Finalize();
    return rv;
}