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

#include <Eddy2D_H.hh>

Inheritance diagram for hp2D::Eddy2D_H:

Collaboration diagram for hp2D::Eddy2D_H:

Public Types
enum	solverType { SUPERLU = 0, BICGSTAB = 1 }
	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.
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.
virtual Real	dissipation ()
	Return dissipation energy.
Real	dissipation ()
	Return dissipation energy.
	Eddy2D_H (concepts::EddyGeometry2D &geom, InputEddy2D_H &input, const uint geomRefAttrib=100)
	Eddy2D_H (concepts::EddyGeometry2D &geom, const concepts::Formula< Real > &H0, const concepts::Formula< Real2d > &curlH0, const concepts::Formula< Real > divgradH0=0, Eddy2D_H_Interior interior=0, const uint geomRefAttrib=100, const Real omega=OMEGA50, const Real mu=MU0, enum solverType type=SUPERLU)
	Constructor.
	Eddy2D_H (concepts::EddyGeometry2D &geom, const Real H0, const Real Omega_i=0, const uint geomRefAttr=100, const Real omega=OMEGA50, const Real mu=MU0)
	Eddy2D_H (concepts::EddyGeometry2D &geom, const Real H0, Eddy2D_H_Interior *interior=0, const uint geomRefAttrib=100, const Real omega=OMEGA50, const Real mu=MU0, enum solverType type=SUPERLU)
	Constructor.
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 (const Real R1, const Real R2)
	Return magnetic energy, when assuming ring coil with constant current.
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)
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.
void	rebuildMesh (concepts::InputAdaptiveModels &input)
const Vector< F > *	solution ()
	Returns solution vector.
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.
void	solve (InputModels &input)
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	~Eddy2D_H ()
virtual	~Eddy2D_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
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.
virtual concepts::Mesh2 &	mesh_ ()
	Mesh.
virtual const std::string	mshAbbr_ ()
	Mesh abbreviation string.
virtual const std::string	mshAbbr_ ()
	Mesh abbreviation string.
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
const uint	geomRefAttr_
	Attribute of vertices or edges for geometric refinement.
const uint	geomRefAttr_
	Attribute of vertices or edges for geometric refinement.
std::string	matrixFile_
	matlab file for matrices and load vector output
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.
std::string	vectorFile_
	file for solution output, which can be reloaded into concepts
Private Member Functions
bool	connectedIdx_ (uint &i)
	Gives connected index at interior border to i, if exists.
bool	connectedIdx_ (uint &i)
	Gives connected index at interior border to i, if exists.
void	constructor_ ()
	Private constructor.
virtual void	dataOut_ ()
	Collect the parameters.
void	identityMatrix_ ()
	Calculate identity matrix.
void	identityMatrix_ ()
	Calculate identity matrix.
void	laplaceMatrix_ ()
	Calculate stiffness matrix.
void	laplaceMatrix_ ()
	Calculate stiffness matrix.
void	linearform_ ()
	Calculate the load vector.
void	linearform_ ()
	Calculate the load vector.
virtual void	matrices_ ()
	Building the matrices.
void	matrices_ ()
	Building the matrices.
virtual hpFull &	prebuild_ ()
	Space Prebuilder.
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.
void	systemMatrix_ ()
	Compress system matrix.
Private Attributes
std::auto_ptr < concepts::SparseMatrix< Real > >	A_
	Stiffness matrix.
std::auto_ptr < concepts::BoundaryConditions >	bc_
	Boundary conditions.
std::auto_ptr < concepts::CellConditions >	cc_
	Cell conditions.
std::auto_ptr< const concepts::PiecewiseFormulaBase < Real2d > >	curlH0_
std::auto_ptr< Real >	dissipation_
	dissipation energy
std::auto_ptr< const concepts::PiecewiseFormulaBase < Real > >	divgradH0_
std::auto_ptr< const concepts::ElementFunction < Cmplx > >	fun_
	Element function, is hp2D::Value + homogenisation.
concepts::EddyGeometry2D &	geom_
	Mesh and material constants (sigma)
std::auto_ptr< const concepts::PiecewiseFormulaBase < Real > >	H0_
	Formula for homogenization.
const Real	H0_
	$h_0$ from coil
std::auto_ptr< Eddy2D_H_Interior >	interior_
	Size and magnetic field of or in the nonconducting imbedded region.
std::auto_ptr < concepts::SparseMatrix< Real > >	M_
	Mass matrix.
std::auto_ptr< Real >	magnEnergy_
	Magnetic energy.
double	matrixtime_
const Real	mu_
	Permeability constant.
const Real	omega_
	Angular frequency.
const Real	Omega_i_
	size of the nonconducting imbedded region
std::auto_ptr < concepts::Vector< Cmplx > >	residual_
	Residual vector.
std::auto_ptr< Real >	residualNorm_
	Euclidian norm of the residual of solving the linear system.
std::auto_ptr < concepts::Vector< Cmplx > >	rhs_
std::auto_ptr < concepts::SparseMatrix < Cmplx > >	S_
	System matrix.
concepts::PiecewiseFormulaFun < Real, Real >	Sigma_Inv_
	Piecewise constant formula $\sigma^{-1}$ .
double	solvetime_
	Time to solve the system, to build the matrices, to rebuild the space.
double	spacetime_
std::auto_ptr< hpAdaptiveSpaceH1 >	spc_
	Space.
enum solverType	type_
	Solver type.

Detailed Description

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

Author:: Kersten Schmidt, 2005

Definition at line 31 of file Eddy2D_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 hp2D::Eddy2D_H::solverType

Type of the solver.

Enumerator:

SUPERLU
BICGSTAB

Definition at line 69 of file Eddy2D_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::Eddy2D_H::Eddy2D_H	(	concepts::EddyGeometry2D &	geom,
		const Real	H0,
		const Real	Omega_i = `0`,
		const uint	geomRefAttr = `100`,
		const Real	omega = `OMEGA50`,
		const Real	mu = `MU0`
	)

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

Definition at line 36 of file Eddy2D_H.hh.

hp2D::Eddy2D_H::Eddy2D_H	(	concepts::EddyGeometry2D &	geom,
		const concepts::Formula< Real > &	H0,
		const concepts::Formula< Real2d > &	curlH0,
		const concepts::Formula< Real > *	divgradH0 = `0`,
		Eddy2D_H_Interior *	interior = `0`,
		const uint	geomRefAttrib = `100`,
		const Real	omega = `OMEGA50`,
		const Real	mu = `MU0`,
		enum solverType	type = `SUPERLU`
	)

Constructor.

Parameters:

geom	geometry, conductivity and source currents
H0	Formula for homogenization.
curlH0	\[{curl}_{2D} h_0\]
divgradH0	pointer to formula \[{div} {grad} h_0\], only needed if offset function has gradient inside the conductive region
omega	angular frequency in 1/s
mu	permeability constant in Ohm*s/m
geomRefAttrib	attrib for geometric refinement
type	solver type

hp2D::Eddy2D_H::Eddy2D_H	(	concepts::EddyGeometry2D &	geom,
		const Real	H0,
		Eddy2D_H_Interior *	interior = `0`,
		const uint	geomRefAttrib = `100`,
		const Real	omega = `OMEGA50`,
		const Real	mu = `MU0`,
		enum solverType	type = `SUPERLU`
	)

Constructor.

Parameters:

geom	geometry, conductivity and source currents
H0	magnetic field, created by the coil
omega	angular frequency in 1/s
mu	permeability constant in Ohm*s/m
geomRefAttrib	attrib for geometric refinement
type	solver type

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

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

Definition at line 102 of file Eddy2D_H.hh.

Member Function Documentation

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

Writes matlab graphics files for each basis function.

bool hp2D::Eddy2D_H::connectedIdx_ ( uint & i ) [private]

Gives connected index at interior border to i, if exists.

bool hp2D::Eddy2D_H::connectedIdx_ ( uint & i ) [private]

Gives connected index at interior border to i, if exists.

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

Private constructor.

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::Eddy2D_H::dataOut_ ( ) [private, virtual]

Collect the parameters.

Implements hp2D::Models< Cmplx >.

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

Return dissipation energy.

Implements concepts::MaxwellModel.

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

Return dissipation energy.

Implements concepts::MaxwellModel.

concepts::ElementFormula<concepts::Cmplx2d>* hp2D::Eddy2D_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::Eddy2D_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::Eddy2D_H::hField ( )

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

Solves the problem, if not yet.

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

Calculate identity matrix.

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

Calculate identity matrix.

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

Returns information in an output stream.

Reimplemented from hp2D::Models< Cmplx >.

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

Returns information in an output stream.

Reimplemented from hp2D::Models< Cmplx >.

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

Calculate stiffness matrix.

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

Calculate stiffness matrix.

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

Calculate the load vector.

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

Calculate the load vector.

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.

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

Return magnetic energy.

Implements concepts::MaxwellModel.

Real hp2D::Eddy2D_H::magnEnergy	(	const Real	R1,
		const Real	R2
	)

Return magnetic energy, when assuming ring coil with constant current.

Parameters:

R1	inner radius of coil
R1	outer radius of coil

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

Building the matrices.

Implements hp2D::Models< Cmplx >.

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

Building the matrices.

Implements hp2D::Models< Cmplx >.

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.

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

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

Mesh.

Implements hp2D::Models< Cmplx >.

Definition at line 54 of file Eddy2D_H.hh.

Here is the call graph for this function:

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

Mesh abbreviation string.

Implements hp2D::Models< Cmplx >.

Definition at line 120 of file Eddy2D_H.hh.

Here is the call graph for this function:

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

Mesh abbreviation string.

Implements hp2D::Models< Cmplx >.

Definition at line 56 of file Eddy2D_H.hh.

Here is the call graph for this function:

virtual hpFull& hp2D::Eddy2D_H::prebuild_ ( ) [inline, private, virtual]

Space Prebuilder.

Implements hp2D::AdaptiveModel< Cmplx >.

Definition at line 125 of file Eddy2D_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 ( InputModels< Cmplx > & input ) [inherited]

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]`