filament.h

Go to the documentation of this file.

#ifndef FILAMENT_H
#define FILAMENT_H
 
#include "vect.h"
#include <vector>
#include "IGBheader.h"
 
#define log_msg(F,L,O,...) fprintf( stderr, __VA_ARGS__ )
 
/*
 * This file is used for testing the local FEM routines
 */
 
typedef enum {face_based=0, vol_based} SingularityDetector_t;
typedef enum {PhaseSingularity=1,FilamentSegment=2} singularity_t;
enum Elem_t {Tetra=0, Hexa, Octa, Pyramid, Prism, Quad, Tri, Line, TetP2, TriP2, TetraBubble, HexaBubble, PrismBubble, PyramidBubble, TriBubble, QuadBubble};
 
typedef struct _PS {
  int type;
  int eIndx;
  Point p0;
} PS;
 
typedef struct _FilSeg {
  int type;
  int eIndx;
  Point p0;
  Point p1;
} FilSeg;
 
typedef union _singularity {
  PS     ps;
  FilSeg fs;
} singularity;
 
 
#define BLOCK_ALLOC 1000
class filament {
    public:
        void add_singularity( singularity *s, int );
        int num_segs(){ return (int)s.size(); }
        int num_pts(){ return numPts; }
        int          numPts=0;
        std::vector<singularity> s;
        void sort();
};
 
typedef struct NodeList {
  Point          *Node;        
  Integer         NNodes;      
}NodeList;
 
struct Elem {
  Elem_t type;
  int    Nnode;
  int   *N;
};
 
typedef struct ElemList {
  Integer        NElems;          
  Elem          *Ele;
}ElemList;
 
// copied from CARP/ioutils.h
typedef struct {
  const char  *dname;         
  FILE        *data;          
  float      **d;             
  int          dsize;         
  int          numSamps;      
  int          nodes;         
  int          h_wdth;        
  int          wdth;          
  int          buffered;      
  float      **dr;            
  int          rwdth;         
  float        rdist;         
  float        t;             
  float        dt_in;         
  bool         filled;        
  IGBheader   *igb;           
} tbuffer;
 
#ifdef __cplusplus
extern "C" {
#endif
 
// prototypes
int  init_filament_IO(const char *meshname, int numT, FILE *fh[]);
void write_filaments(filament *f, FILE *fh[],float t);
int  write_connections(ElemList *elst,NodeList *nlst,filament *f,FILE *filHdls[],float t);
int  write_elems_martin(ElemList *elst, NodeList *nlst, filament *f, FILE *filHdls[],float t);
void close_filament_IO(FILE *fh[]);
 
void interpolateTBuffer(tbuffer *tbf);
void interpolate_vm(int numNodes,float **d,int M,float **rd,float rDist);
bool compute_filament(float *d0, float *d1, float vth, ElemList *elst, 
                      NodeList *nlst, filament *f, SingularityDetector_t meth);
#ifdef __cplusplus
}
#endif
 
#endif