IndependentMHDecision.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 "PMCMCDecision.h"
#include "IndependentGaussianMH.h"
#include "AdaptiveMonteCarloUtils.h"
 
class IndependentMHDecision : public PMCMCDecision
{
public:
  static InputParameters validParams();
 
  IndependentMHDecision(const InputParameters & parameters);
 
protected:
  virtual void computeEvidence(std::vector<Real> & evidence,
                               const DenseMatrix<Real> & input_matrix) override;
 
  virtual void computeTransitionVector(std::vector<Real> & tv,
                                       const std::vector<Real> & evidence) override;
 
  virtual void nextSamples(std::vector<Real> & req_inputs,
                           DenseMatrix<Real> & input_matrix,
                           const std::vector<Real> & tv,
                           const unsigned int & parallel_index) override;
 
  virtual void nextSeeds() override;
 
private:
  const IndependentGaussianMH * const _igmh;
 
  std::vector<Real> & _seed_input;
 
  std::vector<Real> _seed_outputs;
 
  std::vector<Real> _tpm_modified;
 
  std::vector<Real> _outputs_sto;
};