Front-end for the grid manager of the finite element toolbox UG. More...

#include <dune/grid/uggrid.hh>

Public Types
typedef UGGridFamily< dim >	GridFamily
	type of the used GridFamily for this grid

typedef UG::DOUBLE	ctype
	The type used to store coordinates.

Exported constants
enum
	A constant that exports the template parameter dim.

enum
	A constant that exports the template parameter dimworld.

Exported types
typedef Partition< All_Partition >::LevelGridView	LevelGridView
	View types for All_Partition.

typedef Partition< All_Partition >::LeafGridView	LeafGridView

typedef GridFamily::Traits::LeafIntersection	LeafIntersection
	A type that is a model of Dune::Intersection, an intersections of two codimension 1 of two codimension 0 entities in the leaf view.

typedef GridFamily::Traits::LevelIntersection	LevelIntersection
	A type that is a model of Dune::Intersection, an intersections of two codimension 1 of two codimension 0 entities in a level view.

typedef GridFamily::Traits::LeafIntersectionIterator	LeafIntersectionIterator
	A type that is a model of Dune::IntersectionIterator which is an iterator that allows to examine, but not to modify, the intersections of codimension 1 of an leaf element (entity of codimension 0) with other leaf elements.

typedef GridFamily::Traits::LevelIntersectionIterator	LevelIntersectionIterator
	A type that is a model of Dune::IntersectionIterator which is an iterator that allows to examine, but not to modify, the intersections of codimension 1 of an element (entity of codimension 0) with other elements on the same level.

typedef GridFamily::Traits::HierarchicIterator	HierarchicIterator
	A type that is a model of Dune::HierarchicIterator A type of iterator that allows to examine, but not to modify, entities of codimension 0 that result from refinement of an entity of codimension 0.

typedef GridFamily::Traits::LevelIndexSet	LevelIndexSet
	A type that is a model of Dune::IndexSet which provides a consecutive, but non persistent, numbering for entities on a grid level.

typedef GridFamily::Traits::LeafIndexSet	LeafIndexSet
	A type that is a model of Dune::IndexSet which provides a consecutive, but non persistent, numbering for entities in the leaf grid.

typedef GridFamily::Traits::GlobalIdSet	GlobalIdSet
	A type that is a model of Dune::IdSet which provides a unique and persistent numbering for all entities in the grid. The numbering is unique over all processes over which the grid is partitioned. The numbering is not necessarily consecutive.

typedef GridFamily::Traits::LocalIdSet	LocalIdSet
	A type that is a model of Dune::IdSet which provides a unique and persistent numbering for all entities in the grid. The numbering is only unique in a single process and it is not necessarily consecutive.

typedef GridFamily::Traits::CollectiveCommunication	CollectiveCommunication
	A type that is a model of Dune::CollectiveCommunication. It provides a portable way for collective communication on the set of processes used by the grid.

Public Member Functions
	UGGrid ()
	Default constructor. More...

	~UGGrid ()
	Destructor.

int	maxLevel () const

template<int codim>
Traits::template Codim< codim >::LevelIterator	lbegin (int level) const
	Iterator to first entity of given codim on level.

template<int codim>
Traits::template Codim< codim >::LevelIterator	lend (int level) const
	one past the end on this level

template<int codim, PartitionIteratorType PiType>
Traits::template Codim< codim >::template Partition< PiType >::LevelIterator	lbegin (int level) const
	Iterator to first entity of given codim on level.

template<int codim, PartitionIteratorType PiType>
Traits::template Codim< codim >::template Partition< PiType >::LevelIterator	lend (int level) const
	one past the end on this level

template<int codim>
Traits::template Codim< codim >::LeafIterator	leafbegin () const
	Iterator to first leaf entity of given codim.

template<int codim>
Traits::template Codim< codim >::LeafIterator	leafend () const
	one past the end of the sequence of leaf entities

template<int codim, PartitionIteratorType PiType>
Traits::template Codim< codim >::template Partition< PiType >::LeafIterator	leafbegin () const
	Iterator to first leaf entity of given codim.

template<int codim, PartitionIteratorType PiType>
Traits::template Codim< codim >::template Partition< PiType >::LeafIterator	leafend () const
	one past the end of the sequence of leaf entities

template<typename Seed >
Traits::template Codim< Seed::codimension >::EntityPointer	entityPointer (const Seed &seed) const
	Create an EntityPointer from an EntitySeed.

int	size (int level, int codim) const
	Number of grid entities per level and codim.

int	size (int codim) const
	number of leaf entities per codim in this process

int	size (int level, GeometryType type) const
	number of entities per level and geometry type in this process

int	size (GeometryType type) const
	number of leaf entities per geometry type in this process

size_t	numBoundarySegments () const
	Return the number of boundary segments.

const Traits::GlobalIdSet &	globalIdSet () const
	Access to the GlobalIdSet.

const Traits::LocalIdSet &	localIdSet () const
	Access to the LocalIdSet.

const Traits::LevelIndexSet &	levelIndexSet (int level) const
	Access to the LevelIndexSets.

const Traits::LeafIndexSet &	leafIndexSet () const
	Access to the LeafIndexSet.

Traits::template Partition< pitype >::LevelGridView	DUNE_DEPRECATED_MSG ("The method levelView has been renamed to levelGridView.") levelView(int level) const
	View for a grid level.

Traits::template Partition< pitype >::LeafGridView	DUNE_DEPRECATED_MSG ("The method levelView has been renamed to leafGridView.") leafView() const
	View for the leaf grid.

Traits::template Partition< All_Partition >::LevelGridView	levelView (int level) const DUNE_DEPRECATED_MSG("The method levelView has been renamed to levelGridView.")
	View for a grid level for All_Partition.

Traits::template Partition< All_Partition >::LeafGridView	leafView () const DUNE_DEPRECATED_MSG("The method leafView has been renamed to leafGridView.")
	View for the leaf grid for All_Partition.

Traits::template Partition< pitype >::LevelGridView	levelGridView (int level) const
	View for a grid level.

Traits::template Partition< All_Partition >::LevelGridView	levelGridView (int level) const
	View for a grid level for All_Partition.

Traits::template Partition< pitype >::LeafGridView	leafGridView () const
	View for the leaf grid.

Traits::template Partition< All_Partition >::LeafGridView	leafGridView () const
	View for the leaf grid for All_Partition.

bool	mark (int refCount, const typename Traits ::template Codim< 0 >::Entity &e)
	Marks an entity to be refined/coarsened in a subsequent adapt. More...

void	communicate (CommDataHandleIF< DataHandleImp, DataTypeImp > &data, InterfaceType iftype, CommunicationDirection dir, int level) const

void	communicate (CommDataHandleIF< DataHandleImp, DataTypeImp > &data, InterfaceType iftype, CommunicationDirection dir) const

Adaptivity and grid refinement
bool	mark (int refCount, const typename Codim< 0 >::Entity &e)
	Marks an entity to be refined/coarsened in a subsequent adapt. More...

int	getMark (const typename Codim< 0 >::Entity &e) const
	returns adaptation mark for given entity More...

Static Protected Member Functions
static ReturnImplementationType< InterfaceType >::ImplementationType &	getRealImplementation (InterfaceType &i)
	return real implementation of interface class

Grid Refinement Methods
enum	RefinementType { LOCAL , COPY }
	The different forms of grid refinement that UG supports. More...

enum	ClosureType { GREEN , NONE }
	Decide whether to add a green closure to locally refined grid sections or not. More...

bool	mark (int refCount, const typename Traits::template Codim< 0 >::Entity &e)
	Mark element for refinement. More...

bool	mark (const typename Traits::template Codim< 0 >::Entity &e, typename UG_NS< dim >::RefinementRule rule, int side=0)
	Mark method accepting a UG refinement rule. More...

int	getMark (const typename Traits::template Codim< 0 >::Entity &e) const
	Query whether element is marked for refinement.

bool	preAdapt ()
	returns true, if some elements might be coarsend during grid adaption, here always returns true

bool	adapt ()
	Triggers the grid refinement process.

void	postAdapt ()
	Clean up refinement markers.

unsigned int	overlapSize (int codim) const
	Size of the overlap on the leaf level.

unsigned int	ghostSize (int codim) const
	Size of the ghost cell layer on the leaf level.

unsigned int	overlapSize (int level, int codim) const
	Size of the overlap on a given level.

unsigned int	ghostSize (int level, int codim) const
	Size of the ghost cell layer on a given level.

bool	loadBalance ()
	Default load balancing. More...

template<class DataHandle >
bool	loadBalance (DataHandle &dataHandle)
	Distributes the grid and some data over the available nodes in a distributed machine. More...

bool	loadBalance (int strategy, int minlevel, int depth, int maxlevel, int minelement)
	Distributes this grid over the available nodes in a distributed machine. More...

template<class DataHandle >
void	communicate (DataHandle &dataHandle, InterfaceType iftype, CommunicationDirection dir, int level) const
	The communication interface for all codims on a given level. More...

template<class DataHandle >
void	communicate (DataHandle &dataHandle, InterfaceType iftype, CommunicationDirection dir) const
	The communication interface for all codims on the leaf level. More...

const CollectiveCommunication< UGGrid > &	comm () const

void	getChildrenOfSubface (const typename Traits::template Codim< 0 >::EntityPointer &e, int elementSide, int maxl, std::vector< typename Traits::template Codim< 0 >::EntityPointer > &childElements, std::vector< unsigned char > &childElementSides) const
	Rudimentary substitute for a hierarchic iterator on faces. More...

void	setRefinementType (RefinementType type)
	Sets the type of grid refinement.

void	setClosureType (ClosureType type)
	Sets the type of grid refinement closure.

void	setPosition (const typename Traits::template Codim< dim >::EntityPointer &e, const FieldVector< double, dim > &pos)
	Sets a vertex to a new position. More...

void	globalRefine (int n)
	Does uniform refinement. More...

void	saveState (const std::string &filename) const
	Save entire grid hierarchy to disk. More...

void	loadState (const std::string &filename)
	Read entire grid hierarchy from disk. More...

static void	setDefaultHeapSize (unsigned size)
	Sets the default heap size. More...

Parallel data distribution and communication
Codim< EntitySeed::codimension >::EntityPointer	entityPointer (const EntitySeed &seed) const
	obtain EntityPointer from EntitySeed.

GridImp &	asImp ()
	Barton-Nackman trick.

const GridImp &	asImp () const
	Barton-Nackman trick.

Detailed Description

template<int dim>
class Dune::UGGrid< dim >

Front-end for the grid manager of the finite element toolbox UG.

This is the implementation of the grid interface using the UG grid management system (http://www.iwr.uni-heidelberg.de/frame/iwrwikiequipment/software/ug). UG provides conforming grids in two and three space dimensions. The grids can be mixed, i.e. 2d grids can contain triangles and quadrilaterals and 3d grids can contain tetrahedra and hexahedra and also pyramids and prisms. The grid refinement rules are very flexible. Local adaptive red/green refinement is the default, but a special method in the UGGrid class allows you to directly access a number of anisotropic refinements rules. Last but not least, the UG grid manager is completely parallelized, and you can use boundaries parametrized by either analytical expressions or high-resolution piecewise linear surfaces.

To use this module you need the UG library. See the DUNE installation notes on how to obtain and install it.

In your Dune application, you can now instantiate objects of the type UGGrid<2> or UGGrid<3>. You can have more than one, if you choose. It is even possible to have 2d and 3d grids at the same time, even though the original UG system never intended to support this!

See the documentation for the factory class GridFactory<UGGrid<dimworld> > to learn how to create UGGrid objects.

Please send any questions, suggestions, or bug reports to the Dune mailing list dune@.nosp@m.dune.nosp@m.-proj.nosp@m.ect..nosp@m.org

For installation instructions see http://www.dune-project.org/external_libraries/install_ug.html .

Member Enumeration Documentation

◆ ClosureType

template<int dim>

enum Dune::UGGrid::ClosureType

Decide whether to add a green closure to locally refined grid sections or not.

Enumerator
GREEN	Standard red/green refinement.
NONE	No closure, results in nonconforming meshes.

◆ RefinementType

template<int dim>

enum Dune::UGGrid::RefinementType

The different forms of grid refinement that UG supports.

Enumerator
LOCAL	New level consists only of the refined elements and the closure.
COPY	New level consists of the refined elements and the unrefined ones, too.

Constructor & Destructor Documentation

◆ UGGrid()

template<int dim>

Dune::UGGrid< dim >::UGGrid ( )

Default constructor.

Uses the default heapsize, which can be set using the static method setDefaultHeapSize().

Member Function Documentation

◆ comm()

template<int dim>

const CollectiveCommunication< UGGrid > & Dune::UGGrid< dim >::comm ( ) const

inline

the collective communication

◆ communicate() [1/4]

void Dune::GridDefaultImplementation< dim, dimworld, double , UGGridFamily< dim > >::communicate	(	CommDataHandleIF< DataHandleImp, DataTypeImp > &	data,
		InterfaceType	iftype,
		CommunicationDirection	dir
	)		const

inlineinherited

dummy communicate, doing nothing

◆ communicate() [2/4]

void Dune::GridDefaultImplementation< dim, dimworld, double , UGGridFamily< dim > >::communicate	(	CommDataHandleIF< DataHandleImp, DataTypeImp > &	data,
		InterfaceType	iftype,
		CommunicationDirection	dir,
		int	level
	)		const

inlineinherited

dummy communicate, doing nothing

◆ communicate() [3/4]

template<int dim>

template<class DataHandle >

void Dune::UGGrid< dim >::communicate	(	DataHandle &	dataHandle,
		InterfaceType	iftype,
		CommunicationDirection	dir
	)		const

inline

The communication interface for all codims on the leaf level.

Parameters

dataHandle	type used to gather/scatter data in and out of the message buffer
iftype	one of the predifined interface types, throws error if it is not implemented
dir	choose beetween forward and backward communication

Implements a generic communication function sending an object of type P for each entity in the intersection of two processors. P has two methods gather and scatter that implement the protocol. Therefore P is called the "protocol class".

References Dune::className(), DUNE_THROW, and Dune::GridDefaultImplementation< dim, dim, double, UGGridFamily< dim > >::leafView().

◆ communicate() [4/4]

template<int dim>

template<class DataHandle >

void Dune::UGGrid< dim >::communicate	(	DataHandle &	dataHandle,
		InterfaceType	iftype,
		CommunicationDirection	dir,
		int	level
	)		const

inline

The communication interface for all codims on a given level.

Parameters

dataHandle	type used to gather/scatter data in and out of the message buffer
iftype	one of the predifined interface types, throws error if it is not implemented
dir	choose beetween forward and backward communication
level	communicate for entities on the given level

Implements a generic communication function sending an object of type P for each entity in the intersection of two processors. P has two methods gather and scatter that implement the protocol. Therefore P is called the "protocol class".

References Dune::className(), DUNE_THROW, and Dune::GridDefaultImplementation< dim, dim, double, UGGridFamily< dim > >::levelGridView().

◆ getChildrenOfSubface()

template<int dim>

void Dune::UGGrid< dim >::getChildrenOfSubface	(	const typename Traits::template Codim< 0 >::EntityPointer &	e,
		int	elementSide,
		int	maxl,
		std::vector< typename Traits::template Codim< 0 >::EntityPointer > &	childElements,
		std::vector< unsigned char > &	childElementSides
	)		const

Rudimentary substitute for a hierarchic iterator on faces.

Parameters

	e,elementSide	Grid face specified by an element and one of its sides
	maxl	The finest level that should be traversed by the iterator
[out]	childElements	For each subface: element index, elementSide, and level
[out]	childElementSides	Indices for transformation because Dune numbers the faces of several elements differently than UG

◆ getMark()

int Dune::Grid< dim, dimworld, double , UGGridFamily< dim > >::getMark ( const typename Codim< 0 >::Entity & e ) const

inlineinherited

returns adaptation mark for given entity

Parameters

[in] e Entity for which adaptation mark should be determined

Returns: int adaptation mark currently set for given Entity e

◆ globalRefine()

template<int dim>

void Dune::UGGrid< dim >::globalRefine ( int n )

Does uniform refinement.

Parameters

n	Number of uniform refinement steps

◆ loadBalance() [1/3]

template<int dim>

bool Dune::UGGrid< dim >::loadBalance ( )

inline

Default load balancing.

Bug:: The return value is always 'true'

Returns: true if the grid has changed

References Dune::UGGrid< dim >::loadBalance().

Referenced by Dune::UGGrid< dim >::loadBalance().

◆ loadBalance() [2/3]

template<int dim>

template<class DataHandle >

bool Dune::UGGrid< dim >::loadBalance ( DataHandle & dataHandle )

inline

Distributes the grid and some data over the available nodes in a distributed machine.

Template Parameters

DataHandle works like the data handle for the communicate methods.

Returns: True, if grid has changed, false otherwise

References DUNE_THROW, Dune::GridDefaultImplementation< dim, dim, double, UGGridFamily< dim > >::leafView(), and Dune::UGGrid< dim >::loadBalance().

◆ loadBalance() [3/3]

template<int dim>

bool Dune::UGGrid< dim >::loadBalance	(	int	strategy,
		int	minlevel,
		int	depth,
		int	maxlevel,
		int	minelement
	)

Distributes this grid over the available nodes in a distributed machine.

If you want the UG default for the parameters pick

strategy = 0
minlevel = 1
depth = 2
maxlevel = 32
minelement = 1

Bug:: The return value is always 'true'

Parameters

minlevel	The coarsest grid level that gets distributed
maxlevel	does currently get ignored

◆ loadState()

template<int dim>

void Dune::UGGrid< dim >::loadState ( const std::string & filename )

Read entire grid hierarchy from disk.

Test implementation – not working!

◆ mark() [1/4]

template<int dim>

bool Dune::UGGrid< dim >::mark	(	const typename Traits::template Codim< 0 >::Entity &	e,
		typename UG_NS< dim >::RefinementRule	rule,
		int	side = `0`
	)

Mark method accepting a UG refinement rule.

Parameters

e	Pointer to the element to be marked for refinement
rule	One of the UG refinement rules
side	If rule==UG::D2::BLUE (one quadrilateral is split into two rectangles) you can choose the orientation of the cut by setting side==0 or side==1

The available values for RefinementRule are: (see the RefinementRule enum in ug/gm/gm.h)

2D

NO_REFINEMENT
COPY
RED
BLUE
COARSE
BISECTION_1
BISECTION_2_Q
BISECTION_2_T1
BISECTION_2_T2
BISECTION_3

3D

NO_REFINEMENT
COPY
RED
BLUE
COARSE
TETRA_RED_HEX
PRISM_BISECT_1_2
PRISM_QUADSECT
PRISM_BISECT_HEX0
PRISM_BISECT_HEX1
PRISM_BISECT_HEX2
PRISM_ROTATE_LEFT
PRISM_ROTATE_RGHT
PRISM_QUADSECT_HEXPRI0
PRISM_RED_HEX
PRISM_BISECT_0_1
PRISM_BISECT_0_2
PRISM_BISECT_0_3
HEX_BISECT_0_1
HEX_BISECT_0_2
HEX_BISECT_0_3
HEX_TRISECT_0
HEX_TRISECT_5
HEX_QUADSECT_0
HEX_QUADSECT_1
HEX_QUADSECT_2
HEX_BISECT_HEXPRI0
HEX_BISECT_HEXPRI1

◆ mark() [2/4]

bool Dune::Grid< dim, dimworld, double , UGGridFamily< dim > >::mark	(	int	refCount,
		const typename Codim< 0 >::Entity &	e
	)

inlineinherited

Marks an entity to be refined/coarsened in a subsequent adapt.

Parameters

[in]	refCount	Number of subdivisions that should be applied. Negative value means coarsening.
[in]	e	Entity that should be marked

Returns: true if Entity was marked, false otherwise.

◆ mark() [3/4]

bool Dune::GridDefaultImplementation< dim, dimworld, double , UGGridFamily< dim > >::mark	(	int	refCount,
		const typename Traits ::template Codim< 0 >::Entity &	e
	)

inlineinherited

Marks an entity to be refined/coarsened in a subsequent adapt.

Parameters

[in]	refCount	Number of subdivisions that should be applied. Negative value means coarsening.
[in]	e	Entity to Entity that should be refined

Returns: true if Entity was marked, false otherwise.

Note

default implementation is: return false; for grids with no adaptation.
for the grid programmer: this method is implemented as a template method, because the Entity type is not defined when the class is instantiated You won't need this trick in the implementation. In your implementation you should use it as
bool mark( int refCount,

typename Traits::template Codim<0>::Entity & e ).

Dune::UGGrid::mark
bool mark(int refCount, const typename Traits::template Codim< 0 >::Entity &e)
Mark element for refinement.

This template method will vanish due to the inheritance rules.

◆ mark() [4/4]

template<int dim>

bool Dune::UGGrid< dim >::mark	(	int	refCount,
		const typename Traits::template Codim< 0 >::Entity &	e
	)

Mark element for refinement.

Parameters

refCount	1: mark for red refinement -1: mark for coarsening 0: delete a possible refinement mark
e	Element to be marked

Returns

true, if element was marked
false, if nothing changed

◆ maxLevel()

template<int dim>

int Dune::UGGrid< dim >::maxLevel ( ) const

Return maximum level defined in this grid. Levels are numbered 0 ... maxlevel with 0 the coarsest level.

Referenced by Dune::UGGrid< dim >::levelIndexSet().

◆ saveState()

template<int dim>

void Dune::UGGrid< dim >::saveState ( const std::string & filename ) const

Save entire grid hierarchy to disk.

Test implementation – not working!

◆ setDefaultHeapSize()

template<int dim>

static void Dune::UGGrid< dim >::setDefaultHeapSize ( unsigned size )

inlinestatic

Sets the default heap size.

UGGrid keeps an internal heap to allocate memory from, which must be specified on grid creation (at the latest). This sets the default heap size, which is used when no heap size is given to the constructor.

References Dune::UGGrid< dim >::size().

◆ setPosition()

template<int dim>

void Dune::UGGrid< dim >::setPosition	(	const typename Traits::template Codim< dim >::EntityPointer &	e,
		const FieldVector< double, dim > &	pos
	)

Sets a vertex to a new position.

Changing a vertex' position changes its position on all grid levels!

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

dune/grid/uggrid.hh

DUNE Distributed and Unified Numerics Environment

Dune Core Modules (2.3.1)

Public Types

Public Member Functions

Static Protected Member Functions

Grid Refinement Methods

Parallel data distribution and communication

Detailed Description

Member Enumeration Documentation

◆ ClosureType

◆ RefinementType

Constructor & Destructor Documentation

◆ UGGrid()

Member Function Documentation

◆ comm()

◆ communicate() [1/4]

◆ communicate() [2/4]

◆ communicate() [3/4]

◆ communicate() [4/4]

◆ getChildrenOfSubface()

◆ getMark()

◆ globalRefine()

◆ loadBalance() [1/3]

◆ loadBalance() [2/3]

◆ loadBalance() [3/3]

◆ loadState()

◆ mark() [1/4]

2D

3D

◆ mark() [2/4]

◆ mark() [3/4]

◆ mark() [4/4]

◆ maxLevel()

◆ saveState()

◆ setDefaultHeapSize()

◆ setPosition()