#include <AverageSectionValueSampler.h>

Inheritance diagram for AverageSectionValueSampler:

Public Types
typedef DataFileName	DataFileParameterType

Public Member Functions
	AverageSectionValueSampler (const InputParameters &parameters)

virtual void	initialSetup () override

virtual void	meshChanged () override

virtual void	initialize () override

virtual void	execute () override

virtual void	finalize () override

SubProblem &	getSubProblem () const

bool	shouldDuplicateInitialExecution () const

virtual Real	spatialValue (const Point &) const

virtual const std::vector< Point >	spatialPoints () const

void	gatherSum (T &value)

void	gatherMax (T &value)

void	gatherMin (T &value)

void	gatherProxyValueMax (T1 &proxy, T2 &value)

void	gatherProxyValueMin (T1 &proxy, T2 &value)

void	setPrimaryThreadCopy (UserObject *primary)

UserObject *	primaryThreadCopy ()

std::set< UserObjectName >	getDependObjects () const

virtual bool	needThreadedCopy () const

const std::set< std::string > &	getRequestedItems () override

const std::set< std::string > &	getSuppliedItems () override

unsigned int	systemNumber () const

virtual bool	enabled () const

std::shared_ptr< MooseObject >	getSharedPtr ()

std::shared_ptr< const MooseObject >	getSharedPtr () const

MooseApp &	getMooseApp () const

const std::string &	type () const

virtual const std::string &	name () const

std::string	typeAndName () const

std::string	errorPrefix (const std::string &error_type) const

void	callMooseError (std::string msg, const bool with_prefix) const

MooseObjectParameterName	uniqueParameterName (const std::string &parameter_name) const

const InputParameters &	parameters () const

MooseObjectName	uniqueName () const

const T &	getParam (const std::string &name) const

std::vector< std::pair< T1, T2 > >	getParam (const std::string &param1, const std::string &param2) const

const T *	queryParam (const std::string &name) const

const T &	getRenamedParam (const std::string &old_name, const std::string &new_name) const

T	getCheckedPointerParam (const std::string &name, const std::string &error_string="") const

bool	isParamValid (const std::string &name) const

bool	isParamSetByUser (const std::string &nm) const

void	paramError (const std::string &param, Args... args) const

void	paramWarning (const std::string &param, Args... args) const

void	paramInfo (const std::string &param, Args... args) const

void	connectControllableParams (const std::string &parameter, const std::string &object_type, const std::string &object_name, const std::string &object_parameter) const

void	mooseError (Args &&... args) const

void	mooseErrorNonPrefixed (Args &&... args) const

void	mooseDocumentedError (const std::string &repo_name, const unsigned int issue_num, Args &&... args) const

void	mooseWarning (Args &&... args) const

void	mooseWarningNonPrefixed (Args &&... args) const

void	mooseDeprecated (Args &&... args) const

void	mooseInfo (Args &&... args) const

std::string	getDataFileName (const std::string &param) const

std::string	getDataFileNameByName (const std::string &relative_path) const

std::string	getDataFilePath (const std::string &relative_path) const

virtual void	timestepSetup ()

virtual void	jacobianSetup ()

virtual void	residualSetup ()

virtual void	customSetup (const ExecFlagType &)

const ExecFlagEnum &	getExecuteOnEnum () const

UserObjectName	getUserObjectName (const std::string &param_name) const

const T &	getUserObject (const std::string &param_name, bool is_dependency=true) const

const T &	getUserObjectByName (const UserObjectName &object_name, bool is_dependency=true) const

const UserObject &	getUserObjectBase (const std::string &param_name, bool is_dependency=true) const

const UserObject &	getUserObjectBaseByName (const UserObjectName &object_name, bool is_dependency=true) const

const std::vector< MooseVariableScalar *> &	getCoupledMooseScalarVars ()

const std::set< TagID > &	getScalarVariableCoupleableVectorTags () const

const std::set< TagID > &	getScalarVariableCoupleableMatrixTags () const

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialProperty (const std::string &name, const unsigned int state=0)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialProperty (const std::string &name, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialProperty (const std::string &name, MaterialData &material_data, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialProperty (const std::string &name, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialProperty (const std::string &name, const unsigned int state=0)

const ADMaterialProperty< T > &	getADMaterialProperty (const std::string &name, MaterialData &material_data)

const ADMaterialProperty< T > &	getADMaterialProperty (const std::string &name)

const ADMaterialProperty< T > &	getADMaterialProperty (const std::string &name)

const MaterialProperty< T > &	getMaterialPropertyOld (const std::string &name, MaterialData &material_data)

const MaterialProperty< T > &	getMaterialPropertyOld (const std::string &name)

const MaterialProperty< T > &	getMaterialPropertyOld (const std::string &name)

const MaterialProperty< T > &	getMaterialPropertyOlder (const std::string &name, MaterialData &material_data)

const MaterialProperty< T > &	getMaterialPropertyOlder (const std::string &name)

const MaterialProperty< T > &	getMaterialPropertyOlder (const std::string &name)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)

const GenericMaterialProperty< T, is_ad > &	getGenericMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)

const MaterialProperty< T > &	getMaterialPropertyByName (const MaterialPropertyName &name, const unsigned int state=0)

const ADMaterialProperty< T > &	getADMaterialPropertyByName (const MaterialPropertyName &name, MaterialData &material_data)

const ADMaterialProperty< T > &	getADMaterialPropertyByName (const MaterialPropertyName &name)

const ADMaterialProperty< T > &	getADMaterialPropertyByName (const MaterialPropertyName &name)

const MaterialProperty< T > &	getMaterialPropertyOldByName (const MaterialPropertyName &name, MaterialData &material_data)

const MaterialProperty< T > &	getMaterialPropertyOldByName (const MaterialPropertyName &name)

const MaterialProperty< T > &	getMaterialPropertyOldByName (const MaterialPropertyName &name)

const MaterialProperty< T > &	getMaterialPropertyOlderByName (const MaterialPropertyName &name, MaterialData &material_data)

const MaterialProperty< T > &	getMaterialPropertyOlderByName (const MaterialPropertyName &name)

const MaterialProperty< T > &	getMaterialPropertyOlderByName (const MaterialPropertyName &name)

std::pair< const MaterialProperty< T > *, std::set< SubdomainID > >	getBlockMaterialProperty (const MaterialPropertyName &name)

const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty (const std::string &name)

const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialProperty ()

const GenericMaterialProperty< T, is_ad > &	getGenericZeroMaterialPropertyByName (const std::string &prop_name)

const MaterialProperty< T > &	getZeroMaterialProperty (Ts... args)

std::set< SubdomainID >	getMaterialPropertyBlocks (const std::string &name)

std::vector< SubdomainName >	getMaterialPropertyBlockNames (const std::string &name)

std::set< BoundaryID >	getMaterialPropertyBoundaryIDs (const std::string &name)

std::vector< BoundaryName >	getMaterialPropertyBoundaryNames (const std::string &name)

void	checkBlockAndBoundaryCompatibility (std::shared_ptr< MaterialBase > discrete)

std::unordered_map< SubdomainID, std::vector< MaterialBase *> >	buildRequiredMaterials (bool allow_stateful=true)

void	statefulPropertiesAllowed (bool)

bool	getMaterialPropertyCalled () const

virtual const std::unordered_set< unsigned int > &	getMatPropDependencies () const

virtual void	resolveOptionalProperties ()

const GenericMaterialProperty< T, is_ad > &	getPossiblyConstantGenericMaterialPropertyByName (const MaterialPropertyName &prop_name, MaterialData &material_data, const unsigned int state)

bool	isImplicit ()

Moose::StateArg	determineState () const

std::string	PPName () const

bool	containsCompleteHistory () const

bool	isDistributed () const

const std::set< std::string > &	getVectorNames () const

void	buildOutputHideVariableList (std::set< std::string > variable_names)

const std::set< OutputName > &	getOutputs ()

virtual void	threadJoin (const UserObject &) override

virtual void	threadJoin (const UserObject &) override

virtual void	subdomainSetup () override

virtual void	subdomainSetup () override

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObject (const std::string &param_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

bool	hasUserObjectByName (const UserObjectName &object_name) const

const GenericOptionalMaterialProperty< T, is_ad > &	getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)

const GenericOptionalMaterialProperty< T, is_ad > &	getGenericOptionalMaterialProperty (const std::string &name, const unsigned int state=0)

const OptionalMaterialProperty< T > &	getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)

const OptionalMaterialProperty< T > &	getOptionalMaterialProperty (const std::string &name, const unsigned int state=0)

const OptionalADMaterialProperty< T > &	getOptionalADMaterialProperty (const std::string &name)

const OptionalADMaterialProperty< T > &	getOptionalADMaterialProperty (const std::string &name)

const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOld (const std::string &name)

const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOld (const std::string &name)

const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOlder (const std::string &name)

const OptionalMaterialProperty< T > &	getOptionalMaterialPropertyOlder (const std::string &name)

MaterialBase &	getMaterial (const std::string &name)

MaterialBase &	getMaterial (const std::string &name)

MaterialBase &	getMaterialByName (const std::string &name, bool no_warn=false)

MaterialBase &	getMaterialByName (const std::string &name, bool no_warn=false)

bool	hasMaterialProperty (const std::string &name)

bool	hasMaterialProperty (const std::string &name)

bool	hasMaterialPropertyByName (const std::string &name)

bool	hasMaterialPropertyByName (const std::string &name)

bool	hasADMaterialProperty (const std::string &name)

bool	hasADMaterialProperty (const std::string &name)

bool	hasADMaterialPropertyByName (const std::string &name)

bool	hasADMaterialPropertyByName (const std::string &name)

bool	hasGenericMaterialProperty (const std::string &name)

bool	hasGenericMaterialProperty (const std::string &name)

bool	hasGenericMaterialPropertyByName (const std::string &name)

bool	hasGenericMaterialPropertyByName (const std::string &name)

const Function &	getFunction (const std::string &name) const

const Function &	getFunctionByName (const FunctionName &name) const

bool	hasFunction (const std::string &param_name) const

bool	hasFunctionByName (const FunctionName &name) const

bool	isDefaultPostprocessorValue (const std::string &param_name, const unsigned int index=0) const

bool	hasPostprocessor (const std::string &param_name, const unsigned int index=0) const

bool	hasPostprocessorByName (const PostprocessorName &name) const

std::size_t	coupledPostprocessors (const std::string &param_name) const

const PostprocessorName &	getPostprocessorName (const std::string &param_name, const unsigned int index=0) const

const VectorPostprocessorValue &	getVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const

const VectorPostprocessorValue &	getVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const

const VectorPostprocessorValue &	getVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const

const ScatterVectorPostprocessorValue &	getScatterVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const

bool	hasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const

bool	hasVectorPostprocessor (const std::string &param_name) const

bool	hasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const

bool	hasVectorPostprocessorByName (const VectorPostprocessorName &name) const

const VectorPostprocessorName &	getVectorPostprocessorName (const std::string &param_name) const

T &	getSampler (const std::string &name)

Sampler &	getSampler (const std::string &name)

T &	getSamplerByName (const SamplerName &name)

Sampler &	getSamplerByName (const SamplerName &name)

virtual void	meshDisplaced ()

PerfGraph &	perfGraph ()

const PostprocessorValue &	getPostprocessorValue (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValue (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const

const PostprocessorValue &	getPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const

virtual const PostprocessorValue &	getPostprocessorValueByName (const PostprocessorName &name) const

virtual const PostprocessorValue &	getPostprocessorValueByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOldByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOldByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOlderByName (const PostprocessorName &name) const

const PostprocessorValue &	getPostprocessorValueOlderByName (const PostprocessorName &name) const

bool	isVectorPostprocessorDistributed (const std::string &param_name) const

bool	isVectorPostprocessorDistributed (const std::string &param_name) const

bool	isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const

bool	isVectorPostprocessorDistributedByName (const VectorPostprocessorName &name) const

const Distribution &	getDistribution (const std::string &name) const

const T &	getDistribution (const std::string &name) const

const Distribution &	getDistribution (const std::string &name) const

const T &	getDistribution (const std::string &name) const

const Distribution &	getDistributionByName (const DistributionName &name) const

const T &	getDistributionByName (const std::string &name) const

const Distribution &	getDistributionByName (const DistributionName &name) const

const T &	getDistributionByName (const std::string &name) const

const Parallel::Communicator &	comm () const

processor_id_type	n_processors () const

processor_id_type	processor_id () const

Static Public Member Functions
static InputParameters	validParams ()

static void	sort (typename std::vector< T > &vector)

static void	sortDFS (typename std::vector< T > &vector)

static void	cyclicDependencyError (CyclicDependencyException< T2 > &e, const std::string &header)

Public Attributes
const ConsoleStream	_console

Static Public Attributes
static constexpr PropertyValue::id_type	default_property_id

static constexpr PropertyValue::id_type	zero_property_id

static constexpr auto	SYSTEM

static constexpr auto	NAME

Protected Member Functions
virtual void	addPostprocessorDependencyHelper (const PostprocessorName &name) const override

virtual void	addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &name) const override

virtual void	addUserObjectDependencyHelper (const UserObject &uo) const override

void	addReporterDependencyHelper (const ReporterName &reporter_name) override

const ReporterName &	getReporterName (const std::string &param_name) const

T &	declareRestartableData (const std::string &data_name, Args &&... args)

ManagedValue< T >	declareManagedRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)

const T &	getRestartableData (const std::string &data_name) const

T &	declareRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)

T &	declareRecoverableData (const std::string &data_name, Args &&... args)

T &	declareRestartableDataWithObjectName (const std::string &data_name, const std::string &object_name, Args &&... args)

T &	declareRestartableDataWithObjectNameWithContext (const std::string &data_name, const std::string &object_name, void *context, Args &&... args)

std::string	restartableName (const std::string &data_name) const

const T &	getMeshProperty (const std::string &data_name, const std::string &prefix)

const T &	getMeshProperty (const std::string &data_name)

bool	hasMeshProperty (const std::string &data_name, const std::string &prefix) const

bool	hasMeshProperty (const std::string &data_name, const std::string &prefix) const

bool	hasMeshProperty (const std::string &data_name) const

bool	hasMeshProperty (const std::string &data_name) const

std::string	meshPropertyName (const std::string &data_name) const

PerfID	registerTimedSection (const std::string &section_name, const unsigned int level) const

PerfID	registerTimedSection (const std::string &section_name, const unsigned int level, const std::string &live_message, const bool print_dots=true) const

std::string	timedSectionName (const std::string &section_name) const

bool	isCoupledScalar (const std::string &var_name, unsigned int i=0) const

unsigned int	coupledScalarComponents (const std::string &var_name) const

unsigned int	coupledScalar (const std::string &var_name, unsigned int comp=0) const

libMesh::Order	coupledScalarOrder (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarValue (const std::string &var_name, unsigned int comp=0) const

const ADVariableValue &	adCoupledScalarValue (const std::string &var_name, unsigned int comp=0) const

const GenericVariableValue< is_ad > &	coupledGenericScalarValue (const std::string &var_name, unsigned int comp=0) const

const GenericVariableValue< false > &	coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const

const GenericVariableValue< true > &	coupledGenericScalarValue (const std::string &var_name, const unsigned int comp) const

const VariableValue &	coupledVectorTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const

const VariableValue &	coupledMatrixTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const

const VariableValue &	coupledScalarValueOld (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarValueOlder (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDot (const std::string &var_name, unsigned int comp=0) const

const ADVariableValue &	adCoupledScalarDot (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotDot (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotOld (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotDotOld (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotDu (const std::string &var_name, unsigned int comp=0) const

const VariableValue &	coupledScalarDotDotDu (const std::string &var_name, unsigned int comp=0) const

const MooseVariableScalar *	getScalarVar (const std::string &var_name, unsigned int comp) const

virtual void	checkMaterialProperty (const std::string &name, const unsigned int state)

void	markMatPropRequested (const std::string &)

MaterialPropertyName	getMaterialPropertyName (const std::string &name) const

void	checkExecutionStage ()

VectorPostprocessorValue &	declareVector (const std::string &vector_name)

const T &	getReporterValue (const std::string &param_name, const std::size_t time_index=0)

const T &	getReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)

const T &	getReporterValue (const std::string &param_name, const std::size_t time_index=0)

const T &	getReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)

const T &	getReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)

const T &	getReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)

const T &	getReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)

const T &	getReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)

bool	hasReporterValue (const std::string &param_name) const

bool	hasReporterValue (const std::string &param_name) const

bool	hasReporterValue (const std::string &param_name) const

bool	hasReporterValue (const std::string &param_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

bool	hasReporterValueByName (const ReporterName &reporter_name) const

const GenericMaterialProperty< T, is_ad > *	defaultGenericMaterialProperty (const std::string &name)

const GenericMaterialProperty< T, is_ad > *	defaultGenericMaterialProperty (const std::string &name)

const MaterialProperty< T > *	defaultMaterialProperty (const std::string &name)

const MaterialProperty< T > *	defaultMaterialProperty (const std::string &name)

const ADMaterialProperty< T > *	defaultADMaterialProperty (const std::string &name)

const ADMaterialProperty< T > *	defaultADMaterialProperty (const std::string &name)

Static Protected Member Functions
static std::string	meshPropertyName (const std::string &data_name, const std::string &prefix)

Protected Attributes
const std::shared_ptr< MooseMesh > &	_mesh
	Reference to the mesh. More...

std::vector< VariableName >	_variables
	Variables to output. More...

std::vector< unsigned int >	_var_numbers
	Indices of the variables in their systems. More...

std::vector< VectorPostprocessorValue * >	_output_vector
	Vector of outputs, where each entry is the vector of average values for single variable at the selected points along the axis. More...

std::vector< SubdomainID >	_block_ids

const Point	_direction
	Axis direction of the structural component. More...

const Point	_reference_point
	Starting or reference point of the structural component to locate nodes on the cross section. More...

std::vector< Real >	_positions
	Axial positions along the component at which average values are computed. More...

const bool	_have_symmetry_plane
	Whether a symmetry plane has been defined by the user. More...

RealVectorValue	_symmetry_plane
	Vector normal to a symmetry plane, optionally defined if the section has a symmetry plane. More...

const bool	_automatically_locate_positions
	Whether to automatically locate positions along section for averaging field values. More...

const Real	_tolerance
	Tolerance to identify nodes on the user-prescribed cross section. More...

std::vector< unsigned int >	_number_of_nodes
	Number of nodes for computing output (local and global). More...

const Real	_cross_section_maximum_radius
	Tolerance to disambiguate cross section locations in different components within the same block. More...

const bool	_require_equal_node_counts
	Whether to require the number of nodes at each axial location to be equal. More...

bool	_need_mesh_initializations
	Whether node locations need to be identified and nodes at positions need to be counted. More...

SubProblem &	_subproblem

FEProblemBase &	_fe_problem

SystemBase &	_sys

const THREAD_ID	_tid

Assembly &	_assembly

const Moose::CoordinateSystemType &	_coord_sys

const bool	_duplicate_initial_execution

std::set< std::string >	_depend_uo

const bool &	_enabled

MooseApp &	_app

const std::string	_type

const std::string	_name

const InputParameters &	_pars

Factory &	_factory

ActionFactory &	_action_factory

const ExecFlagEnum &	_execute_enum

const ExecFlagType &	_current_execute_flag

MooseApp &	_restartable_app

const std::string	_restartable_system_name

const THREAD_ID	_restartable_tid

const bool	_restartable_read_only

FEProblemBase &	_mci_feproblem

FEProblemBase &	_mdi_feproblem

MooseApp &	_pg_moose_app

const std::string	_prefix

FEProblemBase &	_sc_fe_problem

const THREAD_ID	_sc_tid

const Real &	_real_zero

const VariableValue &	_scalar_zero

const Point &	_point_zero

const InputParameters &	_mi_params

const std::string	_mi_name

const MooseObjectName	_mi_moose_object_name

FEProblemBase &	_mi_feproblem

SubProblem &	_mi_subproblem

const THREAD_ID	_mi_tid

const Moose::MaterialDataType	_material_data_type

MaterialData &	_material_data

bool	_stateful_allowed

bool	_get_material_property_called

std::vector< std::unique_ptr< PropertyValue > >	_default_properties

std::unordered_set< unsigned int >	_material_property_dependencies

const MaterialPropertyName	_get_suffix

const bool	_use_interpolated_state

const InputParameters &	_ti_params

FEProblemBase &	_ti_feproblem

bool	_is_implicit

Real &	_t

const Real &	_t_old

int &	_t_step

Real &	_dt

Real &	_dt_old

bool	_is_transient

const std::string	_vpp_name

FEProblemBase &	_vpp_fe_problem

const MooseEnum &	_parallel_type

const Parallel::Communicator &	_communicator

Static Protected Attributes
static const std::string	_interpolated_old

static const std::string	_interpolated_older

Private Member Functions
Real	axialPosition (const Node &node) const
	Determine axial distance of the point from the component's reference point. More...

void	automaticallyLocatePositions ()
	Automatically identify all axial positions of nodes within the component and store their unique values (within tolerance) in _lengths. More...

Detailed Description

Definition at line 15 of file AverageSectionValueSampler.h.

Constructor & Destructor Documentation

◆ AverageSectionValueSampler()

AverageSectionValueSampler::AverageSectionValueSampler ( const InputParameters & parameters )

Definition at line 67 of file AverageSectionValueSampler.C.

   : GeneralVectorPostprocessor(parameters),
     _mesh(_app.actionWarehouse().mesh()),
     _variables(getParam<std::vector<VariableName>>("variables")),
     _direction(getParam<Point>("axis_direction")),
     _reference_point(getParam<Point>("reference_point")),
     _positions(isParamSetByUser("lengths") ? getParam<std::vector<Real>>("lengths")
                                            : getParam<std::vector<Real>>("positions")),
     _have_symmetry_plane(isParamSetByUser("symmetry_plane")),
     _symmetry_plane(getParam<RealVectorValue>("symmetry_plane")),
     _automatically_locate_positions(_positions.size() == 0),
     _tolerance(getParam<Real>("tolerance")),
     _number_of_nodes(_positions.size()),
     _cross_section_maximum_radius(getParam<Real>("cross_section_maximum_radius")),
     _require_equal_node_counts(getParam<bool>("require_equal_node_counts")),
     _need_mesh_initializations(true)
 {
   if (parameters.isParamSetByUser("lengths") && parameters.isParamSetByUser("positions"))
     paramError("lengths", "The 'lengths' and 'positions' parameters cannot both be set.");
 
   if (parameters.isParamSetByUser("lengths") && _positions.size() == 0)
     paramError("lengths", "If 'lengths' is specified, at least one value must be provided");
 
   if (parameters.isParamSetByUser("positions") && _positions.size() == 0)
     paramError("positions", "If 'positions' is specified, at least one value must be provided");
 
   if (_mesh->dimension() != 3)
     mooseError("The AverageSectionValueSampler postprocessor can only be used with three "
                "dimensional meshes.");
 
   if (!MooseUtils::absoluteFuzzyEqual(_direction.norm_sq(), 1.0))
     paramError("axis_direction", "Axis direction must be a unit vector.");
 
   _output_vector.resize(_variables.size() + 2);
   for (const auto j : make_range(_variables.size()))
     _output_vector[j] = &declareVector(_variables[j]);
   const auto pos_idx = _variables.size();
   _output_vector[pos_idx] = &declareVector("axial_position");
   _output_vector[pos_idx + 1] = &declareVector("node_count");
 
   for (const auto j : make_range(_variables.size()))
   {
     if (!_sys.hasVariable(_variables[j]))
       paramError("variables",
                  "One or more of the specified variables do not exist in this system.");
     const MooseVariable & variable = _sys.getFieldVariable<Real>(_tid, _variables[j]);
     if (!variable.isNodal())
       paramError("variables",
                  "One or more of the specified variables are not defined at the nodes.");
     _var_numbers.push_back(variable.number());
   }
 
   if (_have_symmetry_plane)
   {
     const auto symm_plane_norm = _symmetry_plane.norm();
     if (MooseUtils::absoluteFuzzyEqual(symm_plane_norm, 0.0))
       mooseError("Vector defined in 'symmetry_plane' cannot have a norm of zero");
     _symmetry_plane /= _symmetry_plane.norm();
     if (_variables.size() != 3)
       paramError("variables",
                  "If 'symmetry_plane' is prescribed, three variables must be provided (for the x, "
                  "y, z vector components, in that order)");
     mooseInfo("In AverageSectionValueSampler " + name() +
               ": \nTreating the variables as vector components (" + _variables[0] + ", " +
               _variables[1] + ", " + _variables[2] +
               ") and zeroing out the component normal to the prescribed 'symmetry_plane'\n");
   }
 }

Member Function Documentation

◆ automaticallyLocatePositions()

void AverageSectionValueSampler::automaticallyLocatePositions ( )

private

Automatically identify all axial positions of nodes within the component and store their unique values (within tolerance) in _lengths.

Definition at line 285 of file AverageSectionValueSampler.C.

Referenced by initialize().

 {
   _positions.clear();
 
   // Data structure used to collect the locations with nodes on each processor,
   // and gather in parallel.
   std::vector<Real> pos_vec;
 
   // Data structure used to store parallel-gathered locations
   std::set<Real> pos_set;
 
   _block_ids = _mesh->getSubdomainIDs(getParam<std::vector<SubdomainName>>("block"));
   auto * nodes = _mesh->getLocalNodeRange();
 
   for (const auto node : *nodes)
   {
     const std::set<SubdomainID> & node_blk_ids = _mesh->getNodeBlockIds(*node);
 
     for (const auto i : _block_ids)
     {
       if (node_blk_ids.find(i) != node_blk_ids.end())
       {
         bool found_match = false;
         Real axial_position = axialPosition(*node);
         if (axial_position == std::numeric_limits<Real>::max()) // Node not in this component
           continue;
         for (auto & pos : pos_vec)
         {
           const Real dist_to_plane = std::abs(axial_position - pos);
           if (dist_to_plane <= _tolerance)
           {
             found_match = true;
             break;
           }
         }
         if (!found_match)
           pos_vec.emplace_back(axial_position);
       }
     }
   }
 
   _communicator.allgather(pos_vec);
 
   for (const auto & posv : pos_vec)
   {
     bool found_match = false;
     for (auto & poss : pos_set)
     {
       if (std::abs(posv - poss) < _tolerance)
       {
         found_match = true;
         break;
       }
     }
     if (!found_match)
       pos_set.insert(posv);
   }
   for (const auto & poss : pos_set)
     _positions.emplace_back(poss);
 }

◆ axialPosition()

Real AverageSectionValueSampler::axialPosition ( const Node & node ) const

private

Determine axial distance of the point from the component's reference point.

If the in-plane distance is greater than the input parameter cross section maximum radius, return the maximum value for a Real.

Parameters

node	The node whose distance from a plane is to be considered

Definition at line 266 of file AverageSectionValueSampler.C.

Referenced by automaticallyLocatePositions(), and execute().

 {
   // Compute node location w.r.t. structural component length
   const Point relative_distance{
       node(0) - _reference_point(0), node(1) - _reference_point(1), node(2) - _reference_point(2)};
 
   const Real axial_position = _direction * relative_distance;
   const Real in_plane_distance =
       (relative_distance - relative_distance * _direction * _direction).norm();
 
   // If the in-plane distance is greater than the specified cross-section radius, the point is not
   // in this component
   if (in_plane_distance > _cross_section_maximum_radius)
     return std::numeric_limits<Real>::max();
 
   return axial_position;
 }

◆ execute()

void AverageSectionValueSampler::execute ( )

overridevirtual

Implements GeneralVectorPostprocessor.

Definition at line 209 of file AverageSectionValueSampler.C.

 {
   _block_ids = _mesh->getSubdomainIDs(getParam<std::vector<SubdomainName>>("block"));
 
   auto * active_nodes = _mesh->getActiveSemiLocalNodeRange();
 
   std::vector<std::vector<Real>> output_vector_partial(_variables.size(),
                                                        std::vector<Real>(_positions.size()));
 
   const NumericVector<Number> & _sol = *_sys.currentSolution();
 
   for (const auto & node : *active_nodes)
   {
     const std::set<SubdomainID> & node_blk_ids = _mesh->getNodeBlockIds(*node);
 
     for (const auto li : index_range(_positions))
     {
       for (const auto i : _block_ids)
       {
         if (node_blk_ids.find(i) != node_blk_ids.end())
         {
           // Check if node is within tolerance of user-prescribed plane
           if (std::abs(axialPosition(*node) - _positions[li]) > _tolerance)
             continue;
           if ((*node).processor_id() != processor_id())
             continue;
 
           // Retrieve nodal variables
           std::vector<Real> this_node_vars(_var_numbers.size());
           for (const auto j : make_range(_var_numbers.size()))
             this_node_vars[j] = _sol(node->dof_number(_sys.number(), _var_numbers[j], 0));
 
           if (_have_symmetry_plane)
           {
             RealVectorValue this_node_vec_var(
                 this_node_vars[0], this_node_vars[1], this_node_vars[2]);
             this_node_vec_var -= _symmetry_plane * this_node_vec_var * _symmetry_plane;
             this_node_vars[0] = this_node_vec_var(0);
             this_node_vars[1] = this_node_vec_var(1);
             this_node_vars[2] = this_node_vec_var(2);
           }
 
           for (const auto j : make_range(_var_numbers.size()))
             output_vector_partial[j][li] += this_node_vars[j];
 
           _number_of_nodes[li]++;
           break;
         }
       }
     }
   }
 
   for (const auto j : make_range(_variables.size()))
     *_output_vector[j] = output_vector_partial[j];
 }

◆ finalize()

void AverageSectionValueSampler::finalize ( )

overridevirtual

Reimplemented from GeneralVectorPostprocessor.

Definition at line 168 of file AverageSectionValueSampler.C.

 {
   for (const auto j : make_range(_variables.size()))
     _communicator.sum(*(_output_vector[j]));
 
   _communicator.sum(_number_of_nodes);
 
   for (const auto li : index_range(_positions))
     if (_number_of_nodes[li] < 1)
       mooseError("No nodes were found in AverageSectionValueSampler postprocessor.");
 
   for (const auto li : index_range(_positions))
     for (const auto j : make_range(_variables.size()))
       (*_output_vector[j])[li] /= _number_of_nodes[li];
 
   const auto pos_idx = _variables.size();
   (*_output_vector[pos_idx]) = _positions;
   std::vector<Real> num_nodes_real(_number_of_nodes.begin(), _number_of_nodes.end());
   (*_output_vector[pos_idx + 1]) = num_nodes_real;
 
   if (_require_equal_node_counts)
   {
     for (const auto li : index_range(_number_of_nodes))
     {
       if (_number_of_nodes[li] != _number_of_nodes[0])
       {
         std::ostringstream pos_out;
         for (const auto li2 : index_range(_number_of_nodes))
           pos_out << std::setw(10) << _positions[li2] << std::setw(10) << _number_of_nodes[li2]
                   << "\n";
         mooseError(
             "Node counts do not match for all axial positions. If this behavior "
             "is desired to accommodate nonuniform meshes, set "
             "'require_equal_node_counts=false'\n     Axial      Node\n  Position     Count\n" +
             pos_out.str());
       }
     }
   }
 }

◆ initialize()

void AverageSectionValueSampler::initialize ( )

overridevirtual

Implements GeneralVectorPostprocessor.

Definition at line 148 of file AverageSectionValueSampler.C.

 {
   if (_need_mesh_initializations)
   {
     // Option 1: locate positions first, then check counts
     if (_automatically_locate_positions)
       automaticallyLocatePositions();
     _need_mesh_initializations = false;
   }
 
   for (const auto j : make_range(_variables.size() + 2))
   {
     _output_vector[j]->clear();
     _output_vector[j]->resize(_positions.size());
   }
   _number_of_nodes.clear();
   _number_of_nodes.resize(_positions.size());
 }

◆ initialSetup()

void AverageSectionValueSampler::initialSetup ( )

overridevirtual

Reimplemented from GeneralVectorPostprocessor.

Definition at line 137 of file AverageSectionValueSampler.C.

138 {

139 }

◆ meshChanged()

void AverageSectionValueSampler::meshChanged ( )

overridevirtual

Reimplemented from GeneralVectorPostprocessor.

Definition at line 142 of file AverageSectionValueSampler.C.

 {
   _need_mesh_initializations = true;
 }

◆ validParams()

InputParameters AverageSectionValueSampler::validParams ( )

static

Definition at line 19 of file AverageSectionValueSampler.C.

 {
   InputParameters params = GeneralVectorPostprocessor::validParams();
   params.addClassDescription("Compute the section's variable average in three-dimensions "
                              "given a user-defined definition of the cross section.");
   params.addParam<std::vector<SubdomainName>>(
       "block",
       "The list of blocks in which to search for cross sectional nodes to compute the variable "
       "average.");
   params.addRequiredParam<Point>("axis_direction", "Direction of the structural component's axis");
   params.addRequiredParam<Point>("reference_point",
                                  "Structural component reference starting point from which the "
                                  "input parameter 'lengths' applies.");
   params.addParam<Real>("cross_section_maximum_radius",
                         std::numeric_limits<Real>::max(),
                         "Radial distance with respect to the body axis within which nodes are "
                         "considered to belong to this "
                         "structural component. Used to disambiguate multiple components that share "
                         "the same mesh block.");
   params.addRequiredParam<std::vector<VariableName>>(
       "variables", "Variables for the cross section output. These variables must be nodal.");
   params.addDeprecatedParam<std::vector<Real>>(
       "lengths",
       {},
       "Positions along axis from reference_point at which to compute average values.",
       "Use the 'positions' parameter instead");
   params.addParam<std::vector<Real>>(
       "positions",
       {},
       "Positions along axis from reference_point at which to compute average values.");
   params.addParam<RealVectorValue>(
       "symmetry_plane",
       RealVectorValue(0, 0, 0),
       "Vector normal to a symmetry plane, used if the section has a symmetry plane through it. "
       "Causes the variables to be treated as components of a vector and zeros out the compomnent "
       "in the direction normal to symmetry_plane. Three variables must be defined in 'variables' "
       "to use this option.");
   params.addParam<Real>("tolerance",
                         1.0e-6,
                         "Maximum axial distance of nodes from the specified axial lengths to "
                         "consider them in the cross-section average");
   params.addParam<bool>(
       "require_equal_node_counts",
       true,
       "Whether to require the number of nodes at each axial location to be equal");
   return params;
 }