Merge two codimension-1 surfaces that may be a positive distance apart. More...

#include <dune/grid-glue/merging/contactmerge.hh>

Inheritance diagram for Dune::GridGlue::ContactMerge< dimworld, T >:

Public Types
enum	ProjectionType { OUTER_NORMAL, CLOSEST_POINT }
	Type of the projection, closest point or outer normal projection. More...

typedef T	ctype
	the numeric type used in this interface More...

typedef Dune::FieldVector< T, dimworld >	WorldCoords
	the coordinate type used in this interface More...

typedef Dune::FieldVector< T, dim >	LocalCoords
	the coordinate type used in this interface More...

typedef Merger< T, grid1Dim, grid2Dim, dimworld >::Grid1Coords	Grid1Coords
	Type used for local coordinates on the grid1 side. More...

typedef Merger< T, grid1Dim, grid2Dim, dimworld >::Grid2Coords	Grid2Coords
	Type used for local coordinates on the grid2 side. More...

Public Member Functions
	ContactMerge (const T allowedOverlap=T(0), std::function< WorldCoords(WorldCoords)> domainDirections=nullptr, std::function< WorldCoords(WorldCoords)> targetDirections=nullptr, ProjectionType type=OUTER_NORMAL)
	Construct merger given overlap and possible projection directions. More...

	ContactMerge (const T allowedOverlap, ProjectionType type)
	Construct merger given overlap and type of the projection. More...

void	setSurfaceDirections (std::function< WorldCoords(WorldCoords)> domainDirections, std::function< WorldCoords(WorldCoords)> targetDirections)
	Set surface direction functions. More...

void	setOverlap (T overlap)
	Set the allowed overlap of the surfaces. More...

T	getOverlap () const
	Get the allowed overlap of the surfaces. More...

void	minNormalAngle (T angle)
	set minimum angle in radians between normals at `x` and `Φ(x)` More...

T	minNormalAngle () const
	get minimum angle in radians between normals at `x` and `Φ(x)` More...

unsigned int	nSimplices () const
	get the number of simplices in the merged grid The indices are then in 0..nSimplices()-1 More...

void	clear ()

void	enableFallback (bool fallback)

void	enableBruteForce (bool bruteForce)

unsigned int	parents (unsigned int idx) const

unsigned int	parent (unsigned int idx, unsigned int parId=0) const
	get index of grid-n's parent simplex for given merged grid simplex More...

bool	simplexRefined (unsigned int idx, std::vector< unsigned int > &indices) const
	get the merged grid simplices refining a given grid-n simplex More...

GridTraits< n >::Coords	parentLocal (unsigned int idx, unsigned int corner, unsigned int parId=0) const
	get the grid-n parent's simplex local coordinates for a particular merged grid simplex corner (parent's index can be obtained via "parent<n>") More...

Public Attributes
unsigned int	counter
	Counts the number of times the computeIntersection method has been called. More...

Protected Types
typedef StandardMerge< T, dimworld-1, dimworld-1, dimworld >::RemoteSimplicialIntersection	RemoteSimplicialIntersection

Protected Member Functions
void	build (const std::vector< Dune::FieldVector< T, dimworld > > &grid1Coords, const std::vector< unsigned int > &grid1Elements, const std::vector< Dune::GeometryType > &grid1ElementTypes, const std::vector< Dune::FieldVector< T, dimworld > > &grid2Coords, const std::vector< unsigned int > &grid2Elements, const std::vector< Dune::GeometryType > &grid2ElementTypes)
	builds the merged grid More...

void	computeCyclicOrder (const std::vector< std::array< LocalCoords, 2 > > &polytopeCorners, const LocalCoords &center, std::vector< int > &ordering) const
	Order the corners of the intersection polytope in cyclic order. More...

void	setupNodalDirections (const std::vector< WorldCoords > &coords1, const std::vector< unsigned int > &elements1, const std::vector< Dune::GeometryType > &elementTypes1, const std::vector< WorldCoords > &coords2, const std::vector< unsigned int > &elements2, const std::vector< Dune::GeometryType > &elementTypes2)
	Setup the direction vectors containing the directions for each vertex. More...

void	computeOuterNormalField (const std::vector< WorldCoords > &coords, const std::vector< unsigned int > &elements, const std::vector< Dune::GeometryType > &elementTypes, std::vector< WorldCoords > &normals)
	If no direction field was specified compute the outer normal field. More...

void	removeDoubles (std::vector< std::array< LocalCoords, 2 > > &polytopeCorners)
	Remove all multiples. More...

virtual void	computeIntersections (const Dune::GeometryType &grid1ElementType, const std::vector< Dune::FieldVector< T, dimworld > > &grid1ElementCorners, std::bitset<(1<< grid1Dim)> &neighborIntersects1, unsigned int grid1Index, const Dune::GeometryType &grid2ElementType, const std::vector< Dune::FieldVector< T, dimworld > > &grid2ElementCorners, std::bitset<(1<< grid2Dim)> &neighborIntersects2, unsigned int grid2Index, std::vector< RemoteSimplicialIntersection > &intersections)=0
	Compute the intersection between two overlapping elements. More...

bool	computeIntersection (unsigned int candidate0, unsigned int candidate1, const std::vector< Dune::FieldVector< T, dimworld > > &grid1Coords, const std::vector< Dune::GeometryType > &grid1_element_types, std::bitset<(1<< grid1Dim)> &neighborIntersects1, const std::vector< Dune::FieldVector< T, dimworld > > &grid2Coords, const std::vector< Dune::GeometryType > &grid2_element_types, std::bitset<(1<< grid2Dim)> &neighborIntersects2, bool insert=true)
	Compute the intersection between two overlapping elements. More...

Protected Attributes
bool	valid

std::vector< RemoteSimplicialIntersection >	intersections_
	The computed intersections. More...

std::vector< std::vector< unsigned int > >	grid1ElementCorners_
	Temporary internal data. More...

std::vector< std::vector< unsigned int > >	grid2ElementCorners_

std::vector< std::vector< int > >	elementNeighbors1_

std::vector< std::vector< int > >	elementNeighbors2_

Detailed Description

template<int dimworld, typename T = double>
class Dune::GridGlue::ContactMerge< dimworld, T >

Merge two codimension-1 surfaces that may be a positive distance apart.

Template Parameters

dimworld	Dimension of the world coordinates.
T	Type used for coordinates

Member Typedef Documentation

◆ ctype

template<int dimworld, typename T = double>

typedef T Dune::GridGlue::ContactMerge< dimworld, T >::ctype

the numeric type used in this interface

◆ Grid1Coords

typedef Merger<T,grid1Dim,grid2Dim,dimworld>::Grid1Coords Dune::GridGlue::StandardMerge< T, grid1Dim, grid2Dim, dimworld >::Grid1Coords

inherited

Type used for local coordinates on the grid1 side.

◆ Grid2Coords

typedef Merger<T,grid1Dim,grid2Dim,dimworld>::Grid2Coords Dune::GridGlue::StandardMerge< T, grid1Dim, grid2Dim, dimworld >::Grid2Coords

inherited

Type used for local coordinates on the grid2 side.

◆ LocalCoords

template<int dimworld, typename T = double>

typedef Dune::FieldVector<T, dim> Dune::GridGlue::ContactMerge< dimworld, T >::LocalCoords

the coordinate type used in this interface

◆ RemoteSimplicialIntersection

template<int dimworld, typename T = double>

typedef StandardMerge<T,dimworld-1,dimworld-1,dimworld>::RemoteSimplicialIntersection Dune::GridGlue::ContactMerge< dimworld, T >::RemoteSimplicialIntersection

protected

◆ WorldCoords

template<int dimworld, typename T = double>

typedef Dune::FieldVector<T, dimworld> Dune::GridGlue::ContactMerge< dimworld, T >::WorldCoords

the coordinate type used in this interface

Member Enumeration Documentation

◆ ProjectionType

template<int dimworld, typename T = double>

enum Dune::GridGlue::ContactMerge::ProjectionType

Type of the projection, closest point or outer normal projection.

Enumerator
OUTER_NORMAL
CLOSEST_POINT

Constructor & Destructor Documentation

◆ ContactMerge() [1/2]

template<int dimworld, typename T = double>

Dune::GridGlue::ContactMerge< dimworld, T >::ContactMerge	(	const T	allowedOverlap = `T(0)`,
		std::function< WorldCoords(WorldCoords)>	domainDirections = `nullptr`,
		std::function< WorldCoords(WorldCoords)>	targetDirections = `nullptr`,
		ProjectionType	type = `OUTER_NORMAL`
	)

inline

Construct merger given overlap and possible projection directions.

Parameters

allowedOverlap	Allowed overlap of the surfaces
domainDirections	Projection direction field for the first surface that differ from the defualt normal field
targetDirections	Projection direction field for the second surface that differ from the default normal field

◆ ContactMerge() [2/2]

template<int dimworld, typename T = double>

Dune::GridGlue::ContactMerge< dimworld, T >::ContactMerge	(	const T	allowedOverlap,
		ProjectionType	type
	)

inline

Construct merger given overlap and type of the projection.

Parameters

allowedOverlap	Allowed overlap of the surfacs
type	Type of the projection

Member Function Documentation

◆ build()

template<int dimworld, typename T = double>

void Dune::GridGlue::ContactMerge< dimworld, T >::build	(	const std::vector< Dune::FieldVector< T, dimworld > > &	grid1Coords,
		const std::vector< unsigned int > &	grid1Elements,
		const std::vector< Dune::GeometryType > &	grid1ElementTypes,
		const std::vector< Dune::FieldVector< T, dimworld > > &	grid2Coords,
		const std::vector< unsigned int > &	grid2Elements,
		const std::vector< Dune::GeometryType > &	grid2ElementTypes
	)

inlineprotectedvirtual

builds the merged grid

Note that the indices are used consequently throughout the whole class interface just like they are introduced here.

Parameters

grid1_coords	the grid1 vertices' coordinates ordered like e.g. in 3D x_0 y_0 z_0 x_1 y_1 ... y_(n-1) z_(n-1)
grid1_elements	array with all grid1 elements represented as corner indices into `grid1_coords`
grid1_element_types	array with the GeometryType of the elements listed grid1_elements
grid2_coords	the grid2 vertices' coordinates ordered like e.g. in 3D x_0 y_0 z_0 x_1 y_1 ... y_(n-1) z_(n-1)
grid2_elements	just like with the grid1_elements and grid1_coords
grid2_element_types	array with the GeometryType of the elements listed grid2_elements

Reimplemented from Dune::GridGlue::StandardMerge< T, dimworld-1, dimworld-1, dimworld >.

◆ clear()

void Dune::GridGlue::StandardMerge< T, grid1Dim, grid2Dim, dimworld >::clear ( )

inlinevirtualinherited

Implements Dune::GridGlue::Merger< T, grid1Dim, grid2Dim, dimworld >.

◆ computeCyclicOrder()

template<int dimworld, typename T >

void Dune::GridGlue::ContactMerge< dimworld, T >::computeCyclicOrder	(	const std::vector< std::array< LocalCoords, 2 > > &	polytopeCorners,
		const LocalCoords &	center,
		std::vector< int > &	ordering
	)		const

protected

Order the corners of the intersection polytope in cyclic order.

◆ computeIntersection()

bool Dune::GridGlue::StandardMerge< T, grid1Dim, grid2Dim, dimworld >::computeIntersection	(	unsigned int	candidate0,
		unsigned int	candidate1,
		const std::vector< Dune::FieldVector< T, dimworld > > &	grid1Coords,
		const std::vector< Dune::GeometryType > &	grid1_element_types,
		std::bitset<(1<< grid1Dim)> &	neighborIntersects1,
		const std::vector< Dune::FieldVector< T, dimworld > > &	grid2Coords,
		const std::vector< Dune::GeometryType > &	grid2_element_types,
		std::bitset<(1<< grid2Dim)> &	neighborIntersects2,
		bool	insert = `true`
	)

protectedinherited

Compute the intersection between two overlapping elements.

Returns: true if at least one intersection point was found

◆ computeIntersections()

virtual void Dune::GridGlue::StandardMerge< T, grid1Dim, grid2Dim, dimworld >::computeIntersections	(	const Dune::GeometryType &	grid1ElementType,
		const std::vector< Dune::FieldVector< T, dimworld > > &	grid1ElementCorners,
		std::bitset<(1<< grid1Dim)> &	neighborIntersects1,
		unsigned int	grid1Index,
		const Dune::GeometryType &	grid2ElementType,
		const std::vector< Dune::FieldVector< T, dimworld > > &	grid2ElementCorners,
		std::bitset<(1<< grid2Dim)> &	neighborIntersects2,
		unsigned int	grid2Index,
		std::vector< RemoteSimplicialIntersection > &	intersections
	)

protectedpure virtualinherited

Compute the intersection between two overlapping elements.

The result is a set of simplices stored in the vector intersections.

◆ computeOuterNormalField()

template<int dimworld, typename T >

void Dune::GridGlue::ContactMerge< dimworld, T >::computeOuterNormalField	(	const std::vector< WorldCoords > &	coords,
		const std::vector< unsigned int > &	elements,
		const std::vector< Dune::GeometryType > &	elementTypes,
		std::vector< WorldCoords > &	normals
	)

protected

If no direction field was specified compute the outer normal field.

◆ enableBruteForce()

void Dune::GridGlue::StandardMerge< T, grid1Dim, grid2Dim, dimworld >::enableBruteForce ( bool bruteForce )

inlineinherited

◆ enableFallback()

void Dune::GridGlue::StandardMerge< T, grid1Dim, grid2Dim, dimworld >::enableFallback ( bool fallback )

inlineinherited

◆ getOverlap()

template<int dimworld, typename T = double>

T Dune::GridGlue::ContactMerge< dimworld, T >::getOverlap ( ) const

inline

Get the allowed overlap of the surfaces.

◆ minNormalAngle() [1/2]

template<int dimworld, typename T = double>

void Dune::GridGlue::ContactMerge< dimworld, T >::minNormalAngle ( T angle )

inline

set minimum angle in radians between normals at x and Φ(x)

◆ minNormalAngle() [2/2]

template<int dimworld, typename T = double>

T Dune::GridGlue::ContactMerge< dimworld, T >::minNormalAngle ( ) const

inline

get minimum angle in radians between normals at x and Φ(x)

◆ nSimplices()

unsigned int Dune::GridGlue::StandardMerge< T, grid1Dim, grid2Dim, dimworld >::nSimplices ( ) const

inlinevirtualinherited

get the number of simplices in the merged grid The indices are then in 0..nSimplices()-1

Implements Dune::GridGlue::Merger< T, grid1Dim, grid2Dim, dimworld >.

◆ parent()

unsigned int Dune::GridGlue::Merger< T , grid1Dim, grid2Dim, dimworld >::parent	(	unsigned int	idx,
		unsigned int	parId = `0`
	)		const

inlineinherited

get index of grid-n's parent simplex for given merged grid simplex

Template Parameters

n	specify which grid

Parameters

idx	index of the merged grid simplex

Returns: index of the parent simplex

◆ parentLocal()

GridTraits<n>::Coords Dune::GridGlue::Merger< T , grid1Dim, grid2Dim, dimworld >::parentLocal	(	unsigned int	idx,
		unsigned int	corner,
		unsigned int	parId = `0`
	)		const

inlineinherited

get the grid-n parent's simplex local coordinates for a particular merged grid simplex corner (parent's index can be obtained via "parent<n>")

Template Parameters

n	specify which grid

Parameters

idx	the index of the merged grid simplex
corner	the index of the simplex' corner

Returns: local coordinates in grid-n grid1

◆ parents()

unsigned int Dune::GridGlue::Merger< T , grid1Dim, grid2Dim, dimworld >::parents ( unsigned int idx ) const

inlineinherited

doc me

◆ removeDoubles()

template<int dimworld, typename T >

void Dune::GridGlue::ContactMerge< dimworld, T >::removeDoubles ( std::vector< std::array< LocalCoords, 2 > > & polytopeCorners )

protected

Remove all multiples.

◆ setOverlap()

template<int dimworld, typename T = double>

void Dune::GridGlue::ContactMerge< dimworld, T >::setOverlap ( T overlap )

inline

Set the allowed overlap of the surfaces.

◆ setSurfaceDirections()

template<int dimworld, typename T = double>

void Dune::GridGlue::ContactMerge< dimworld, T >::setSurfaceDirections	(	std::function< WorldCoords(WorldCoords)>	domainDirections,
		std::function< WorldCoords(WorldCoords)>	targetDirections
	)

inline

Set surface direction functions.

The matching of the geometries offers the possibility to specify a function for the exact evaluation of domain surface normals. If no such function is specified (default) normals are interpolated.

Parameters

value the new function (or nullptr to unset the function)

◆ setupNodalDirections()

template<int dimworld, typename T >

void Dune::GridGlue::ContactMerge< dimworld, T >::setupNodalDirections	(	const std::vector< WorldCoords > &	coords1,
		const std::vector< unsigned int > &	elements1,
		const std::vector< Dune::GeometryType > &	elementTypes1,
		const std::vector< WorldCoords > &	coords2,
		const std::vector< unsigned int > &	elements2,
		const std::vector< Dune::GeometryType > &	elementTypes2
	)

protected

Setup the direction vectors containing the directions for each vertex.

◆ simplexRefined()

bool Dune::GridGlue::Merger< T , grid1Dim, grid2Dim, dimworld >::simplexRefined	(	unsigned int	idx,
		std::vector< unsigned int > &	indices
	)		const

inlineinherited

get the merged grid simplices refining a given grid-n simplex

Template Parameters

n	specify which grid (grid1/grid2: 0/1)

Parameters

idx	index of grid-n simplex
indices	will be resized first and then filled with the refining simplices

Returns: TRUE <=> given simplex could be matched and is part of the merged grid

Member Data Documentation

◆ counter

unsigned int Dune::GridGlue::Merger< T , grid1Dim, grid2Dim, dimworld >::counter

inherited

Counts the number of times the computeIntersection method has been called.

Used temporarily to speed up the implementation

◆ elementNeighbors1_

std::vector<std::vector<int> > Dune::GridGlue::StandardMerge< T, grid1Dim, grid2Dim, dimworld >::elementNeighbors1_

protectedinherited

◆ elementNeighbors2_

std::vector<std::vector<int> > Dune::GridGlue::StandardMerge< T, grid1Dim, grid2Dim, dimworld >::elementNeighbors2_

protectedinherited

◆ grid1ElementCorners_

std::vector<std::vector<unsigned int> > Dune::GridGlue::StandardMerge< T, grid1Dim, grid2Dim, dimworld >::grid1ElementCorners_

protectedinherited

Temporary internal data.

◆ grid2ElementCorners_

std::vector<std::vector<unsigned int> > Dune::GridGlue::StandardMerge< T, grid1Dim, grid2Dim, dimworld >::grid2ElementCorners_

protectedinherited

◆ intersections_

std::vector<RemoteSimplicialIntersection> Dune::GridGlue::StandardMerge< T, grid1Dim, grid2Dim, dimworld >::intersections_

protectedinherited

The computed intersections.

◆ valid

bool Dune::GridGlue::StandardMerge< T, grid1Dim, grid2Dim, dimworld >::valid

protectedinherited

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

Public Types

Public Member Functions

Public Attributes

Protected Types

Protected Member Functions

Protected Attributes

Detailed Description

template<int dimworld, typename T = double> class Dune::GridGlue::ContactMerge< dimworld, T >

Member Typedef Documentation

◆ ctype

◆ Grid1Coords

◆ Grid2Coords

◆ LocalCoords

◆ RemoteSimplicialIntersection

◆ WorldCoords

Member Enumeration Documentation

◆ ProjectionType

Constructor & Destructor Documentation

◆ ContactMerge() [1/2]

◆ ContactMerge() [2/2]

Member Function Documentation

◆ build()

◆ clear()

◆ computeCyclicOrder()

◆ computeIntersection()

◆ computeIntersections()

◆ computeOuterNormalField()

◆ enableBruteForce()

◆ enableFallback()

◆ getOverlap()

◆ minNormalAngle() [1/2]

◆ minNormalAngle() [2/2]

◆ nSimplices()

◆ parent()

◆ parentLocal()

◆ parents()

◆ removeDoubles()

◆ setOverlap()

◆ setSurfaceDirections()

◆ setupNodalDirections()

◆ simplexRefined()

Member Data Documentation

◆ counter

◆ elementNeighbors1_

◆ elementNeighbors2_

◆ grid1ElementCorners_

◆ grid2ElementCorners_

◆ intersections_

◆ valid

template<int dimworld, typename T = double>
class Dune::GridGlue::ContactMerge< dimworld, T >