PiecewiseConstantByBlockMaterial.C

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
 
#include "PiecewiseConstantByBlockMaterial.h"
 
registerMooseObject("MooseApp", PiecewiseConstantByBlockMaterial);
registerMooseObject("MooseApp", ADPiecewiseConstantByBlockMaterial);
 
template <bool is_ad>
InputParameters
PiecewiseConstantByBlockMaterialTempl<is_ad>::validParams()
{
  auto params = Material::validParams();
  params.addClassDescription("Computes a property value on a per-subdomain basis");
  // Somehow min gcc doesn't know the type of params here
  params.template addRequiredParam<MaterialPropertyName>("prop_name",
                                                         "The name of the property to declare");
  params.template addRequiredParam<std::map<std::string, Real>>(
      "subdomain_to_prop_value", "Map from subdomain to property value");
  return params;
}
 
template <bool is_ad>
PiecewiseConstantByBlockMaterialTempl<is_ad>::PiecewiseConstantByBlockMaterialTempl(
    const InputParameters & params)
  : Material(params), _prop(declareGenericProperty<Real, is_ad>("prop_name"))
{
  for (const auto & map_pr : getParam<std::map<std::string, Real>>("subdomain_to_prop_value"))
    _sub_id_to_prop.emplace(std::make_pair(_mesh.getSubdomainID(map_pr.first), map_pr.second));
}
 
template <bool is_ad>
void
PiecewiseConstantByBlockMaterialTempl<is_ad>::computeQpProperties()
{
  mooseAssert(_current_elem,
              "We should be on a block which means we should definitely have a current element");
  auto it = _sub_id_to_prop.find(_current_elem->subdomain_id());
  mooseAssert(it != _sub_id_to_prop.end(),
              "Block restriction must match the subdomain names passed in the "
              "subdomain_to_prop_value parameter");
  _prop[_qp] = it->second;
}
 
template class PiecewiseConstantByBlockMaterialTempl<false>;
template class PiecewiseConstantByBlockMaterialTempl<true>;