MercuryDPM  Trunk
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Base_XY_XZ_YZ.cc
Go to the documentation of this file.
1 //Copyright (c) 2013-2020, The MercuryDPM Developers Team. All rights reserved.
2 //For the list of developers, see <http://www.MercuryDPM.org/Team>.
3 //
4 //Redistribution and use in source and binary forms, with or without
5 //modification, are permitted provided that the following conditions are met:
6 // * Redistributions of source code must retain the above copyright
7 // notice, this list of conditions and the following disclaimer.
8 // * Redistributions in binary form must reproduce the above copyright
9 // notice, this list of conditions and the following disclaimer in the
10 // documentation and/or other materials provided with the distribution.
11 // * Neither the name MercuryDPM nor the
12 // names of its contributors may be used to endorse or promote products
13 // derived from this software without specific prior written permission.
14 //
15 //THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND
16 //ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
17 //WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
18 //DISCLAIMED. IN NO EVENT SHALL THE MERCURYDPM DEVELOPERS TEAM BE LIABLE FOR ANY
19 //DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
20 //(INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
21 //LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
22 //ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
23 //(INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
24 //SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
25 
26 #include "Base_XY_XZ_YZ.h"
27 #include "Particles/BaseParticle.h"
28 #include "DPMBase.h"
29 
30 using namespace CGCoordinates;
31 
33 {
34  //Wolfram alpha: integrate(x*exp(-x^2/(2w^2)),{x,0,c})/integrate(x*exp(-x^2/(2w^2)),{x,0,inf})=1-e^(-c^2/(2 w^2))
35  Mdouble prefactor = 1.0 / (constants::sqrt_2 * constants::sqrt_pi * width);
36  return mathsFunc::square(prefactor) / (1.0 - exp(-0.5 * mathsFunc::square(cutoff / width)));
37 }
38 
40 {
41  Mdouble widthSqrt2 = width * constants::sqrt_2;
42  Mdouble a = -cutoff;
43  Mdouble b = cutoff + distance;
44  //1D prefactor
45  Mdouble prefactor_ = 1.0 / (widthSqrt2 * constants::sqrt_pi);
46  prefactor_ /= erf(cutoff / (widthSqrt2));
47  return prefactor_ * 0.5 / (
48  +erf(b / widthSqrt2) * b
49  + widthSqrt2 / constants::sqrt_pi * exp(-mathsFunc::square(b / widthSqrt2))
50  - erf(a / widthSqrt2) * a
51  - widthSqrt2 / constants::sqrt_pi * exp(-mathsFunc::square(a / widthSqrt2))
52  );
53 }
54 
60 void Base_XY_XZ_YZ::normalisePolynomialCoefficients(std::vector<Mdouble>& coefficients, Mdouble cutoff)
61 {
62  Mdouble volume = 0.0;
63  for (std::size_t i = 0; i < coefficients.size(); i++)
64  volume += coefficients[i] / static_cast<Mdouble>(i + 2);
65  volume *= 2.0 * constants::pi * mathsFunc::square(cutoff);
66  for (double& coefficient : coefficients)
67  coefficient /= volume;
68 }
69 
71 {
72  return 2;
73 }
double Mdouble
Definition: GeneralDefine.h:34
Mdouble exp(Mdouble Exponent)
Definition: ExtendedMath.cc:84
static Mdouble getGaussIntegralPrefactor(Mdouble distance, Mdouble width, Mdouble cutoff)
Computes the prefactor of the Gauss line integral, which is dependent on the number of non-averaged d...
const std::complex< Mdouble > i
Definition: ExtendedMath.h:50
const Mdouble sqrt_pi
Definition: ExtendedMath.h:46
static Mdouble getGaussPrefactor(Mdouble width, Mdouble cutoff)
Computes the prefactor of the Gauss CGFunction, which is dependent on the number of non-averaged dime...
const Mdouble pi
Definition: ExtendedMath.h:45
static void normalisePolynomialCoefficients(std::vector< Mdouble > &coefficients, Mdouble cutoff)
Normalises the coefficients of Polynomial CGFunction such that the integral over all non-averaged dim...
const Mdouble sqrt_2
Definition: ExtendedMath.h:48
static const unsigned countVariables()
T square(const T val)
squares a number
Definition: ExtendedMath.h:104