#include <FMGPoissonSolver.h>

Inheritance diagram for FMGPoissonSolver:

Collaboration diagram for FMGPoissonSolver:

[legend]

Public Member Functions
	FMGPoissonSolver (AmrBoxLib *itsAmrObject_p)

void	solve (AmrScalarFieldContainer_t &rho, AmrScalarFieldContainer_t &phi, AmrVectorFieldContainer_t &efield, unsigned short baseLevel, unsigned short finestLevel, bool prevAsGuess=true)

double	getXRangeMin (unsigned short level=0)

double	getXRangeMax (unsigned short level=0)

double	getYRangeMin (unsigned short level=0)

double	getYRangeMax (unsigned short level=0)

double	getZRangeMin (unsigned short level=0)

double	getZRangeMax (unsigned short level=0)

Inform &	print (Inform &os) const

Public Member Functions inherited from AmrPoissonSolver< AmrBoxLib >
	AmrPoissonSolver (AmrBoxLib *itsAmrObject_p)

virtual	~AmrPoissonSolver ()

void	computePotential (Field_t &, Vector_t)

void	computePotential (Field_t &, Vector_t, double)

void	test (PartBunchBase< double, 3 > *)

void	hasToRegrid ()

Public Member Functions inherited from PoissonSolver
virtual void	computePotential (Field_t &rho, Vector_t hr)=0

virtual void	solve (AmrScalarFieldContainer_t &, AmrScalarFieldContainer_t &, AmrVectorFieldContainer_t &, unsigned short, unsigned short, bool=true)

virtual void	hasToRegrid ()

virtual void	computePotential (Field_t &rho, Vector_t hr, double zshift)=0

virtual double	getXRangeMin (unsigned short level=0)=0

virtual double	getXRangeMax (unsigned short level=0)=0

virtual double	getYRangeMin (unsigned short level=0)=0

virtual double	getYRangeMax (unsigned short level=0)=0

virtual double	getZRangeMin (unsigned short level=0)=0

virtual double	getZRangeMax (unsigned short level=0)=0

virtual void	test (PartBunchBase< double, 3 > *bunch)=0

virtual	~PoissonSolver ()

virtual void	resizeMesh (Vector_t &, Vector_t &, const Vector_t &, const Vector_t &, double)

virtual double	getinteractionRadius () const

virtual void	calculatePairForces (PartBunchBase< double, 3 > *, double)

Private Types
typedef AmrBoxLib::AmrGeomContainer_t	GeomContainer_t

typedef amrex::Vector< AmrBoxLib::AmrField_t * >	AmrFieldContainer_pt

typedef AmrBoxLib::AmrGeometry_t	AmrGeometry_t

typedef AmrBoxLib::AmrGrid_t	AmrGrid_t

typedef AmrBoxLib::AmrProcMap_t	AmrProcMap_t

typedef AmrBoxLib::AmrScalarFieldContainer_t	AmrScalarFieldContainer_t

typedef AmrBoxLib::AmrVectorFieldContainer_t	AmrVectorFieldContainer_t

Private Member Functions
void	initParameters_m ()

double	solveWithF90_m (const AmrFieldContainer_pt &rho, const AmrFieldContainer_pt &phi, const amrex::Vector< AmrFieldContainer_pt > &grad_phi_edge, const GeomContainer_t &geom, int baseLevel, int finestLevel)

Private Attributes
int	bc_m [2 *AMREX_SPACEDIM]
	Boundary conditions. More...

double	reltol_m
	Relative tolearance for solver. More...

double	abstol_m
	Absolute tolerance for solver. More...

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

Protected Attributes inherited from AmrPoissonSolver< AmrBoxLib >
AmrBoxLib *	itsAmrObject_mp

bool	regrid_m
	is set to true by itsAmrObject_mp and reset to false by solver More...

Detailed Description

Definition at line 33 of file FMGPoissonSolver.h.

Member Typedef Documentation

◆ AmrFieldContainer_pt

typedef amrex::Vector<AmrBoxLib::AmrField_t*> FMGPoissonSolver::AmrFieldContainer_pt

private

Definition at line 37 of file FMGPoissonSolver.h.

◆ AmrGeometry_t

typedef AmrBoxLib::AmrGeometry_t FMGPoissonSolver::AmrGeometry_t

private

Definition at line 38 of file FMGPoissonSolver.h.

◆ AmrGrid_t

typedef AmrBoxLib::AmrGrid_t FMGPoissonSolver::AmrGrid_t

private

Definition at line 39 of file FMGPoissonSolver.h.

◆ AmrProcMap_t

typedef AmrBoxLib::AmrProcMap_t FMGPoissonSolver::AmrProcMap_t

private

Definition at line 40 of file FMGPoissonSolver.h.

◆ AmrScalarFieldContainer_t

typedef AmrBoxLib::AmrScalarFieldContainer_t FMGPoissonSolver::AmrScalarFieldContainer_t

private

Definition at line 41 of file FMGPoissonSolver.h.

◆ AmrVectorFieldContainer_t

typedef AmrBoxLib::AmrVectorFieldContainer_t FMGPoissonSolver::AmrVectorFieldContainer_t

private

Definition at line 42 of file FMGPoissonSolver.h.

◆ GeomContainer_t

typedef AmrBoxLib::AmrGeomContainer_t FMGPoissonSolver::GeomContainer_t

private

Definition at line 36 of file FMGPoissonSolver.h.

Constructor & Destructor Documentation

◆ FMGPoissonSolver()

FMGPoissonSolver::FMGPoissonSolver ( AmrBoxLib * itsAmrObject_p )

This solver only works with AmrBoxLib. In order the solver to work the cells need to be of cubic shape. Otherwise the solver stops with the error message that it did not converge (AMReX internal).

Parameters

itsAmrObject_p has information about refinemen ratios, etc.

Definition at line 32 of file FMGPoissonSolver.cpp.

References bc_m, and initParameters_m().

Here is the call graph for this function:

Member Function Documentation

◆ getXRangeMax()

double FMGPoissonSolver::getXRangeMax ( unsigned short level = 0 )

virtual

Implements PoissonSolver.

Definition at line 130 of file FMGPoissonSolver.cpp.

References AmrPoissonSolver< AmrBoxLib >::itsAmrObject_mp.

◆ getXRangeMin()

double FMGPoissonSolver::getXRangeMin ( unsigned short level = 0 )

virtual

Implements PoissonSolver.

Definition at line 125 of file FMGPoissonSolver.cpp.

References AmrPoissonSolver< AmrBoxLib >::itsAmrObject_mp.

◆ getYRangeMax()

double FMGPoissonSolver::getYRangeMax ( unsigned short level = 0 )

virtual

Implements PoissonSolver.

Definition at line 140 of file FMGPoissonSolver.cpp.

References AmrPoissonSolver< AmrBoxLib >::itsAmrObject_mp.

◆ getYRangeMin()

double FMGPoissonSolver::getYRangeMin ( unsigned short level = 0 )

virtual

Implements PoissonSolver.

Definition at line 135 of file FMGPoissonSolver.cpp.

References AmrPoissonSolver< AmrBoxLib >::itsAmrObject_mp.

◆ getZRangeMax()

double FMGPoissonSolver::getZRangeMax ( unsigned short level = 0 )

virtual

Implements PoissonSolver.

Definition at line 150 of file FMGPoissonSolver.cpp.

References AmrPoissonSolver< AmrBoxLib >::itsAmrObject_mp.

◆ getZRangeMin()

double FMGPoissonSolver::getZRangeMin ( unsigned short level = 0 )

virtual

Implements PoissonSolver.

Definition at line 145 of file FMGPoissonSolver.cpp.

References AmrPoissonSolver< AmrBoxLib >::itsAmrObject_mp.

◆ initParameters_m()

void FMGPoissonSolver::initParameters_m ( )

private

Initialize additional parameters for the multigrid solver that are given by ParmParse. The function is called in the constructory only.

Definition at line 164 of file FMGPoissonSolver.cpp.

Referenced by FMGPoissonSolver().

◆ print()

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

Print information abour tolerances.

Parameters

os	output stream where to write to

Definition at line 155 of file FMGPoissonSolver.cpp.

References abstol_m, endl(), and reltol_m.

Here is the call graph for this function:

◆ solve()

void FMGPoissonSolver::solve	(	AmrScalarFieldContainer_t &	rho,
		AmrScalarFieldContainer_t &	phi,
		AmrVectorFieldContainer_t &	efield,
		unsigned short	baseLevel,
		unsigned short	finestLevel,
		bool	prevAsGuess = `true`
	)

virtual

Multigrid solve based on AMReX FMultiGrid solver. The relative tolerance is set to 1.0e-9 and the absolute tolerance to 0.0.

Parameters

rho	right-hand side charge density on grid [C / m]
phi	electrostatic potential (unknown) [V]
efield	electric field [V / m]
baseLevel	for solve
finestLevel	for solve
prevAsGuess	use of previous solution as initial guess

Reimplemented from PoissonSolver.

Definition at line 46 of file FMGPoissonSolver.cpp.

References abstol_m, bc_m, AmrPoissonSolver< AmrBoxLib >::itsAmrObject_mp, Hypervolume::n, PAssert, and solveWithF90_m().

Here is the call graph for this function:

◆ solveWithF90_m()

double FMGPoissonSolver::solveWithF90_m	(	const AmrFieldContainer_pt &	rho,
		const AmrFieldContainer_pt &	phi,
		const amrex::Vector< AmrFieldContainer_pt > &	grad_phi_edge,
		const GeomContainer_t &	geom,
		int	baseLevel,
		int	finestLevel
	)

private

Does the actual solve. It calls the FMultiGrid solver of AMReX. It uses an approximation order of 3 at Dirichlet boundaries.

Parameters

rho	charge density on grids [C / m]
phi	electrostatic potential on grid [V]
grad_phi_edge	gradient of the potential (values at faces)
geom	geometry of the problem, i.e. physical domain boundaries
baseLevel	for solve
finestLevel	for solve

Returns: the residuum norm

Definition at line 264 of file FMGPoissonSolver.cpp.

References abstol_m, bc_m, AmrPoissonSolver< AmrBoxLib >::itsAmrObject_mp, and reltol_m.

Referenced by solve().

Member Data Documentation

◆ abstol_m

double FMGPoissonSolver::abstol_m

private

Absolute tolerance for solver.

Definition at line 128 of file FMGPoissonSolver.h.

Referenced by print(), solve(), and solveWithF90_m().

◆ bc_m

int FMGPoissonSolver::bc_m[2 *AMREX_SPACEDIM]

private

Boundary conditions.

In case of reusing the previous solution the grids might have changed, thus, we interpolate from the coarse grids to the fine grids

Definition at line 126 of file FMGPoissonSolver.h.

Referenced by FMGPoissonSolver(), solve(), and solveWithF90_m().

◆ reltol_m

double FMGPoissonSolver::reltol_m

private

Relative tolearance for solver.

Definition at line 127 of file FMGPoissonSolver.h.

Referenced by print(), and solveWithF90_m().

The documentation for this class was generated from the following files:

/Users/gsell/src/OPAL/src/src/Solvers/BoxLibSolvers/FMGPoissonSolver.h
/Users/gsell/src/OPAL/src/src/Solvers/BoxLibSolvers/FMGPoissonSolver.cpp

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ AmrFieldContainer_pt

◆ AmrGeometry_t

◆ AmrGrid_t

◆ AmrProcMap_t

◆ AmrScalarFieldContainer_t

◆ AmrVectorFieldContainer_t

◆ GeomContainer_t

Constructor & Destructor Documentation

◆ FMGPoissonSolver()

Member Function Documentation

◆ getXRangeMax()

◆ getXRangeMin()

◆ getYRangeMax()

◆ getYRangeMin()

◆ getZRangeMax()

◆ getZRangeMin()

◆ initParameters_m()

◆ print()

◆ solve()

◆ solveWithF90_m()

Member Data Documentation

◆ abstol_m

◆ bc_m

◆ reltol_m