/usr/include/openturns/RandomWalkMetropolisHastings.hxx is in libopenturns-dev 1.9-5.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 | // -*- C++ -*-
/**
* @brief RandomWalkMetropolisHastings updates the chain
*
* Copyright 2005-2017 Airbus-EDF-IMACS-Phimeca
*
* This library is free software: you can redistribute it and/or modify
* it under the terms of the GNU Lesser General Public License as published by
* the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* This library is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU Lesser General Public License for more details.
*
* You should have received a copy of the GNU Lesser General Public
* along with this library. If not, see <http://www.gnu.org/licenses/>.
*
*/
#ifndef OPENTURNS_RANDOMWALKMETROPOLISHASTINGS_HXX
#define OPENTURNS_RANDOMWALKMETROPOLISHASTINGS_HXX
#include "openturns/OTprivate.hxx"
#include "openturns/MCMC.hxx"
#include "openturns/CalibrationStrategy.hxx"
#include "openturns/Interval.hxx"
#include "openturns/ResourceMap.hxx"
BEGIN_NAMESPACE_OPENTURNS
/**
* @class RandomWalkMetropolisHastings
*
* @brief The class implements the concept of comparison operator defined in Sampler.
*
* @see Sampler
*/
class OT_API RandomWalkMetropolisHastings
: public MCMC
{
CLASSNAME;
public:
typedef Collection<Distribution> DistributionCollection;
typedef PersistentCollection<Distribution> DistributionPersistentCollection;
typedef Collection<CalibrationStrategy> CalibrationStrategyCollection;
typedef PersistentCollection<CalibrationStrategy> CalibrationStrategyPersistentCollection;
/** Default constructor */
RandomWalkMetropolisHastings();
/** Constructor with parameters*/
RandomWalkMetropolisHastings(const Distribution & prior,
const Distribution & conditional,
const Sample & observations,
const Point & initialState,
const DistributionCollection & proposal);
/** Constructor with parameters*/
RandomWalkMetropolisHastings(const Distribution & prior,
const Distribution & conditional,
const Function & model,
const Sample & parameters,
const Sample & observations,
const Point & initialState,
const DistributionCollection & proposal);
/** String converter */
virtual String __repr__() const;
/* Here is the interface that all derived class must implement */
/** Virtual constructor */
virtual RandomWalkMetropolisHastings * clone() const;
/** @copydoc Sampler::getRealization() const */
virtual Point getRealization() const;
/** Calibration strategy accessor */
void setCalibrationStrategy(const CalibrationStrategy & calibrationStrategy);
void setCalibrationStrategyPerComponent(const CalibrationStrategyCollection & calibrationStrategy);
CalibrationStrategyCollection getCalibrationStrategyPerComponent() const;
/** Proposal accessor */
void setProposal(const DistributionCollection & proposal);
DistributionCollection getProposal() const;
/** Acceptance rate accessor*/
Point getAcceptanceRate() const;
/** Method save() stores the object through the StorageManager */
void save(Advocate & adv) const;
/** Method load() reloads the object from the StorageManager */
void load(Advocate & adv);
private:
/// proposal densities of the markov chain
DistributionPersistentCollection proposal_;
/// how to calibrate the each chain component
CalibrationStrategyPersistentCollection calibrationStrategy_;
/// number of samples
mutable UnsignedInteger samplesNumber_;
/// number of samples accepted
mutable Indices acceptedNumber_;
mutable Scalar currentLogLikelihood_;
}; /* class RandomWalkMetropolisHastings */
END_NAMESPACE_OPENTURNS
#endif /* OPENTURNS_RANDOMWALKMETROPOLISHASTINGS_HXX */
|