PMCMCBase.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 "Sampler.h"
#include "TransientInterface.h"
#include "Distribution.h"

class PMCMCBase : public Sampler, public TransientInterface
{
public:
  static InputParameters validParams();

  PMCMCBase(const InputParameters & parameters);

  dof_id_type getNumberOfConfigValues() const { return _confg_values[0].size(); }

  dof_id_type getNumberOfConfigParams() const { return _confg_values.size(); }

  dof_id_type getNumParallelProposals() const { return _num_parallel_proposals; }

  const std::vector<Real> & getRandomNumbers() const;

  const std::vector<Real> & getVarSamples() const;

  const std::vector<const Distribution *> getPriors() const;

  const Distribution * getVarPrior() const;

  virtual int decisionStep() const { return 1; }

protected:
  virtual void proposeSamples(const unsigned int seed_value);

  // See Sampler.h for description
  virtual void sampleSetUp(const Sampler::SampleMode mode) override;

  // See Sampler.h for description
  virtual Real computeSample(dof_id_type row_index, dof_id_type col_index) override;

  void randomIndex(const unsigned int & upper_bound,
                   const unsigned int & exclude,
                   const unsigned int & seed,
                   unsigned int & req_index);

  void randomIndexPair(const unsigned int & upper_bound,
                       const unsigned int & exclude,
                       const unsigned int & seed,
                       unsigned int & req_index1,
                       unsigned int & req_index2);

  const unsigned int _num_parallel_proposals;

  std::vector<const Distribution *> _priors;

  const Distribution * _var_prior;

  const std::vector<Real> * _lower_bound;

  const std::vector<Real> * _upper_bound;

  int _check_step;

  const std::vector<Real> & _initial_values;

  std::vector<std::vector<Real>> _new_samples;

  std::vector<Real> _new_var_samples;

  std::vector<Real> _rnd_vec;

private:
  void combineWithExperimentalConfig();

  const unsigned int _num_random_seeds;

  std::vector<std::vector<Real>> _confg_values;

  std::vector<std::vector<Real>> _new_samples_confg;
};