AdaptiveMonteCarloDecision.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html
#pragma once
 
#include "GeneralReporter.h"
#include "AdaptiveImportanceSampler.h"
#include "ParallelSubsetSimulation.h"
 
class AdaptiveMonteCarloDecision : public GeneralReporter
{
public:
  static InputParameters validParams();
  AdaptiveMonteCarloDecision(const InputParameters & parameters);
  virtual void initialize() override {}
  virtual void finalize() override {}
  virtual void execute() override;
 
protected:
  const std::vector<Real> & _output_value;
 
  std::vector<Real> & _output_required;
 
  std::vector<std::vector<Real>> & _inputs;
 
private:
  void reinitChain();
 
  Sampler & _sampler;
 
  const AdaptiveImportanceSampler * const _ais;
 
  const ParallelSubsetSimulation * const _pss;
 
  int _check_step;
 
  libMesh::Parallel::Communicator & _local_comm;
 
  std::vector<std::vector<Real>> _prev_val;
 
  std::vector<Real> _prev_val_out;
 
  std::vector<std::vector<Real>> _inputs_sto;
 
  std::vector<std::vector<Real>> _inputs_sorted;
 
  std::vector<Real> _outputs_sto;
 
  std::vector<Real> _output_sorted;
 
  Real _output_limit;
 
  const bool _gp_used;
 
  const int * const _gp_training_samples;
};