BiCrystal.h

Go to the documentation of this file.

/* This file is part of gbLAB.
 *
 * gbLAB is distributed without any warranty under the MIT License.
 */
 
 
#ifndef gbLAB_BiCrystal_h_
#define gbLAB_BiCrystal_h_
 
#include "../Math/LLL.h"
#include "../Math/RLLL.h"
#include "../Math/RationalMatrix.h"
#include "../Math/SmithDecomposition.h"
#include "../Utilities/Rotation.h"
#include "LatticeModule.h"
#include <unordered_set>
 
namespace oILAB {
    template <int dim>
    class BiCrystal : public RationalMatrix<dim>
    /*             */,public SmithDecomposition<dim>
    {
        typedef typename LatticeCore<dim>::IntScalarType IntScalarType;
        typedef typename LatticeCore<dim>::VectorDimD VectorDimD;
        typedef typename LatticeCore<dim>::MatrixDimD MatrixDimD;
        typedef typename LatticeCore<dim>::VectorDimI VectorDimI;
        typedef typename LatticeCore<dim>::MatrixDimI MatrixDimI;
 
        
        static MatrixDimI getM(const RationalMatrix<dim>& rm, const SmithDecomposition<dim>& sd);
        static MatrixDimI getN(const RationalMatrix<dim>& rm, const SmithDecomposition<dim>& sd);
        static MatrixDimI getLambdaA(const MatrixDimI& M, const MatrixDimI& N);
        static MatrixDimI getLambdaB(const MatrixDimI& M, const MatrixDimI& N);
        static MatrixDimD getCSLBasis(const Lattice<dim>& A,
                                      const Lattice<dim>& B,
                                      const SmithDecomposition<dim>& sd,
                                      const MatrixDimI& M,
                                      const MatrixDimI& N);
        static MatrixDimD getDSCLBasis(const Lattice<dim>& A,
                                       const Lattice<dim>& B,
                                       const SmithDecomposition<dim>& sd,
                                       const MatrixDimI& M,
                                       const MatrixDimI& N);
 
    public:
 
        const Lattice<dim>& A;
        const Lattice<dim>& B;
 
        const MatrixDimI M;
 
        const MatrixDimI N;
 
        const IntScalarType sigmaA;
 
        const IntScalarType sigmaB;
 
        const int sigma;
 
    private:
        const Lattice<dim> cslp;
 
        const Lattice<dim> dsclp;
 
        const MatrixDimI dscl2dsclp;
 
        const MatrixDimI csl2cslp;
 
    public:
        const Lattice<dim> csl;
 
        const Lattice<dim> dscl;
 
        const Lattice<dim> Ap;
 
        const Lattice<dim> Bp;
 
        const MatrixDimI LambdaA;
 
        const MatrixDimI LambdaB;
 
 
        LatticeVector<dim> shiftTensorA(const LatticeVector<dim>& d) const;
        LatticeVector<dim> shiftTensorB(const LatticeVector<dim>& d) const;
 
        /**********************************************************************/
        BiCrystal(const Lattice<dim>& A,
                  const Lattice<dim>& B,
                  const bool& useRLLL=false);
 
        LatticeVector<dim> getLatticeVectorInA(const LatticeVector<dim>& v) const;
        LatticeVector<dim> getLatticeVectorInB(const LatticeVector<dim>& v) const;
        LatticeVector<dim> getLatticeVectorInD(const LatticeVector<dim>& v) const;
        LatticeDirection<dim> getLatticeDirectionInC(const LatticeVector<dim>& v) const;
        LatticeDirection<dim> getLatticeDirectionInD(const LatticeVector<dim>& v) const;
 
        ReciprocalLatticeDirection<dim> getReciprocalLatticeDirectionInA(const ReciprocalLatticeVector<dim>& v) const;
        ReciprocalLatticeDirection<dim> getReciprocalLatticeDirectionInB(const ReciprocalLatticeVector<dim>& v) const;
        ReciprocalLatticeDirection<dim> getReciprocalLatticeDirectionInC(const ReciprocalLatticeVector<dim>& v) const;
        ReciprocalLatticeDirection<dim> getReciprocalLatticeDirectionInD(const ReciprocalLatticeVector<dim>& v) const;
 
        template<int dm=dim>
        typename std::enable_if<dm==2 || dm==3,std::map<IntScalarType,Gb<dm>>>::type
        generateGrainBoundaries(const LatticeDirection<dim>& d, int div=30) const;
 
 
        template<int dm=dim>
        typename std::enable_if<dm==2 || dm==3,std::vector<LatticeVector<dim>>>::type
        box(std::vector<LatticeVector<dim>>& boxVectors, 
                const double& orthogonality, 
                const int& dsclFactor,
                std::string filename= "", 
                bool orient=false) const;
    };
    
    
} // end namespace
#endif