TetGenMeshLoader.h

// =======================================================================
// @(#)Domain.C        1.23 05/05/00
// 
// Class:       TTetGenMeshLoader
// Purpose:     creates domain with tetgen from a .smesh file
//
// Author:      Andreas Hahn 26.04.2010
//
// History:
//
// =======================================================================

#ifndef __TETGENMESHLOADER__
#define __TETGENMESHLOADER__

#define __TETGEN_14X__ 

#include <tetgen.h>

#include <Constants.h>
#include <MooNMD_Io.h>

class TBoundPart;
class TBoundComp3D;
class TVertex;
class TBaseCell;
class TJoint;
class TDomain;

class TTetGenMeshLoader
{
  protected:
    char *mFileName;
    
    char *mBoundFile;
    
    TBoundComp3D **mAllBoundComps;
    
    TVertex **mAllVertices;
    
    TJoint **mAllJoints;
    
    int **mTrifaceHash;
    
    tetgenio mTetIn;
    tetgenio mTetOut;
#ifdef __TETGEN_14X__
    tetgenbehavior mTetBeh;
    tetgenio mTetAddIn;
#endif
    
    bool plc;
    bool reconstruct;
    bool insertpoints;
    
  protected:
    int ReadBoundFile(int N_BoundComps);
    int ReadBoundComp(std::ifstream *in, int &CurrBoundComp, int &Range);
    int ReadBoundType(std::ifstream *in, char *buff);
    int ReadBoundParams(std::ifstream *in, char *type, int CurrBoundComp, int Range);
    int ReadBoundParamsCylinder(std::ifstream *in, double &radius,
                double &px, double &py, double &pz,
                double &ax, double &ay, double &az,
                double &nx, double &ny, double &nz);
    int ReadBoundParamsSphere(std::ifstream *in, double &radius, 
                  double &mx, double &my, double &mz);
                
    int ExtractValue(const char *buff, int increment, double &val);
    int ExtractVector(const char *buff, int increment,
              double &val1, double &val2, double &val3);
              
    int SetBdPlaneParams(int N_BoundComps);
    
    void Normalize3D(double *vec);
    void Cross3D(double *vec1, double *vec2, double *res);
    
    int AllocVertices(double &StartX, double &StartY, double &StartZ,
              double &BoundX, double &BoundY, double &BoundZ);
    int AllocRootCells(TBaseCell **&CellTree, int &N_RootCells);
    int AllocJoints(TBaseCell **CellTree);
    int DistributeJoints(TBaseCell **CellTree);
    
    int  FindTriface(int a, int b, int c);
    void HashTrifaces();
    int  CreateAdjacency();
    
    void MakeOptionString(std::string &opts);
    
    void MakeBoundaryLayer_smesh();
    void MakeBoundaryLayer_msh();
    void FindBoundaryPoints_msh(int *pointlist, double *normallist);
    
    int Load_msh(const char *filename);    
    
    int BuildBoundary  (TBoundPart **&BdParts, int &N_BoundParts,
                int &N_BoundComps, int *&StartBdCompID, int *&Interfaces);
              
    int BuildBoundary2 (TBoundPart **&BdParts, int &N_BoundParts,
                int &N_BoundComps, int *&StartBdCompID, int *&Interfaces);
          
    int Tetgen();
    
    int BuildMooNMDMesh(TBaseCell **&CellTree, int &N_RootCells,
            double &StartX, double &StartY, double &StartZ,
            double &BoundX, double &BoundY, double &BoundZ);
  
  public:
    typedef struct
    {
      TBoundPart **BdParts;
      int N_BoundParts;
      int N_BoundComps;
      int *StartBdCompID;
      int *Interfaces;
      TBaseCell **CellTree;
      int N_RootCells;
      double StartX;
      double StartY;
      double StartZ;
      double BoundX;
      double BoundY;
      double BoundZ;
    } TDummyDomain;
            
  public:
    // CTOR
//     TTetGenMeshLoader() : mFileName(0), mAllBoundComps(0), mBoundFile(0),
//            mAllVertices(0), mAllJoints(0) {};
    TTetGenMeshLoader();
    TTetGenMeshLoader(const char *FileName);
    
    // DTOR
    ~TTetGenMeshLoader();
    
    int Generate(TDomain *Domain);
    int Generate(int N_Points, double *Points, int N_Facets, int *Facets,
         int N_Regions, double *Regions, TDummyDomain *Domain);
};
#endif