ParallelSubsetSimulation.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 "ReporterInterface.h"

class ParallelSubsetSimulation : public Sampler
{
public:
  static InputParameters validParams();

  ParallelSubsetSimulation(const InputParameters & parameters);

  const unsigned int & getNumSamplesSub() const;

  const bool & getUseAbsoluteValue() const;

  const Real & getSubsetProbability() const;

  virtual bool isAdaptiveSamplingCompleted() const override { return _is_sampling_completed; }

protected:
  virtual void sampleSetUp(const Sampler::SampleMode mode) override;
  virtual Real computeSample(dof_id_type row_index, dof_id_type col_index) override;

  const unsigned int & _num_samplessub;

  const unsigned int & _num_subsets;

  const bool & _use_absolute_value;

  const Real & _subset_probability;

  const unsigned int & _num_random_seeds;

  const std::vector<Real> & _outputs;

  const std::vector<std::vector<Real>> & _inputs;

  const int & _step;

  const unsigned int _count_max;

  int _check_step;

  unsigned int _subset;

  std::vector<Distribution const *> _distributions;

  bool _is_sampling_completed;

private:
  std::vector<std::vector<Real>> _inputs_sto;

  std::vector<Real> _outputs_sto;

  std::vector<std::vector<Real>> _inputs_sorted;

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