Home Research Papers Posters Manuals Talks Download People Help Internal

libs/ck-libs/fem/fem_adapt_algs.h

Go to the documentation of this file.

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

// This module implements high level mesh adaptivity algorithms that make use 
// of the primitive mesh adaptivity operations provided by fem_adapt(_new).
// Ask: TLW
#ifndef __CHARM_FEM_ADAPT_ALGS_H
#define __CHARM_FEM_ADAPT_ALGS_H

#include "charm++.h"
#include "tcharm.h"
#include "charm-api.h"
#include "ckvector3d.h"
#include "fem.h"
#include "fem_mesh.h"
#include "vector2d.h"  // **CW** tentative

#define SLIVERAREA 1.0e-18
#define REFINE_TOL 1.1  // Refine elements with average edge length > 
                        // REFINE_TOL*desiredEdgeLength
#define COARSEN_TOL 0.8 // Coarsen element with average edge length <
                        // COARSEN_TOL*desiredEdgeLength
#define QUALITY_MIN 0.6

class FEM_Adapt_Algs;
CtvExtern(FEM_Adapt_Algs *, _adaptAlgs);

class femMeshModify;
class FEM_Adapt;
class FEM_AdaptL;

class FEM_Adapt_Algs {
  friend class FEM_AdaptL;
  friend class FEM_Adapt;
  friend class femMeshModify;
  friend class FEM_Interpolate;
  friend class FEM_MUtil;

 public:
  int coord_attr;
  int bc_attr;

 protected: 
  FEM_Mesh *theMesh;
  femMeshModify *theMod;
  //FEM_Adapt *theAdaptor;
  FEM_AdaptL *theAdaptor;
  int numNodes, numElements, dim;
  // These are for element sorting
  typedef struct {
    int elID;
    double len;
  } elemHeap;
  elemHeap *coarsenElements;
  elemHeap *refineElements;
  elemHeap *refineStack;
  int refineTop, refineHeapSize, coarsenHeapSize;

 public:
  FEM_Adapt_Algs() {
    theMesh = NULL; theMod = NULL; theAdaptor = NULL;
  }
  FEM_Adapt_Algs(FEM_Mesh *m, femMeshModify *fm, int dimension);
  void FEM_Adapt_Algs_Init(int coord_at, int bc_at) {
    coord_attr = coord_at;
    bc_attr = bc_at;
  }

  void FEM_Refine(int qm, int method, double factor, double *sizes);

  void FEM_Coarsen(int qm, int method, double factor, double *sizes);

  void FEM_AdaptMesh(int qm, int method, double factor, double *sizes);
  void FEM_Smooth(int qm, int method);

  // FEM_Mesh_mooth
  //    Inputs  : meshP - a pointer to the FEM_Mesh object to smooth
  //        : nodes - an array of local node numbers to be smoothed.  Send
  //              NULL pointer to smooth all nodes.
  //        : nNodes - the size of the nodes array
  //        : attrNo - the attribute number where the coords are registered
  //    Shifts nodes around to improve mesh quality.  FEM_BOUNDARY attribute
  //    and interpolator function must be registered by user to maintain 
  //    boundary information.
  void FEM_mesh_smooth(FEM_Mesh *meshP, int *nodes, int nNodes, int attrNo);

  void FEM_Repair(int qm);

  void FEM_Remesh(int qm, int method, double factor, double *sizes);
  
  void SetReferenceMesh();
  void GradateMesh(double smoothness);
 private:
  // Helper methods
  int Refine(int qm, int method, double factor, double *sizes);
  int Coarsen(int qm, int method, double factor, double *sizes);
  void SetMeshSize(int method, double factor, double *sizes);
  void Insert(int elID, double len, int cFlag);
  int Delete_Min(int cFlag);
 public:

  virtual int refine_element_leb(int e);
  virtual void refine_flip_element_leb(int e, int p, int n1, int n2, 
                       double le);

  int simple_refine(double targetA, double xmin=0.0, double ymin=0.0, double xmax=1.0, double ymax=1.0);
  int simple_coarsen(double targetA, double xmin=0.0, double ymin=0.0, double xmax=1.0, double ymax=1.0);
  double length(int n1, int n2);
  double getArea(int n1, int n2, int n3);
  double length(double *n1_coord, double *n2_coord);
  double getArea(double *n1_coord, double *n2_coord, double *n3_coord);
  int getCoord(int n1, double *crds);
  int getShortestEdge(int n1, int n2, int n3, int* shortestEdge);
  double getAreaQuality(int elem);
  bool didItFlip(int n1, int n2, int n3, double *n4_coord);
  bool didItFlip(double *n1_coord, double *n2_coord, double *n3_coord, double *n4_coord);
  double getSignedArea(double *n1_coord, double *n2_coord, double *n3_coord);
  void tests(void);
};


#endif

---