This UserObject computes a volumes and centers of grains. More...

#include <ComputeGrainCenterUserObject.h>

Inheritance diagram for ComputeGrainCenterUserObject:

Public Member Functions
	ComputeGrainCenterUserObject (const InputParameters &parameters)

virtual void	initialize ()

virtual void	execute ()

virtual void	threadJoin (const UserObject &y)

virtual void	finalize ()

const std::vector< Real > &	getGrainVolumes () const

const std::vector< Point > &	getGrainCenters () const

Protected Attributes
unsigned int	_qp

unsigned int	_ncrys

std::vector< const VariableValue * >	_vals

unsigned int	_ncomp


std::vector< Real >	_grain_data
	storing volumes and centers of all the grains More...

std::vector< Real >	_grain_volumes

std::vector< Point >	_grain_centers

Detailed Description

This UserObject computes a volumes and centers of grains.

Definition at line 23 of file ComputeGrainCenterUserObject.h.

Constructor & Destructor Documentation

◆ ComputeGrainCenterUserObject()

ComputeGrainCenterUserObject::ComputeGrainCenterUserObject ( const InputParameters & parameters )

Definition at line 24 of file ComputeGrainCenterUserObject.C.

   : ElementUserObject(parameters),
     _ncrys(coupledComponents("etas")), // determine number of grains from the number of names passed
                                        // in.  Note this is the actual number -1
     _vals(_ncrys),                     // Size variable arrays
     _ncomp(4 * _ncrys),
     _grain_data(_ncomp),
     _grain_volumes(_ncrys),
     _grain_centers(_ncrys)
 {
   for (unsigned int i = 0; i < _ncrys; ++i)
     _vals[i] = &coupledValue("etas", i);
 }

Member Function Documentation

◆ execute()

void ComputeGrainCenterUserObject::execute ( )

virtual

Definition at line 46 of file ComputeGrainCenterUserObject.C.

 {
   for (unsigned int i = 0; i < _ncrys; ++i)
     for (_qp = 0; _qp < _qrule->n_points(); ++_qp)
     {
       _grain_data[4 * i + 0] += _JxW[_qp] * _coord[_qp] * (*_vals[i])[_qp];
       _grain_data[4 * i + 1] += _JxW[_qp] * _coord[_qp] * _q_point[_qp](0) * (*_vals[i])[_qp];
       _grain_data[4 * i + 2] += _JxW[_qp] * _coord[_qp] * _q_point[_qp](1) * (*_vals[i])[_qp];
       _grain_data[4 * i + 3] += _JxW[_qp] * _coord[_qp] * _q_point[_qp](2) * (*_vals[i])[_qp];
     }
 }

◆ finalize()

void ComputeGrainCenterUserObject::finalize ( )

virtual

Definition at line 59 of file ComputeGrainCenterUserObject.C.

 {
   gatherSum(_grain_data);
  
   for (unsigned int i = 0; i < _ncrys; ++i)
   {
     _grain_volumes[i] = _grain_data[4 * i + 0];
     _grain_centers[i](0) = _grain_data[4 * i + 1] / _grain_volumes[i];
     _grain_centers[i](1) = _grain_data[4 * i + 2] / _grain_volumes[i];
     _grain_centers[i](2) = _grain_data[4 * i + 3] / _grain_volumes[i];
   }
 }

◆ getGrainCenters()

const std::vector< Point > & ComputeGrainCenterUserObject::getGrainCenters ( ) const

Definition at line 87 of file ComputeGrainCenterUserObject.C.

 {
   return _grain_centers;
 }

◆ getGrainVolumes()

const std::vector< Real > & ComputeGrainCenterUserObject::getGrainVolumes ( ) const

Definition at line 81 of file ComputeGrainCenterUserObject.C.

 {
   return _grain_volumes;
 }

◆ initialize()

void ComputeGrainCenterUserObject::initialize ( )

virtual

Definition at line 39 of file ComputeGrainCenterUserObject.C.

 {
   for (unsigned int i = 0; i < _ncomp; ++i)
     _grain_data[i] = 0;
 }

◆ threadJoin()

void ComputeGrainCenterUserObject::threadJoin ( const UserObject & y )

virtual

Definition at line 73 of file ComputeGrainCenterUserObject.C.

 {
   const ComputeGrainCenterUserObject & pps = static_cast<const ComputeGrainCenterUserObject &>(y);
   for (unsigned int i = 0; i < _ncomp; ++i)
     _grain_data[i] += pps._grain_data[i];
 }