Mesh-based parameter optimization. More...

#include <ParameterMeshOptimization.h>

Inheritance diagram for ParameterMeshOptimization:

Public Types
typedef DataFileName	DataFileParameterType

Public Member Functions
	ParameterMeshOptimization (const InputParameters &parameters)

virtual Real	computeObjective () override
	Function to compute objective.

virtual void	computeGradient (libMesh::PetscVector< Number > &gradient) const override
	Function to compute gradient.

virtual void	initialize () override

virtual void	execute () override

virtual void	finalize () override

void	setInitialCondition (libMesh::PetscVector< Number > &param)
	Function to initialize petsc vectors from vpp data.

virtual void	setMisfitToSimulatedValues ()
	Function to override misfit values with the simulated values from the matrix free hessian forward solve.

Real	getUpperBound (dof_id_type i) const
	Upper and lower bounds for each parameter being controlled.

Real	getLowerBound (dof_id_type i) const

virtual void	computeEqualityConstraints (libMesh::PetscVector< Number > &eqs_constraints) const
	Function to compute the equality constraints.

virtual void	computeInequalityConstraints (libMesh::PetscVector< Number > &ineqs_constraints) const
	Function to compute the inequality constraints.

virtual void	computeEqualityGradient (libMesh::PetscMatrix< Number > &gradient) const
	Function to compute the gradient of the equality constraints/ This is the last call of the equality constraint gradient routine.

virtual void	computeInequalityGradient (libMesh::PetscMatrix< Number > &gradient) const
	Function to compute the gradient of the inequality constraints/ This is the last call of the inequality constraint gradient routine.

dof_id_type	getNumEqCons () const
	Function to get the total number of equalities.

dof_id_type	getNumInEqCons () const
	Function to get the total number of inequalities.

void	threadJoin (const UserObject &) final

bool	shouldStore () const override final

virtual Real	spatialValue (const Point &) const

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

void	setPrimaryThreadCopy (UserObject *primary)

UserObject *	primaryThreadCopy ()

SubProblem &	getSubProblem () const

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

std::set< UserObjectName >	getDependObjects () const

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

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

unsigned int	systemNumber () const

virtual bool	needThreadedCopy () const

virtual bool	enabled () const

std::shared_ptr< MooseObject >	getSharedPtr ()

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

bool	isKokkosObject () const

MooseApp &	getMooseApp () const

const std::string &	type () const

const std::string &	name () const

std::string	typeAndName () const

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

MooseObjectName	uniqueName () const

const InputParameters &	parameters () const

const hit::Node *	getHitNode () const

bool	hasBase () const

const std::string &	getBase () 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 &name) const

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

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

void	paramWarning (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

std::string	messagePrefix (const bool hit_prefix=true) const

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

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

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

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

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

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

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

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

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

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

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

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

void	callMooseError (std::string msg, const bool with_prefix, const hit::Node *node=nullptr, const bool show_trace=true) 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	initialSetup ()

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 UserObjectBase &	getUserObjectBase (const std::string &param_name, bool is_dependency=true) const

const UserObjectBase &	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)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyByName (const std::string &prop_name_in)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOldByName (const std::string &prop_name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOlderByName (const std::string &prop_name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialProperty (const std::string &name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOld (const std::string &name)

Moose::Kokkos::MaterialProperty< T, dimension >	getKokkosMaterialPropertyOlder (const std::string &name)

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

std::pair< Moose::Kokkos::MaterialProperty< T, dimension >, std::set< SubdomainID > >	getKokkosBlockMaterialProperty (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)

virtual 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	store (nlohmann::json &json) const

virtual void	declareLateValues ()

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

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

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	hasKokkosMaterialProperty (const std::string &name)

bool	hasKokkosMaterialProperty (const std::string &name)

bool	hasKokkosMaterialPropertyByName (const std::string &name)

bool	hasKokkosMaterialPropertyByName (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

Moose::Kokkos::Function	getKokkosFunction (const std::string &name) const

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

Moose::Kokkos::Function	getKokkosFunctionByName (const FunctionName &name) const

const T &	getKokkosFunctionByName (const FunctionName &name) const

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

bool	hasKokkosFunctionByName (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

Static Public Member Functions
static InputParameters	validParams ()

static void	callMooseError (MooseApp const app, const InputParameters &params, std::string msg, const bool with_prefix, const hit::Node node, const bool show_trace=true)

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, NameFunc &&name_func)

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

Public Attributes
	usingCombinedWarningSolutionWarnings

const ConsoleStream	_console

Static Public Attributes
static const std::string	type_param

static const std::string	name_param

static const std::string	unique_name_param

static const std::string	app_param

static const std::string	moose_base_param

static const std::string	kokkos_object_param

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	setICsandBounds () override
	Sets the initial conditions and bounds right before it is needed.

virtual dof_id_type	getNumParams () const override
	Function to get the total number of parameters.

virtual void	updateParameters (const libMesh::PetscVector< Number > &x)
	Function to set parameters.

std::vector< Real >	parseInputData (std::string type, Real default_value, unsigned int param_id) const
	Function to to parse bounds and initial conditions from input file.

virtual void	addPostprocessorDependencyHelper (const PostprocessorName &name) const override

virtual void	addVectorPostprocessorDependencyHelper (const VectorPostprocessorName &name) const override

virtual void	addUserObjectDependencyHelper (const UserObjectBase &uo) const override

void	addReporterDependencyHelper (const ReporterName &reporter_name) override

void	flagInvalidSolutionInternal (const InvalidSolutionID invalid_solution_id) const

InvalidSolutionID	registerInvalidSolutionInternal (const std::string &message, const bool warning) const

const ReporterContextBase &	getReporterContextBaseByName (const ReporterName &reporter_name) const

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)

virtual void	getKokkosMaterialPropertyHook (const std::string &, const unsigned int)

void	markMatPropRequested (const std::string &)

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

void	checkExecutionStage ()

T &	declareUnusedValue (Args &&... args)

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)

T &	declareValue (const std::string &param_name, Args &&... args)

T &	declareValue (const std::string &param_name, ReporterMode mode, Args &&... args)

T &	declareValue (const std::string &param_name, Args &&... args)

T &	declareValue (const std::string &param_name, ReporterMode mode, Args &&... args)

T &	declareValue (const std::string &param_name, Args &&... args)

T &	declareValue (const std::string &param_name, ReporterMode mode, Args &&... args)

T &	declareValue (const std::string &param_name, Args &&... args)

T &	declareValue (const std::string &param_name, ReporterMode mode, Args &&... args)

T &	declareValueByName (const ReporterValueName &value_name, Args &&... args)

T &	declareValueByName (const ReporterValueName &value_name, ReporterMode mode, Args &&... args)

T &	declareValueByName (const ReporterValueName &value_name, Args &&... args)

T &	declareValueByName (const ReporterValueName &value_name, ReporterMode mode, Args &&... args)

T &	declareValueByName (const ReporterValueName &value_name, Args &&... args)

T &	declareValueByName (const ReporterValueName &value_name, ReporterMode mode, Args &&... args)

T &	declareValueByName (const ReporterValueName &value_name, Args &&... args)

T &	declareValueByName (const ReporterValueName &value_name, ReporterMode mode, Args &&... args)

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

Protected Attributes
Real &	_objective_val
	Reporter that will hold the objective value.

std::vector< dof_id_type > *	_num_values_reporter

const std::vector< ReporterValueName > &	_parameter_names
	Parameter names.

const unsigned int	_nparams
	Number of parameter vectors.

std::vector< std::vector< Real > * >	_parameters
	Parameter values declared as reporter data.

std::vector< std::vector< Real > * >	_gradients
	Gradient values declared as reporter data.

const Real	_tikhonov_coeff
	Tikhonov Coefficient for regularization.

const std::vector< ReporterValueName > *	_equality_names
	Equality constraint names.

const unsigned int	_n_eq_cons
	Number of equality constraint names.

std::vector< std::vector< Real > * >	_eq_constraints
	Equality values declared as reporter data.

std::vector< std::vector< Real > * >	_eq_gradients
	Gradient values declared as reporter data.

const std::vector< ReporterValueName > *	_inequality_names
	Inequality constraint names.

const unsigned int	_n_ineq_cons
	Number of inequality constraint names.

std::vector< std::vector< Real > * >	_ineq_constraints
	Inequality values declared as reporter data.

std::vector< std::vector< Real > * >	_ineq_gradients
	Gradient values declared as reporter data.

std::vector< Real >	_lower_bounds
	Bounds of the parameters.

std::vector< Real >	_upper_bounds

std::vector< dof_id_type >	_nvalues
	Number of values for each parameter.

dof_id_type	_ndof
	Total number of parameters.

const Moose::CoordinateSystemType &	_coord_sys

const THREAD_ID	_tid

SubProblem &	_subproblem

FEProblemBase &	_fe_problem

SystemBase &	_sys

Assembly &	_assembly

const bool	_duplicate_initial_execution

std::set< std::string >	_depend_uo

const bool &	_enabled

MooseApp &	_app

Factory &	_factory

ActionFactory &	_action_factory

const std::string &	_type

const std::string &	_name

const InputParameters &	_pars

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 bool	_is_kokkos_object

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
std::vector< Real >	parseExodusData (const FEType fetype, const FileName mesh_file_name, const std::vector< unsigned int > &exodus_timestep, const std::string &mesh_var_name) const
	Read initialization data off of parameter mesh and error check.

const Moose::FunctionBase *	getKokkosFunctionByNameHelper (const FunctionName &name) const

const UserObjectBase &	getUserObjectFromFEProblem (const UserObjectName &object_name, const THREAD_ID tid=0) const

const T &	castUserObject (const UserObjectBase &uo_base, const std::string &param_name="") const

void	mooseObjectError (const std::string &param_name, std::stringstream &oss) const

const std::string &	userObjectType (const UserObjectBase &uo) const

const std::string &	userObjectName (const UserObjectBase &uo) const

const PostprocessorName &	getPostprocessorNameInternal (const std::string &param_name, const unsigned int index, const bool allow_default_value=true) const

bool	isDefaultPostprocessorValueByName (const PostprocessorName &name) const

PostprocessorValue	getDefaultPostprocessorValueByName (const PostprocessorName &name) const

void	checkParam (const std::string &param_name, const unsigned int index=std::numeric_limits< unsigned int >::max()) const

bool	postprocessorsAdded () const

const VectorPostprocessorValue &	getVectorPostprocessorByNameHelper (const VectorPostprocessorName &name, const std::string &vector_name, bool broadcast, std::size_t t_index) const

const VectorPostprocessorContext< VectorPostprocessorValue > &	getVectorPostprocessorContextByNameHelper (const VectorPostprocessorName &name, const std::string &vector_name) const

bool	vectorPostprocessorsAdded () const

bool	reportersAdded () const

void	possiblyCheckHasReporter (const ReporterName &reporter_name, const std::string &param_name="") const

RestartableDataValue &	registerRestartableDataOnApp (std::unique_ptr< RestartableDataValue > data, THREAD_ID tid) const

void	registerRestartableNameWithFilterOnApp (const std::string &name, Moose::RESTARTABLE_FILTER filter)

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

virtual std::string	meshPropertyPrefix (const std::string &data_name) const

const RestartableDataValue &	getMeshPropertyInternal (const std::string &data_name, const std::string &prefix) const

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

const VariableValue *	getDefaultValue (const std::string &var_name) const

const ADVariableValue *	getADDefaultValue (const std::string &var_name) const

void	checkVar (const std::string &var_name) const

void	validateExecutionerType (const std::string &name, const std::string &fn_name) const

Moose::MaterialDataType	getMaterialDataType (const std::set< BoundaryID > &boundary_ids) const

unsigned int	getMaxQps () const

void	addConsumedPropertyName (const MooseObjectName &obj_name, const std::string &prop_name)

const ReporterValueName &	getReporterValueName (const std::string &param_name) const

const PostprocessorValue &	getPostprocessorValueInternal (const std::string &param_name, unsigned int index, std::size_t t_index) const

const PostprocessorValue &	getPostprocessorValueInternal (const std::string &param_name, unsigned int index, std::size_t t_index) const

const PostprocessorValue &	getPostprocessorValueByNameInternal (const PostprocessorName &name, std::size_t t_index) const

const PostprocessorValue &	getPostprocessorValueByNameInternal (const PostprocessorName &name, std::size_t t_index) const

void	possiblyCheckHasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const

void	possiblyCheckHasVectorPostprocessor (const std::string &param_name, const std::string &vector_name) const

void	possiblyCheckHasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const

void	possiblyCheckHasVectorPostprocessorByName (const VectorPostprocessorName &name, const std::string &vector_name) const

Static Private Member Functions
static const hit::Node *	getHitNode (const InputParameters &params)

static std::string	messagePrefix (const InputParameters &params, const bool hit_prefix)

Private Attributes
std::vector< std::unique_ptr< ParameterMesh > >	_parameter_meshes
	Store parameter meshes for regularization computation.

const std::vector< Real >	_regularization_coeffs
	Vector of regularization coefficients corresponding to each type.

const std::vector< ParameterMesh::RegularizationType >	_regularization_types
	Regularization types to apply.

const bool	_always_store

UserObject *	_primary_thread_copy

std::set< std::string >	_supplied_uo

const ParallelParamObject &	_parent

const MooseBase &	_si_moose_base

const FEProblemBase *	_si_problem

ExecFlagEnum	_empty_execute_enum

const MooseObject &	_fni_object

const InputParameters &	_fni_params

FEProblemBase &	_fni_feproblem

const THREAD_ID	_fni_tid

const MooseObject &	_uoi_moose_object

const FEProblemBase &	_uoi_feproblem

const THREAD_ID	_uoi_tid

const MooseObject &	_ppi_moose_object

const InputParameters &	_ppi_params

const FEProblemBase &	_ppi_feproblem

std::map< PostprocessorName, std::unique_ptr< PostprocessorValue > >	_default_values

const bool	_broadcast_by_default

const MooseObject &	_vpi_moose_object

const FEProblemBase &	_vpi_feproblem

const THREAD_ID	_vpi_tid

const InputParameters &	_ri_params

FEProblemBase &	_ri_fe_problem_base

const ReporterData &	_ri_reporter_data

const MooseObject &	_ri_moose_object

const InputParameters &	_dni_params

FEProblemBase &	_dni_feproblem

const MooseObject *const	_dni_moose_object_ptr

const InputParameters &	_si_params

FEProblemBase &	_si_feproblem

THREAD_ID	_si_tid

const RestartableDataMapName	_metaname

std::string	_restartable_name

MooseApp &	_meta_data_app

const MooseObject *const	_meta_data_object

const InputParameters &	_sc_parameters

const std::string &	_sc_name

const bool	_sc_is_implicit

std::unordered_map< std::string, std::vector< MooseVariableScalar * > >	_coupled_scalar_vars

std::unordered_map< std::string, std::unique_ptr< VariableValue > >	_default_value

std::unordered_map< std::string, std::unique_ptr< ADVariableValue > >	_dual_default_value

std::vector< MooseVariableScalar * >	_coupled_moose_scalar_vars

std::unordered_map< std::string, std::vector< MooseVariableFieldBase * > >	_sc_coupled_vars

std::set< TagID >	_sc_coupleable_vector_tags

std::set< TagID >	_sc_coupleable_matrix_tags

const MooseObject &	_mi_moose_object

const bool	_mi_boundary_restricted

const std::set< SubdomainID > &	_mi_block_ids

const std::set< BoundaryID > &	_mi_boundary_ids

std::vector< std::unique_ptr< OptionalMaterialPropertyProxyBase< MaterialPropertyInterface > > >	_optional_property_proxies

const std::string	_ti_name

const MooseObject &	_reporter_moose_object

const InputParameters &	_reporter_params

const std::string &	_reporter_name

FEProblemBase &	_reporter_fe_problem

ReporterData &	_reporter_data

std::vector< std::unique_ptr< UnusedWrapperBase > >	_unused_values

MooseApp &	_oi_moose_app

OutputWarehouse &	_oi_output_warehouse

std::set< OutputName >	_oi_outputs

Detailed Description

Mesh-based parameter optimization.

Definition at line 18 of file ParameterMeshOptimization.h.

Constructor & Destructor Documentation

◆ ParameterMeshOptimization()

ParameterMeshOptimization::ParameterMeshOptimization ( const InputParameters & parameters )

Definition at line 90 of file ParameterMeshOptimization.C.

  : GeneralOptimization(parameters),
    _regularization_coeffs(getParam<std::vector<Real>>("regularization_coeffs")),
    _regularization_types(getParam<MultiMooseEnum>("regularization_types")
                              .getSetValueIDs<ParameterMesh::RegularizationType>())
{
  // Validate that regularization coefficients match types
  if (_regularization_coeffs.size() != _regularization_types.size())
    paramError("regularization_coeffs",
               "Number of regularization coefficients (",
               _regularization_coeffs.size(),
               ") must match number of regularization types (",
               _regularization_types.size(),
               ")");
}

Member Function Documentation

◆ computeEqualityConstraints()

void OptimizationReporterBase::computeEqualityConstraints ( libMesh::PetscVector< Number > & eqs_constraints ) const

virtualinherited

Function to compute the equality constraints.

This is the last call of the equality function routine.

Reimplemented in QuadraticMinimizeConstrained.

Definition at line 159 of file OptimizationReporterBase.C.

{
  OptUtils::copyReporterIntoPetscVector(_eq_constraints, eqs_constraints);
}

Referenced by OptimizeSolve::equalityFunctionWrapper().

◆ computeEqualityGradient()

void OptimizationReporterBase::computeEqualityGradient ( libMesh::PetscMatrix< Number > & gradient ) const

virtualinherited

Function to compute the gradient of the equality constraints/ This is the last call of the equality constraint gradient routine.

Reimplemented in QuadraticMinimizeConstrained.

Definition at line 173 of file OptimizationReporterBase.C.

{
  for (const auto & p : make_range(_n_eq_cons))
    if (_eq_gradients[p]->size() != _ndof)
      mooseError("The equality jacobian for parameter ",
                 _parameter_names[p],
                 " has changed, expected ",
                 _ndof,
                 " versus ",
                 _eq_gradients[p]->size(),
                 ".");
  OptUtils::copyReporterIntoPetscMatrix(_eq_gradients, jacobian);
}

Referenced by OptimizeSolve::equalityGradientFunctionWrapper().

◆ computeGradient()

void ParameterMeshOptimization::computeGradient ( libMesh::PetscVector< Number > & gradient ) const

overridevirtual

Function to compute gradient.

This is the last call of the gradient routine.

Reimplemented from OptimizationReporterBase.

Definition at line 288 of file ParameterMeshOptimization.C.

{
  // Add regularization gradient contributions to the reporter gradients before base computation
  for (const auto reg_idx : index_range(_regularization_types))
  {
    if (_regularization_coeffs[reg_idx] > 0.0)
    {
      // Convert MultiMooseEnum to RegularizationType using get() method
      ParameterMesh::RegularizationType reg_type = _regularization_types[reg_idx];
 
      for (const auto & param_id : make_range(_nparams))
      {
        // Get current parameter values for this group
        const auto & param_values = *_parameters[param_id];
        auto grad_values = _gradients[param_id];
 
        // Compute regularization gradient for this type
        std::vector<Real> reg_grad =
            _parameter_meshes[param_id]->computeRegularizationGradient(param_values, reg_type);
 
        // Add to gradient with coefficient
        for (unsigned int i = 0; i < param_values.size(); ++i)
          (*grad_values)[i] += _regularization_coeffs[reg_idx] * reg_grad[i];
      }
    }
  }
 
  // Now call base class method which includes Tikhonov and copies to PETSc vector
  OptimizationReporterBase::computeGradient(gradient);
}

◆ computeInequalityConstraints()

void OptimizationReporterBase::computeInequalityConstraints ( libMesh::PetscVector< Number > & ineqs_constraints ) const

virtualinherited

Function to compute the inequality constraints.

This is the last call of the inequality function routine.

Definition at line 166 of file OptimizationReporterBase.C.

{
  OptUtils::copyReporterIntoPetscVector(_ineq_constraints, ineqs_constraints);
}

Referenced by OptimizeSolve::inequalityFunctionWrapper().

◆ computeInequalityGradient()

void OptimizationReporterBase::computeInequalityGradient ( libMesh::PetscMatrix< Number > & gradient ) const

virtualinherited

Function to compute the gradient of the inequality constraints/ This is the last call of the inequality constraint gradient routine.

Definition at line 188 of file OptimizationReporterBase.C.

{
  for (const auto & p : make_range(_n_ineq_cons))
    if (_ineq_gradients[p]->size() != _ndof)
      mooseError("The inequality jacobian for parameter ",
                 _parameter_names[p],
                 " has changed, expected ",
                 _ndof,
                 " versus ",
                 _ineq_gradients[p]->size(),
                 ".");
  OptUtils::copyReporterIntoPetscMatrix(_ineq_gradients, jacobian);
}

Referenced by OptimizeSolve::inequalityGradientFunctionWrapper().

◆ computeObjective()

Real ParameterMeshOptimization::computeObjective ( )

overridevirtual

Function to compute objective.

This is the last function called in objective routine

Reimplemented from GeneralOptimization.

Definition at line 256 of file ParameterMeshOptimization.C.

{
  Real val = GeneralOptimization::computeObjective();
 
  // Apply each regularization type with its coefficient
  for (const auto reg_idx : index_range(_regularization_types))
  {
    if (_regularization_coeffs[reg_idx] > 0.0)
    {
      Real regularization_value = 0.0;
 
      // Convert MultiMooseEnum to RegularizationType using get() method
      ParameterMesh::RegularizationType reg_type = _regularization_types[reg_idx];
 
      for (const auto & param_id : make_range(_nparams))
      {
        // Get current parameter values for this group
        const auto & param_values = *_parameters[param_id];
 
        // Compute regularization objective for this type
        regularization_value +=
            _parameter_meshes[param_id]->computeRegularizationObjective(param_values, reg_type);
      }
 
      val += _regularization_coeffs[reg_idx] * regularization_value;
    }
  }
 
  return val;
}

◆ execute()

virtual void OptimizationReporterBase::execute ( )

inlineoverridevirtualinherited

Implements GeneralReporter.

Reimplemented in OptimizationDataTempl< OptimizationReporterBase >.

Definition at line 32 of file OptimizationReporterBase.h.

32{}

◆ finalize()

virtual void OptimizationReporterBase::finalize ( )

inlineoverridevirtualinherited

Implements GeneralReporter.

Reimplemented in OptimizationDataTempl< OptimizationReporterBase >.

Definition at line 33 of file OptimizationReporterBase.h.

33{}

◆ getLowerBound()

Real OptimizationReporterBase::getLowerBound ( dof_id_type i ) const

inherited

Definition at line 112 of file OptimizationReporterBase.C.

{
  return _lower_bounds[i];
}

Referenced by OptimizationReporterTest::initialSetup(), and OptimizeSolve::variableBounds().

◆ getNumEqCons()

dof_id_type OptimizationReporterBase::getNumEqCons ( ) const

inlineinherited

Function to get the total number of equalities.

Returns: total number of equality constraints

Definition at line 100 of file OptimizationReporterBase.h.

100{ return _n_eq_cons; }

Referenced by OptimizeSolve::taoALCreate(), and OptimizeSolve::taoALDestroy().

◆ getNumInEqCons()

dof_id_type OptimizationReporterBase::getNumInEqCons ( ) const

inlineinherited

Function to get the total number of inequalities.

Returns: total number of inequalities constraints

Definition at line 106 of file OptimizationReporterBase.h.

106{ return _n_ineq_cons; }

Referenced by OptimizeSolve::taoALCreate(), and OptimizeSolve::taoALDestroy().

◆ getNumParams()

dof_id_type GeneralOptimization::getNumParams ( ) const

overrideprotectedvirtualinherited

Function to get the total number of parameters.

Returns: total number of parameters

Reimplemented from OptimizationReporterBase.

Definition at line 82 of file GeneralOptimization.C.

{
  if (_ndof == 0)
    mooseError(
        "The number of parameters you have is zero and this shouldn't happen. Make sure you are "
        "running your forward problem on \'INITIAL\' if you are using a reporter transfer to "
        "supply "
        "that information.");
 
  return _ndof;
}

◆ getUpperBound()

Real OptimizationReporterBase::getUpperBound ( dof_id_type i ) const

inherited

Upper and lower bounds for each parameter being controlled.

Parameters

i	Parameter index

Returns: The upper/lower bound for parameter i

Definition at line 118 of file OptimizationReporterBase.C.

{
  return _upper_bounds[i];
}

Referenced by OptimizationReporterTest::initialSetup(), and OptimizeSolve::variableBounds().

◆ initialize()

virtual void OptimizationReporterBase::initialize ( )

inlineoverridevirtualinherited

Implements GeneralReporter.

Reimplemented in OptimizationDataTempl< OptimizationReporterBase >.

Definition at line 31 of file OptimizationReporterBase.h.

31{}

◆ parseExodusData()

std::vector< Real > ParameterMeshOptimization::parseExodusData	(	const FEType	fetype,
		const FileName	mesh_file_name,
		const std::vector< unsigned int > &	exodus_timestep,
		const std::string &	mesh_var_name
	)		const

private

Read initialization data off of parameter mesh and error check.

Returns: values read from mesh

Definition at line 107 of file ParameterMeshOptimization.C.

{
  // read data off Exodus mesh
  ReadExodusMeshVars data_mesh(fetype, mesh_file_name, mesh_var_name);
  std::vector<Real> parsed_data;
  // read from mesh
  for (auto const & step : exodus_timestep)
  {
    std::vector<Real> data = data_mesh.getParameterValues(step);
    parsed_data.insert(parsed_data.end(), data.begin(), data.end());
  }
 
  return parsed_data;
}

Referenced by setICsandBounds().

◆ parseInputData()

std::vector< Real > OptimizationReporterBase::parseInputData	(	std::string	type,
		Real	default_value,
		unsigned int	param_id
	)		const

protectedinherited

Function to to parse bounds and initial conditions from input file.

Definition at line 124 of file OptimizationReporterBase.C.

{
  // fill with default values
  std::vector<Real> parsed_data_id(_nvalues[param_id], default_value);
  if (isParamValid(type))
  {
    std::vector<std::vector<Real>> parsed_data(getParam<std::vector<std::vector<Real>>>(type));
    parsed_data_id.assign(parsed_data[param_id].begin(), parsed_data[param_id].end());
    if (parsed_data.size() != _nvalues.size())
      paramError(type,
                 "There must be a vector of ",
                 type,
                 " per parameter group.  The ",
                 type,
                 " input format is std::vector<std::vector<Real>> so each vector should be "
                 "seperated by \";\" even if it is a single value per group for a constant ",
                 type,
                 ".");
    // The case when the initial condition is constant for each parameter group
    if (parsed_data[param_id].size() == 1)
      parsed_data_id.assign(_nvalues[param_id], parsed_data[param_id][0]);
    else if (parsed_data[param_id].size() != _nvalues[param_id])
      paramError(type,
                 "When ",
                 type,
                 " are given in input file, there must either be a single value per parameter "
                 "group or a value for every parameter in the group.");
  }
 
  return parsed_data_id;
}

Referenced by GeneralOptimization::setICsandBounds(), OptimizationReporter::setICsandBounds(), and setICsandBounds().

◆ setICsandBounds()

void ParameterMeshOptimization::setICsandBounds ( )

overrideprotectedvirtual

Sets the initial conditions and bounds right before it is needed.

Reimplemented from GeneralOptimization.

Definition at line 126 of file ParameterMeshOptimization.C.

{
  if ((isParamValid("num_values_name") || isParamValid("num_values")))
    paramError("num_values_name or num_values should not be used with ParameterMeshOptimization. "
               "Instead the number of dofs is set by the parameter meshes.");
 
  _nvalues.resize(_nparams, 0);
  // Fill the mesh information
  const auto & meshes = getParam<std::vector<FileName>>("parameter_meshes");
  const auto & families = getParam<MultiMooseEnum>("parameter_families");
  const auto & orders = getParam<MultiMooseEnum>("parameter_orders");
  const auto & ntimes = getParam<unsigned int>("num_parameter_times");
 
  // Fill exodus parameter bounds and IC information
  std::vector<std::string> initial_condition_mesh_variable;
  std::vector<std::string> lower_bound_mesh_variable;
  std::vector<std::string> upper_bound_mesh_variable;
  if (isParamValid("initial_condition_mesh_variable"))
    initial_condition_mesh_variable =
        getParam<std::vector<std::string>>("initial_condition_mesh_variable");
  if (isParamValid("lower_bound_mesh_variable"))
    lower_bound_mesh_variable = getParam<std::vector<std::string>>("lower_bound_mesh_variable");
  if (isParamValid("upper_bound_mesh_variable"))
    upper_bound_mesh_variable = getParam<std::vector<std::string>>("upper_bound_mesh_variable");
 
  std::vector<unsigned int> exodus_timestep;
  if (isParamValid("exodus_timesteps_for_parameter_mesh_variable"))
    exodus_timestep =
        getParam<std::vector<unsigned int>>("exodus_timesteps_for_parameter_mesh_variable");
  else // get last timestep in file
    exodus_timestep = {std::numeric_limits<unsigned int>::max()};
 
  // now do a bunch of error checking
  // Size checks for data
  if (meshes.size() != _nparams)
    paramError("parameter_meshes",
               "There must be a mesh associated with each group of parameters.");
  if (families.size() > 1 && families.size() != _nparams)
    paramError("parameter_families",
               "There must be a family associated with each group of parameters.");
  if (orders.size() > 1 && orders.size() != _nparams)
    paramError("parameter_orders",
               "There must be an order associated with each group of parameters.");
 
  // error checking that initial conditions and bounds are only read from a single location
  if (isParamValid("initial_condition_mesh_variable") && isParamValid("initial_condition"))
    paramError("initial_condition_mesh_variable",
               "Initial conditions for all parameter groups can only be defined by "
               "initial_condition_mesh_variable or "
               "initial_condition but not both.");
  else if (isParamValid("lower_bound_mesh_variable") && isParamValid("lower_bounds"))
    paramError(
        "lower_bound_mesh_variable",
        "Lower bounds for all parameter groups can only be defined by lower_bound_mesh_variable or "
        "lower_bounds but not both.");
  else if (isParamValid("upper_bound_mesh_variable") && isParamValid("upper_bounds"))
    paramError(
        "upper_bound_mesh_variable",
        "Upper bounds for all parameter groups can only be defined by upper_bound_mesh_variable or "
        "upper_bounds but not both.");
 
  // Make sure they did not specify too many timesteps
  if (isParamValid("exodus_timesteps_for_parameter_mesh_variable") &&
      (!isParamValid("lower_bound_mesh_variable") + !isParamValid("upper_bound_mesh_variable") +
           !isParamValid("initial_condition_mesh_variable") ==
       3))
    paramError("\"exodus_timesteps_for_parameter_mesh_variable\" should only be specified if "
               "reading values from a mesh.");
  else if (exodus_timestep.size() != ntimes && exodus_timestep.size() != 1)
    paramError("exodus_timesteps_for_parameter_mesh_variable",
               "Number of timesteps to read mesh data specified by "
               "\"exodus_timesteps_for_parameter_mesh_variable\" incorrect. "
               "\"exodus_timesteps_for_parameter_mesh_variable\" can specify a single timestep or "
               "\"num_parameter_times\" timesteps.");
 
  _ndof = 0;
  _parameter_meshes.resize(_nparams);
  for (const auto & param_id : make_range(_nparams))
  {
    const std::string family = families.size() > 1 ? families[param_id] : families[0];
    const std::string order = orders.size() > 1 ? orders[param_id] : orders[0];
    const FEType fetype(Utility::string_to_enum<Order>(order),
                        Utility::string_to_enum<FEFamily>(family));
 
    _parameter_meshes[param_id] = std::make_unique<ParameterMesh>(fetype, meshes[param_id]);
    _nvalues[param_id] = _parameter_meshes[param_id]->size() * ntimes;
    _ndof += _nvalues[param_id];
 
    // read and assign initial conditions
    {
      std::vector<Real> initial_condition;
      if (isParamValid("initial_condition_mesh_variable"))
        initial_condition = parseExodusData(
            fetype, meshes[param_id], exodus_timestep, initial_condition_mesh_variable[param_id]);
      else
        initial_condition = parseInputData("initial_condition", 0, param_id);
 
      _parameters[param_id]->assign(initial_condition.begin(), initial_condition.end());
    }
 
    // read and assign lower bound
    {
      std::vector<Real> lower_bound;
      if (isParamValid("lower_bound_mesh_variable"))
        lower_bound = parseExodusData(
            fetype, meshes[param_id], exodus_timestep, lower_bound_mesh_variable[param_id]);
      else
        lower_bound = parseInputData("lower_bounds", std::numeric_limits<Real>::lowest(), param_id);
 
      _lower_bounds.insert(_lower_bounds.end(), lower_bound.begin(), lower_bound.end());
    }
 
    // read and assign upper bound
    {
      std::vector<Real> upper_bound;
      if (isParamValid("upper_bound_mesh_variable"))
        upper_bound = parseExodusData(
            fetype, meshes[param_id], exodus_timestep, upper_bound_mesh_variable[param_id]);
      else
        upper_bound = parseInputData("upper_bounds", std::numeric_limits<Real>::max(), param_id);
 
      _upper_bounds.insert(_upper_bounds.end(), upper_bound.begin(), upper_bound.end());
    }
 
    // resize gradient vector to be filled later
    _gradients[param_id]->resize(_nvalues[param_id]);
  }
}

◆ setInitialCondition()

void OptimizationReporterBase::setInitialCondition ( libMesh::PetscVector< Number > & param )

inherited

Function to initialize petsc vectors from vpp data.

Definition at line 98 of file OptimizationReporterBase.C.

{
  setICsandBounds();
  x.init(getNumParams());
  OptUtils::copyReporterIntoPetscVector(_parameters, x);
}

Referenced by OptimizationReporterTest::initialSetup(), and OptimizeSolve::solve().

◆ setMisfitToSimulatedValues()

virtual void OptimizationReporterBase::setMisfitToSimulatedValues ( )

inlinevirtualinherited

Function to override misfit values with the simulated values from the matrix free hessian forward solve.

Reimplemented in OptimizationReporter.

Definition at line 55 of file OptimizationReporterBase.h.

55{}

Referenced by OptimizeSolve::applyHessian().

◆ updateParameters()

void OptimizationReporterBase::updateParameters ( const libMesh::PetscVector< Number > & x )

protectedvirtualinherited

Function to set parameters.

This is the first function called in objective/gradient/hessian routine

Definition at line 106 of file OptimizationReporterBase.C.

{
  OptUtils::copyPetscVectorIntoReporter(x, _parameters);
}

Referenced by OptimizeSolve::applyHessian(), OptimizationReporterTest::execute(), OptimizeSolve::gradientFunction(), and OptimizeSolve::objectiveFunction().

◆ validParams()

InputParameters ParameterMeshOptimization::validParams ( )

static

Definition at line 24 of file ParameterMeshOptimization.C.

{
  InputParameters params = GeneralOptimization::validParams();
 
  params.addClassDescription(
      "Computes objective function, gradient and contains reporters for communicating between "
      "optimizeSolve and subapps using mesh-based parameter definition.");
 
  params.addRequiredParam<std::vector<FileName>>(
      "parameter_meshes", "Exodus file containing meshes describing parameters.");
 
  const auto family = AddVariableAction::getNonlinearVariableFamilies();
  MultiMooseEnum families(family.getRawNames(), "LAGRANGE");
  params.addParam<MultiMooseEnum>(
      "parameter_families",
      families,
      "Specifies the family of FE shape functions for each group of parameters. If a single value "
      "is "
      "specified, then that value is used for all groups of parameters.");
  const auto order = AddVariableAction::getNonlinearVariableOrders();
  MultiMooseEnum orders(order.getRawNames(), "FIRST");
  params.addParam<MultiMooseEnum>(
      "parameter_orders",
      orders,
      "Specifies the order of FE shape functions for each group of parameters. If a single value "
      "is "
      "specified, then that value is used for all groups of parameters.");
 
  params.addParam<unsigned int>(
      "num_parameter_times", 1, "The number of time points the parameters represent.");
 
  params.addParam<std::vector<std::string>>(
      "initial_condition_mesh_variable",
      "Name of variable on parameter mesh to use as initial condition.");
  params.addParam<std::vector<std::string>>(
      "lower_bound_mesh_variable", "Name of variable on parameter mesh to use as lower bound.");
  params.addParam<std::vector<std::string>>(
      "upper_bound_mesh_variable", "Name of variable on parameter mesh to use as upper bound.");
  params.addParam<std::vector<unsigned int>>(
      "exodus_timesteps_for_parameter_mesh_variable",
      "Timesteps to read all parameter group bounds and initial conditions from Exodus mesh.  The "
      "options are to give no timestep, a single timestep or \"num_parameter_times\" timesteps.  "
      "No timestep results in the final timestep from the mesh being used.  A single timestep "
      "results in values at that timestep being used for all timesteps.  \"num_parameter_times\" "
      "timesteps results in values from the mesh at those steps being used.  The same timesteps "
      "are used for all parameter groups and all meshes, the capability to define different "
      "timesteps for different meshes is not supported.");
 
  // New parameters for multiple regularization types
  MultiMooseEnum reg_types("L2_GRADIENT");
  params.addParam<MultiMooseEnum>(
      "regularization_types",
      reg_types,
      "Types of regularization to apply. Multiple types can be specified.");
 
  params.addParam<std::vector<Real>>("regularization_coeffs",
                                     {},
                                     "Coefficients for each regularization type. Must match the "
                                     "number of regularization_types specified.");
 
  params.addParamNamesToGroup("tikhonov_coeff regularization_types regularization_coeffs",
                              "Regularization");
 
  return params;
}

Member Data Documentation

◆ _eq_constraints

std::vector<std::vector<Real> *> OptimizationReporterBase::_eq_constraints

protectedinherited

Equality values declared as reporter data.

Definition at line 142 of file OptimizationReporterBase.h.

Referenced by OptimizationReporterBase::computeEqualityConstraints(), and OptimizationReporterBase::OptimizationReporterBase().

◆ _eq_gradients

std::vector<std::vector<Real> *> OptimizationReporterBase::_eq_gradients

protectedinherited

Gradient values declared as reporter data.

Definition at line 144 of file OptimizationReporterBase.h.

Referenced by OptimizationReporterBase::computeEqualityGradient(), and OptimizationReporterBase::OptimizationReporterBase().

◆ _equality_names

const std::vector<ReporterValueName>* OptimizationReporterBase::_equality_names

protectedinherited

Equality constraint names.

Definition at line 138 of file OptimizationReporterBase.h.

Referenced by OptimizationReporterBase::OptimizationReporterBase().

◆ _gradients

std::vector<std::vector<Real> *> OptimizationReporterBase::_gradients

protectedinherited

Gradient values declared as reporter data.

Definition at line 132 of file OptimizationReporterBase.h.

Referenced by OptimizationReporterBase::computeGradient(), computeGradient(), OptimizationReporterBase::OptimizationReporterBase(), GeneralOptimization::setICsandBounds(), OptimizationReporter::setICsandBounds(), and setICsandBounds().

◆ _ineq_constraints

std::vector<std::vector<Real> *> OptimizationReporterBase::_ineq_constraints

protectedinherited

Inequality values declared as reporter data.

Definition at line 151 of file OptimizationReporterBase.h.

Referenced by OptimizationReporterBase::computeInequalityConstraints(), and OptimizationReporterBase::OptimizationReporterBase().

◆ _ineq_gradients

std::vector<std::vector<Real> *> OptimizationReporterBase::_ineq_gradients

protectedinherited

Gradient values declared as reporter data.

Definition at line 153 of file OptimizationReporterBase.h.

Referenced by OptimizationReporterBase::computeInequalityGradient(), and OptimizationReporterBase::OptimizationReporterBase().

◆ _inequality_names

const std::vector<ReporterValueName>* OptimizationReporterBase::_inequality_names

protectedinherited

Inequality constraint names.

Definition at line 147 of file OptimizationReporterBase.h.

Referenced by OptimizationReporterBase::OptimizationReporterBase().

◆ _lower_bounds

std::vector<Real> OptimizationReporterBase::_lower_bounds

protectedinherited

Bounds of the parameters.

Definition at line 156 of file OptimizationReporterBase.h.

Referenced by OptimizationReporterBase::getLowerBound(), GeneralOptimization::setICsandBounds(), OptimizationReporter::setICsandBounds(), and setICsandBounds().

◆ _n_eq_cons

const unsigned int OptimizationReporterBase::_n_eq_cons

protectedinherited

Number of equality constraint names.

Definition at line 140 of file OptimizationReporterBase.h.

Referenced by OptimizationReporterBase::computeEqualityGradient(), QuadraticMinimizeConstrained::computeEqualityGradient(), OptimizationReporterBase::getNumEqCons(), and OptimizationReporterBase::OptimizationReporterBase().

◆ _n_ineq_cons

const unsigned int OptimizationReporterBase::_n_ineq_cons

protectedinherited

Number of inequality constraint names.

Definition at line 149 of file OptimizationReporterBase.h.

Referenced by OptimizationReporterBase::computeInequalityGradient(), OptimizationReporterBase::getNumInEqCons(), and OptimizationReporterBase::OptimizationReporterBase().

◆ _ndof

dof_id_type OptimizationReporterBase::_ndof

protectedinherited

Total number of parameters.

Definition at line 162 of file OptimizationReporterBase.h.

Referenced by OptimizationReporterBase::computeEqualityGradient(), OptimizationReporterBase::computeInequalityGradient(), OptimizationReporterBase::getNumParams(), GeneralOptimization::getNumParams(), QuadraticMinimize::QuadraticMinimize(), GeneralOptimization::setICsandBounds(), OptimizationReporter::setICsandBounds(), and setICsandBounds().

◆ _nparams

const unsigned int OptimizationReporterBase::_nparams

protectedinherited

Number of parameter vectors.

Definition at line 127 of file OptimizationReporterBase.h.

Referenced by OptimizationReporterBase::computeGradient(), computeGradient(), computeObjective(), OptimizationReporterBase::OptimizationReporterBase(), GeneralOptimization::setICsandBounds(), OptimizationReporter::setICsandBounds(), and setICsandBounds().

◆ _num_values_reporter

std::vector<dof_id_type>* GeneralOptimization::_num_values_reporter

protectedinherited

Definition at line 32 of file GeneralOptimization.h.

Referenced by GeneralOptimization::setICsandBounds().

◆ _nvalues

std::vector<dof_id_type> OptimizationReporterBase::_nvalues

protectedinherited

Number of values for each parameter.

Definition at line 160 of file OptimizationReporterBase.h.

Referenced by OptimizationReporterBase::computeGradient(), OptimizationReporterBase::parseInputData(), GeneralOptimization::setICsandBounds(), OptimizationReporter::setICsandBounds(), and setICsandBounds().

◆ _objective_val

Real& GeneralOptimization::_objective_val

protectedinherited

Reporter that will hold the objective value.

Definition at line 30 of file GeneralOptimization.h.

Referenced by GeneralOptimization::computeObjective().

◆ _parameter_meshes

std::vector<std::unique_ptr<ParameterMesh> > ParameterMeshOptimization::_parameter_meshes

private

Store parameter meshes for regularization computation.

Definition at line 42 of file ParameterMeshOptimization.h.

Referenced by computeGradient(), computeObjective(), and setICsandBounds().

◆ _parameter_names

const std::vector<ReporterValueName>& OptimizationReporterBase::_parameter_names

protectedinherited

Parameter names.

Definition at line 125 of file OptimizationReporterBase.h.

Referenced by OptimizationReporterBase::computeEqualityGradient(), OptimizationReporterBase::computeGradient(), OptimizationReporterBase::computeInequalityGradient(), OptimizationReporterBase::OptimizationReporterBase(), GeneralOptimization::setICsandBounds(), and OptimizationReporter::setICsandBounds().

◆ _parameters

std::vector<std::vector<Real> *> OptimizationReporterBase::_parameters

protectedinherited

Parameter values declared as reporter data.

Definition at line 130 of file OptimizationReporterBase.h.

Referenced by QuadraticMinimizeConstrained::computeEqualityConstraints(), QuadraticMinimizeConstrained::computeEqualityGradient(), OptimizationReporterBase::computeGradient(), computeGradient(), QuadraticMinimize::computeGradient(), QuadraticMinimizeConstrained::computeGradient(), GeneralOptimization::computeObjective(), OptimizationReporter::computeObjective(), computeObjective(), QuadraticMinimize::computeObjective(), QuadraticMinimizeConstrained::computeObjective(), OptimizationReporterBase::OptimizationReporterBase(), GeneralOptimization::setICsandBounds(), OptimizationReporter::setICsandBounds(), setICsandBounds(), OptimizationReporterBase::setInitialCondition(), and OptimizationReporterBase::updateParameters().

◆ _regularization_coeffs

const std::vector<Real> ParameterMeshOptimization::_regularization_coeffs

private

Vector of regularization coefficients corresponding to each type.

Definition at line 45 of file ParameterMeshOptimization.h.

Referenced by computeGradient(), computeObjective(), and ParameterMeshOptimization().

◆ _regularization_types

const std::vector<ParameterMesh::RegularizationType> ParameterMeshOptimization::_regularization_types

private

Regularization types to apply.

Definition at line 48 of file ParameterMeshOptimization.h.

Referenced by computeGradient(), computeObjective(), and ParameterMeshOptimization().

◆ _tikhonov_coeff

const Real OptimizationReporterBase::_tikhonov_coeff

protectedinherited

Tikhonov Coefficient for regularization.

Definition at line 135 of file OptimizationReporterBase.h.

Referenced by OptimizationReporterBase::computeGradient(), GeneralOptimization::computeObjective(), and OptimizationReporter::computeObjective().

◆ _upper_bounds

std::vector<Real> OptimizationReporterBase::_upper_bounds

protectedinherited

Definition at line 157 of file OptimizationReporterBase.h.

Referenced by OptimizationReporterBase::getUpperBound(), GeneralOptimization::setICsandBounds(), OptimizationReporter::setICsandBounds(), and setICsandBounds().

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

optimization/include/optimizationreporters/ParameterMeshOptimization.h
optimization/src/optimizationreporters/ParameterMeshOptimization.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

Static Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ ParameterMeshOptimization()

Member Function Documentation

◆ computeEqualityConstraints()

◆ computeEqualityGradient()

◆ computeGradient()

◆ computeInequalityConstraints()

◆ computeInequalityGradient()

◆ computeObjective()

◆ execute()

◆ finalize()

◆ getLowerBound()

◆ getNumEqCons()

◆ getNumInEqCons()

◆ getNumParams()

◆ getUpperBound()

◆ initialize()

◆ parseExodusData()

◆ parseInputData()

◆ setICsandBounds()

◆ setInitialCondition()

◆ setMisfitToSimulatedValues()

◆ updateParameters()

◆ validParams()

Member Data Documentation

◆ _eq_constraints

◆ _eq_gradients

◆ _equality_names

◆ _gradients

◆ _ineq_constraints

◆ _ineq_gradients

◆ _inequality_names

◆ _lower_bounds

◆ _n_eq_cons

◆ _n_ineq_cons

◆ _ndof

◆ _nparams

◆ _num_values_reporter

◆ _nvalues

◆ _objective_val

◆ _parameter_meshes

◆ _parameter_names

◆ _parameters

◆ _regularization_coeffs

◆ _regularization_types

◆ _tikhonov_coeff

◆ _upper_bounds