OPAL (Object Oriented Parallel Accelerator Library)  2021.1.99
OPAL
Public Types | Public Member Functions | Private Attributes | List of all members
FFTBoxPoissonSolver Class Reference

#include <FFTBoxPoissonSolver.h>

Inheritance diagram for FFTBoxPoissonSolver:
Inheritance graph
[legend]
Collaboration diagram for FFTBoxPoissonSolver:
Collaboration graph
[legend]

Public Types

typedef FFT< SineTransform, 3, double > SINE_t
 

Public Member Functions

 FFTBoxPoissonSolver (PartBunch &bunch, std::string greensFuntion)
 
 FFTBoxPoissonSolver (Mesh_t *mesh, FieldLayout_t *fl, std::string greensFunction, double boxSize)
 
 ~FFTBoxPoissonSolver ()
 
void computePotential (Field_t &rho, Vector_t hr, double zshift)
 
void computePotential (Field_t &rho, Vector_t hr)
 
void greensFunction ()
 
void integratedGreensFunction ()
 compute the integrated Green function as described in Three-dimensional quasistatic model for high brightness beam dynamics simulation by Qiang et al. More...
 
void shiftedIntGreensFunction (double zshift)
 compute the shifted integrated Green function as described in Three-dimensional quasistatic model for high brightness beam dynamics simulation by Qiang et al. More...
 
double getXRangeMin (unsigned short)
 
double getXRangeMax (unsigned short)
 
double getYRangeMin (unsigned short)
 
double getYRangeMax (unsigned short)
 
double getZRangeMin (unsigned short)
 
double getZRangeMax (unsigned short)
 
void test (PartBunchBase< double, 3 > *)
 
Inform & print (Inform &os) const
 
- Public Member Functions inherited from PoissonSolver
virtual void solve (AmrScalarFieldContainer_t &, AmrScalarFieldContainer_t &, AmrVectorFieldContainer_t &, unsigned short, unsigned short, bool=true)
 
virtual void hasToRegrid ()
 
virtual ~PoissonSolver ()
 
virtual void resizeMesh (Vector_t &, Vector_t &, const Vector_t &, const Vector_t &, double)
 

Private Attributes

Field_t rho2_m
 
Field_t grntr_m
 
IField_t grnIField_m [3]
 
SINE_t * sine_m
 
Mesh_t * mesh_m
 
FieldLayout_t * layout_m
 
Mesh_t * mesh2_m
 
FieldLayout_t * layout2_m
 
Field_t tmpgreen
 
NDIndex< 3 > domain_m
 
NDIndex< 3 > domain2_m
 
Vector_t hr_m
 
Vektor< int, 3 > nr_m
 
std::string greensFunction_m
 
double a_m
 
IpplTimings::TimerRef GreensFunctionTimer_m
 
IpplTimings::TimerRef IntGreensFunctionTimer1_m
 
IpplTimings::TimerRef IntGreensFunctionTimer2_m
 
IpplTimings::TimerRef IntGreensFunctionTimer3_m
 
IpplTimings::TimerRef IntGreensFunctionTimer4_m
 
IpplTimings::TimerRef ShIntGreensFunctionTimer1_m
 
IpplTimings::TimerRef ShIntGreensFunctionTimer2_m
 
IpplTimings::TimerRef ShIntGreensFunctionTimer3_m
 
IpplTimings::TimerRef ShIntGreensFunctionTimer4_m
 
IpplTimings::TimerRef GreensFunctionTimer1_m
 
IpplTimings::TimerRef GreensFunctionTimer4_m
 

Additional Inherited Members

- Protected Types inherited from PoissonSolver
typedef Field< int, 3, Mesh_t, Center_t > IField_t
 
typedef Field< std::complex< double >, 3, Mesh_t, Center_t > CxField_t
 

Detailed Description

Definition at line 32 of file FFTBoxPoissonSolver.h.

Member Typedef Documentation

◆ SINE_t

typedef FFT<SineTransform, 3, double> FFTBoxPoissonSolver::SINE_t

Definition at line 34 of file FFTBoxPoissonSolver.h.

Constructor & Destructor Documentation

◆ FFTBoxPoissonSolver() [1/2]

FFTBoxPoissonSolver::FFTBoxPoissonSolver ( PartBunch &  bunch,
std::string  greensFuntion 
)

Definition at line 110 of file FFTBoxPoissonSolver.cpp.

References domain_m, UniformCartesian< Dim, MFLOAT >::get_meshSpacing(), BareField< T, Dim >::getDomain(), FieldLayout< Dim >::getDomain(), IpplTimings::getTimer(), greensFunction_m, GreensFunctionTimer1_m, GreensFunctionTimer4_m, GreensFunctionTimer_m, grnIField_m, grntr_m, hr_m, Field< T, Dim, M, C >::initialize(), IntGreensFunctionTimer1_m, IntGreensFunctionTimer2_m, IntGreensFunctionTimer3_m, IntGreensFunctionTimer4_m, layout_m, cmp_diff::lt(), mesh_m, nr_m, ShIntGreensFunctionTimer1_m, ShIntGreensFunctionTimer2_m, ShIntGreensFunctionTimer3_m, ShIntGreensFunctionTimer4_m, sine_m, tmpgreen, and where().

Here is the call graph for this function:

◆ FFTBoxPoissonSolver() [2/2]

FFTBoxPoissonSolver::FFTBoxPoissonSolver ( Mesh_t *  mesh,
FieldLayout_t *  fl,
std::string  greensFunction,
double  boxSize 
)

Definition at line 37 of file FFTBoxPoissonSolver.cpp.

References domain_m, UniformCartesian< Dim, MFLOAT >::get_meshSpacing(), BareField< T, Dim >::getDomain(), FieldLayout< Dim >::getDomain(), IpplTimings::getTimer(), greensFunction_m, GreensFunctionTimer1_m, GreensFunctionTimer4_m, GreensFunctionTimer_m, grnIField_m, grntr_m, hr_m, Field< T, Dim, M, C >::initialize(), IntGreensFunctionTimer1_m, IntGreensFunctionTimer2_m, IntGreensFunctionTimer3_m, IntGreensFunctionTimer4_m, layout_m, cmp_diff::lt(), mesh_m, nr_m, ShIntGreensFunctionTimer1_m, ShIntGreensFunctionTimer2_m, ShIntGreensFunctionTimer3_m, ShIntGreensFunctionTimer4_m, sine_m, tmpgreen, and where().

Here is the call graph for this function:

◆ ~FFTBoxPoissonSolver()

FFTBoxPoissonSolver::~FFTBoxPoissonSolver ( )

Definition at line 178 of file FFTBoxPoissonSolver.cpp.

References sine_m.

Member Function Documentation

◆ computePotential() [1/2]

void FFTBoxPoissonSolver::computePotential ( Field_t &  rho,
Vector_t  hr 
)
virtual

Implements PoissonSolver.

Definition at line 232 of file FFTBoxPoissonSolver.cpp.

References greensFunction(), greensFunction_m, GreensFunctionTimer_m, grntr_m, hr_m, integratedGreensFunction(), sine_m, IpplTimings::startTimer(), and IpplTimings::stopTimer().

Here is the call graph for this function:

◆ computePotential() [2/2]

void FFTBoxPoissonSolver::computePotential ( Field_t &  rho,
Vector_t  hr,
double  zshift 
)
virtual

Implements PoissonSolver.

Definition at line 188 of file FFTBoxPoissonSolver.cpp.

References GreensFunctionTimer_m, grntr_m, hr_m, nr_m, shiftedIntGreensFunction(), sine_m, IpplTimings::startTimer(), and IpplTimings::stopTimer().

Here is the call graph for this function:

◆ getXRangeMax()

double FFTBoxPoissonSolver::getXRangeMax ( unsigned short  )
inlinevirtual

Implements PoissonSolver.

Definition at line 63 of file FFTBoxPoissonSolver.h.

References a_m.

◆ getXRangeMin()

double FFTBoxPoissonSolver::getXRangeMin ( unsigned short  )
inlinevirtual

Implements PoissonSolver.

Definition at line 62 of file FFTBoxPoissonSolver.h.

References a_m.

◆ getYRangeMax()

double FFTBoxPoissonSolver::getYRangeMax ( unsigned short  )
inlinevirtual

Implements PoissonSolver.

Definition at line 65 of file FFTBoxPoissonSolver.h.

References a_m.

◆ getYRangeMin()

double FFTBoxPoissonSolver::getYRangeMin ( unsigned short  )
inlinevirtual

Implements PoissonSolver.

Definition at line 64 of file FFTBoxPoissonSolver.h.

References a_m.

◆ getZRangeMax()

double FFTBoxPoissonSolver::getZRangeMax ( unsigned short  )
inlinevirtual

Implements PoissonSolver.

Definition at line 67 of file FFTBoxPoissonSolver.h.

References a_m.

◆ getZRangeMin()

double FFTBoxPoissonSolver::getZRangeMin ( unsigned short  )
inlinevirtual

Implements PoissonSolver.

Definition at line 66 of file FFTBoxPoissonSolver.h.

References a_m.

◆ greensFunction()

void FFTBoxPoissonSolver::greensFunction ( )

Definition at line 270 of file FFTBoxPoissonSolver.cpp.

References GreensFunctionTimer1_m, GreensFunctionTimer4_m, grnIField_m, grntr_m, hr_m, sine_m, IpplTimings::startTimer(), and IpplTimings::stopTimer().

Referenced by computePotential().

Here is the call graph for this function:

◆ integratedGreensFunction()

void FFTBoxPoissonSolver::integratedGreensFunction ( )

compute the integrated Green function as described in Three-dimensional quasistatic model for high brightness beam dynamics simulation by Qiang et al.

If the beam has a longitudinal size >> transverse size the direct Green function at each mesh point is not efficient (needs a lot of mesh points along the transverse size to get a good resolution)

If the charge density function is uniform within each cell the following Green's function can be defined:

\[ \overline{G}(x_i - x_{i'}, y_j - y_{j'}, z_k - z_{k'} cout << I << endl; cout << J << endl; cout << K << endl; cout << IE << endl; cout << JE << endl; cout << KE << endl; ) = \int_{x_{i'} - h_x/2}^{x_{i'} + h_x/2} dx' \int_{y_{j'} - h_y/2}^{y_{j'} + h_y/2} dy' \int_{z_{k'} - h_z/2}^{z_{k'} + h_z/2} dz' G(x_i - x_{i'}, y_j - y_{j'}, z_k - z_{k'}). \]

This integral can be calculated analytically in a closed from:

Definition at line 305 of file FFTBoxPoissonSolver.cpp.

References atan(), FieldLayout< Dim >::getLocalNDIndex(), grntr_m, hr_m, IntGreensFunctionTimer1_m, IntGreensFunctionTimer2_m, IntGreensFunctionTimer3_m, IntGreensFunctionTimer4_m, layout2_m, log(), min(), nr_m, sine_m, sqrt(), IpplTimings::startTimer(), IpplTimings::stopTimer(), and tmpgreen.

Referenced by computePotential().

Here is the call graph for this function:

◆ print()

Inform & FFTBoxPoissonSolver::print ( Inform &  os) const

Definition at line 512 of file FFTBoxPoissonSolver.cpp.

References endl(), and hr_m.

Here is the call graph for this function:

◆ shiftedIntGreensFunction()

void FFTBoxPoissonSolver::shiftedIntGreensFunction ( double  zshift)

compute the shifted integrated Green function as described in Three-dimensional quasistatic model for high brightness beam dynamics simulation by Qiang et al.

(x[0:nr_m[0]-1]^2 + y[0:nr_m[1]-1]^2 + (z_c + z[0:nr_m[2]-1])^2)^{-0.5} (x[nr_m[0]:1]^2 + y[0:nr_m[1]-1]^2 + (z_c + z[0:nr_m[2]-1])^2)^{-0.5} (x[0:nr_m[0]-1]^2 + y[nr_m[1]:1]^2 + (z_c + z[0:nr_m[2]-1])^2)^{-0.5} (x[nr_m[0]:1]^2 + y[nr_m[1]:1]^2 + (z_c + z[0:nr_m[2]-1])^2)^{-0.5}

(x[0:nr_m[0]-1]^2 + y[0:nr_m[1]-1]^2 + (z_c - z[nr_m[2]:1])^2)^{-0.5} (x[nr_m[0]:1]^2 + y[0:nr_m[1]-1]^2 + (z_c - z[nr_m[2]:1])^2)^{-0.5} (x[0:nr_m[0]-1]^2 + y[nr_m[1]:1]^2 + (z_c - z[nr_m[2]:1])^2)^{-0.5} (x[nr_m[0]:1]^2 + y[nr_m[1]:1]^2 + (z_c - z[nr_m[2]:1])^2)^{-0.5}

Definition at line 387 of file FFTBoxPoissonSolver.cpp.

References atan(), FieldLayout< Dim >::getLocalNDIndex(), grntr_m, hr_m, Field< T, Dim, M, C >::initialize(), layout2_m, log(), mesh2_m, min(), nr_m, ShIntGreensFunctionTimer1_m, ShIntGreensFunctionTimer2_m, ShIntGreensFunctionTimer3_m, ShIntGreensFunctionTimer4_m, sine_m, sqrt(), IpplTimings::startTimer(), IpplTimings::stopTimer(), and tmpgreen.

Referenced by computePotential().

Here is the call graph for this function:

◆ test()

void FFTBoxPoissonSolver::test ( PartBunchBase< double, 3 > *  )
inlinevirtual

Implements PoissonSolver.

Definition at line 68 of file FFTBoxPoissonSolver.h.

Member Data Documentation

◆ a_m

double FFTBoxPoissonSolver::a_m
private

Definition at line 110 of file FFTBoxPoissonSolver.h.

Referenced by getXRangeMax(), getXRangeMin(), getYRangeMax(), getYRangeMin(), getZRangeMax(), and getZRangeMin().

◆ domain2_m

NDIndex<3> FFTBoxPoissonSolver::domain2_m
private

Definition at line 102 of file FFTBoxPoissonSolver.h.

◆ domain_m

NDIndex<3> FFTBoxPoissonSolver::domain_m
private

Definition at line 99 of file FFTBoxPoissonSolver.h.

Referenced by FFTBoxPoissonSolver().

◆ greensFunction_m

std::string FFTBoxPoissonSolver::greensFunction_m
private

Definition at line 108 of file FFTBoxPoissonSolver.h.

Referenced by computePotential(), and FFTBoxPoissonSolver().

◆ GreensFunctionTimer1_m

IpplTimings::TimerRef FFTBoxPoissonSolver::GreensFunctionTimer1_m
private

Definition at line 124 of file FFTBoxPoissonSolver.h.

Referenced by FFTBoxPoissonSolver(), and greensFunction().

◆ GreensFunctionTimer4_m

IpplTimings::TimerRef FFTBoxPoissonSolver::GreensFunctionTimer4_m
private

Definition at line 125 of file FFTBoxPoissonSolver.h.

Referenced by FFTBoxPoissonSolver(), and greensFunction().

◆ GreensFunctionTimer_m

IpplTimings::TimerRef FFTBoxPoissonSolver::GreensFunctionTimer_m
private

Definition at line 112 of file FFTBoxPoissonSolver.h.

Referenced by computePotential(), and FFTBoxPoissonSolver().

◆ grnIField_m

IField_t FFTBoxPoissonSolver::grnIField_m[3]
private

Definition at line 83 of file FFTBoxPoissonSolver.h.

Referenced by FFTBoxPoissonSolver(), and greensFunction().

◆ grntr_m

Field_t FFTBoxPoissonSolver::grntr_m
private

Definition at line 79 of file FFTBoxPoissonSolver.h.

Referenced by computePotential(), FFTBoxPoissonSolver(), greensFunction(), integratedGreensFunction(), and shiftedIntGreensFunction().

◆ hr_m

Vector_t FFTBoxPoissonSolver::hr_m
private

Definition at line 105 of file FFTBoxPoissonSolver.h.

Referenced by computePotential(), FFTBoxPoissonSolver(), greensFunction(), integratedGreensFunction(), print(), and shiftedIntGreensFunction().

◆ IntGreensFunctionTimer1_m

IpplTimings::TimerRef FFTBoxPoissonSolver::IntGreensFunctionTimer1_m
private

Definition at line 114 of file FFTBoxPoissonSolver.h.

Referenced by FFTBoxPoissonSolver(), and integratedGreensFunction().

◆ IntGreensFunctionTimer2_m

IpplTimings::TimerRef FFTBoxPoissonSolver::IntGreensFunctionTimer2_m
private

Definition at line 115 of file FFTBoxPoissonSolver.h.

Referenced by FFTBoxPoissonSolver(), and integratedGreensFunction().

◆ IntGreensFunctionTimer3_m

IpplTimings::TimerRef FFTBoxPoissonSolver::IntGreensFunctionTimer3_m
private

Definition at line 116 of file FFTBoxPoissonSolver.h.

Referenced by FFTBoxPoissonSolver(), and integratedGreensFunction().

◆ IntGreensFunctionTimer4_m

IpplTimings::TimerRef FFTBoxPoissonSolver::IntGreensFunctionTimer4_m
private

Definition at line 117 of file FFTBoxPoissonSolver.h.

Referenced by FFTBoxPoissonSolver(), and integratedGreensFunction().

◆ layout2_m

FieldLayout_t* FFTBoxPoissonSolver::layout2_m
private

Definition at line 93 of file FFTBoxPoissonSolver.h.

Referenced by integratedGreensFunction(), and shiftedIntGreensFunction().

◆ layout_m

FieldLayout_t* FFTBoxPoissonSolver::layout_m
private

Definition at line 89 of file FFTBoxPoissonSolver.h.

Referenced by FFTBoxPoissonSolver().

◆ mesh2_m

Mesh_t* FFTBoxPoissonSolver::mesh2_m
private

Definition at line 92 of file FFTBoxPoissonSolver.h.

Referenced by shiftedIntGreensFunction().

◆ mesh_m

Mesh_t* FFTBoxPoissonSolver::mesh_m
private

Definition at line 88 of file FFTBoxPoissonSolver.h.

Referenced by FFTBoxPoissonSolver().

◆ nr_m

Vektor<int, 3> FFTBoxPoissonSolver::nr_m
private

Definition at line 106 of file FFTBoxPoissonSolver.h.

Referenced by computePotential(), FFTBoxPoissonSolver(), integratedGreensFunction(), and shiftedIntGreensFunction().

◆ rho2_m

Field_t FFTBoxPoissonSolver::rho2_m
private

Definition at line 76 of file FFTBoxPoissonSolver.h.

◆ ShIntGreensFunctionTimer1_m

IpplTimings::TimerRef FFTBoxPoissonSolver::ShIntGreensFunctionTimer1_m
private

Definition at line 119 of file FFTBoxPoissonSolver.h.

Referenced by FFTBoxPoissonSolver(), and shiftedIntGreensFunction().

◆ ShIntGreensFunctionTimer2_m

IpplTimings::TimerRef FFTBoxPoissonSolver::ShIntGreensFunctionTimer2_m
private

Definition at line 120 of file FFTBoxPoissonSolver.h.

Referenced by FFTBoxPoissonSolver(), and shiftedIntGreensFunction().

◆ ShIntGreensFunctionTimer3_m

IpplTimings::TimerRef FFTBoxPoissonSolver::ShIntGreensFunctionTimer3_m
private

Definition at line 121 of file FFTBoxPoissonSolver.h.

Referenced by FFTBoxPoissonSolver(), and shiftedIntGreensFunction().

◆ ShIntGreensFunctionTimer4_m

IpplTimings::TimerRef FFTBoxPoissonSolver::ShIntGreensFunctionTimer4_m
private

Definition at line 122 of file FFTBoxPoissonSolver.h.

Referenced by FFTBoxPoissonSolver(), and shiftedIntGreensFunction().

◆ sine_m

SINE_t* FFTBoxPoissonSolver::sine_m
private

Definition at line 85 of file FFTBoxPoissonSolver.h.

Referenced by computePotential(), FFTBoxPoissonSolver(), greensFunction(), integratedGreensFunction(), shiftedIntGreensFunction(), and ~FFTBoxPoissonSolver().

◆ tmpgreen

Field_t FFTBoxPoissonSolver::tmpgreen
private

Definition at line 96 of file FFTBoxPoissonSolver.h.

Referenced by FFTBoxPoissonSolver(), integratedGreensFunction(), and shiftedIntGreensFunction().

The documentation for this class was generated from the following files:
Generated on Wed Aug 25 2021 16:41:11 for OPAL (Object Oriented Parallel Accelerator Library) by doxygen 1.9.1