ESPResSo
Extensible Simulation Package for Research on Soft Matter Systems
Loading...
Searching...
No Matches
ljcos2.cpp
Go to the documentation of this file.
1
/*
2
* Copyright (C) 2010-2022 The ESPResSo project
3
* Copyright (C) 2002,2003,2004,2005,2006,2007,2008,2009,2010
4
* Max-Planck-Institute for Polymer Research, Theory Group
5
*
6
* This file is part of ESPResSo.
7
*
8
* ESPResSo is free software: you can redistribute it and/or modify
9
* it under the terms of the GNU General Public License as published by
10
* the Free Software Foundation, either version 3 of the License, or
11
* (at your option) any later version.
12
*
13
* ESPResSo is distributed in the hope that it will be useful,
14
* but WITHOUT ANY WARRANTY; without even the implied warranty of
15
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16
* GNU General Public License for more details.
17
*
18
* You should have received a copy of the GNU General Public License
19
* along with this program. If not, see <http://www.gnu.org/licenses/>.
20
*/
21
/** \file
22
*
23
* Implementation of \ref ljcos2.hpp
24
*/
25
#include "
ljcos2.hpp
"
26
27
#ifdef LJCOS2
28
#include "
nonbonded_interaction_data.hpp
"
29
30
#include <cmath>
31
#include <stdexcept>
32
33
LJcos2_Parameters::LJcos2_Parameters
(
double
epsilon,
double
sigma,
34
double
offset,
double
width)
35
: eps{epsilon}, sig{sigma}, offset{offset}, w{width} {
36
if
(epsilon < 0.) {
37
throw
std::domain_error(
"LJcos2 parameter 'epsilon' has to be >= 0"
);
38
}
39
if
(sigma < 0.) {
40
throw
std::domain_error(
"LJcos2 parameter 'sigma' has to be >= 0"
);
41
}
42
if
(width < 0.) {
43
throw
std::domain_error(
"LJcos2 parameter 'width' has to be >= 0"
);
44
}
45
if
(sigma != 0.) {
46
rchange
= std::pow(2., 1. / 6.) * sigma;
47
cut
= width +
rchange
;
48
}
49
}
50
51
#endif
/* ifdef LJCOS2 */
ljcos2.hpp
Routines to calculate the Lennard-Jones with cosine tail potential between particle pairs.
nonbonded_interaction_data.hpp
Various procedures concerning interactions between particles.
LJcos2_Parameters::cut
double cut
Definition
nonbonded_interaction_data.hpp:223
LJcos2_Parameters::LJcos2_Parameters
LJcos2_Parameters()=default
LJcos2_Parameters::rchange
double rchange
Definition
nonbonded_interaction_data.hpp:226
src
core
nonbonded_interactions
ljcos2.cpp
Generated on Fri Nov 22 2024 02:13:54 for ESPResSo by
1.9.8