Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : #include "OutputEulerAngles.h"
      11             : #include "GrainTracker.h"
      12             : #include "EulerAngleProvider.h"
      13             : 
      14             : registerMooseObject("PhaseFieldApp", OutputEulerAngles);
      15             : 
      16             : InputParameters
      17         253 : OutputEulerAngles::validParams()
      18             : {
      19         253 :   InputParameters params = AuxKernel::validParams();
      20         253 :   params.addClassDescription("Output Euler angles from user object to an AuxVariable.");
      21         506 :   params.addRequiredParam<UserObjectName>("euler_angle_provider",
      22             :                                           "Name of Euler angle provider user object");
      23         506 :   params.addRequiredParam<UserObjectName>("grain_tracker",
      24             :                                           "The GrainTracker UserObject to get values from.");
      25         506 :   MooseEnum euler_angles("phi1 Phi phi2");
      26         506 :   params.addRequiredParam<MooseEnum>("output_euler_angle", euler_angles, "Euler angle to output");
      27         253 :   return params;
      28         253 : }
      29             : 
      30         132 : OutputEulerAngles::OutputEulerAngles(const InputParameters & parameters)
      31             :   : AuxKernel(parameters),
      32         132 :     _euler(getUserObject<EulerAngleProvider>("euler_angle_provider")),
      33         132 :     _grain_tracker(getUserObject<GrainTracker>("grain_tracker")),
      34         396 :     _output_euler_angle(getParam<MooseEnum>("output_euler_angle"))
      35             : {
      36         132 : }
      37             : 
      38             : void
      39      949842 : OutputEulerAngles::precalculateValue()
      40             : {
      41             :   // ID of unique grain at current point
      42             :   const auto grain_id =
      43      949842 :       _grain_tracker.getEntityValue((isNodal() ? _current_node->id() : _current_elem->id()),
      44             :                                     FeatureFloodCount::FieldType::UNIQUE_REGION,
      45             :                                     0);
      46             : 
      47             :   // Recover euler angles for current grain
      48             :   RealVectorValue angles;
      49      949842 :   if (grain_id >= 0)
      50      923188 :     angles = _euler.getEulerAngles(grain_id);
      51             : 
      52             :   // Return specific euler angle
      53      949842 :   _value = angles(_output_euler_angle);
      54      949842 : }
      55             : 
      56             : Real
      57     3870712 : OutputEulerAngles::computeValue()
      58             : {
      59     3870712 :   return _value;
      60             : }