https://mooseframework.inl.gov
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 | List of all members
GenericActiveLearnerTempl< SamplerType > Class Template Reference

A generic reporter to support parallel active learning: re-trains GP and picks the next best batch. More...

#include <GenericActiveLearner.h>

Inheritance diagram for GenericActiveLearnerTempl< SamplerType >:
[legend]

Public Types

typedef DataFileName DataFileParameterType
 

Public Member Functions

 GenericActiveLearnerTempl (const InputParameters &parameters)
 
virtual void initialize () override
 
virtual void finalize () override
 
virtual void execute () override
 
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)
 
UserObjectprimaryThreadCopy ()
 
SubProblemgetSubProblem () 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< MooseObjectgetSharedPtr ()
 
std::shared_ptr< const MooseObjectgetSharedPtr () const
 
bool isKokkosObject () const
 
MooseAppgetMooseApp () 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 InputParametersparameters () const
 
const hit::Node * getHitNode () const
 
bool hasBase () const
 
const std::string & getBase () const
 
const TgetParam (const std::string &name) const
 
std::vector< std::pair< T1, T2 > > getParam (const std::string &param1, const std::string &param2) const
 
const TqueryParam (const std::string &name) const
 
const TgetRenamedParam (const std::string &old_name, const std::string &new_name) const
 
T getCheckedPointerParam (const std::string &name, const std::string &error_string="") const
 
bool haveParameter (const std::string &name) 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 ExecFlagEnumgetExecuteOnEnum () const
 
UserObjectName getUserObjectName (const std::string &param_name) const
 
const TgetUserObject (const std::string &param_name, bool is_dependency=true) const
 
const TgetUserObjectByName (const UserObjectName &object_name, bool is_dependency=true) const
 
const UserObjectBasegetUserObjectBase (const std::string &param_name, bool is_dependency=true) const
 
const UserObjectBasegetUserObjectBaseByName (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< SubdomainIDgetMaterialPropertyBlocks (const std::string &name)
 
std::vector< SubdomainName > getMaterialPropertyBlockNames (const std::string &name)
 
std::set< BoundaryIDgetMaterialPropertyBoundaryIDs (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)
 
MaterialBasegetMaterial (const std::string &name)
 
MaterialBasegetMaterial (const std::string &name)
 
MaterialBasegetMaterialByName (const std::string &name, bool no_warn=false)
 
MaterialBasegetMaterialByName (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 FunctiongetFunction (const std::string &name) const
 
const FunctiongetFunctionByName (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 TgetKokkosFunction (const std::string &name) const
 
Moose::Kokkos::Function getKokkosFunctionByName (const FunctionName &name) const
 
const TgetKokkosFunctionByName (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 VectorPostprocessorValuegetVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name, bool needs_broadcast) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const VectorPostprocessorValuegetVectorPostprocessorValueOldByName (const VectorPostprocessorName &name, const std::string &vector_name, bool needs_broadcast) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValue (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueByName (const VectorPostprocessorName &name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueOld (const std::string &param_name, const std::string &vector_name) const
 
const ScatterVectorPostprocessorValuegetScatterVectorPostprocessorValueOldByName (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
 
TgetSampler (const std::string &name)
 
SamplergetSampler (const std::string &name)
 
TgetSamplerByName (const SamplerName &name)
 
SamplergetSamplerByName (const SamplerName &name)
 
virtual void meshChanged ()
 
virtual void meshDisplaced ()
 
PerfGraphperfGraph ()
 
const PostprocessorValuegetPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValue (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOld (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
 
const PostprocessorValuegetPostprocessorValueOlder (const std::string &param_name, const unsigned int index=0) const
 
virtual const PostprocessorValuegetPostprocessorValueByName (const PostprocessorName &name) const
 
virtual const PostprocessorValuegetPostprocessorValueByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOldByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOldByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOlderByName (const PostprocessorName &name) const
 
const PostprocessorValuegetPostprocessorValueOlderByName (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 DistributiongetDistribution (const std::string &name) const
 
const TgetDistribution (const std::string &name) const
 
const DistributiongetDistribution (const std::string &name) const
 
const TgetDistribution (const std::string &name) const
 
const DistributiongetDistributionByName (const DistributionName &name) const
 
const TgetDistributionByName (const std::string &name) const
 
const DistributiongetDistributionByName (const DistributionName &name) const
 
const TgetDistributionByName (const std::string &name) const
 
const Parallel::Communicator & comm () const
 
processor_id_type n_processors () const
 
processor_id_type processor_id () const
 
template<>
SurrogateModelgetSurrogateModel (const std::string &name) const
 
template<>
SurrogateTrainerBasegetSurrogateTrainer (const std::string &name) const
 
template<>
SurrogateModelgetSurrogateModelByName (const UserObjectName &name) const
 
template<>
SurrogateTrainerBasegetSurrogateTrainerByName (const UserObjectName &name) const
 
template<typename T = SurrogateModel>
TgetSurrogateModel (const std::string &name) const
 Get a SurrogateModel/Trainer with a given name. More...
 
template<typename T = SurrogateTrainerBase>
TgetSurrogateTrainer (const std::string &name) const
 
template<typename T = SurrogateModel>
TgetSurrogateModelByName (const UserObjectName &name) const
 Get a sampler with a given name. More...
 
template<typename T = SurrogateTrainerBase>
TgetSurrogateTrainerByName (const UserObjectName &name) 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 setupGPData (const std::vector< Real > &data_out, const DenseMatrix< Real > &data_in)
 Sets up the training data for the GP model. More...
 
virtual void computeGPOutput (std::vector< Real > &eval_outputs)
 Computes the outputs of the trained GP model. More...
 
virtual Real computeConvergenceValue ()
 Computes the convergence value during active learning. More...
 
virtual void evaluateGPTest ()
 Evaluate the GP on all the test samples sent by the Sampler. More...
 
virtual void setupGeneric ()
 Setup the generic variable for acquisition computation (depends on the objective: optimization, UQ, etc.) More...
 
virtual void includeAdditionalInputs ()
 Include additional inputs before evaluating the acquisition function. More...
 
virtual void getAcquisition (std::vector< Real > &acq_new, std::vector< unsigned int > &indices)
 Output the acquisition function values and ordering of the indices. More...
 
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 ReporterContextBasegetReporterContextBaseByName (const ReporterName &reporter_name) const
 
const ReporterNamegetReporterName (const std::string &param_name) const
 
TdeclareRestartableData (const std::string &data_name, Args &&... args)
 
ManagedValue< TdeclareManagedRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
 
const TgetRestartableData (const std::string &data_name) const
 
TdeclareRestartableDataWithContext (const std::string &data_name, void *context, Args &&... args)
 
TdeclareRecoverableData (const std::string &data_name, Args &&... args)
 
TdeclareRestartableDataWithObjectName (const std::string &data_name, const std::string &object_name, Args &&... args)
 
TdeclareRestartableDataWithObjectNameWithContext (const std::string &data_name, const std::string &object_name, void *context, Args &&... args)
 
std::string restartableName (const std::string &data_name) const
 
const TgetMeshProperty (const std::string &data_name, const std::string &prefix)
 
const TgetMeshProperty (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 VariableValuecoupledScalarValue (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledScalarValue (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 VariableValuecoupledVectorTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
 
const VariableValuecoupledMatrixTagScalarValue (const std::string &var_name, TagID tag, unsigned int comp=0) const
 
const VariableValuecoupledScalarValueOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarValueOlder (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDot (const std::string &var_name, unsigned int comp=0) const
 
const ADVariableValueadCoupledScalarDot (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDot (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDotOld (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDu (const std::string &var_name, unsigned int comp=0) const
 
const VariableValuecoupledScalarDotDotDu (const std::string &var_name, unsigned int comp=0) const
 
const MooseVariableScalargetScalarVar (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 ()
 
TdeclareUnusedValue (Args &&... args)
 
const TgetReporterValue (const std::string &param_name, const std::size_t time_index=0)
 
const TgetReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)
 
const TgetReporterValue (const std::string &param_name, const std::size_t time_index=0)
 
const TgetReporterValue (const std::string &param_name, ReporterMode mode, const std::size_t time_index=0)
 
const TgetReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)
 
const TgetReporterValueByName (const ReporterName &reporter_name, ReporterMode mode, const std::size_t time_index=0)
 
const TgetReporterValueByName (const ReporterName &reporter_name, const std::size_t time_index=0)
 
const TgetReporterValueByName (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)
 
TdeclareValue (const std::string &param_name, Args &&... args)
 
TdeclareValue (const std::string &param_name, ReporterMode mode, Args &&... args)
 
TdeclareValue (const std::string &param_name, Args &&... args)
 
TdeclareValue (const std::string &param_name, ReporterMode mode, Args &&... args)
 
TdeclareValue (const std::string &param_name, Args &&... args)
 
TdeclareValue (const std::string &param_name, ReporterMode mode, Args &&... args)
 
TdeclareValue (const std::string &param_name, Args &&... args)
 
TdeclareValue (const std::string &param_name, ReporterMode mode, Args &&... args)
 
TdeclareValueByName (const ReporterValueName &value_name, Args &&... args)
 
TdeclareValueByName (const ReporterValueName &value_name, ReporterMode mode, Args &&... args)
 
TdeclareValueByName (const ReporterValueName &value_name, Args &&... args)
 
TdeclareValueByName (const ReporterValueName &value_name, ReporterMode mode, Args &&... args)
 
TdeclareValueByName (const ReporterValueName &value_name, Args &&... args)
 
TdeclareValueByName (const ReporterValueName &value_name, ReporterMode mode, Args &&... args)
 
TdeclareValueByName (const ReporterValueName &value_name, Args &&... args)
 
TdeclareValueByName (const ReporterValueName &value_name, ReporterMode mode, Args &&... args)
 
ParallelAcquisitionFunctionBasegetParallelAcquisitionFunctionByName (const UserObjectName &name) const
 Lookup a ParallelAcquisitionFunction object by name and return pointer. More...
 

Static Protected Member Functions

static std::string meshPropertyName (const std::string &data_name, const std::string &prefix)
 

Protected Attributes

SamplerType & _al_sampler
 The base sampler. More...
 
unsigned int _n_dim
 The input dimension for GP, equal to Sampler columns. More...
 
dof_id_type _props
 Storage for the number of parallel proposals. More...
 
const std::vector< std::vector< Real > > & _inputs_test
 Storage for all the proposed samples to test the GP model. More...
 
const std::vector< Real > & _output_value
 Model output value from SubApp. More...
 
std::vector< Real > & _output_comm
 Modified value of model output by this reporter class. More...
 
std::vector< unsigned int > & _sorted_indices
 The selected sample indices to evaluate the subApp. More...
 
const ActiveLearningGaussianProcess_al_gp
 The active learning GP trainer that permits re-training. More...
 
const SurrogateModel_gp_eval
 The GP evaluator object that permits re-evaluations. More...
 
ParallelAcquisitionFunctionBase_acquisition_obj
 Storage for the parallel acquisition object to be utilized. More...
 
std::vector< Real > & _acquisition_value
 The acquistion function values in the current iteration. More...
 
Real_convergence_value
 For monitoring convergence of active learning. More...
 
std::vector< std::vector< Real > > _inputs_test_modified
 Storage for all the modified proposed samples to test the GP model. More...
 
std::vector< std::vector< Real > > & _inputs_required
 Transmit the required inputs to the json file. More...
 
const bool & _penalize_acquisition
 Penalize acquisition to prevent clustering when operating in parallel. More...
 
int _check_step
 Ensure that the MCMC algorithm proceeds in a sequential fashion. More...
 
std::vector< std::vector< Real > > _gp_inputs
 Storage for the GP re-training inputs. More...
 
std::vector< Real_gp_outputs
 Storage for the GP re-training outputs. More...
 
std::vector< Real_gp_outputs_test
 Outputs of GP model for the test samples. More...
 
std::vector< Real_gp_std_test
 Outputs of GP model standard deviation for the test samples. More...
 
std::vector< Real_length_scales
 Storage for the length scales after the GP training. More...
 
std::vector< Real_generic
 A generic parameter to be passed to the acquisition function. More...
 
std::vector< Real_eval_outputs_current
 The GP outputs from the current iteration before re-training (to evaluate convergence) More...
 
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
 

Detailed Description

template<typename SamplerType>
class GenericActiveLearnerTempl< SamplerType >

A generic reporter to support parallel active learning: re-trains GP and picks the next best batch.

Definition at line 24 of file GenericActiveLearner.h.

Constructor & Destructor Documentation

◆ GenericActiveLearnerTempl()

template<typename SamplerType >
GenericActiveLearnerTempl< SamplerType >::GenericActiveLearnerTempl ( const InputParameters parameters)

Definition at line 197 of file GenericActiveLearner.h.

202  _al_sampler(getSampler<SamplerType>("sampler")),
203  _n_dim(_al_sampler.getNumberOfCols()),
204  _props(_al_sampler.getNumParallelProposals()),
205  _inputs_test(_al_sampler.getSampleTries()),
206  _output_value(getReporterValue<std::vector<Real>>("output_value", REPORTER_MODE_DISTRIBUTED)),
207  _output_comm(declareValue<std::vector<Real>>("outputs_required")),
208  _sorted_indices(declareValue<std::vector<unsigned int>>("sorted_indices")),
209  _al_gp(getUserObject<ActiveLearningGaussianProcess>("al_gp")),
210  _gp_eval(getSurrogateModel<GaussianProcessSurrogate>("gp_evaluator")),
211  _acquisition_obj(getParallelAcquisitionFunctionByName(getParam<UserObjectName>("acquisition"))),
212  _acquisition_value(declareValue<std::vector<Real>>("acquisition_function")),
213  _convergence_value(declareValue<Real>("convergence_value")),
214  _inputs_required(declareValue<std::vector<std::vector<Real>>>("inputs")),
215  _penalize_acquisition(getParam<bool>("penalize_acquisition")),
216  _check_step(std::numeric_limits<int>::max())
217 {
218  // Setting up the variable sizes to facilitate active learning.
219  _gp_outputs_test.resize(_inputs_test.size());
220  _gp_std_test.resize(_inputs_test.size());
221  _acquisition_value.resize(_props);
222  _length_scales.resize(_n_dim);
224  _generic.resize(1);
225  _inputs_required.resize(_props, std::vector<Real>(_n_dim, 0.0));
226  _sorted_indices.resize(_props, 1u);
227 }
const std::vector< Real > & _output_value
Model output value from SubApp.
std::vector< Real > _eval_outputs_current
The GP outputs from the current iteration before re-training (to evaluate convergence) ...
const InputParameters & parameters() const
std::vector< Real > & _acquisition_value
The acquistion function values in the current iteration.
ParallelAcquisitionFunctionBase & getParallelAcquisitionFunctionByName(const UserObjectName &name) const
Lookup a ParallelAcquisitionFunction object by name and return pointer.
const SurrogateModel & _gp_eval
The GP evaluator object that permits re-evaluations.
std::vector< Real > _length_scales
Storage for the length scales after the GP training.
dof_id_type _props
Storage for the number of parallel proposals.
unsigned int _n_dim
The input dimension for GP, equal to Sampler columns.
const T & getReporterValue(const std::string &param_name, const std::size_t time_index=0)
std::vector< unsigned int > & _sorted_indices
The selected sample indices to evaluate the subApp.
const std::vector< std::vector< Real > > & _inputs_test
Storage for all the proposed samples to test the GP model.
const bool & _penalize_acquisition
Penalize acquisition to prevent clustering when operating in parallel.
Real & _convergence_value
For monitoring convergence of active learning.
std::vector< Real > & _output_comm
Modified value of model output by this reporter class.
GeneralReporter(const InputParameters &parameters)
ParallelAcquisitionFunctionBase & _acquisition_obj
Storage for the parallel acquisition object to be utilized.
T & declareValue(const std::string &param_name, Args &&... args)
ParallelAcquisitionInterface(const InputParameters &parameters)
const ActiveLearningGaussianProcess & _al_gp
The active learning GP trainer that permits re-training.
std::vector< Real > _gp_std_test
Outputs of GP model standard deviation for the test samples.
SamplerType & _al_sampler
The base sampler.
std::vector< std::vector< Real > > & _inputs_required
Transmit the required inputs to the json file.
std::vector< Real > _generic
A generic parameter to be passed to the acquisition function.
int _check_step
Ensure that the MCMC algorithm proceeds in a sequential fashion.
SurrogateModelInterface(const MooseObject *moose_object)
std::vector< Real > _gp_outputs_test
Outputs of GP model for the test samples.

Member Function Documentation

◆ computeConvergenceValue()

template<typename SamplerType >
Real GenericActiveLearnerTempl< SamplerType >::computeConvergenceValue ( )
protectedvirtual

Computes the convergence value during active learning.

Reimplemented in BayesianActiveLearner.

Definition at line 283 of file GenericActiveLearner.h.

284 {
285  Real convergence_value = 0.0;
286  for (unsigned int ii = 0; ii < _output_comm.size(); ++ii)
287  convergence_value += Utility::pow<2>(_output_comm[ii] - _eval_outputs_current[ii]);
288  convergence_value = std::sqrt(convergence_value) / _output_comm.size();
289  return convergence_value;
290 }
std::vector< Real > _eval_outputs_current
The GP outputs from the current iteration before re-training (to evaluate convergence) ...
std::vector< Real > & _output_comm
Modified value of model output by this reporter class.
DIE A HORRIBLE DEATH HERE typedef LIBMESH_DEFAULT_SCALAR_TYPE Real

◆ computeGPOutput()

template<typename SamplerType >
void GenericActiveLearnerTempl< SamplerType >::computeGPOutput ( std::vector< Real > &  eval_outputs)
protectedvirtual

Computes the outputs of the trained GP model.

Parameters
eval_outputsThe outputs predicted by the GP model

Definition at line 245 of file GenericActiveLearner.h.

246 {
247  for (unsigned int i = 0; i < eval_outputs.size(); ++i)
248  eval_outputs[i] = _gp_eval.evaluate(_gp_inputs[i]);
249 }
std::vector< std::vector< Real > > _gp_inputs
Storage for the GP re-training inputs.
const SurrogateModel & _gp_eval
The GP evaluator object that permits re-evaluations.
virtual Real evaluate(const std::vector< Real > &x) const
Evaluate surrogate model given a row of parameters.

◆ evaluateGPTest()

template<typename SamplerType >
void GenericActiveLearnerTempl< SamplerType >::evaluateGPTest ( )
protectedvirtual

Evaluate the GP on all the test samples sent by the Sampler.

Reimplemented in BayesianActiveLearner.

Definition at line 294 of file GenericActiveLearner.h.

295 {
296  for (unsigned int i = 0; i < _gp_outputs_test.size(); ++i)
298 }
const SurrogateModel & _gp_eval
The GP evaluator object that permits re-evaluations.
const std::vector< std::vector< Real > > & _inputs_test
Storage for all the proposed samples to test the GP model.
virtual Real evaluate(const std::vector< Real > &x) const
Evaluate surrogate model given a row of parameters.
std::vector< Real > _gp_std_test
Outputs of GP model standard deviation for the test samples.
std::vector< Real > _gp_outputs_test
Outputs of GP model for the test samples.

◆ execute()

template<typename SamplerType >
void GenericActiveLearnerTempl< SamplerType >::execute ( )
overridevirtual

Implements GeneralReporter.

Definition at line 302 of file GenericActiveLearner.h.

303 {
304  if (_al_sampler.getNumberOfLocalRows() == 0 || _check_step == _t_step)
305  {
307  return;
308  }
309 
310  DenseMatrix<Real> data_in(_al_sampler.getNumberOfRows(), _al_sampler.getNumberOfCols());
311  for (dof_id_type ss = _al_sampler.getLocalRowBegin(); ss < _al_sampler.getLocalRowEnd(); ++ss)
312  {
313  const auto data = _al_sampler.getNextLocalRow();
314  for (unsigned int j = 0; j < _al_sampler.getNumberOfCols(); ++j)
315  data_in(ss, j) = data[j];
316  }
317  _communicator.sum(data_in.get_values());
320 
321  if (_t_step > 1)
322  {
323  // Setup the GP training data
324  setupGPData(_output_comm, data_in);
325 
326  // Compute the convergence value before re-training the GP
327  if (_t_step > 2)
328  {
331  }
332 
333  // Retrain the GP and get the length scales
336 
337  // Evaluate the GP on all the test samples sent by the Sampler
338  evaluateGPTest();
339 
340  // Setup the generic variable for acquisition computation (depends on the objective:
341  // optimization, UQ, etc.)
342  setupGeneric();
343 
344  // Get the acquisition function values and ordering of indices as per the acquisition
345  std::vector<Real> acq_new;
346  std::vector<unsigned int> indices;
347  indices.resize(_inputs_test.size());
348  getAcquisition(acq_new, indices);
349 
350  // Output the acquisition function values and the best ordering of the indices
351  std::copy_n(indices.begin(), _props, _sorted_indices.begin());
352  std::copy_n(acq_new.begin(), _props, _acquisition_value.begin());
353  }
354  else
355  std::iota(_sorted_indices.begin(), _sorted_indices.end(), 0);
356 
357  // Track the current step
359 }
void allgather(const T &send_data, std::vector< T, A > &recv_data) const
std::vector< std::vector< Real > > _gp_inputs
Storage for the GP re-training inputs.
const std::vector< Real > & _output_value
Model output value from SubApp.
virtual void reTrain(const std::vector< std::vector< Real >> &inputs, const std::vector< Real > &outputs) const final
virtual void computeGPOutput(std::vector< Real > &eval_outputs)
Computes the outputs of the trained GP model.
std::vector< Real > _eval_outputs_current
The GP outputs from the current iteration before re-training (to evaluate convergence) ...
std::vector< Real > & _acquisition_value
The acquistion function values in the current iteration.
std::vector< Real > _length_scales
Storage for the length scales after the GP training.
const Parallel::Communicator & _communicator
const std::vector< Real > & getLengthScales() const
Return the current length scales from GP training.
dof_id_type _props
Storage for the number of parallel proposals.
std::vector< unsigned int > & _sorted_indices
The selected sample indices to evaluate the subApp.
const std::vector< std::vector< Real > > & _inputs_test
Storage for all the proposed samples to test the GP model.
virtual void evaluateGPTest()
Evaluate the GP on all the test samples sent by the Sampler.
Real & _convergence_value
For monitoring convergence of active learning.
std::vector< Real > & _output_comm
Modified value of model output by this reporter class.
virtual void setupGeneric()
Setup the generic variable for acquisition computation (depends on the objective: optimization...
virtual Real computeConvergenceValue()
Computes the convergence value during active learning.
virtual void setupGPData(const std::vector< Real > &data_out, const DenseMatrix< Real > &data_in)
Sets up the training data for the GP model.
virtual void getAcquisition(std::vector< Real > &acq_new, std::vector< unsigned int > &indices)
Output the acquisition function values and ordering of the indices.
const ActiveLearningGaussianProcess & _al_gp
The active learning GP trainer that permits re-training.
std::vector< Real > _gp_outputs
Storage for the GP re-training outputs.
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
SamplerType & _al_sampler
The base sampler.
int _check_step
Ensure that the MCMC algorithm proceeds in a sequential fashion.
uint8_t dof_id_type

◆ finalize()

template<typename SamplerType>
virtual void GenericActiveLearnerTempl< SamplerType >::finalize ( )
inlineoverridevirtual

Implements GeneralReporter.

Definition at line 42 of file GenericActiveLearner.h.

42 {}

◆ getAcquisition()

template<typename SamplerType >
void GenericActiveLearnerTempl< SamplerType >::getAcquisition ( std::vector< Real > &  acq_new,
std::vector< unsigned int > &  indices 
)
protectedvirtual

Output the acquisition function values and ordering of the indices.

Parameters
acq_newThe computed values of the acquisition function
indicesThe indices ordered according to the acqusition values to be sent to Sampler

Definition at line 267 of file GenericActiveLearner.h.

269 {
270  std::vector<Real> acq;
271  acq.resize(_inputs_test.size());
275  acq_new = acq;
278  acq_new, indices, acq, _length_scales, _inputs_test_modified);
279 }
std::vector< std::vector< Real > > _inputs_test_modified
Storage for all the modified proposed samples to test the GP model.
std::vector< std::vector< Real > > _gp_inputs
Storage for the GP re-training inputs.
std::vector< Real > _length_scales
Storage for the length scales after the GP training.
void penalizeAcquisition(std::vector< Real > &modified_acq, std::vector< unsigned int > &sorted_indices, const std::vector< Real > &acq, const std::vector< Real > &length_scales, const std::vector< std::vector< Real >> &inputs)
Return the modified acquisition function values and sorted indices considering local penalization (in...
const std::vector< std::vector< Real > > & _inputs_test
Storage for all the proposed samples to test the GP model.
const bool & _penalize_acquisition
Penalize acquisition to prevent clustering when operating in parallel.
virtual void includeAdditionalInputs()
Include additional inputs before evaluating the acquisition function.
ParallelAcquisitionFunctionBase & _acquisition_obj
Storage for the parallel acquisition object to be utilized.
std::vector< Real > _gp_std_test
Outputs of GP model standard deviation for the test samples.
std::vector< Real > _generic
A generic parameter to be passed to the acquisition function.
std::vector< Real > _gp_outputs_test
Outputs of GP model for the test samples.
void computeAcquisition(std::vector< Real > &acq, const std::vector< Real > &gp_mean, const std::vector< Real > &gp_std, const std::vector< std::vector< Real >> &test_inputs, const std::vector< std::vector< Real >> &train_inputs, const std::vector< Real > &generic) const
Compute the acquisition function values.

◆ getParallelAcquisitionFunctionByName()

ParallelAcquisitionFunctionBase & ParallelAcquisitionInterface::getParallelAcquisitionFunctionByName ( const UserObjectName &  name) const
protectedinherited

Lookup a ParallelAcquisitionFunction object by name and return pointer.

Definition at line 24 of file ParallelAcquisitionInterface.C.

26 {
27  std::vector<ParallelAcquisitionFunctionBase *> models;
29  .query()
30  .condition<AttribName>(name)
31  .condition<AttribSystem>("ParallelAcquisitionFunctionBase")
32  .queryInto(models);
33 
34  if (models.empty())
35  mooseError("Unable to find a ParallelAcquisitionFunction object with the name '" + name + "'");
36 
37  return *models.front();
38 }
void mooseError(Args &&... args)
FEProblemBase & _parallelacquisition_feproblem
Reference to FEProblemBase instance.
TheWarehouse & theWarehouse() const
const std::string name
Definition: Setup.h:21
Query query()

◆ getSurrogateModel() [1/2]

template<>
SurrogateModel& SurrogateModelInterface::getSurrogateModel ( const std::string &  name) const
inherited

Definition at line 46 of file SurrogateModelInterface.C.

47 {
48  return getSurrogateModelByName<SurrogateModel>(_smi_params.get<UserObjectName>(name));
49 }
const InputParameters & _smi_params
Parameters of the object with this interface.
std::vector< std::pair< R1, R2 > > get(const std::string &param1, const std::string &param2) const
const std::string name
Definition: Setup.h:21

◆ getSurrogateModel() [2/2]

template<typename T >
T & SurrogateModelInterface::getSurrogateModel ( const std::string &  name) const
inherited

Get a SurrogateModel/Trainer with a given name.

Parameters
nameThe name of the parameter key of the sampler to retrieve
Returns
The sampler with name associated with the parameter 'name'

Definition at line 81 of file SurrogateModelInterface.h.

Referenced by SurrogateTrainer::initialize().

82 {
83  return getSurrogateModelByName<T>(_smi_params.get<UserObjectName>(name));
84 }
const InputParameters & _smi_params
Parameters of the object with this interface.
std::vector< std::pair< R1, R2 > > get(const std::string &param1, const std::string &param2) const
const std::string name
Definition: Setup.h:21

◆ getSurrogateModelByName() [1/2]

template<>
SurrogateModel& SurrogateModelInterface::getSurrogateModelByName ( const UserObjectName &  name) const
inherited

Definition at line 31 of file SurrogateModelInterface.C.

32 {
33  std::vector<SurrogateModel *> models;
35  .query()
36  .condition<AttribName>(name)
37  .condition<AttribSystem>("SurrogateModel")
38  .queryInto(models);
39  if (models.empty())
40  mooseError("Unable to find a SurrogateModel object with the name '" + name + "'");
41  return *(models[0]);
42 }
void mooseError(Args &&... args)
FEProblemBase & _smi_feproblem
Reference to FEProblemBase instance.
TheWarehouse & theWarehouse() const
const std::string name
Definition: Setup.h:21
Query query()

◆ getSurrogateModelByName() [2/2]

template<typename T >
T & SurrogateModelInterface::getSurrogateModelByName ( const UserObjectName &  name) const
inherited

Get a sampler with a given name.

Parameters
nameThe name of the sampler to retrieve
Returns
The sampler with name 'name'

Definition at line 88 of file SurrogateModelInterface.h.

Referenced by CrossValidationScores::CrossValidationScores(), EvaluateSurrogate::EvaluateSurrogate(), and InverseMapping::initialSetup().

89 {
90  std::vector<T *> models;
92  .query()
93  .condition<AttribName>(name)
94  .condition<AttribSystem>("SurrogateModel")
95  .queryInto(models);
96  if (models.empty())
97  mooseError("Unable to find a SurrogateModel object of type " + std::string(typeid(T).name()) +
98  " with the name '" + name + "'");
99  return *(models[0]);
100 }
const double T
void mooseError(Args &&... args)
FEProblemBase & _smi_feproblem
Reference to FEProblemBase instance.
TheWarehouse & theWarehouse() const
const std::string name
Definition: Setup.h:21
Query query()

◆ getSurrogateTrainer() [1/2]

template<typename T >
T & SurrogateModelInterface::getSurrogateTrainer ( const std::string &  name) const
inherited

Definition at line 104 of file SurrogateModelInterface.h.

105 {
106  return getSurrogateTrainerByName<T>(_smi_params.get<UserObjectName>(name));
107 }
const InputParameters & _smi_params
Parameters of the object with this interface.
std::vector< std::pair< R1, R2 > > get(const std::string &param1, const std::string &param2) const
const std::string name
Definition: Setup.h:21

◆ getSurrogateTrainer() [2/2]

template<>
SurrogateTrainerBase& SurrogateModelInterface::getSurrogateTrainer ( const std::string &  name) const
inherited

Definition at line 60 of file SurrogateModelInterface.C.

61 {
62  return getSurrogateTrainerByName<SurrogateTrainerBase>(_smi_params.get<UserObjectName>(name));
63 }
const InputParameters & _smi_params
Parameters of the object with this interface.
std::vector< std::pair< R1, R2 > > get(const std::string &param1, const std::string &param2) const
const std::string name
Definition: Setup.h:21

◆ getSurrogateTrainerByName() [1/2]

template<>
SurrogateTrainerBase& SurrogateModelInterface::getSurrogateTrainerByName ( const UserObjectName &  name) const
inherited

Definition at line 53 of file SurrogateModelInterface.C.

54 {
56 }
T & getUserObject(const std::string &name, unsigned int tid=0) const
FEProblemBase & _smi_feproblem
Reference to FEProblemBase instance.
const std::string name
Definition: Setup.h:21
This is the base trainer class whose main functionality is the API for declaring model data...

◆ getSurrogateTrainerByName() [2/2]

template<typename T >
T & SurrogateModelInterface::getSurrogateTrainerByName ( const UserObjectName &  name) const
inherited

Definition at line 111 of file SurrogateModelInterface.h.

Referenced by SurrogateTrainerOutput::output().

112 {
113  SurrogateTrainerBase * base_ptr =
115  T * obj_ptr = dynamic_cast<T *>(base_ptr);
116  if (!obj_ptr)
117  mooseError("Failed to find a SurrogateTrainer object of type " + std::string(typeid(T).name()) +
118  " with the name '",
119  name,
120  "' for the desired type.");
121  return *obj_ptr;
122 }
T & getUserObject(const std::string &name, unsigned int tid=0) const
const double T
void mooseError(Args &&... args)
FEProblemBase & _smi_feproblem
Reference to FEProblemBase instance.
const std::string name
Definition: Setup.h:21
const THREAD_ID _smi_tid
Thread ID.
This is the base trainer class whose main functionality is the API for declaring model data...

◆ includeAdditionalInputs()

template<typename SamplerType >
void GenericActiveLearnerTempl< SamplerType >::includeAdditionalInputs ( )
protectedvirtual

Include additional inputs before evaluating the acquisition function.

Has trivial function in base, but can be modified in derived if necessary depending upon the objective of active learning (i.e., forward UQ, inverse UQ, optimization, etc.)

Reimplemented in BayesianActiveLearner.

Definition at line 260 of file GenericActiveLearner.h.

261 {
263 }
std::vector< std::vector< Real > > _inputs_test_modified
Storage for all the modified proposed samples to test the GP model.
const std::vector< std::vector< Real > > & _inputs_test
Storage for all the proposed samples to test the GP model.

◆ initialize()

template<typename SamplerType>
virtual void GenericActiveLearnerTempl< SamplerType >::initialize ( )
inlineoverridevirtual

Implements GeneralReporter.

Definition at line 41 of file GenericActiveLearner.h.

41 {}

◆ setupGeneric()

template<typename SamplerType >
void GenericActiveLearnerTempl< SamplerType >::setupGeneric ( )
protectedvirtual

Setup the generic variable for acquisition computation (depends on the objective: optimization, UQ, etc.)

Definition at line 253 of file GenericActiveLearner.h.

254 {
256 }
std::vector< Real > _gp_outputs
Storage for the GP re-training outputs.
std::vector< Real > _generic
A generic parameter to be passed to the acquisition function.

◆ setupGPData()

template<typename SamplerType >
void GenericActiveLearnerTempl< SamplerType >::setupGPData ( const std::vector< Real > &  data_out,
const DenseMatrix< Real > &  data_in 
)
protectedvirtual

Sets up the training data for the GP model.

Parameters
data_outThe data vector containing the outputs to train the GP
data_inThe data matrix containing the inputs to train the GP

Reimplemented in BayesianActiveLearner.

Definition at line 231 of file GenericActiveLearner.h.

233 {
234  for (unsigned int i = 0; i < data_out.size(); ++i)
235  {
236  for (unsigned int j = 0; j < _n_dim; ++j)
237  _inputs_required[i][j] = data_in(i, j);
238  _gp_inputs.push_back(_inputs_required[i]);
239  _gp_outputs.push_back(data_out[i]);
240  }
241 }
std::vector< std::vector< Real > > _gp_inputs
Storage for the GP re-training inputs.
unsigned int _n_dim
The input dimension for GP, equal to Sampler columns.
std::vector< Real > _gp_outputs
Storage for the GP re-training outputs.
static const std::complex< double > j(0, 1)
Complex number "j" (also known as "i")
std::vector< std::vector< Real > > & _inputs_required
Transmit the required inputs to the json file.

◆ validParams()

template<typename SamplerType >
InputParameters GenericActiveLearnerTempl< SamplerType >::validParams ( )
static

Definition at line 161 of file GenericActiveLearner.h.

Referenced by BayesianActiveLearner::validParams().

162 {
165  params.addClassDescription("A generic reporter to support parallel active learning: re-trains GP "
166  "and picks the next best batch.");
167  params.addRequiredParam<ReporterName>("output_value",
168  "Value of the model output from the SubApp.");
169  params.addParam<ReporterValueName>(
170  "outputs_required",
171  "outputs_required",
172  "Modified value of the model output from this reporter class.");
173  params.addRequiredParam<SamplerName>("sampler", "The sampler object.");
174  params.addRequiredParam<UserObjectName>("al_gp", "Active learning GP trainer.");
175  params.addRequiredParam<UserObjectName>("gp_evaluator", "Evaluator for the trained GP.");
176  params.addParam<ReporterValueName>(
177  "sorted_indices",
178  "sorted_indices",
179  "The sorted sample indices in order of importance to evaluate the subApp.");
180  params.addParam<ReporterValueName>(
181  "acquisition_function",
182  "acquisition_function",
183  "The values of the acquistion function in the current iteration.");
184  params.addParam<ReporterValueName>(
185  "convergence_value", "convergence_value", "Value to measure convergence of active learning.");
186  params.addParam<ReporterValueName>(
187  "inputs", "inputs", "Modified value of the model inputs from this reporter class.");
188  params.addRequiredParam<UserObjectName>("acquisition", "Name of the acquisition function.");
189  params.addParam<bool>(
190  "penalize_acquisition",
191  true,
192  "Set true to prevent clustering of the best batch inputs when operating in parallel.");
193  return params;
194 }
void addParam(const std::string &name, const std::initializer_list< typename T::value_type > &value, const std::string &doc_string)
static InputParameters validParams()
void addRequiredParam(const std::string &name, const std::string &doc_string)
void addClassDescription(const std::string &doc_string)

Member Data Documentation

◆ _acquisition_obj

template<typename SamplerType>
ParallelAcquisitionFunctionBase& GenericActiveLearnerTempl< SamplerType >::_acquisition_obj
protected

Storage for the parallel acquisition object to be utilized.

Definition at line 117 of file GenericActiveLearner.h.

◆ _acquisition_value

template<typename SamplerType>
std::vector<Real>& GenericActiveLearnerTempl< SamplerType >::_acquisition_value
protected

The acquistion function values in the current iteration.

Definition at line 120 of file GenericActiveLearner.h.

Referenced by GenericActiveLearnerTempl< BayesianActiveLearningSampler >::GenericActiveLearnerTempl().

◆ _al_gp

template<typename SamplerType>
const ActiveLearningGaussianProcess& GenericActiveLearnerTempl< SamplerType >::_al_gp
protected

The active learning GP trainer that permits re-training.

Definition at line 111 of file GenericActiveLearner.h.

◆ _al_sampler

template<typename SamplerType>
SamplerType& GenericActiveLearnerTempl< SamplerType >::_al_sampler
protected

The base sampler.

Definition at line 90 of file GenericActiveLearner.h.

◆ _check_step

template<typename SamplerType>
int GenericActiveLearnerTempl< SamplerType >::_check_step
protected

Ensure that the MCMC algorithm proceeds in a sequential fashion.

Definition at line 135 of file GenericActiveLearner.h.

◆ _convergence_value

template<typename SamplerType>
Real& GenericActiveLearnerTempl< SamplerType >::_convergence_value
protected

For monitoring convergence of active learning.

Definition at line 123 of file GenericActiveLearner.h.

◆ _eval_outputs_current

template<typename SamplerType>
std::vector<Real> GenericActiveLearnerTempl< SamplerType >::_eval_outputs_current
protected

The GP outputs from the current iteration before re-training (to evaluate convergence)

Definition at line 156 of file GenericActiveLearner.h.

Referenced by GenericActiveLearnerTempl< BayesianActiveLearningSampler >::GenericActiveLearnerTempl().

◆ _generic

template<typename SamplerType>
std::vector<Real> GenericActiveLearnerTempl< SamplerType >::_generic
protected

A generic parameter to be passed to the acquisition function.

Definition at line 153 of file GenericActiveLearner.h.

Referenced by GenericActiveLearnerTempl< BayesianActiveLearningSampler >::GenericActiveLearnerTempl().

◆ _gp_eval

template<typename SamplerType>
const SurrogateModel& GenericActiveLearnerTempl< SamplerType >::_gp_eval
protected

The GP evaluator object that permits re-evaluations.

Definition at line 114 of file GenericActiveLearner.h.

◆ _gp_inputs

template<typename SamplerType>
std::vector<std::vector<Real> > GenericActiveLearnerTempl< SamplerType >::_gp_inputs
protected

Storage for the GP re-training inputs.

Definition at line 138 of file GenericActiveLearner.h.

◆ _gp_outputs

template<typename SamplerType>
std::vector<Real> GenericActiveLearnerTempl< SamplerType >::_gp_outputs
protected

Storage for the GP re-training outputs.

Definition at line 141 of file GenericActiveLearner.h.

◆ _gp_outputs_test

template<typename SamplerType>
std::vector<Real> GenericActiveLearnerTempl< SamplerType >::_gp_outputs_test
protected

Outputs of GP model for the test samples.

Definition at line 144 of file GenericActiveLearner.h.

Referenced by GenericActiveLearnerTempl< BayesianActiveLearningSampler >::GenericActiveLearnerTempl().

◆ _gp_std_test

template<typename SamplerType>
std::vector<Real> GenericActiveLearnerTempl< SamplerType >::_gp_std_test
protected

Outputs of GP model standard deviation for the test samples.

Definition at line 147 of file GenericActiveLearner.h.

Referenced by GenericActiveLearnerTempl< BayesianActiveLearningSampler >::GenericActiveLearnerTempl().

◆ _inputs_required

template<typename SamplerType>
std::vector<std::vector<Real> >& GenericActiveLearnerTempl< SamplerType >::_inputs_required
protected

Transmit the required inputs to the json file.

Definition at line 129 of file GenericActiveLearner.h.

Referenced by GenericActiveLearnerTempl< BayesianActiveLearningSampler >::GenericActiveLearnerTempl().

◆ _inputs_test

template<typename SamplerType>
const std::vector<std::vector<Real> >& GenericActiveLearnerTempl< SamplerType >::_inputs_test
protected

Storage for all the proposed samples to test the GP model.

Definition at line 99 of file GenericActiveLearner.h.

Referenced by GenericActiveLearnerTempl< BayesianActiveLearningSampler >::GenericActiveLearnerTempl().

◆ _inputs_test_modified

template<typename SamplerType>
std::vector<std::vector<Real> > GenericActiveLearnerTempl< SamplerType >::_inputs_test_modified
protected

Storage for all the modified proposed samples to test the GP model.

Definition at line 126 of file GenericActiveLearner.h.

◆ _length_scales

template<typename SamplerType>
std::vector<Real> GenericActiveLearnerTempl< SamplerType >::_length_scales
protected

Storage for the length scales after the GP training.

Definition at line 150 of file GenericActiveLearner.h.

Referenced by GenericActiveLearnerTempl< BayesianActiveLearningSampler >::GenericActiveLearnerTempl().

◆ _n_dim

template<typename SamplerType>
unsigned int GenericActiveLearnerTempl< SamplerType >::_n_dim
protected

The input dimension for GP, equal to Sampler columns.

Definition at line 93 of file GenericActiveLearner.h.

Referenced by GenericActiveLearnerTempl< BayesianActiveLearningSampler >::GenericActiveLearnerTempl().

◆ _output_comm

template<typename SamplerType>
std::vector<Real>& GenericActiveLearnerTempl< SamplerType >::_output_comm
protected

Modified value of model output by this reporter class.

Definition at line 105 of file GenericActiveLearner.h.

◆ _output_value

template<typename SamplerType>
const std::vector<Real>& GenericActiveLearnerTempl< SamplerType >::_output_value
protected

Model output value from SubApp.

Definition at line 102 of file GenericActiveLearner.h.

◆ _penalize_acquisition

template<typename SamplerType>
const bool& GenericActiveLearnerTempl< SamplerType >::_penalize_acquisition
protected

Penalize acquisition to prevent clustering when operating in parallel.

Definition at line 132 of file GenericActiveLearner.h.

◆ _props

template<typename SamplerType>
dof_id_type GenericActiveLearnerTempl< SamplerType >::_props
protected

Storage for the number of parallel proposals.

Definition at line 96 of file GenericActiveLearner.h.

Referenced by GenericActiveLearnerTempl< BayesianActiveLearningSampler >::GenericActiveLearnerTempl().

◆ _sorted_indices

template<typename SamplerType>
std::vector<unsigned int>& GenericActiveLearnerTempl< SamplerType >::_sorted_indices
protected

The selected sample indices to evaluate the subApp.

Definition at line 108 of file GenericActiveLearner.h.

Referenced by GenericActiveLearnerTempl< BayesianActiveLearningSampler >::GenericActiveLearnerTempl().


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