MercuryDPM  Trunk
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Dipole.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 #include "Dipole.h"
26 #include "Multipole.h"
27 #include "Math/ExtendedMath.h"
28 #include "Math/NumericalVector.h"
29 #include <complex>
30 #include <vector>
31 
32 Dipole::Dipole(int p, NumericalVector<>* squaredFactorials, Vec3D location, Vec3D velocity, Mdouble strength) :
33  Multipole(p, squaredFactorials, location),
34  velocity_(velocity),
35  strength_(strength)
36 {
37 }
38 
40 {
41  size_t nTerms = (p_ + 1) * (p_ + 1);
42  NumericalVector<std::complex<Mdouble>> multipoleExpansionCoefficients(nTerms);
43 
44  //Calculate dipole coefficients for spherical harmonics
48 
49  multipoleExpansionCoefficients[1] = 1.0 / std::sqrt(2.0) * (-s1 + s2 / constants::i);
50  multipoleExpansionCoefficients[2] = s3;
51  multipoleExpansionCoefficients[3] = -1.0 / (sqrt(2.0)) * (s1 + s2 / constants::i);
52 
53  multipoleExpansionCoefficients_ = multipoleExpansionCoefficients;
54 }
55 
Vec3D velocity_
Definition: Dipole.h:47
Mdouble strength_
Definition: Dipole.h:48
double Mdouble
Definition: GeneralDefine.h:34
const std::complex< Mdouble > i
Definition: ExtendedMath.h:50
int p_
Definition: Multipole.h:71
Dipole(int p, NumericalVector<> *squaredFactorials, Vec3D location, Vec3D velocity, Mdouble strength)
Definition: Dipole.cc:32
NumericalVector< std::complex< Mdouble > > multipoleExpansionCoefficients_
Definition: Multipole.h:74
Mdouble getComponent(int index) const
Returns the requested component of this Vec3D.
Definition: Vector.cc:194
void computeMultipoleExpansion() final
Definition: Dipole.cc:39
Definition: Vector.h:49