LatticeGCD.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_LatticeGCD_cpp_
#define model_LatticeGCD_cpp_
 
#include <LatticeGCD.h>
#include <Eigen/Dense>
 
namespace model
{
 
    template <int dim>
    long int LatticeGCD<dim>::gcd(const long int &a, const long int &b)
    {
        const long int absA(abs(a));
        const long int absB(abs(b));
        return absB > 0 ? gcd(absB, absA % absB) : (absA > 0 ? absA : 1);
    }
 
    template <int dim>
    long int LatticeGCD<dim>::gcd(const long int &a, const long int &b, const long int &c)
    {
        return gcd(a, gcd(b, c));
    }
 
    template <int dim>
    long int LatticeGCD<dim>::gcd(const VectorDimI &v)
    {
        return gcd(v(0), v(1), v(2));
    }
 
    template struct LatticeGCD<3>;
} // end namespace
#endif