Line data    Source code

       1             : /**********************************************************************/
       2             : /*                     DO NOT MODIFY THIS HEADER                      */
       3             : /* MAGPIE - Mesoscale Atomistic Glue Program for Integrated Execution */
       4             : /*                                                                    */
       5             : /*            Copyright 2017 Battelle Energy Alliance, LLC            */
       6             : /*                        ALL RIGHTS RESERVED                         */
       7             : /**********************************************************************/
       8             : 
       9             : #include "RadialGreensAux.h"
      10             : #include "RadialGreensConvolution.h"
      11             : 
      12             : registerMooseObject("MagpieApp", RadialGreensAux);
      13             : 
      14             : InputParameters
      15          18 : RadialGreensAux::validParams()
      16             : {
      17          18 :   InputParameters params = AuxKernel::validParams();
      18          18 :   params.addClassDescription("Visualize data generated in a RadialGreensConvolution user object");
      19          36 :   params.addRequiredParam<UserObjectName>("convolution", "RadialGreensConvolution user object");
      20          18 :   return params;
      21           0 : }
      22             : 
      23          10 : RadialGreensAux::RadialGreensAux(const InputParameters & parameters)
      24             :   : AuxKernel(parameters),
      25          10 :     _convolution(getUserObject<RadialGreensConvolution>("convolution").getConvolution())
      26             : {
      27          10 :   if (isNodal())
      28           0 :     paramError("variable", "RadialGreensAux must be applied to an elemental AuxVariable");
      29          10 : }
      30             : 
      31             : Real
      32      179200 : RadialGreensAux::computeValue()
      33             : {
      34      179200 :   auto it = _convolution.find(_current_elem->id());
      35      179200 :   if (it != _convolution.end())
      36      179200 :     return it->second[_qp] / (_JxW[_qp] * _coord[_qp]);
      37             : 
      38             :   return 0.0;
      39             : }