Function checkAbsoluteValueGroup

test/cpp/test_SymmetryFunctionGroups.cpp:203–232

}

void checkAbsoluteValueGroup(Structure&   s,
                             ElementMap&  em,
                             size_t const type,
                             string const setupLine,
                             double const value)
{
    SymFnc* sf = setupSymmetryFunction(em, type, setupLine);

    SymGrp* sfg = setupSymmetryFunctionGroup(em, *sf);

    // Allocate symmetry function arrays.
    s.atoms.at(0).numSymmetryFunctions = 1;
    s.atoms.at(0).numSymmetryFunctionDerivatives = vector<size_t>(em.size(), 1);
#ifndef NOSFCACHE
    s.atoms.at(0).cacheSizePerElement = vector<size_t>(maxElements,
                                                       maxCacheSize);
#endif
    s.atoms.at(0).allocate(true);

    // Calculate symmetry function for atom 0.
    recalculateSymmetryFunctionGroup(s, s.atoms.at(0), *sfg);

    BOOST_TEST_INFO(string("Symmetry function values, type ")
                    << type << "\n");
    BOOST_REQUIRE_SMALL(s.atoms.at(0).G.at(0) - value, accuracy);

    delete sf;

    return;
}

BOOST_AUTO_TEST_SUITE(IntegrationTests)