ComputeGBMisorientationType.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#pragma once

#include "Material.h"
#include "GrainTracker.h"
#include "EBSDReader.h"

class ComputeGBMisorientationType : public Material
{
public:
  static InputParameters validParams();

  ComputeGBMisorientationType(const InputParameters & parameters);

protected:
  virtual void computeQpProperties() override;

  virtual unsigned int getTotalLineNum() const;
  virtual unsigned int getLineNum(unsigned int grain_i, unsigned int grain_j);
  virtual Real getTripleJunctionType();
  void rotationSymmetryToQuaternion(const double O[3][3], Eigen::Quaternion<Real> & q);
  void defineSymmetryOperator();
  double getMisorientationFromQuaternion(const Eigen::Quaternion<Real> qi,
                                         const Eigen::Quaternion<Real> qj);
  void getMisorientationAngles();

  const GrainTracker & _grain_tracker;

  const EBSDReader & _ebsd_reader;

  std::vector<double> _misorientation_angles;

  std::vector<unsigned int> _gb_pairs;
  std::vector<Real> _gb_op_pairs;

  const unsigned int _op_num;
  const std::vector<const VariableValue *> _vals;

  const Real _angle_threshold;

  std::vector<Eigen::Quaternion<Real>> _sym_quat;
  int _o_sym = 24;
  std::vector<EulerAngles> _euler_angle;
  std::vector<Eigen::Quaternion<Real>> _quat_angle;

  ADMaterialProperty<Real> & _gb_type;
};