Home Research Papers Posters Manuals Talks Download People Help Internal

libs/ck-libs/TMRC2D/refine.h

Go to the documentation of this file.

/*Charm++ Mesh Refinement Framework:
  C++ implementation file
  
  Orion Sky Lawlor, olawlor@acm.org, 4/8/2002
  Modified by Terry Wilmarth, wilmarth@cse.uiuc.edu, 4/16/2002
*/
#ifndef _UIUC_CHARM_REFINE_IMPL_H
#define _UIUC_CHARM_REFINE_IMPL_H

#ifdef __cplusplus
extern "C" {
#endif
  
  // Create a new refinement object for this virtual processor.
  // Must be called exactly once at startup.
  void REFINE2D_Init(void);
  
  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);
  
     class FEM_Refine_Operation_Data;
  void REFINE2D_Split(int nNode, double *coord, int nEl, double *desiredArea,FEM_Refine_Operation_Data *refine_data);
 
    typedef struct{
    int tri,A,B,C,D, _new,origEdgeB, newEdge1B, newEdge2B; 
    double frac;
    int flag;
    } refineData;

 
  // Get the number of split triangles.
  int REFINE2D_Get_Split_Length(void);

  void REFINE2D_Get_Split(int splitNo, refineData *d);


  // Data structures to store the changes caused by coarsening.  
  // Store the data for 1 collapse operation
  typedef struct {
    int elemID, nodeToKeep, nodeToDelete;
    double newX, newY;
    int flag;
    int boundaryFlag;
    double frac;
  } collapseData;
  
  // Store the data for 1 nodeUpdate operation
  typedef struct {
    int nodeID, boundaryFlag;
    double newX, newY;
  } updateData;
  
  // Store the data for 1 nodeReplaceDelete operation
  typedef struct {
    int elemID, relnodeID, oldNodeID, newNodeID;
  } replacedeleteData;
  
  // Store the data for a change caused by coarsening, might be of 3 types: 
  // collapse, updatenode or replace and delete a node
  enum {INVALID=0,COLLAPSE,UPDATE,REPLACE};
  typedef struct {
    int type;
    union {
      collapseData cdata;
      updateData udata;
      replacedeleteData rddata;
    } data;
  } coarsenData;
  
  class FEM_Operation_Data;
  void REFINE2D_Coarsen(int nNode, double *coord, int nEl, double *desiredArea,
            FEM_Operation_Data *data);
  
  /*void REFINE2D_Get_Collapse(int i,int *conn,int *tri,int *nodeToThrow,int *nodeToKeep,double *nx,double *ny,int *flag,int idxbase);*/
  
  void REFINE2D_Get_Collapse(int i, coarsenData *output);
  int REFINE2D_Get_Collapse_Length();
  
  // Check to make sure our connectivity and the refine connectivity agree.
  void REFINE2D_Check(int nEle, const int *conn, int nNodes);
  
#ifdef __cplusplus
};
#endif

#endif

---