Home Research Papers Posters Manuals Talks Download People Help Internal

libs/ck-libs/ParFUM/ParFUM_Adapt.h

Go to the documentation of this file.

/* This module contains the class which presents the 
 * edge_bisect, edge_contract, edge_flip functions.
 *
 * File: ParFUM_Adapt.h
 * Authors: Terry Wilmarth, Nilesh Choudhury
 */

#ifndef __ParFUM_Adapt_H
#define __ParFUM_Adapt_H

class femMeshModify;


class FEM_Adapt {
 protected:
  FEM_Mesh *theMesh;
  femMeshModify *theMod;
  // Helper methods: see bottom of this file
  int check_orientation(int e1, int e3, int n, int n1, int n2);

 public:
  FEM_Adapt() {
    theMesh = NULL; theMod = NULL;
  }
  FEM_Adapt(FEM_Mesh *m, femMeshModify *fm) { 
    theMesh = m; theMod = fm; 
  }
  FEM_Adapt(femMeshModify *fm) { theMesh = NULL; theMod = fm; }
  void FEM_AdaptSetMesh(FEM_Mesh *m) {theMesh = m;}
  void pup(PUP::er &p) {
  }

  int get_edge_index(int local_node1, int local_node2);
  int find_local_node_index(int e, int n);

  int findAdjData(int n1, int n2, int *e1, int *e2, int *e1n1, int *e1n2, 
           int *e1n3, int *e2n1, int *e2n2, int *e2n3, int *n3, 
           int *n4);
  int e2n_getNot(int e, int n1, int n2);
  int n2e_exists(int n, int e);
  int findElementWithNodes(int n1, int n2, int n3);

  int getSharedNodeIdxl(int n, int chk);
  int getGhostNodeIdxl(int n, int chk);
  int getGhostElementIdxl(int e, int chk);

  void printAdjacencies(int *nodes, int numNodes, int *elems, int numElems);

  bool isFixedNode(int n1);
  bool isCorner(int n1);
  bool isEdgeBoundary(int n1, int n2);


  int edge_flip_help(int e1, int e2, int n1, int n2, int e1_n1, 
                 int e1_n2, int e1_n3, int n3, int n4,int *locknodes);
  

  int edge_bisect_help(int e1, int e2, int n1, int n2, int e1_n1, 
                   int e1_n2, int e1_n3, int e2_n1, int e2_n2, 
                   int e2_n3, int n3, int n4);


  int vertex_remove_help(int e1, int e2, int n1, int n2, int e1_n1, 
                 int e1_n2, int e1_n3, int e2_n1, int e2_n2, 
                 int e2_n3, int n3, int n4, int n5);


  int vertex_split(int n, int n1, int n2);
  int vertex_split_help(int n, int n1, int n2, int e1, int e3);
};



class FEM_AdaptL : public FEM_Adapt {
 public:
  FEM_AdaptL() {
    theMesh = NULL; theMod = NULL;
  }
  FEM_AdaptL(FEM_Mesh *m, femMeshModify *fm) { theMesh = m; theMod = fm; }
  FEM_AdaptL(femMeshModify *fm) { theMesh = NULL; theMod = fm; }
  void FEM_AdaptLSetMesh(FEM_Mesh *m) {theMesh = m;}
  void pup(PUP::er &p) {
  }

  int lockNodes(int *gotlocks, int *lockrnodes, int numRNodes, int *lockwnodes, int numWNodes);
  int unlockNodes(int *gotlocks, int *lockrnodes, int numRNodes, int *lockwnodes, int numWNodes);

  int edge_flip(int n1, int n2);
  int edge_bisect(int n1, int n2);
  int vertex_remove(int n1, int n2);

  int edge_contraction(int n1, int n2);
  int edge_contraction_help(int *e1P, int *e2P, int n1, int n2, int e1_n1, 
                    int e1_n2, int e1_n3, int e2_n1, int e2_n2,
                    int e2_n3, int n3, int n4);

  int eatIntoElement(int e, bool aggressive_node_removal=false);
  void residualLockTest();
  void unlockAll();
  void structureTest();
};

#endif

---