libs/ck-libs/tmr/refine.h File Reference

Go to the source code of this file.

Functions
void	REFINE2D_Init (void)
	Create a new refinement object for this virtual processor.
void	REFINE2D_NewMesh (int nEl, int nGhost, const int *conn, const int *gid)
	Push a new mesh into the refinement system.
void	REFINE2D_Split (int nNode, double *coord, int nEl, double *desiredArea)
	Refine the mesh so each element has, at most, the given desired area.
int	REFINE2D_Get_Split_Length (void)
	Get the number of split triangles.
void	REFINE2D_Get_Split (int splitNo, const int *conn, int *tri, int *A, int *B, int *C, double *frac, int *flags)
	Return one split triangle.
void	REFINE2D_Check (int nEle, const int *conn, int nNodes)
	Check to make sure our connectivity and the refine connectivity agree.

---

Function Documentation

void REFINE2D_Init

(

void

)

Create a new refinement object for this virtual processor.

Must be called exactly once at startup.

Definition at line 18 of file refine.C.

void REFINE2D_NewMesh	(	int	nEl,
		int	nGhost,
		const int *	conn,
		const int *	gid
	)

Push a new mesh into the refinement system.

This is the first call user programs make into the refinement system. This call need *not* be repeated when the refinement system changes the mesh; only when the user changes the mesh (e.g., to coarsen it).

Conn is row-major, and maps an element number to three node numbers. Hence conn[i*3+j] gives the local number of node j of element i. Because of shared nodes, every node of every local element will be local. Ghost elements may not have the complete set of nodes-- some of their nodes may have the invalid number -1.

Elements with numbers between 0 and nEl-1 (inclusive) are local. Elements with numbers between nEl and nGhost-1 (inclusive) are "ghosts"-- elements that are actually local on another processor. There are guaranteed to be enough ghosts that every local element's non-boundary edge will face a ghost element.

gid maps an element number to a chunk number and local number on that chunk. These are stored at gid[i*2+0] (chunk number) and gid[i*2+1] (local number).

Definition at line 48 of file refine.C.

References CkWaitQD(), and MPI_Barrier.

Referenced by FEM_REFINE2D_Newmesh().

Here is the call graph for this function:

Here is the caller graph for this function:

void REFINE2D_Split	(	int	nNode,
		double *	coord,
		int	nEl,
		double *	desiredArea
	)

Refine the mesh so each element has, at most, the given desired area.

Coord gives the (x,y) coordinates of nodes 0..nNode-1. coord[2*j+0] is the x coordinate of node j coord[2*j+1] is the y coordinate of node j

desiredArea[i] gives the desired area of element i. nEl must be equal to the "nEl" value passed in via REFINE_NewMesh plus the number of "split" calls received during earlier refinements.

Definition at line 147 of file refine.C.

References C, CkWaitQD(), chunk::multipleRefine(), chunk::refineResultsStorage, and chunk::updateNodeCoords().

Referenced by FEM_REFINE2D_Split(), and FTN_NAME().

Here is the call graph for this function:

Here is the caller graph for this function:

int REFINE2D_Get_Split_Length

(

void

)

Get the number of split triangles.

Definition at line 176 of file refine.C.

void REFINE2D_Get_Split	(	int	splitNo,
		const int *	conn,
		int *	tri,
		int *	A,
		int *	B,
		int *	C,
		double *	frac,
		int *	flags
	)

Return one split triangle.

A and B are the nodes along the splitting edge:

C C / \ /|\ / \ / | \ / \ => / | \ / tri \ / | \ / \ /tri | new\ B --------- A B --- D --- A

The original triangle's node A should be replaced by D; while a new triangle should be inserted with nodes CAD.

The new node D's location should equal A*(1-frac)+B*frac. For a simple splitter, frac will always be 0.5.

If nodes A and B are shared with some other processor, that processor will also receive a "split" call for the same edge. If nodes A and B are shared by some other local triangle, that triangle will immediately receive a "split" call for the same edge.

Parameters: -splitNo is the number of this split. Pass splitNo in increasing order from 0 to REFINE2D_Get_Split_Length()-1. -conn is the triangle connectivity array, used to look up the node numbers A, B, and C. This array is *not* modified. -tri returns the number of the old triangle being split; as labelled above. -A,B,C returns the numbers of the nodes in the above diagram. -frac returns the weighting for D between A and B; for now, this is always 0.5.

Client's responsibilities: -Add the new node D. Since both sides of a shared local edge will receive a "split" call, you must ensure the node is not added twice; you can do this by checking this split's nodes A and B against the previous split. -Update connectivity for source triangle -Add new triangle.

Definition at line 187 of file refine.C.

References refineResults::extract(), getResults(), and r.

Referenced by FEM_REFINE2D_Split().

Here is the call graph for this function:

Here is the caller graph for this function:

void REFINE2D_Check	(	int	nEle,
		const int *	conn,
		int	nNodes
	)

Check to make sure our connectivity and the refine connectivity agree.

Definition at line 256 of file refine.C.