Home Research Papers Posters Manuals Talks Download People Help Internal

libs/ck-libs/fem/fem_adapt_new.h

Go to the documentation of this file.

/* File: fem_adapt_new.h
 * Authors: Terry Wilmarth, Nilesh Choudhury
 *
 */


// This will eventually take the place of the previous version when Nilesh and 
// Isaac's fem_mesh_modify module is complete.
#ifndef __CHARM_FEM_ADAPT_H
#define __CHARM_FEM_ADAPT_H

#include "charm-api.h"
#include "ckvector3d.h"
#include "fem.h"
#include "fem_mesh.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:

  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 eConn[3];
    theMesh->e2n_getAll(e, eConn);
    for (int i=0; i<3; i++) 
      if ((eConn[i] != n2) && (eConn[i] != n1)) return eConn[i];
    return -1; //should never come here
  }
  int n2e_exists(int n, int e) {
    int *nConn, nSz;
    theMesh->n2e_getAll(n, &nConn, &nSz);
    for (int i=0; i<nSz; i++) {
      if (nConn[i] == e) {
    if(nSz!=0) free(nConn);
    return 1;
      }
    }
    if(nSz!=0) free(nConn);
    return 0;
  }
  int findElementWithNodes(int n1, int n2, int n3) {
    int *nConn, nSz;
    int ret = -1;
    theMesh->n2e_getAll(n1, &nConn, &nSz);
    for (int i=0; i<nSz; i++) {
      if ((n2e_exists(n2, nConn[i])) && (n2e_exists(n3, nConn[i]))) {
    ret = nConn[i];
    break;
      }
    }
    if(nSz!=0) free(nConn);
    return ret; //should never come here
  }
  int getSharedNodeIdxl(int n, int chk);
  int getGhostNodeIdxl(int n, int chk);
  int getGhostElementIdxl(int e, int chk);

  FEM_Adapt() {
    theMesh = NULL; theMod = NULL;
  }

  FEM_Adapt(FEM_Mesh *m, femMeshModify *fm) { theMesh = m; theMod = fm; }

  //virtual int edge_flip(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);
  

  //virtual int edge_bisect(int n1, int n2);
  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);


  //virtual int vertex_remove(int n1, int n2);
  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);


  //virtual int edge_contraction(int n1, int n2);
  //virtual 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);


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

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

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


#endif

---