24 void reset(std::vector<Atom>
const& atoms,
size_t const valuesPerPair);
34 double getf(
size_t const atomIndex,
35 size_t const neighIndex,
36 size_t const valIndex,
41 std::vector<std::vector<double>>
f;
Helper class to store previously calculated values of erfc() that are needed during the charge equili...
std::vector< std::vector< double > > f
2d vector to store already calculated results.
size_t numValuesPerPair
Typically one needs erfc(a_i * rij), where the number of a_i's correspond to numValuesPerPair.
double getf(size_t const atomIndex, size_t const neighIndex, size_t const valIndex, double const x)
Either returns already stored value erfc(x) or calculates it if not yet stored.
void reset(std::vector< Atom > const &atoms, size_t const valuesPerPair)
Resizes and resets the storage array to fit the current structure.