Gb.h

Go to the documentation of this file.

//
// Created by Nikhil Chandra Admal on 11/5/22.
//
 
#ifndef OILAB_GB_H
#define OILAB_GB_H
 
#include "../Utilities/Rotation.h"
#include "BiCrystal.h"
#include "ReciprocalLatticeDirection.h"
 
namespace oILAB {
template <int dim> class Gb {
  using VectorDimI = typename LatticeCore<dim>::VectorDimI;
  using VectorDimD = typename LatticeCore<dim>::VectorDimD;
  using MatrixDimD = typename LatticeCore<dim>::MatrixDimD;
  using MatrixDimI = typename LatticeCore<dim>::MatrixDimI;
  using IntScalarType = typename LatticeCore<dim>::IntScalarType;
 
private:
  MatrixDimI getBasisT(const BiCrystal<dim> &bc,
                       const ReciprocalLatticeDirection<dim> &n);
 
public:
  const BiCrystal<dim> &bc;
  const ReciprocalLatticeDirection<dim> nA;
  const ReciprocalLatticeDirection<dim> nB;
 
  const Lattice<dim> T;
 
  const MatrixDimI basisT;
 
  LatticeVector<dim> getLatticeVectorInT(const LatticeVector<dim> &v) const;
 
  ReciprocalLatticeVector<dim>
  getReciprocalLatticeVectorInT(const ReciprocalLatticeVector<dim> &v) const;
 
  ReciprocalLatticeDirection<dim>
  getReciprocalLatticeDirectionInT(const ReciprocalLatticeVector<dim> &v) const;
 
  double stepHeightA(const LatticeVector<dim> &d) const;
  double stepHeightB(const LatticeVector<dim> &d) const;
 
  double stepHeight(const LatticeVector<dim> &d) const;
 
  Gb(const BiCrystal<dim> &bc, const ReciprocalLatticeDirection<dim> &n);
 
  template <int dm = dim>
  typename std::enable_if<dm == 2, LatticeVector<dim>>::type
  getPeriodVector(const ReciprocalLatticeVector<dim> &axis) const;
 
  template <int dm = dim>
  typename std::enable_if<dm == 3, LatticeVector<dim>>::type
  getPeriodVector(const ReciprocalLatticeVector<dm> &axis) 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;
 
    };
 
    } // namespace oILAB
 
#endif //OILAB_GB_H