ESPResSo
Extensible Simulation Package for Research on Soft Matter Systems
Loading...
Searching...
No Matches
AS_erfc_part.hpp
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
22
#ifndef UTILS_MATH_AC_ERF_PART_HPP
23
#define UTILS_MATH_AC_ERF_PART_HPP
24
25
namespace
Utils
{
26
/** Approximate \f$ \exp(d^2) \mathrm{erfc}(d)\f$ by applying a formula from
27
* @cite abramowitz65a chapter 7.
28
*/
29
template
<
typename
T>
constexpr
T
AS_erfc_part
(T d) {
30
T
const
constexpr
a1 = 0.254829592;
31
T
const
constexpr
a2 = -0.284496736;
32
T
const
constexpr
a3 = 1.421413741;
33
T
const
constexpr
a4 = -1.453152027;
34
T
const
constexpr
a5 = 1.061405429;
35
T
const
constexpr
p = 0.3275911;
36
37
auto
const
t = 1.0 / (1.0 + p * d);
38
39
return
t * (a1 + t * (a2 + t * (a3 + t * (a4 + t * a5))));
40
}
41
}
// namespace Utils
42
43
#endif
Utils
Definition
Variant.hpp:40
Utils::AS_erfc_part
constexpr T AS_erfc_part(T d)
Approximate by applying a formula from chapter 7.
Definition
AS_erfc_part.hpp:29
src
utils
include
utils
math
AS_erfc_part.hpp
Generated on Fri Nov 8 2024 02:12:53 for ESPResSo by
1.9.8