Home Research Papers Posters Manuals Talks Download People Help Internal

libs/ck-libs/ParFUM/bulk_adapt_ops.h

Go to the documentation of this file.

#ifndef __BULK_ADAPT_OPS_H
#define __BULK_ADAPT_OPS_H

#include "charm++.h"
#include "ParFUM.h"
#include "ParFUM_internals.h"
#include "ParFUM_SA.h"


class BulkAdapt {
 private:
  int partitionID;
  int meshID;
  FEM_Mesh *meshPtr;
  CProxy_ParFUMShadowArray shadowProxy;
  ParFUMShadowArray *localShadow;

  adaptAdj *elemPairs[10];
  int freeTable[10];
  int numGathered[10];
  int firstFree;
 public:
  BulkAdapt(int meshid, FEM_Mesh *mPtr, int partID, CProxy_ParFUMShadowArray sa_proxy);
  ~BulkAdapt();

  void pup(PUP::er &p);

  /* BASIC DATA ACCESS FUNCTIONS */
  
  int getPartition() { return partitionID; }
  FEM_Mesh *getMeshPtr() { return meshPtr; }

  int getTableID() { 
    int x = firstFree; 
    freeTable[x]=0; 
    numGathered[x] = 0;
    firstFree++; 
    while (!freeTable[firstFree]) {
      firstFree++;
      if (firstFree == 10) {
    CkPrintf("ERROR: elemPairs table is full!\n");
    break;
      }
    }
    return x; 
  }
  void freeTableID(int x) { 
    freeTable[x]=1; 
    free(elemPairs[x]);
    if (x<firstFree) firstFree = x;
  }

  /* BULK MESH OPERATIONS: These are all called locally, but may invoke
     remote operations. */


    bool edge_bisect(int elemID, int elemType, int edgeID, int dim, RegionID lockRegionID);
  bool edge_bisect_2D(int elemID, int elemType, int edgeID);
  bool edge_bisect_3D(int elemID, int elemType, int edgeID, RegionID lockRegionID);

  //TODO: add elemType to the prototype of all the following mesh-modification functions

  int edge_flip(int elemID, int edgeID);


  int flip_23(int elemID, int faceID);


  int flip_32(int elemID, int edgeID);


  int edge_collapse(int elemID, int edgeID);

  void one_side_split_2D(adaptAdj &startElem, adaptAdj &splitElem, int edgeID,
             int *node1idx, int *node2idx, int *newNodeID,
             bool startSide);

  /* COMMUNICATION HELPERS FOR BULK ADAPTIVITY OPERATIONS */
  adaptAdj remote_edge_bisect_2D(adaptAdj nbrElem, adaptAdj splitElem, 
                 int new_idxl, int n1_idxl, int n2_idxl, 
                 int remotePartID);

  void remote_adaptAdj_replace(adaptAdj elem, adaptAdj oldElem, 
                   adaptAdj newElem);
  void remote_edgeAdj_replace(int remotePartID, adaptAdj adj, adaptAdj elem, 
                  adaptAdj splitElem, int n1_idxl, int n2_idxl);
  void remote_edgeAdj_add(int remotePartID, adaptAdj adj, adaptAdj splitElem,
              int n1_idxl, int n2_idxl);

  adaptAdj remote_edge_bisect_3D(adaptAdj nbrElem, adaptAdj splitElem, 
                 int new_idxl, int n1_idxl, int n2_idxl, 
                 int remotePartID);

  void handle_split_3D(int remotePartID, int pos, int tableID, adaptAdj elem, RegionID lockRegionID,
               int n1_idxl, int n2_idxl, int n5_idxl);
  void recv_split_3D(int pos, int tableID, adaptAdj elem, adaptAdj splitElem);
  bool all_splits_received(int tableID, int expectedSplits);
  void update_asterisk_3D(int remotePartID, int i, adaptAdj elem, 
              int numElemPairs, adaptAdj *elemPairs, RegionID lockRegionID,
              int n1_idxl, int n2_idxl, int n5_idxl);

  bool isLongest(int elem, int elemType, double len);

  /* LOCAL HELPERS FOR BULK ADAPTIVITY OPERATIONS */

  int lock_3D_region(int elemID, int elemType, int edgeID, double prio, RegionID *lockRegionID);
  void unlock_3D_region(RegionID lockRegionID);
  void unpend_3D_region(RegionID lockRegionID);
    
  int add_element(int elemType,int nodesPerElem,int *conn,double sizing);
  
  void update_element_conn(int elemType,int elemID,int nodesPerElem,int *conn);
  
  int add_node(int dim,double *coords);
  
  void update_node_coord(int nodeID,int dim,double *coords);
  
  void make_node_shared(int nodeID,int numSharedChunks,int *sharedChunks);
  
  int get_idxl_for_node(int nodeID, int partID);
  int get_node_from_idxl(int node_idxl, int partID);
  bool is_node_in_idxl(int node_idxl, int partID);

  void get_elemsToLock(adaptAdj startElem, adaptAdj **elemsToLock, int edgeID, int *count);

  adaptAdj *local_split_3D(const adaptAdj elem, int n1, int n2, int n5);

  void local_update_asterisk_3D(int i, adaptAdj elem, int numElemPairs, 
                adaptAdj *elemPairs, int n1, 
                int n2, int n5);
  void update_local_face_adj(const adaptAdj elem, const adaptAdj splitElem, int n1, int n2, int n5);
  void update_local_edge_adj(const adaptAdj elem, const adaptAdj splitElem, int n1, int n2, int n5);
  double length(int n1, int n2, int dim);
  double length(double *n1, double *n2, int dim);
  void dumpConn();
};

// GENERAL HELPER FUNCTIONS

void midpoint(double *n1, double *n2, int dim, double *result);

int getRelNode(int nodeIdx, int *conn, int nodesPerElem);
void getRelNodes(int edgeID, int nodesPerElem, int *r1, int *r2);
int getEdgeID(int node1, int node2, int nodePerElem, int dim);
int getFaceID(int node1, int node2, int node3, int nodesPerElem);

void fillNodes(int *relNode, int *nodeIDs, int *conn);
void fillNodes(int *relNode, int n1, int n2, int *conn);

#endif

---