Lattices/Lattice.h

Go to the documentation of this file.

/* This file is part of PIL, the Particle Interaction Library.
 *
 *
 * MODEL is distributed without any warranty under the
 * GNU General Public License (GPL) v2 <http://www.gnu.org/licenses/>.
 */
 
#ifndef model_Lattice_h_
#define model_Lattice_h_
 
#include <iomanip> // FLT_EPSILON
#include <cfloat> // FLT_EPSILON
#include <Eigen/Dense>
#include <StaticID.h>
#include <LatticeVector.h>
#include <ReciprocalLatticeVector.h>
#include <LatticeDirection.h>
#include <ReciprocalLatticeDirection.h>
#include <BestRationalApproximation.h>
 
#include <RationalLatticeDirection.h>
 
namespace model
{
 
    template <int dim>
    class Lattice : public StaticID<Lattice<dim>>
    {
        
    public:
        static constexpr double roundTol=FLT_EPSILON;
        typedef Eigen::Matrix<  double,dim,1> VectorDimD;
        typedef Eigen::Matrix<long int,dim,1> VectorDimI;
        typedef Eigen::Matrix<double,dim,dim> MatrixDimD;
        typedef LatticeVector<dim> LatticeVectorType;
        typedef LatticeDirection<dim> LatticeDirectionType;
        typedef RationalLatticeDirection<dim> RationalLatticeDirectionType;
        typedef ReciprocalLatticeVector<dim> ReciprocalLatticeVectorType;
        typedef ReciprocalLatticeDirection<dim> ReciprocalLatticeDirectionType;
 
        static MatrixDimD getLatticeBasis(const MatrixDimD& A,const MatrixDimD& Q);
 
        const MatrixDimD    latticeBasis;
        const MatrixDimD reciprocalBasis;
        const MatrixDimD C2G;
        
        Lattice(const MatrixDimD& A,const MatrixDimD& Q) ;
        static Eigen::Matrix<long int,dim,1> rationalApproximation(VectorDimD nd);
        LatticeVectorType snapToLattice(const VectorDimD& d) const;
        LatticeDirectionType latticeDirection(const VectorDimD& d) const;
        ReciprocalLatticeDirectionType reciprocalLatticeDirection(const VectorDimD& d) const;
        RationalLatticeDirectionType rationalLatticeDirection(const VectorDimD& d,
                                                              const typename BestRationalApproximation::LongIntType& maxDen=1000) const;
        LatticeVectorType latticeVector(const VectorDimD& p) const;
        ReciprocalLatticeVectorType reciprocalLatticeVector(const VectorDimD& p) const;
    };
}
#endif