#include <PODReducedBasisTrainer.h>

Inheritance diagram for PODReducedBasisTrainer:

Public Types
typedef DataFileName	DataFileParameterType

Public Member Functions
	PODReducedBasisTrainer (const InputParameters &parameters)

virtual void	initialSetup () override

virtual void	initialize () override

virtual void	execute () override

virtual void	finalize () override

void	initReducedOperators ()
	Initializing the reduced operators. More...

void	addSnapshot (unsigned int var_i, unsigned int glob_i, const std::shared_ptr< DenseVector< Real >> &snapshot)
	Adding a snapshot for a variable. More...

void	addToReducedOperator (unsigned int base_i, unsigned int tag_i, std::vector< DenseVector< Real >> &residual)
	Adding the contribution of a residual to the reduced operators. More...

const std::vector< std::string > &	getVarNames () const

const std::vector< std::string > &	getTagNames () const

const std::vector< std::string > &	getTagTypes () const

unsigned int	getSnapsSize (unsigned int var_i) const
	Getting the snapshot size across all of the processors for a given variable with index var_i. More...

unsigned int	getBaseSize (unsigned int var_i) const
	Getting the base size for variable with index v_ind. More...

unsigned int	getSumBaseSize () const
	Getting the overall base size, which is the sum of the individual bases. More...

const DenseVector< Real > &	getBasisVector (unsigned int var_i, unsigned int base_i) const
	Getting a basis vector for a given variable. More...

const DenseVector< Real > &	getBasisVector (unsigned int glob_i) const
	Getting basis vector based on its global index. More...

unsigned int	getVariableIndex (unsigned int glob_i) const
	Getting appropriate variable index for a global base index. More...

virtual void	threadJoin (const UserObject &) final

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

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	meshChanged ()

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

const std::string &	modelMetaDataName () const
	Accessor for the name of the model meta data. More...

const FileName &	getModelDataFileName () const
	Get the associated filename. More...

bool	hasModelData () const
	Check if we need to load model data (if the filename parameter is used) More...


template<typename T , typename... Args>
T &	declareModelData (const std::string &data_name, Args &&... args)
	Declare model data for loading from file as well as restart. More...


template<typename T , typename... Args>
const T &	getModelData (const std::string &data_name, Args &&... args) const
	Retrieve model data from the interface. More...

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
void	computeCorrelationMatrix ()
	Computes the correlation matrices using the snapshots. More...

void	computeEigenDecomposition ()
	Computes the eigen-decomposition of the stored correlation matrices. More...

void	computeBasisVectors ()
	Generates the basis vectors using the snapshots together with the eigen-decomposition of the correlation matrices. More...

void	printEigenvalues ()
	Prints the eigenvalues of the correlation matrix for each variable. More...

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 ()

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
std::vector< std::string > &	_var_names
	Vector containing the names of the variables we want to use for constructing the surrogates. More...

std::vector< Real >	_error_res
	Energy limits that define how many basis functions will be kept for each variable. More...

std::vector< std::string > &	_tag_names
	Names of the tags that should be used to fetch residuals from the MultiApp. More...

std::vector< std::string > &	_tag_types
	list of bools describing which tag is indepedent of the solution. More...

std::vector< DistributedSnapshots >	_snapshots
	Distributed container for snapshots per variable. More...

std::vector< DenseMatrix< Real > >	_corr_mx
	The correlation matrices for the variables. More...

std::vector< DenseVector< Real > >	_eigenvalues
	The eigenvalues of the correalation matrix for each variable. More...

std::vector< DenseMatrix< Real > >	_eigenvectors
	The eigenvectors of the correalation matrix for each variable. More...

std::vector< std::vector< DenseVector< Real > > > &	_base
	The reduced basis for the variables. More...

std::vector< DenseMatrix< Real > > &	_red_operators
	The reduced operators that should be transferred to the surrogate. More...

bool	_base_completed
	Switch that tells if the object has already computed the necessary basis vectors. More...

bool	_empty_operators
	Flag describing if the reduced operators are empty or not. 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 Parallel::Communicator &	_communicator

Static Protected Attributes
static const std::string	_interpolated_old

static const std::string	_interpolated_older

Private Member Functions
unsigned int	determineNumberOfModes (Real error, const std::vector< Real > &inp_vec) const
	Computes the number of bases necessary for a given error indicator. More...

void	receiveObjects (ReplicatedMesh &mesh, std::unordered_map< unsigned int, std::vector< std::shared_ptr< DenseVector< Real >>>> &received_vectors, std::unordered_map< unsigned int, std::vector< std::shared_ptr< DenseVector< Real >>>> &local_vectors, processor_id_type, const std::vector< std::tuple< unsigned int, unsigned int, std::shared_ptr< DenseVector< Real >>>> &vectors)
	Function that manipulates the received objects and computes the correlation matrices on the fly. More...

Detailed Description

Definition at line 27 of file PODReducedBasisTrainer.h.

Constructor & Destructor Documentation

◆ PODReducedBasisTrainer()

PODReducedBasisTrainer::PODReducedBasisTrainer ( const InputParameters & parameters )

Definition at line 44 of file PODReducedBasisTrainer.C.

   : SurrogateTrainerBase(parameters),
     _var_names(declareModelData<std::vector<std::string>>(
         "_var_names", getParam<std::vector<std::string>>("var_names"))),
     _error_res(getParam<std::vector<Real>>("error_res")),
     _tag_names(declareModelData<std::vector<std::string>>(
         "_tag_names", getParam<std::vector<std::string>>("tag_names"))),
     _tag_types(declareModelData<std::vector<std::string>>(
         "_tag_types", getParam<std::vector<std::string>>("tag_types"))),
     _base(declareModelData<std::vector<std::vector<DenseVector<Real>>>>("_base")),
     _red_operators(declareModelData<std::vector<DenseMatrix<Real>>>("_red_operators")),
     _base_completed(false),
     _empty_operators(true)
 {
   if (_error_res.size() != _var_names.size())
     paramError("error_res",
                "The number of elements is not equal to the number"
                " of elements in 'var_names'!");
 
   if (_tag_names.size() != _tag_types.size())
     paramError("tag_types",
                "The number of elements is not equal to the number"
                " of elements in 'tag_names'!");
 
   std::vector<std::string> available_names{"op", "src", "op_dir", "src_dir"};
   for (auto tag_type : _tag_types)
   {
     auto it = std::find(available_names.begin(), available_names.end(), tag_type);
     if (it == available_names.end())
       paramError("tag_types",
                  "Tag type '",
                  tag_type,
                  "' is not valid, available names are:",
                  " op, op_dir, src, src_dir");
   }
 }

Member Function Documentation

◆ addSnapshot()

void PODReducedBasisTrainer::addSnapshot	(	unsigned int	var_i,
		unsigned int	glob_i,
		const std::shared_ptr< DenseVector< Real >> &	snapshot
	)

Adding a snapshot for a variable.

Definition at line 134 of file PODReducedBasisTrainer.C.

Referenced by PODSamplerSolutionTransfer::execute(), and PODSamplerSolutionTransfer::executeFromMultiapp().

 {
   _snapshots[var_i].addNewEntry(glob_i, snapshot);
 }

◆ addToReducedOperator()

void PODReducedBasisTrainer::addToReducedOperator	(	unsigned int	base_i,
		unsigned int	tag_i,
		std::vector< DenseVector< Real >> &	residual
	)

Adding the contribution of a residual to the reduced operators.

Definition at line 500 of file PODReducedBasisTrainer.C.

Referenced by PODResidualTransfer::transferResidual().

 {
   // Computing the elements of the reduced operator using Galerkin projection
   // on the residual. This is done in parallel and the local contributions are
   // gathered in finalize().
   unsigned int counter = 0;
   for (unsigned int var_i = 0; var_i < _var_names.size(); ++var_i)
   {
     for (unsigned int base_j = 0; base_j < _base[var_i].size(); ++base_j)
     {
       _red_operators[tag_i](counter, base_i) = residual[var_i].dot(_base[var_i][base_j]);
       counter++;
     }
   }
 
   _empty_operators = false;
 }

◆ computeBasisVectors()

void PODReducedBasisTrainer::computeBasisVectors ( )

protected

Generates the basis vectors using the snapshots together with the eigen-decomposition of the correlation matrices.

Definition at line 443 of file PODReducedBasisTrainer.C.

Referenced by execute().

 {
   // Looping over all the variables.
   for (unsigned int var_i = 0; var_i < _eigenvectors.size(); ++var_i)
   {
     unsigned int no_bases = _eigenvalues[var_i].size();
     unsigned int no_snaps = _snapshots[var_i].getNumberOfLocalEntries();
 
     _base[var_i].resize(no_bases);
 
     // Filling the containers using the local snapshots and the eigenvalues and
     // eigenvectors of the correlation matrices.
     for (unsigned int base_i = 0; base_i < no_bases; ++base_i)
     {
       _base[var_i][base_i].resize(_snapshots[var_i].getLocalEntry(0)->size());
 
       for (unsigned int loc_i = 0; loc_i < no_snaps; ++loc_i)
       {
         unsigned int glob_i = _snapshots[var_i].getGlobalIndex(loc_i);
         const DenseVector<Real> & snapshot = *_snapshots[var_i].getLocalEntry(loc_i);
 
         for (unsigned int i = 0; i < _base[var_i][base_i].size(); ++i)
         {
           _base[var_i][base_i](i) += _eigenvectors[var_i](glob_i, base_i) * snapshot(i);
         }
       }
 
       // Gathering and summing the local contributions over all of the processes.
       // This makes sure that every process sees all of the basis functions.
       gatherSum(_base[var_i][base_i].get_values());
 
       // Normalizing the basis functions to make sure they are orthonormal.
       _base[var_i][base_i].scale(1.0 / sqrt(_eigenvalues[var_i](base_i)));
     }
   }
 }

◆ computeCorrelationMatrix()

void PODReducedBasisTrainer::computeCorrelationMatrix ( )

protected

Computes the correlation matrices using the snapshots.

Definition at line 142 of file PODReducedBasisTrainer.C.

Referenced by execute().

 {
   // Getting the number of snapshots. It is assumed that every variable has the
   // same number of snapshots. This assumption is used at multiple locations in
   // this source file.
   const auto no_snaps = getSnapsSize(0);
 
   // Initializing the correlation matrices for each variable.
   for (unsigned int var_i = 0; var_i < _snapshots.size(); ++var_i)
     _corr_mx[var_i] = DenseMatrix<Real>(no_snaps, no_snaps);
 
   /*
   Creating a simple 2D map that distributes the elements in the correlation
   matrices to each processor. A square mesh with elements corresponding to
   the entries in the correlation matrix is used, since optimal partitioning
   routines are already implemented for it. In this case, the centroids of the
   elements are equivalent to the indices of the correlation matrix. This method
   may seem too convoluted, but is necessary to ensure that the communication
   stays balanced and scales well for runs with a large number of processors.
   */
   ReplicatedMesh mesh(_communicator, 2);
   MeshTools::Generation::build_square(
       mesh, no_snaps, no_snaps, -0.5, no_snaps - 0.5, -0.5, no_snaps - 0.5);
 
   // A flag is added to each element saying if the result has already been compuuted
   // or not.
   mesh.add_elem_integer("filled");
 
   // Since the correlation matrix is symmetric, it is enough to compute the
   // entries on the diagonal in addition to the elements above the diagonal.
   // Therefore, the elements in the mesh below the diagonal are deleted.
   for (Elem * elem : mesh.active_element_ptr_range())
   {
     const auto centroid = elem->vertex_average();
     if (centroid(0) > centroid(1))
       mesh.delete_elem(elem);
   }
 
   // The mesh is distributed among the processors.
   mesh.prepare_for_use();
 
   /*
   This step restructures the snapshots into an unordered map to make it easier
   to do cross-products in the later stages. Since _snapshots only contains
   pointers to the data, this should not include a considerable amount of copy
   operations.
   */
   std::unordered_map<unsigned int, std::vector<std::shared_ptr<DenseVector<Real>>>> local_vectors;
   for (unsigned int loc_vec_i = 0; loc_vec_i < _snapshots[0].getNumberOfLocalEntries(); ++loc_vec_i)
   {
     const unsigned int glob_vec_i = _snapshots[0].getGlobalIndex(loc_vec_i);
     auto & entry = local_vectors[glob_vec_i];
 
     for (unsigned int v_ind = 0; v_ind < _snapshots.size(); ++v_ind)
       entry.push_back(_snapshots[v_ind].getLocalEntry(loc_vec_i));
   }
 
   /*
   Initializing containers for the objects that will be sent to other processors.
   send_vectors is just a temporary object that ensures that the same snapshot is
   sent to other processors only once. send_map, on the other hand, will be used
   by the communicator and contains the required snapshots for each processor
   which is not the the current rank.
   std::tuple<unsigned int, unsigned int, std::shared_ptr<DenseVector<Real>>> type
   is a container that uses (global snapshot index, variable index, snapshot) to
   identify and send/receive snapshots during the communication.
   */
 
   std::unordered_map<processor_id_type, std::set<unsigned int>> send_vectors;
   std::unordered_map<
       processor_id_type,
       std::vector<std::tuple<unsigned int, unsigned int, std::shared_ptr<DenseVector<Real>>>>>
       send_map;
 
   // Fill the send map with snapshots we need to send for each processor. First,
   // we loop over the matrix entries (elements in the mesh)
   for (Elem * elem : mesh.active_element_ptr_range())
     if (elem->processor_id() != processor_id())
     {
       // The centroids in the mesh correspond to the 2D corrdinates of the elements
       // in the mesh.
       const auto centroid = elem->vertex_average();
       const unsigned int i = centroid(0);
       const unsigned int j = centroid(1);
 
       /*
       Checking if the current processor has the required snapshot.
       We assume that every variable has the same number of snapshots with the
       same distribution among processors. Therefore, it is enough to test
       the first variable only.
       */
       if (_snapshots[0].hasGlobalEntry(i))
       {
         // Continue loop if the snapshot is already being sent.
         if (send_vectors[elem->processor_id()].count(i))
           continue;
 
         // Add another entry to the map if another processor needs the owned
         // snapshot.
         send_vectors[elem->processor_id()].insert(i);
         for (unsigned int v_ind = 0; v_ind < _snapshots.size(); ++v_ind)
           send_map[elem->processor_id()].emplace_back(
               i, v_ind, _snapshots[v_ind].getGlobalEntry(i));
       }
       else if (_snapshots[0].hasGlobalEntry(j))
       {
         // Continue loop if the snapshot is already being sent.
         if (send_vectors[elem->processor_id()].count(j))
           continue;
 
         // Add another entry to the map if another processor needs the owned
         // snapshot.
         send_vectors[elem->processor_id()].insert(j);
         for (unsigned int v_ind = 0; v_ind < _snapshots.size(); ++v_ind)
           send_map[elem->processor_id()].emplace_back(
               j, v_ind, _snapshots[v_ind].getGlobalEntry(j));
       }
     }
     else
       // Initializing the flag value with 0 (not computed) for each element.
       elem->set_extra_integer(0, 0);
 
   // Creating container for received data. In this case the map contains the snapshots
   // for each variable for each required global snapshot index.
   std::unordered_map<unsigned int, std::vector<std::shared_ptr<DenseVector<Real>>>>
       received_vectors;
 
   // Converting function to functor to be able to pass to push packed range.
   auto functor =
       [this, &mesh, &received_vectors, &local_vectors](
           processor_id_type pid,
           const std::vector<
               std::tuple<unsigned int, unsigned int, std::shared_ptr<DenseVector<Real>>>> & vectors)
   { PODReducedBasisTrainer::receiveObjects(mesh, received_vectors, local_vectors, pid, vectors); };
 
   Parallel::push_parallel_packed_range(_communicator, send_map, (void *)nullptr, functor);
 
   // This extra call is necessary in case a processor has all the elements it needs
   // (hence doesn't receive any).
   functor(
       processor_id(),
       std::vector<std::tuple<unsigned int, unsigned int, std::shared_ptr<DenseVector<Real>>>>());
 
   // Now, the correlation matrices aregathered and  summed to make sure every processor
   // sees them.
   gatherSum(_corr_mx[0].get_values());
 
   // The lower triangle of the matrices are then filled using symmetry.
   for (auto & corr_mx : _corr_mx)
     for (unsigned int row_i = 1; row_i < corr_mx.m(); ++row_i)
       for (unsigned int col_i = 0; col_i < row_i; ++col_i)
         corr_mx(row_i, col_i) = corr_mx(col_i, row_i);
 }

◆ computeEigenDecomposition()

void PODReducedBasisTrainer::computeEigenDecomposition ( )

protected

Computes the eigen-decomposition of the stored correlation matrices.

Definition at line 377 of file PODReducedBasisTrainer.C.

Referenced by execute().

 {
   for (unsigned int v_ind = 0; v_ind < _corr_mx.size(); ++v_ind)
   {
     unsigned int no_snaps = _corr_mx[v_ind].n();
 
     // Initializing temporary objects for the eigenvalues and eigenvectors since
     // evd_left() returns an unordered vector of eigenvalues.
     DenseVector<Real> eigenvalues(no_snaps);
     DenseMatrix<Real> eigenvectors(no_snaps, no_snaps);
 
     // Creating a temporary placeholder for the imaginary parts of the eigenvalues
     DenseVector<Real> eigenvalues_imag(no_snaps);
 
     // Performing the eigenvalue decomposition
     _corr_mx[v_ind].evd_left(eigenvalues, eigenvalues_imag, eigenvectors);
 
     // Sorting the eigenvectors and eigenvalues based on the magnitude of
     // the eigenvalues
     std::vector<unsigned int> idx(eigenvalues.size());
     std::iota(idx.begin(), idx.end(), 0);
     std::vector<Real> & v = eigenvalues.get_values();
 
     // Getting the indices to be able to copy the corresponding eigenvector too.
     std::stable_sort(
         idx.begin(), idx.end(), [&v](unsigned int i, unsigned int j) { return v[i] > v[j]; });
 
     // Getting a cutoff for the number of modes. The function requires a sorted list,
     // thus the temporary vector is sorted.
     std::stable_sort(v.begin(), v.end(), std::greater<Real>());
     unsigned int cutoff = determineNumberOfModes(_error_res[v_ind], v);
 
     // Initializing the actual containers for the eigenvectors and eigenvalues.
     _eigenvalues[v_ind] = DenseVector<Real>(cutoff);
     _eigenvectors[v_ind] = DenseMatrix<Real>(eigenvectors.m(), cutoff);
 
     // Copying the kept eigenvalues and eigenvectors in a sorted container.
     for (unsigned int j = 0; j < cutoff; ++j)
     {
       _eigenvalues[v_ind](j) = eigenvalues(j);
       for (unsigned int k = 0; k < _eigenvectors[v_ind].m(); ++k)
       {
         _eigenvectors[v_ind](k, j) = eigenvectors(k, idx[j]);
       }
     }
 
     printEigenvalues();
   }
 }

◆ declareModelData()

template<typename T , typename... Args>

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

inherited

Declare model data for loading from file as well as restart.

Definition at line 78 of file RestartableModelInterface.h.

 {
   return _model_restartable.declareRestartableData<T>(data_name, std::forward<Args>(args)...);
 }

◆ determineNumberOfModes()

unsigned int PODReducedBasisTrainer::determineNumberOfModes	(	Real	error,
		const std::vector< Real > &	inp_vec
	)		const

private

Computes the number of bases necessary for a given error indicator.

This needs a sorted vector as input.

Definition at line 428 of file PODReducedBasisTrainer.C.

Referenced by computeEigenDecomposition().

 {
   Real sum = std::accumulate(inp_vec.begin(), inp_vec.end(), 0.0);
 
   Real part_sum = 0.0;
   for (unsigned int i = 0; i < inp_vec.size(); ++i)
   {
     part_sum += inp_vec[i];
     if (part_sum / sum > 1 - error)
       return (i + 1);
   }
   return (inp_vec.size());
 }

◆ execute()

void PODReducedBasisTrainer::execute ( )

overridevirtual

Implements GeneralUserObject.

Definition at line 104 of file PODReducedBasisTrainer.C.

 {
   // If the base is not ready yet, create it by performing the POD on the
   // snapshot matrices. Also, initialize the containers for the reduced
   // operators.
   if (!_base_completed)
   {
     computeCorrelationMatrix();
     computeEigenDecomposition();
     computeBasisVectors();
     initReducedOperators();
     _base_completed = true;
     _empty_operators = true;
   }
 }

◆ finalize()

void PODReducedBasisTrainer::finalize ( )

overridevirtual

Reimplemented from SurrogateTrainerBase.

Definition at line 121 of file PODReducedBasisTrainer.C.

 {
   // If the operators are already filled on each processor, gather and sum them
   // together. This way every processor has access to every complete reduced
   // operator.
   if (!_empty_operators)
     for (unsigned int tag_i = 0; tag_i < _red_operators.size(); ++tag_i)
     {
       gatherSum(_red_operators[tag_i].get_values());
     }
 }

◆ getBaseSize()

unsigned int PODReducedBasisTrainer::getBaseSize ( unsigned int var_i ) const

inline

Getting the base size for variable with index v_ind.

Definition at line 66 of file PODReducedBasisTrainer.h.

Referenced by PODSamplerSolutionTransfer::execute().

66 { return _base[var_i].size(); }

PODReducedBasisTrainer::_base

std::vector< std::vector< DenseVector< Real > > > & _base

The reduced basis for the variables.

Definition: PODReducedBasisTrainer.h:120

◆ getBasisVector() [1/2]

const DenseVector< Real > & PODReducedBasisTrainer::getBasisVector	(	unsigned int	var_i,
		unsigned int	base_i
	)		const

Getting a basis vector for a given variable.

Definition at line 538 of file PODReducedBasisTrainer.C.

Referenced by PODSamplerSolutionTransfer::execute(), and PODSamplerSolutionTransfer::executeToMultiapp().

 {
   return _base[var_i][base_i];
 }

◆ getBasisVector() [2/2]

const DenseVector< Real > & PODReducedBasisTrainer::getBasisVector ( unsigned int glob_i ) const

Getting basis vector based on its global index.

Definition at line 544 of file PODReducedBasisTrainer.C.

 {
   unsigned int counter = 0;
 
   for (unsigned int var_i = 0; var_i < _var_names.size(); ++var_i)
     for (unsigned int base_i = 0; base_i < _base[var_i].size(); ++base_i)
     {
       if (glob_i == counter)
         return _base[var_i][base_i];
 
       counter += 1;
     }
 
   mooseError("The basis vector with global index ", glob_i, "is not available!");
   return _base[0][0];
 }

◆ getModelData()

template<typename T , typename... Args>

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

inherited

Retrieve model data from the interface.

Definition at line 85 of file RestartableModelInterface.h.

 {
   return _model_restartable.getRestartableData<T>(data_name, std::forward<Args>(args)...);
 }

◆ getModelDataFileName()

const FileName & RestartableModelInterface::getModelDataFileName ( ) const

inherited

Get the associated filename.

Definition at line 33 of file RestartableModelInterface.C.

 {
   return _model_object.getParam<FileName>("filename");
 }

◆ getSnapsSize()

unsigned int PODReducedBasisTrainer::getSnapsSize ( unsigned int var_i ) const

Getting the snapshot size across all of the processors for a given variable with index var_i.

Definition at line 521 of file PODReducedBasisTrainer.C.

Referenced by computeCorrelationMatrix().

 {
   return _snapshots[var_i].getNumberOfGlobalEntries();
 }

◆ getSumBaseSize()

unsigned int PODReducedBasisTrainer::getSumBaseSize ( ) const

Getting the overall base size, which is the sum of the individual bases.

Definition at line 527 of file PODReducedBasisTrainer.C.

Referenced by PODFullSolveMultiApp::computeResidualBatch(), PODResidualTransfer::execute(), initReducedOperators(), and PODFullSolveMultiApp::preTransfer().

 {
   unsigned int sum = 0;
 
   for (unsigned int i = 0; i < _base.size(); ++i)
     sum += _base[i].size();
 
   return sum;
 }

◆ getTagNames()

const std::vector<std::string>& PODReducedBasisTrainer::getTagNames ( ) const

inline

Definition at line 57 of file PODReducedBasisTrainer.h.

Referenced by PODFullSolveMultiApp::computeResidual(), and PODResidualTransfer::transferResidual().

57 { return _tag_names; }

PODReducedBasisTrainer::_tag_names

std::vector< std::string > & _tag_names

Names of the tags that should be used to fetch residuals from the MultiApp.

Definition: PODReducedBasisTrainer.h:102

◆ getTagTypes()

const std::vector<std::string>& PODReducedBasisTrainer::getTagTypes ( ) const

inline

Definition at line 59 of file PODReducedBasisTrainer.h.

Referenced by PODResidualTransfer::transferResidual().

59 { return _tag_types; }

PODReducedBasisTrainer::_tag_types

std::vector< std::string > & _tag_types

list of bools describing which tag is indepedent of the solution.

Definition: PODReducedBasisTrainer.h:105

◆ getVariableIndex()

unsigned int PODReducedBasisTrainer::getVariableIndex ( unsigned int glob_i ) const

Getting appropriate variable index for a global base index.

Definition at line 562 of file PODReducedBasisTrainer.C.

Referenced by PODSamplerSolutionTransfer::executeToMultiapp().

 {
   unsigned int counter = 0;
 
   for (unsigned int var_i = 0; var_i < _var_names.size(); ++var_i)
     for (unsigned int base_i = 0; base_i < _base[var_i].size(); ++base_i)
     {
       if (glob_i == counter)
         return var_i;
 
       counter += 1;
     }
 
   mooseError("Variable with global base index ", glob_i, "is not available!");
   return 0;
 }

◆ getVarNames()

const std::vector<std::string>& PODReducedBasisTrainer::getVarNames ( ) const

inline

Definition at line 55 of file PODReducedBasisTrainer.h.

Referenced by PODSamplerSolutionTransfer::execute(), PODSamplerSolutionTransfer::executeFromMultiapp(), PODSamplerSolutionTransfer::executeToMultiapp(), PODSamplerSolutionTransfer::initialSetup(), and PODResidualTransfer::transferResidual().

55 { return _var_names; }

PODReducedBasisTrainer::_var_names

std::vector< std::string > & _var_names

Vector containing the names of the variables we want to use for constructing the surrogates.

Definition: PODReducedBasisTrainer.h:96

◆ hasModelData()

bool RestartableModelInterface::hasModelData ( ) const

inherited

Check if we need to load model data (if the filename parameter is used)

Definition at line 39 of file RestartableModelInterface.C.

 {
   return _model_object.isParamValid("filename");
 }

◆ initialize()

virtual void PODReducedBasisTrainer::initialize ( )

inlineoverridevirtual

Reimplemented from SurrogateTrainerBase.

Definition at line 36 of file PODReducedBasisTrainer.h.

36 {};

◆ initialSetup()

void PODReducedBasisTrainer::initialSetup ( )

overridevirtual

Reimplemented from GeneralUserObject.

Definition at line 82 of file PODReducedBasisTrainer.C.

 {
   // Initializing the containers for the essential data to construct the
   // reduced operators.
 
   _snapshots.clear();
   _snapshots.resize(_var_names.size(), DistributedSnapshots(_communicator));
 
   _base.clear();
   _base.resize(_var_names.size());
 
   _corr_mx.clear();
   _corr_mx.resize(_var_names.size());
 
   _eigenvalues.clear();
   _eigenvalues.resize(_var_names.size());
 
   _eigenvectors.clear();
   _eigenvectors.resize(_var_names.size());
 }

◆ initReducedOperators()

void PODReducedBasisTrainer::initReducedOperators ( )

Initializing the reduced operators.

Definition at line 481 of file PODReducedBasisTrainer.C.

Referenced by execute().

 {
   _red_operators.resize(_tag_names.size());
 
   // Getting the sum of the ranks of the different bases. Every reduced operator
   // is resized with this number. This way the construction can be more general.
   unsigned int base_num = getSumBaseSize();
 
   // Initializing each operator (each operator with a tag).
   for (unsigned int tag_i = 0; tag_i < _red_operators.size(); ++tag_i)
   {
     if (_tag_types[tag_i] == "src" || _tag_types[tag_i] == "src_dir")
       _red_operators[tag_i].resize(base_num, 1);
     else
       _red_operators[tag_i].resize(base_num, base_num);
   }
 }

◆ modelMetaDataName()

const std::string& RestartableModelInterface::modelMetaDataName ( ) const

inlineinherited

Accessor for the name of the model meta data.

Definition at line 47 of file RestartableModelInterface.h.

Referenced by SurrogateTrainerOutput::output(), and MappingOutput::output().

47 { return _model_meta_data_name; }

RestartableModelInterface::_model_meta_data_name

const std::string _model_meta_data_name

The model meta data name.

Definition: RestartableModelInterface.h:61

◆ printEigenvalues()

void PODReducedBasisTrainer::printEigenvalues ( )

protected

Prints the eigenvalues of the correlation matrix for each variable.

Definition at line 580 of file PODReducedBasisTrainer.C.

Referenced by computeEigenDecomposition().

 {
   if (processor_id() == 0 && _tid == 0 && isParamValid("filenames"))
   {
     std::vector<std::string> filenames = getParam<std::vector<std::string>>("filenames");
 
     if (filenames.size() != _var_names.size())
       paramError("filenames",
                  "The number of file names is not equal to the number of variable names!");
 
     for (unsigned int var_i = 0; var_i < _var_names.size(); ++var_i)
     {
       std::filebuf fb;
       fb.open(filenames[var_i], std::ios::out);
       std::ostream os(&fb);
       os << "evs" << std::endl;
       _eigenvalues[var_i].print_scientific(os);
     }
   }
 }

◆ receiveObjects()

void PODReducedBasisTrainer::receiveObjects	(	ReplicatedMesh &	mesh,
		std::unordered_map< unsigned int, std::vector< std::shared_ptr< DenseVector< Real >>>> &	received_vectors,
		std::unordered_map< unsigned int, std::vector< std::shared_ptr< DenseVector< Real >>>> &	local_vectors,
		processor_id_type	,
		const std::vector< std::tuple< unsigned int, unsigned int, std::shared_ptr< DenseVector< Real >>>> &	vectors
	)

private

Function that manipulates the received objects and computes the correlation matrices on the fly.

Definition at line 297 of file PODReducedBasisTrainer.C.

Referenced by computeCorrelationMatrix().

 {
   for (auto & tuple : vectors)
   {
     const auto glob_vec_i = std::get<0>(tuple);
     const auto var_i = std::get<1>(tuple);
     const auto & vector = std::get<2>(tuple);
 
     // The entry that will be filled with all of the variable solutions for a vector
     std::vector<std::shared_ptr<DenseVector<Real>>> & entry = received_vectors[glob_vec_i];
 
     // Size it to the number of variables in case we haven't already
     entry.resize(_snapshots.size(), nullptr);
 
     // Add this varaible's contribution - this is shared_ptr so we are claiming partial
     // ownership of this vector and do not have to do a copy
     entry[var_i] = vector;
   }
 
   // Looping over the locally owned entries in the matrix and computing the
   // values.
   for (Elem * elem : mesh.active_local_element_ptr_range())
   {
     // If the matrix entry is already filled, the loop skips this element.
     if (elem->get_extra_integer(0))
       continue;
 
     // Getting pointers to the necessary snapshots for the matrix entry.
     // This points to a vector of size (number of variables).
     const unsigned int i = elem->vertex_average()(0);
     const unsigned int j = elem->vertex_average()(1);
     std::vector<std::shared_ptr<DenseVector<Real>>> * i_vec = nullptr;
     std::vector<std::shared_ptr<DenseVector<Real>>> * j_vec = nullptr;
 
     const auto find_i_local = local_vectors.find(i);
     if (find_i_local != local_vectors.end())
       i_vec = &find_i_local->second;
     else
     {
       const auto find_i_received = received_vectors.find(i);
       if (find_i_received != received_vectors.end())
       {
         i_vec = &find_i_received->second;
       }
       else
         continue;
     }
 
     const auto find_j_local = local_vectors.find(j);
     if (find_j_local != local_vectors.end())
       j_vec = &find_j_local->second;
     else
     {
       const auto find_j_received = received_vectors.find(j);
       if (find_j_received != received_vectors.end())
       {
         j_vec = &find_j_received->second;
       }
       else
         continue;
     }
 
     // Coputing the available matrix entries for every variable.
     for (unsigned int v_ind = 0; v_ind < _snapshots.size(); ++v_ind)
       _corr_mx[v_ind](i, j) = (*i_vec)[v_ind]->dot(*(*j_vec)[v_ind]);
 
     // Set the 'filled' flag to 1 (true) to make sure it is not recomputed.
     elem->set_extra_integer(0, 1);
   }
 }

◆ threadJoin()

virtual void SurrogateTrainerBase::threadJoin ( const UserObject & )

inlinefinalvirtualinherited

Reimplemented from GeneralUserObject.

Definition at line 40 of file SurrogateTrainer.h.

40 {} // GeneralUserObjects are not threaded

◆ validParams()

InputParameters PODReducedBasisTrainer::validParams ( )

static

Definition at line 21 of file PODReducedBasisTrainer.C.

 {
   InputParameters params = SurrogateTrainerBase::validParams();
 
   params.addClassDescription("Computes the reduced subspace plus the reduced operators for "
                              "POD-RB surrogate.");
   params.addRequiredParam<std::vector<std::string>>("var_names",
                                                     "Names of variables we want to "
                                                     "extract from solution vectors.");
   params.addRequiredParam<std::vector<Real>>("error_res",
                                              "The errors allowed in the snapshot reconstruction.");
   params.addRequiredParam<std::vector<std::string>>("tag_names",
                                                     "Names of tags for the reduced operators.");
   params.addParam<std::vector<std::string>>(
       "filenames", "Files where the eigenvalues are printed for each variable (if given).");
 
   params.addRequiredParam<std::vector<std::string>>(
       "tag_types",
       "List of keywords describing if the tags"
       " correspond to independent operators or not. (op/op_dir/src/src_dir)");
   return params;
 }

Member Data Documentation

◆ _base

std::vector<std::vector<DenseVector<Real> > >& PODReducedBasisTrainer::_base

protected

The reduced basis for the variables.

Definition at line 120 of file PODReducedBasisTrainer.h.

Referenced by addToReducedOperator(), computeBasisVectors(), getBaseSize(), getBasisVector(), getSumBaseSize(), getVariableIndex(), and initialSetup().

◆ _base_completed

bool PODReducedBasisTrainer::_base_completed

protected

Switch that tells if the object has already computed the necessary basis vectors.

This switch is used in execute() to determine if we want to compute basis vectors or do something else.

Definition at line 128 of file PODReducedBasisTrainer.h.

Referenced by execute().

◆ _corr_mx

std::vector<DenseMatrix<Real> > PODReducedBasisTrainer::_corr_mx

protected

The correlation matrices for the variables.

Definition at line 111 of file PODReducedBasisTrainer.h.

Referenced by computeCorrelationMatrix(), computeEigenDecomposition(), and initialSetup().

◆ _eigenvalues

std::vector<DenseVector<Real> > PODReducedBasisTrainer::_eigenvalues

protected

The eigenvalues of the correalation matrix for each variable.

Definition at line 114 of file PODReducedBasisTrainer.h.

Referenced by computeBasisVectors(), computeEigenDecomposition(), initialSetup(), and printEigenvalues().

◆ _eigenvectors

std::vector<DenseMatrix<Real> > PODReducedBasisTrainer::_eigenvectors

protected

The eigenvectors of the correalation matrix for each variable.

Definition at line 117 of file PODReducedBasisTrainer.h.

Referenced by computeBasisVectors(), computeEigenDecomposition(), and initialSetup().

◆ _empty_operators

bool PODReducedBasisTrainer::_empty_operators

protected

Flag describing if the reduced operators are empty or not.

This is necessary to be able to do both the base and reduced operator generation in the same object.

Definition at line 133 of file PODReducedBasisTrainer.h.

Referenced by addToReducedOperator(), execute(), and finalize().

◆ _error_res

std::vector<Real> PODReducedBasisTrainer::_error_res

protected

Energy limits that define how many basis functions will be kept for each variable.

Definition at line 99 of file PODReducedBasisTrainer.h.

Referenced by computeEigenDecomposition(), and PODReducedBasisTrainer().

◆ _red_operators

std::vector<DenseMatrix<Real> >& PODReducedBasisTrainer::_red_operators

protected

The reduced operators that should be transferred to the surrogate.

Definition at line 123 of file PODReducedBasisTrainer.h.

Referenced by addToReducedOperator(), finalize(), and initReducedOperators().

◆ _snapshots

std::vector<DistributedSnapshots> PODReducedBasisTrainer::_snapshots

protected

Distributed container for snapshots per variable.

Definition at line 108 of file PODReducedBasisTrainer.h.

Referenced by addSnapshot(), computeBasisVectors(), computeCorrelationMatrix(), getSnapsSize(), initialSetup(), and receiveObjects().

◆ _tag_names

std::vector<std::string>& PODReducedBasisTrainer::_tag_names

protected

Names of the tags that should be used to fetch residuals from the MultiApp.

Definition at line 102 of file PODReducedBasisTrainer.h.

Referenced by getTagNames(), initReducedOperators(), and PODReducedBasisTrainer().

◆ _tag_types

std::vector<std::string>& PODReducedBasisTrainer::_tag_types

protected

list of bools describing which tag is indepedent of the solution.

Definition at line 105 of file PODReducedBasisTrainer.h.

Referenced by getTagTypes(), initReducedOperators(), and PODReducedBasisTrainer().

◆ _var_names

std::vector<std::string>& PODReducedBasisTrainer::_var_names

protected

Vector containing the names of the variables we want to use for constructing the surrogates.

Definition at line 96 of file PODReducedBasisTrainer.h.

Referenced by addToReducedOperator(), getBasisVector(), getVariableIndex(), getVarNames(), initialSetup(), PODReducedBasisTrainer(), and printEigenvalues().

The documentation for this class was generated from the following files:

stochastic_tools/include/trainers/PODReducedBasisTrainer.h
stochastic_tools/src/trainers/PODReducedBasisTrainer.C

Public Types

Public Member Functions

Static Public Member Functions

Public Attributes

Static Public Attributes

Protected Member Functions

Static Protected Member Functions

Protected Attributes

Static Protected Attributes

Private Member Functions

Detailed Description

Constructor & Destructor Documentation

◆ PODReducedBasisTrainer()

Member Function Documentation

◆ addSnapshot()

◆ addToReducedOperator()

◆ computeBasisVectors()

◆ computeCorrelationMatrix()

◆ computeEigenDecomposition()

◆ declareModelData()

◆ determineNumberOfModes()

◆ execute()

◆ finalize()

◆ getBaseSize()

◆ getBasisVector() [1/2]

◆ getBasisVector() [2/2]

◆ getModelData()

◆ getModelDataFileName()

◆ getSnapsSize()

◆ getSumBaseSize()

◆ getTagNames()

◆ getTagTypes()

◆ getVariableIndex()

◆ getVarNames()

◆ hasModelData()

◆ initialize()

◆ initialSetup()

◆ initReducedOperators()

◆ modelMetaDataName()

◆ printEigenvalues()

◆ receiveObjects()

◆ threadJoin()

◆ validParams()

Member Data Documentation

◆ _base

◆ _base_completed

◆ _corr_mx

◆ _eigenvalues

◆ _eigenvectors

◆ _empty_operators

◆ _error_res

◆ _red_operators

◆ _snapshots

◆ _tag_names

◆ _tag_types

◆ _var_names