TFIMapping.hpp

Go to the documentation of this file.

//-----------------------------------C++-------------------------------------//
// File: src/algs/TFIMapping.hpp
// Wednesday February 11 10:50 2011
// Brief: TFIMapping class definition: generate the all-quad mesh by
//       TFI mapping. Before using this class, four edges should have already
//       been meshed by EdgeMesher
//---------------------------------------------------------------------------//


#ifndef MESHKIT_TFIMAPPING_HPP
#define MESHKIT_TFIMAPPING_HPP

#include <stdlib.h>
#include <stdio.h>
#include <assert.h>
#include <string>
#include <iostream>
#include <fstream>
#include <string.h>
#include <limits.h>

#include <iGeom.h>
#include "Global.hpp"
#include "meshkit/MKCore.hpp"
#include "meshkit/SizingFunction.hpp"
#include "moab/ReadUtilIface.hpp"
#include <iMesh.h>
#include <set>
#include <iRel.h>
#include <vector>

#include "meshkit/MeshScheme.hpp"

using namespace std;

namespace MeshKit {
//===========================================================================//
//===========================================================================//
class TFIMapping: public MeshScheme {
public:

  TFIMapping(MKCore *mk_core, const MEntVector &me_vec);
  ~TFIMapping();

  virtual void setup_this();
  virtual void execute_this();

  static const char* name()
  {
    return "TFIMapping";
  }

  static bool can_mesh(iBase_EntityType dim)
  {
    return iBase_FACE == dim;
  }

  static bool can_mesh(ModelEnt *me)
  {
    return canmesh_region(me);
  }

  static const moab::EntityType* output_types();

  virtual const moab::EntityType* mesh_types_arr() const
  {
    return output_types();
  }
  void skip_improve()
  {
    _shapeImprove = false;// by default it is true, but sometimes we need to skip it
  }

private:

  bool _shapeImprove;
  //implement the transfinite interpolation between (pt_0s, pt_1s) and (pt_r0, pt_r1)
  void parametricTFI3D(double *pts, double r, double s, double pt_0s[3], double pt_1s[3], double pt_r0[3], double pt_r1[3]);

  //generate the mesh on the linking surface
  int SurfMapping(ModelEnt *ent);

 //generate the mesh on the cylinder-like linking surface
  int cylinderSurfMapping(ModelEnt *ent);

  // transform from source line to target line, such that
  // target= A * ( source - 2*sc + tc) + sc
  void computeTransformation(Vector3D & A, Vector3D & B, Vector3D & C, Vector3D & D,
      Matrix3D & transMatrix);

  void SurfImprove(iBase_EntityHandle surface, iBase_EntitySetHandle surfMesh, iBase_EntityType entity_type);

  //smooth the quadrilateral mesh on the linking surface
  void SmoothWinslow(std::vector<iBase_EntityHandle> &List_i, std::vector<iBase_EntityHandle> &List_ii, std::vector<
      iBase_EntityHandle> &List_j, std::vector<iBase_EntityHandle> &List_jj, std::vector<iBase_EntityHandle> &InteriorNodes,
      std::vector<iBase_EntityHandle> &quads, iBase_TagHandle &taghandle, ModelEnt *ent);

};

}

#endif