Lattices/LatticeDirection.cpp

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_LatticeDirection_cpp_
#define model_LatticeDirection_cpp_
 
#include <LatticeModule.h>
namespace model
{
 
    template <int dim>
    LatticeDirection<dim>::LatticeDirection(const LatticeVectorType& v) :
    /* delegate */ LatticeDirection<dim>(v.base(),v.lattice)
    {
    }
 
    template <int dim>
    LatticeDirection<dim>::LatticeDirection(const VectorDimI& v,const Lattice<dim>& lat):
    /* base init */ LatticeVectorType(((v.squaredNorm()==0)? v : (v/LatticeGCD<dim>::gcd(v)).eval()),lat)
    {
        
    }
 
    template <int dim>
    LatticeDirection<dim>::LatticeDirection(const ReciprocalLatticeVectorType& r1,const ReciprocalLatticeVectorType& r2) :
    /* delegating */ LatticeDirection(LatticeVectorType(r1.cross(r2)))
    {
    }
 
    template <int dim>
    typename LatticeDirection<dim>::LatticeVectorType LatticeDirection<dim>::snapToLattice(const VectorDimD &dP) const
    { 
        const VectorDimD dc(this->cartesian());
        return LatticeVectorType((round(dP.dot(dc) / dc.squaredNorm()) * dc).eval(),
                                 this->lattice);
    }
 
    template <int dim>
    typename LatticeDirection<dim>::VectorDimD LatticeDirection<dim>::snapToDirection(const VectorDimD &dP) const
    {
        VectorDimD dc(this->cartesian());
        const double dNorm(dc.norm());
        assert(dNorm > FLT_EPSILON);
        dc /= dNorm;
        return dP.dot(dc) * dc;
    }
 
    template struct LatticeDirection<3>;
 
} // end namespace
#endif