slopes.h

#ifndef __SLOPES_H__
#define __SLOPES_H__

#include <iostream>
#include <fstream>
#include <typeinfo>
#include <algorithm>
#include <map>

#include "precs.h"
#include "../bm/coords.h"
#include "../bm/blocksel.h"
#include "limits"
#include "slopedef.h"
using namespace std;


namespace delphos{ 


typedef struct relcoordindex 
    {
        long idx;
        long idy;
        long idz; 
    } RelativeCoordinatesIndex; 

typedef struct coordindex 
    {
        BlockIndexType idx;
        BlockIndexType idy;
        BlockIndexType idz; 
    } CoordinatesIndex; 
    
class SlopePrecedence : public Precedence 
{
private:
    vector<SlopeDefinition> _slope_defs;
    
    static vector<Coordinates>* RosettaPattern(const delphos::RosettaSlope& rosetta,const unsigned int & levels,const double & dx,const double & dy,const double & dz);
    
    static vector<Coordinates>* RosettaLightPattern(vector<delphos::Coordinates>* pattern1,vector<delphos::Coordinates>* pattern, const unsigned int & levels,const double & dz);
    
    
    map<double,vector<Coordinates>*>* RosettaLightPatterns(const map<double,delphos::RosettaSlope>& rosettas, const unsigned int & levels,const double & dx,const double & dy,const double & dz);

    map<double,vector<Coordinates>*>* RosettaPatterns(const map<double,delphos::RosettaSlope>& rosettas, const unsigned int & levels,const double & dx,const double & dy,const double & dz);
    
    
public:
    SlopePrecedence(delphos::BlockModel & bm);

    virtual void CreateArcs(const double & angle_degrees, const double & max_z_diff);

    virtual void CreateAdditionalArcs(const double & angle_degrees, const double & max_z_diff);


    virtual void CreateArcs(const double & angle_degrees, const double & max_z_diff,const delphos::BlockSelection& bs);
    
    virtual void CreateAdditionalArcs(const double & angle_degrees, const double & max_z_diff,const delphos::BlockSelection& bs);
    
    
    virtual void CreateArcs(
        const double & angle_degrees, 
        const unsigned int & levels,
        const double & dx,
        const double & dy,
        const double & dz);
    
    virtual void CreateAdditionalArcs(
        const double & angle_degrees, 
        const unsigned int & levels,
        const double & dx,
        const double & dy,
        const double & dz); 
    
    virtual void CreateAdditionalArcs(
        const double & angle_degrees, 
        const unsigned int & levels,
        const double & dx,
        const double & dy,
        const double & dz,
        const delphos::BlockSelection& bs);
    
    virtual void CreateArcs(
        const double & angle_degrees, 
        const unsigned int & levels,
        const double & dx,
        const double & dy,
        const double & dz,
        const delphos::BlockSelection& bs);
    
    virtual void CreateArcs(const delphos::SlopeDefinition& slope_def);
    
    virtual void CreateAdditionalArcs(const delphos::SlopeDefinition& slope_def);
    
    #ifndef DOXYGEN 
    virtual void CreateArcs(
        vector<Coordinates>& pattern1,
        vector<Coordinates>& pattern,
        const double & dx,
        const double & dy,
        const double & dz); 
    
    
    static vector<delphos::Coordinates> * SlopePattern(
        const double & angle_degrees, 
        const unsigned int & levels,
        const double & dx,
        const double & dy,
        const double & dz);
    
    static vector<delphos::Coordinates> * SlopeLightPattern(
        vector<delphos::Coordinates>* pattern1,
        vector<delphos::Coordinates>* pattern, 
        const unsigned int & levels,
        const double & dz);
        
    vector<delphos::RelativeCoordinatesIndex>* ToIndex(const double & dx,const double & dy,const double & dz,vector<delphos::Coordinates>* vcoords) const;  
    #endif  
};

}
#endif