PPL Logo

libs/ck-libs/TMRC2D/refine.C File Reference

Go to the source code of this file.

Data Structures

class  refineResults
class  resultsRefineClient
class  coarsenResults
class  resultsCoarsenClient

Functions

CDECL void REFINE2D_Init (void)
 Create a new refinement object for this virtual processor.
FDECL void FTN_NAME (REFINE2D_INIT, refine2d_init)
CDECL void REFINE2D_NewMesh (int meshID, int nEl, int nGhost, int nnodes, const int *conn, const int *gid, const int *boundaries, const int *edgeBounds, const int *edgeConn, int nEdges)
 Push a new mesh into the refinement system.
void FEM_Modify_IDXL (FEM_Refine_Operation_Data *data, refineData &d)
void FEM_Coarsen_Operation (FEM_Operation_Data *coarsen_data, coarsenData &operation)
CDECL void REFINE2D_Split (int nNode, double *coord, int nEl, double *desiredArea, FEM_Refine_Operation_Data *refine_data)
CDECL void REFINE2D_Coarsen (int nNode, double *coord, int nEl, double *desiredArea, FEM_Operation_Data *data)
FDECL void FTN_NAME (REFINE2D_SPLIT, refine2d_split)
static refineResultsgetResults (void)
CDECL int REFINE2D_Get_Split_Length (void)
 Get the number of split triangles.
FDECL int FTN_NAME (REFINE2D_GET_SPLIT_LENGTH, refine2d_get_split_length)
CDECL void REFINE2D_Get_Split (int splitNo, refineData *d)
 Return one split triangle.
FDECL void FTN_NAME (REFINE2D_GET_SPLIT, refine2d_get_split)
static coarsenResultsgetCoarsenResults (void)
CDECL int REFINE2D_Get_Collapse_Length ()
CDECL void REFINE2D_Get_Collapse (int i, coarsenData *output)
static int checkElement (chunk *C, const element &e, const int *uc, int idxBase)
static void checkConn (int nEl, const int *conn, int idxBase, int nNode)
CDECL void REFINE2D_Check (int nEl, const int *conn, int nNodes)
 Check to make sure our connectivity and the refine connectivity agree.
FDECL void FTN_NAME (REFINE2D_CHECK, refine2d_check)


Function Documentation

CDECL void REFINE2D_Init ( void   ) 

Create a new refinement object for this virtual processor.

Must be called exactly once at startup.

Definition at line 16 of file refine.C.

References CkArrayOptions::bindTo(), TCharm::get(), TCharm::getNumElements(), TCharm::getProxy(), mesh, MPI_Bcast(), MPI_BYTE, MPI_Comm_rank(), MPI_COMM_WORLD, chunkMsg::myThreads, chunkMsg::nChunks, opts, and TCharm::suspend().

FDECL void FTN_NAME ( REFINE2D_INIT  ,
refine2d_init   
)

Definition at line 40 of file refine.C.

References REFINE2D_Init().

CDECL void REFINE2D_NewMesh ( int  meshID,
int  nEl,
int  nGhost,
int  nnodes,
const int conn,
const int gid,
const int boundaries,
const int edgeBounds,
const int edgeConn,
int  nEdges 
)

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).

boundaries is the boundary flags for nodes. boundaries may be null, if the user doesnt specify boundary flags. nnodes specifies the number of nodes

Definition at line 46 of file refine.C.

References CkWaitQD(), MPI_Barrier(), and MPI_COMM_WORLD.

void FEM_Modify_IDXL ( FEM_Refine_Operation_Data data,
refineData d 
)

void FEM_Coarsen_Operation ( FEM_Operation_Data coarsen_data,
coarsenData operation 
)

CDECL void REFINE2D_Split ( int  nNode,
double coord,
int  nEl,
double desiredArea,
FEM_Refine_Operation_Data refine_data 
)

CDECL void REFINE2D_Coarsen ( int  nNode,
double coord,
int  nEl,
double desiredArea,
FEM_Operation_Data data 
)

FDECL void FTN_NAME ( REFINE2D_SPLIT  ,
refine2d_split   
)

Definition at line 272 of file refine.C.

References data, and REFINE2D_Split().

static refineResults* getResults ( void   )  [static]

Definition at line 278 of file refine.C.

References C, and chunk::refineResultsStorage.

CDECL int REFINE2D_Get_Split_Length ( void   ) 

Get the number of split triangles.

Definition at line 288 of file refine.C.

References refineResults::countResults(), and getResults().

FDECL int FTN_NAME ( REFINE2D_GET_SPLIT_LENGTH  ,
refine2d_get_split_length   
)

Definition at line 293 of file refine.C.

References REFINE2D_Get_Split_Length().

CDECL void REFINE2D_Get_Split ( int  splitNo,
refineData d 
)

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 299 of file refine.C.

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

FDECL void FTN_NAME ( REFINE2D_GET_SPLIT  ,
refine2d_get_split   
)

Definition at line 305 of file refine.C.

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

static coarsenResults* getCoarsenResults ( void   )  [static]

Definition at line 313 of file refine.C.

References C, and chunk::coarsenResultsStorage.

Referenced by REFINE2D_Get_Collapse(), and REFINE2D_Get_Collapse_Length().

CDECL int REFINE2D_Get_Collapse_Length (  ) 

Definition at line 324 of file refine.C.

References coarsenResults::countResults(), and getCoarsenResults().

Referenced by FEM_REFINE2D_Coarsen().

CDECL void REFINE2D_Get_Collapse ( int  i,
coarsenData output 
)

Definition at line 332 of file refine.C.

References coarsenResults::extract(), and getCoarsenResults().

static int checkElement ( chunk C,
const element e,
const int uc,
int  idxBase 
) [static]

Definition at line 337 of file refine.C.

static void checkConn ( int  nEl,
const int conn,
int  idxBase,
int  nNode 
) [static]

Definition at line 354 of file refine.C.

References C, checkElement(), chunk::cid, chunk::numElements, chunk::numNodes, and chunk::theElements.

CDECL void REFINE2D_Check ( int  nEl,
const int conn,
int  nNodes 
)

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

Definition at line 388 of file refine.C.

References checkConn().

FDECL void FTN_NAME ( REFINE2D_CHECK  ,
refine2d_check   
)

Definition at line 391 of file refine.C.

References checkConn().


Generated on Tue May 22 07:44:31 2012 for Charm++ by  doxygen 1.5.5