structures.h

Go to the documentation of this file.

#ifndef _polyStruct_h_
#define _polyStruct_h_
#include <vector>
#include <cmath>


/**************************************************************************************/
/**************************************************************************************/
struct Poly{
    int numberOfNodes;
      
    int familyNum;
    
    unsigned int groupNum;

    float area;

    double xradius; 
    
    double yradius;        

    float aspectRatio; 

    double translation[3];

    double normal[3]; 

    double boundingBox[6];

    double aperture;  

    double *vertices;

    double permeability; 

    bool faces[6]; // Touching boundary faces of the fractures group, not neccesarily the faces of the fracture

    bool XYPlane;

    bool truncated; 

    std::vector<unsigned int> intersectionIndex; 
    
    // Constructor
    Poly();
};



/**************************************************************************************/
/**************************************************************************************/
struct Point{
    double x;
    double y;
    double z;

    Point();
    Point(double _x, double _y, double _z);
};



/**************************************************************************************/
/**************************************************************************************/
struct IntPoints {
    long int fract1; 

    long int fract2;

    double x1;
    double y1;
    double z1;
    double x2;
    double y2;
    double z2;
    std::vector<unsigned int> triplePointsIdx; 
    
    bool intersectionShortened;

    IntPoints(); // Constructor. initializes fract1, fract2 = -1
};



/**************************************************************************************/
/**************************************************************************************/
struct TriplePtTempData {
    Point triplePoint;
    std::vector<int> intIndex;
};



/**************************************************************************************/
/**************************************************************************************/
struct FractureGroups {
    unsigned long long int groupNum;
    std::vector<unsigned int> polyList;
    FractureGroups();
};



/**************************************************************************************/
/**************************************************************************************/
// GroupData works in conjunction with 'FractureGroups'
// GroupData keeps track of which group numbers ('groupNum' in above struct) 
// Connect to which boundaries. When a cluster of fractures bridges two 
// groups/clusters, One of the groups 'valid' bool is set to 0 (not valid). 
// The fractures whos valid bool becomes 0 are moved to the other group 
// In the 'GroupData' structure array, FractureGroups.groupNum-1 is the index
// to the corresponding GroupData struct in the GroupData array. 
struct GroupData {
    unsigned int size;
    bool valid;
    bool faces[6];

    // Constructor sets size to zero an sets 'valid' bit and 'faces' elements 
    //    to false.
    GroupData();
};



/**************************************************************************************/
/**************************************************************************************/
struct RejectionReasons {
    unsigned long long int shortIntersection;
    unsigned long long int closeToNode;
    unsigned long long int closeToEdge;
    /* Counter no longer in use. */
    unsigned long long int closePointToEdge;
    unsigned long long int outside;
    unsigned long long int triple;
    unsigned long long int interCloseToInter;
    RejectionReasons();
};


/**************************************************************************************/
/**************************************************************************************/
// TODO: Make singleton
struct Stats { 

    int *acceptedFromFam;

    int *rejectedFromFam;

    int *expectedFromFam;

    unsigned int acceptedPolyCount;

    unsigned long long int rejectedPolyCount;

    unsigned int retranslatedPolyCount;

    unsigned int truncated;

 //   unsigned int numIntPoints;

    double areaBeforeRemoval;

    double areaAfterRemoval;

    double volBeforeRemoval; 
 
    double volAfterRemoval;

    unsigned long long int nextGroupNum;
    
    struct RejectionReasons rejectionReasons;
    
    unsigned int intersectionsShortened;

    double originalLength; // Length of all intersections if none were shortened

    double discardedLength; 
    
    unsigned int intersectionNodeCount;

    unsigned int tripleNodeCount;
 
    std::vector<unsigned int> rejectsPerAttempt;   

    std::vector<struct FractureGroups> fractGroup;
    std::vector<struct GroupData> groupData;
    // Constructor 
    Stats();
};



/**************************************************************************************/
/**************************************************************************************/
struct Shape {

    short shapeFamily; 

    short distributionType; 

    short numPoints;   

    float *thetaList;     

    unsigned int radiiIdx;      

    std::vector<double> radiiList;  

    short layer;  

    float aspectRatio;

    float p32Target; 

    float currentP32; 

    bool angleOption;

    bool betaDistribution;

    float beta;
    
    double theta; 

    double phi; 

    double kappa;
    

    /**************** Distribution Variables *********************/
    /*************************************************************/
    double minDistInput; 

    double maxDistInput; 
    
    float expMean;

    float expLambda;

    float expMin;

    float expMax;
    
    float mean;
    
    float sd;

    float logMin;
    float logMax;

    float constRadi;
    
    float min; 

    float max; 

    float alpha;
    
    Shape(); // Constructor
};


void printPolyData(struct Poly &Poly); 
void printStats(struct stats *obj);    

#endif