EBSDReader.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 "EulerAngleProvider.h"
#include "EBSDAccessFunctors.h"
 
class EBSDReader : public EulerAngleProvider, public EBSDAccessFunctors
{
public:
  static InputParameters validParams();
 
  EBSDReader(const InputParameters & params);
  virtual ~EBSDReader();
 
  virtual void readFile();
 
  virtual void initialize() {}
  virtual void execute() {}
  virtual void finalize() {}
 
  const EBSDPointData & getData(const Point & p) const;
 
  const EBSDAvgData & getAvgData(unsigned int i) const;
 
  const EBSDAvgData & getAvgData(unsigned int phase, unsigned int local_id) const;
 
  virtual const EulerAngles & getEulerAngles(unsigned int) const;
 
  virtual unsigned int getGrainNum() const;
 
  virtual unsigned int getPhaseNum() const { return _global_id.size(); }
 
  virtual unsigned int getGrainNum(unsigned int phase) const;
 
  unsigned int getFeatureID(unsigned int phase, unsigned int local_id) const
  {
    return _avg_data[_global_id[phase][local_id]]._feature_id;
  }
  unsigned int getFeatureID(unsigned int global_id) const
  {
    return _avg_data[global_id]._feature_id;
  }
 
  virtual unsigned int getGlobalID(unsigned int phase, unsigned int local_id) const
  {
    return _global_id[phase][local_id];
  }
  virtual unsigned int getGlobalID(unsigned int feature_id) const;
 
  MooseSharedPointer<EBSDPointDataFunctor>
  getPointDataAccessFunctor(const MooseEnum & field_name) const;
  MooseSharedPointer<EBSDAvgDataFunctor>
  getAvgDataAccessFunctor(const MooseEnum & field_name) const;
 
  const std::map<dof_id_type, std::vector<Real>> & getNodeToGrainWeightMap() const;
 
  const std::map<dof_id_type, std::vector<Real>> & getNodeToPhaseWeightMap() const;
 
  void meshChanged();
 
protected:
  MooseMesh & _mesh;
  NonlinearSystemBase & _nl;
 
  unsigned int _grain_num;
  Point _bottom_left;
  Point _top_right;
  Point _range;
 
  unsigned int _custom_columns;
 
  std::vector<EBSDPointData> _data;
 
  std::vector<EBSDAvgData> _avg_data;
 
  std::vector<EulerAngles> _avg_angles;
 
  std::map<unsigned int, unsigned int> _global_id_map;
 
  std::vector<std::vector<unsigned int>> _global_id;
 
  std::map<dof_id_type, std::vector<Real>> _node_to_grain_weight_map;
 
  std::map<dof_id_type, std::vector<Real>> _node_to_phase_weight_map;
 
  const int & _time_step;
 
  unsigned int _mesh_dimension;
 
  unsigned int _bins;
 
  unsigned int _L_norm;
 
  unsigned _nx, _ny, _nz;
 
  Real _dx, _dy, _dz;
 
  Real _minx, _miny, _minz;
 
  Real _maxx, _maxy, _maxz;
 
  unsigned indexFromPoint(const Point & p) const;
 
  unsigned indexFromIndex(unsigned int var) const;
 
  void buildNodeWeightMaps();
};