hexagonal_packing.h

Go to the documentation of this file.

#ifndef __HEXAGONALPACKING_H__
#define __HEXAGONALPACKING_H__

#include <cstddef>
#include <cmath>


class HexagonalPacking2D {
public:
    HexagonalPacking2D(double xmin_, double xmax_, double ymin_, double ymax_, double h_r_);
    
    void getParameters(size_t& m0_, size_t& m1_, 
                       size_t& n0_odd_, size_t& n1_odd_, 
                       size_t& n0_even_, size_t& n1_even_) {
        m0_ = m0; 
        m1_ = m1;
        n0_odd_ = n0_odd; 
        n1_odd_ = n1_odd;
        n0_even_ = n0_even;
        n1_even_ = n1_even;
    }
    
    double computeX(int tag, size_t k) {
        if(tag==0) return xmin+2.0*h_r+k*2.0*h_r; //odd-numbered rows
        else return  xmin+h_r+k*2.0*h_r; //even-numbered rows
    }
    
    double computeY(size_t j) {
        return ymin+h_r+j*sqrt(3.0)*h_r;
    }

private:
    double xmin, xmax, ymin, ymax, h_r;
    size_t m0, m1, n0_odd, n1_odd, n0_even, n1_even;        
    //HexagonalPacking2D() {}
    //HexagonalPacking2D(const HexagonalPacking2D& other) {}
    //HexagonalPacking2D& operator=(const HexagonalPacking2D& other) {return *this;}
};










class HexagonalPacking3D {
public:
    HexagonalPacking3D(double xmin_, double xmax_, double ymin_, double ymax_, double zmin_, double zmax_, double h_r_);
    
    void getParameters(size_t& l0_, size_t& l1_, 
                       size_t& m0_odd_, size_t& m1_odd_, 
                       size_t& m0_even_, size_t& m1_even_,
                       size_t& n0_odd_, size_t& n1_odd_,
                       size_t& n0_even_, size_t& n1_even_,
                       size_t& nn0_odd_, size_t& nn1_odd_,
                       size_t& nn0_even_, size_t& nn1_even_) {
        l0_ = l0; 
        l1_ = l1; 
        m0_odd_ = m0_odd; 
        m1_odd_ = m1_odd; 
        m0_even_ = m0_even; 
        m1_even_ = m1_even;
        n0_odd_ = n0_odd; 
        n1_odd_ = n1_odd;
        n0_even_ = n0_even; 
        n1_even_ = n1_even;
        nn0_odd_ = nn0_odd; 
        nn1_odd_ = nn1_odd;
        nn0_even_ = nn0_even; 
        nn1_even_ = nn1_even;
    }
    

    double computeX(int tag, size_t k) {
        if(tag==0) return xmin+2.0*h_r+k*2.0*h_r; //(odd layers && odd rows) || (even layers && even rows)
        else return  xmin+h_r+k*2.0*h_r; // (odd layers && even rows) || (even layers && odd rows)
    }

    double computeY(int tag, size_t j) {
        if(tag==0) return ymin+h_r+j*sqrt(3.0)*h_r; //odd-numbered layers
        else return ymin+sqrt(3.0)/3.0*h_r+j*sqrt(3.0)*h_r; //even-numbered layers
    }


    double computeZ(size_t i) {
        return zmin+h_r+i*2.0*sqrt(6.0)/3.0*h_r;
    }
    
private:
    
    double xmin, xmax, ymin, ymax, zmin, zmax, h_r;
    size_t l0,l1;
    size_t m0_odd, m1_odd, m0_even, m1_even, n0_odd, n1_odd, n0_even, n1_even;
    size_t nn0_odd, nn1_odd, nn0_even, nn1_even;        
    //HexagonalPacking3D() {}
    //HexagonalPacking3D(const HexagonalPacking3D& other) {}
    //HexagonalPacking3D& operator=(const HexagonalPacking3D& other) {return *this;}
};


#endif // __HEXAGONALPACKING_H__