Class for calculating Eddy current problem with Maxwell modell in h formulation. More...

#include <Maxwell2D_H.hh>

Inheritance diagram for hp2D::Maxwell2D_H:

Collaboration diagram for hp2D::Maxwell2D_H:

Public Types
enum	boundaryType { PMC = 0, PEC = 1, MAX_TYPE }
	Boundary type. More...
enum	solverType { SUPERLU = 0, SUPERLU2 = 1, BICGSTAB = 2, BICGSTAB2 = 3, BICGSTABSUPERLU = 4 }
	Type of the solver. More...
typedef F	type
Public Member Functions
void	basisFunctions (const std::string path="")
	Writes matlab graphics files for each basis function.
boundaryType &	bType ()
const boundaryType	bType () const
	Returns boundary type.
const std::string	bTypeStr () const
	Returns name of boundary type as string.
virtual void	dataOut (concepts::InOutParameters &parameters, const bool store=false, const std::string path="")
	Append the parameters of the calculation to `parameters` and writes it in matlab format if `path` is given.
Real	dissipation ()
	Return dissipation power loss.
virtual Real	dissipation ()
	Return dissipation power loss.
concepts::ElementFormula < concepts::Cmplx2d > *	eField ()
	Returns a pointer to the e-Field.
virtual void	graphicsOut (const uint numPoints, const std::string path="")
	Writes matlab graphics files.
concepts::ElementFormula< Cmplx > *	hField ()
	Returns a pointer to the h-Field (solution + h0).
virtual void	loadSolution (const uint l=0, const uint p=1, const uint g=0, const std::string path="")
	Loads solution vector from standard file in `path` with use of `l`, `p` and `g`.
virtual void	loadSolution (const std::string vectorFile)
	Loads solution vector from `vectorFile`.
Real	magnEnergy ()
	Return magnetic energy.
virtual Real	magnEnergy ()
	Return magnetic energy.
void	matrixBuild (const uint l=0, const uint p=1, const uint g=0, const std::string matrixFile="")
	Build the space, matrices and the load vector.
void	matrixBuild (InputModels &input)
	Maxwell2D_H (concepts::EddyGeometry2D &geom, InputMaxwell2D_H &input, const uint geomRefAttrib=100)
	Maxwell2D_H (concepts::EddyGeometry2D &geom, InputMaxwell2D_H &input)
	Maxwell2D_H (concepts::EddyGeometry2D &geom, enum solverType type=SUPERLU, bool diagPrecond=true, bool afterIter=false, const Real eps=EPS0, const Real omega=OMEGA50, const Real mu=MU0)
	Maxwell2D_H (concepts::EddyGeometry2D &geom, enum concepts::MaxwellBoundary::boundaryType bType=PMC, enum solverType type=SUPERLU, bool diagPrecond=true, bool afterIter=false, const Real eps=EPS0, const Real omega=OMEGA50, const Real mu=MU0, const uint geomRefAttrib=100)
	Constructor.
void	rebuildMesh (concepts::InputAdaptiveModels &input)
void	rebuildMesh (const uint l=0, const uint p=1, const uint g=0, const uint subdiv=X\|Y)
	Rebuilds only the mesh and sets the polynomial degrees.
const Vector< F > *	solution ()
	Returns solution vector.
void	solve (InputModels &input)
void	solve (const uint l=0, const uint p=1, const uint g=0, const std::string vectorFile="", const std::string matrixFile="")
	Build the space, matrices and the load vector and solve the system.
virtual hpAdaptiveSpaceH1 &	space () const
	Returns the space.
void	spaceBuild (const uint l=0, const uint p=1, const uint g=0)
	Builds only the space.
virtual	~Maxwell2D_H ()
virtual	~Maxwell2D_H ()
Protected Types
enum	subdivTypes
Protected Member Functions
std::string	filename_ (const std::string varname, const std::string subpath, const bool p=false)
	Gives a file name, which consists of p, l, g.
void	formulaOut_ (const concepts::PiecewiseFormulaBase< G > &frm, const std::string frmstr, const std::string path="")
	graphical output of formulas
concepts::SparseMatrix< Real > *	identityMatrix_ (concepts::Space< Real > &spc, concepts::SparseMatrix< Cmplx > *S)
	Calculate identity matrix and add's it to system matrix `S`.
virtual std::ostream &	info (std::ostream &os) const
	Returns information in an output stream.
virtual std::ostream &	info (std::ostream &os) const
	Returns information in an output stream.
concepts::SparseMatrix< Cmplx > *	laplaceMatrix_ (concepts::Space< Real > &spc, concepts::SparseMatrix< Cmplx > *S)
	Calculate stiffness matrix and add's it to system matrix `S`.
virtual concepts::Mesh2 &	mesh_ ()
	Mesh.
virtual const std::string	mshAbbr_ ()
	Mesh abbreviation string.
virtual hpFull &	prebuild_ ()
	Space Prebuilder.
virtual void	rebuildMesh_ (const uint l=0, const uint p=1, const uint g=0, const uint subdiv=X\|Y)
	Rebuilds only the mesh and sets the polynomial degrees.
Protected Attributes
std::auto_ptr < concepts::BoundaryConditions >	bc_
	Boundary conditions.
std::auto_ptr < concepts::BoundaryConditions >	bc_
	Boundary conditions.
const uint	geomRefAttr_
	Attribute of vertices or edges for geometric refinement.
const uint	geomRefAttr_
	Attribute of vertices or edges for geometric refinement.
uint	iterations_
	Number of iterations for iterative solver.
std::string	matrixFile_
	matlab file for matrices and load vector output
double	matrixtime_
concepts::InOutParameters	parameters_
	Parameters of the calculation.
const std::string	problemName_
	Name of the problem.
const std::string	problemName_
	Name of the problem.
std::auto_ptr< Vector< F > >	sol_
	Solution vector.
std::auto_ptr < concepts::Vector< Cmplx > >	sol_
	Solution vector.
double	solvetime_
	Time to solve the system, to build the matrices, to rebuild the space.
double	spacetime_
std::string	vectorFile_
	file for solution output, which can be reloaded into concepts
Private Member Functions
void	constructor_ ()
	private constructor
void	constructSpace_ ()
	Constructs the space.
virtual void	dataOut_ ()
	Collect the parameters.
void	identityMatrix_ ()
	Calculate identity matrix.
void	laplaceMatrix_ ()
	Calculate stiffness matrix.
void	linearform_ ()
	Calculate the load vector, assumes sigma = 0 inside the coil.
void	matrices_ ()
	Building the matrices.
virtual void	matrices_ ()
	Building the matrices.
virtual void	solve_ ()
	Method for solving, throws exception when it wasn't successfull.
virtual void	solve_ ()
	Method for solving, throws exception when it wasn't successfull.
void	systemMatrix_ ()
	Compress system matrix.
Private Attributes
std::auto_ptr < concepts::SparseMatrix < Cmplx > >	A_
	Stiffness matrix.
bool	afterIter_
	After iterations.
bool	diagPrecond_
	Using diagonal preconditioning.
std::auto_ptr< Real >	dissipation_
	Dissipation power loss.
const Real	eps_
	Dielectricity constant.
concepts::EddyGeometry2D &	geom_
	Mesh and material constants (sigma, j0)
concepts::PiecewiseFormulaFun < Cmplx, Real >	iOmegaEps_plus_Sigma_Inv_
	Piecewise constant formula $(i \omega \varepsilon_0 + \sigma)^-1$ .
std::auto_ptr < concepts::SparseMatrix< Real > >	M_
	Mass matrix.
std::auto_ptr< Real >	magnEnergy_
	Magnetic energy.
const Real	mu_
	Permeability constant.
const Real	omega_
	Angular frequency.
std::auto_ptr < concepts::Vector< Cmplx > >	residual_
	Residual of solving the linear system.
std::auto_ptr< Real >	residualNorm_
	Euclidian norm of the residual of solving the linear system.
std::auto_ptr < concepts::Vector< Cmplx > >	rhs_
concepts::Vector< Cmplx > *	rhsScaled_
	Pointer to scaled load vector if diagonal preconditioning is choosed.
std::auto_ptr < concepts::SparseMatrix < Cmplx > >	S_
	System matrix.
std::auto_ptr< hpAdaptiveSpaceH1 >	spc_
	Space.
bool	statusAfterIter_
	Internal status of after iteration false : before after iteration true : at the after iteration or it's done.
enum solverType	type_
	Solver type.
Friends
class	concepts::ModelControl< Maxwell2D_H >
class	concepts::ModelControl< Maxwell2D_H_Base >

Detailed Description

Class for calculating Eddy current problem with Maxwell modell in h formulation.

Author:: Kersten Schmidt, 2005

Definition at line 36 of file Maxwell2D_H.hh.

Member Typedef Documentation

template<class F>

typedef F concepts::Model< F >::type [inherited]

Definition at line 58 of file models.hh.

Member Enumeration Documentation

enum concepts::MaxwellBoundary::boundaryType [inherited]

Boundary type.

Either perfect magnetic conductor (PMC) or perfect electric conductor (PEC). Wether it's dirichlet or neumann boundary is dependent from the formulation.

Enumerator:

PMC
PEC
MAX_TYPE

Definition at line 25 of file maxwell.hh.

enum hp2D::Maxwell2D_H::solverType

Type of the solver.

Enumerator:

SUPERLU
SUPERLU2
BICGSTAB
BICGSTAB2
BICGSTABSUPERLU

Reimplemented from hp2D::Maxwell2D_H_Base.

Definition at line 38 of file Maxwell2D_H.hh.

enum hp2D::AdaptiveModel::subdivTypes [protected, inherited]

Reimplemented from concepts::AdaptiveModel< F, dim >.

Definition at line 82 of file adaptiveModels.hh.

Constructor & Destructor Documentation

hp2D::Maxwell2D_H::Maxwell2D_H	(	concepts::EddyGeometry2D &	geom,
		enum solverType	type = `SUPERLU`,
		bool	diagPrecond = `true`,
		bool	afterIter = `false`,
		const Real	eps = `EPS0`,
		const Real	omega = `OMEGA50`,
		const Real	mu = `MU0`
	)

hp2D::Maxwell2D_H::Maxwell2D_H	(	concepts::EddyGeometry2D &	geom,
		InputMaxwell2D_H &	input
	)

virtual hp2D::Maxwell2D_H::~Maxwell2D_H ( ) [virtual]

hp2D::Maxwell2D_H::Maxwell2D_H	(	concepts::EddyGeometry2D &	geom,
		enum concepts::MaxwellBoundary::boundaryType	bType = `PMC`,
		enum solverType	type = `SUPERLU`,
		bool	diagPrecond = `true`,
		bool	afterIter = `false`,
		const Real	eps = `EPS0`,
		const Real	omega = `OMEGA50`,
		const Real	mu = `MU0`,
		const uint	geomRefAttrib = `100`
	)

Constructor.

Parameters:

geom	geometry, conductivity and source currents
bType	type of boundary condition
type	solver type
diagPrecond	flag for preconditioning with diagonal matrix
aterIter	flag for after iterations
eps	dielectricity constant
omega	angular frequency in 1/s
mu	permeability constant in Ohm*s/m
geomRefAttrib	attrib for geometric refinement

hp2D::Maxwell2D_H::Maxwell2D_H	(	concepts::EddyGeometry2D &	geom,
		InputMaxwell2D_H &	input,
		const uint	geomRefAttrib = `100`
	)

virtual hp2D::Maxwell2D_H::~Maxwell2D_H ( ) [inline, virtual]

Definition at line 144 of file Maxwell2D_H.hh.

Member Function Documentation

void hp2D::Models< Cmplx >::basisFunctions ( const std::string path = "" ) [inherited]

Writes matlab graphics files for each basis function.

const boundaryType concepts::MaxwellBoundary::bType ( ) const [inline, inherited]

Returns boundary type.

Definition at line 31 of file maxwell.hh.

boundaryType& concepts::MaxwellBoundary::bType ( ) [inline, inherited]

Definition at line 32 of file maxwell.hh.

const std::string concepts::MaxwellBoundary::bTypeStr ( ) const [inline, inherited]

Returns name of boundary type as string.

Definition at line 34 of file maxwell.hh.

void hp2D::Maxwell2D_H::constructor_ ( ) [private]

private constructor

Reimplemented from hp2D::Maxwell2D_H_Base.

void hp2D::Maxwell2D_H::constructSpace_ ( ) [private]

Constructs the space.

virtual void hp2D::Models< Cmplx >::dataOut	(	concepts::InOutParameters &	parameters,
		const bool	store = `false`,
		const std::string	path = `""`
	)		`[virtual, inherited]`

Append the parameters of the calculation to parameters and writes it in matlab format if path is given.

virtual void hp2D::Maxwell2D_H::dataOut_ ( ) [private, virtual]

Collect the parameters.

Implements hp2D::Models< Cmplx >.

virtual Real hp2D::Maxwell2D_H::dissipation ( ) [virtual]

Return dissipation power loss.

Implements concepts::MaxwellModel.

Real hp2D::Maxwell2D_H::dissipation ( ) [virtual]

Return dissipation power loss.

Implements concepts::MaxwellModel.

concepts::ElementFormula<concepts::Cmplx2d>* hp2D::Maxwell2D_H::eField ( )

Returns a pointer to the e-Field.

Solves the problem, if not yet.

std::string hp2D::Models< Cmplx >::filename_	(	const std::string	varname,
		const std::string	subpath,
		const bool	p = `false`
	)		`[protected, inherited]`

Gives a file name, which consists of p, l, g.

void hp2D::Models< Cmplx >::formulaOut_	(	const concepts::PiecewiseFormulaBase< G > &	frm,
		const std::string	frmstr,
		const std::string	path = `""`
	)		`[protected, inherited]`

graphical output of formulas

virtual void hp2D::Maxwell2D_H::graphicsOut	(	const uint	numPoints,
		const std::string	path = `""`
	)		`[virtual]`

Writes matlab graphics files.

Parameters:

numPoints number of points per cell in one direction

Implements hp2D::Models< Cmplx >.

concepts::ElementFormula<Cmplx>* hp2D::Maxwell2D_H::hField ( )

Returns a pointer to the h-Field (solution + h0).

Solves the problem, if not yet.

void hp2D::Maxwell2D_H::identityMatrix_ ( ) [private]

Calculate identity matrix.

concepts::SparseMatrix<Real>* hp2D::Maxwell2D_H_Base::identityMatrix_	(	concepts::Space< Real > &	spc,
		concepts::SparseMatrix< Cmplx > *	S
	)		`[protected, inherited]`

Calculate identity matrix and add's it to system matrix S.

virtual std::ostream& hp2D::Maxwell2D_H::info ( std::ostream & os ) const [protected, virtual]

Returns information in an output stream.

Reimplemented from hp2D::Maxwell2D_H_Base.

virtual std::ostream& hp2D::Maxwell2D_H::info ( std::ostream & os ) const [protected, virtual]

Returns information in an output stream.

Reimplemented from hp2D::Maxwell2D_H_Base.

void hp2D::Maxwell2D_H::laplaceMatrix_ ( ) [private]

Calculate stiffness matrix.

concepts::SparseMatrix<Cmplx>* hp2D::Maxwell2D_H_Base::laplaceMatrix_	(	concepts::Space< Real > &	spc,
		concepts::SparseMatrix< Cmplx > *	S
	)		`[protected, inherited]`

Calculate stiffness matrix and add's it to system matrix S.

void hp2D::Maxwell2D_H::linearform_ ( ) [private]

Calculate the load vector, assumes sigma = 0 inside the coil.

Reimplemented from hp2D::Maxwell2D_H_Base.

virtual void hp2D::Models< Cmplx >::loadSolution ( const std::string vectorFile ) [virtual, inherited]

Loads solution vector from vectorFile.

virtual void hp2D::Models< Cmplx >::loadSolution	(	const uint	l = `0`,
		const uint	p = `1`,
		const uint	g = `0`,
		const std::string	path = `""`
	)		`[virtual, inherited]`

Loads solution vector from standard file in path with use of l, p and g.

Real hp2D::Maxwell2D_H::magnEnergy ( ) [virtual]

Return magnetic energy.

Implements concepts::MaxwellModel.

virtual Real hp2D::Maxwell2D_H::magnEnergy ( ) [virtual]

Return magnetic energy.

Implements concepts::MaxwellModel.

virtual void hp2D::Maxwell2D_H::matrices_ ( ) [private, virtual]

Building the matrices.

Implements hp2D::Models< Cmplx >.

void hp2D::Maxwell2D_H::matrices_ ( ) [private, virtual]

Building the matrices.

Implements hp2D::Models< Cmplx >.

void hp2D::Models< Cmplx >::matrixBuild ( InputModels< Cmplx > & input ) [inherited]

void hp2D::Models< Cmplx >::matrixBuild	(	const uint	l = `0`,
		const uint	p = `1`,
		const uint	g = `0`,
		const std::string	matrixFile = `""`
	)		`[inherited]`

Build the space, matrices and the load vector.

virtual concepts::Mesh2& hp2D::Maxwell2D_H::mesh_ ( ) [inline, protected, virtual]

Mesh.

Implements hp2D::Models< Cmplx >.

Definition at line 61 of file Maxwell2D_H.hh.

Here is the call graph for this function:

virtual const std::string hp2D::Maxwell2D_H::mshAbbr_ ( ) [inline, protected, virtual]

Mesh abbreviation string.

Reimplemented from hp2D::Maxwell2D_H_Base.

Definition at line 63 of file Maxwell2D_H.hh.

Here is the call graph for this function:

virtual hpFull& hp2D::Maxwell2D_H::prebuild_ ( ) [inline, protected, virtual]

Space Prebuilder.

Implements hp2D::AdaptiveModel< Cmplx >.

Definition at line 162 of file Maxwell2D_H.hh.

template<class F, uint dim>

void concepts::AdaptiveModel< F, dim >::rebuildMesh ( concepts::InputAdaptiveModels & input ) [inherited]

template<class F, uint dim>

void concepts::AdaptiveModel< F, dim >::rebuildMesh	(	const uint	l = `0`,
		const uint	p = `1`,
		const uint	g = `0`,
		const uint	subdiv = `X\|Y`
	)		`[inherited]`

Rebuilds only the mesh and sets the polynomial degrees.

Parameters:

l	number of uniform refinements
p	number of polynomial enlargements
g	number of geometric refinements
subdiv	possibility to restrict subdivision strategy for geometric refinement

virtual void hp2D::AdaptiveModel< Cmplx >::rebuildMesh_	(	const uint	l = `0`,
		const uint	p = `1`,
		const uint	g = `0`,
		const uint	subdiv = `X \| Y`
	)		`[protected, virtual, inherited]`

Rebuilds only the mesh and sets the polynomial degrees.

Implements concepts::AdaptiveModel< F, dim >.

template<class F>

const Vector<F>* concepts::Model< F >::solution ( ) [inline, inherited]

Returns solution vector.

Definition at line 66 of file models.hh.

void hp2D::Models< Cmplx >::solve	(	const uint	l = `0`,
		const uint	p = `1`,
		const uint	g = `0`,
		const std::string	vectorFile = `""`,
		const std::string	matrixFile = `""`
	)		`[inherited]`