8 #ifndef IMPCORE_MONTE_CARLO_H
9 #define IMPCORE_MONTE_CARLO_H
11 #include <IMP/core/core_config.h>
16 #include <IMP/kernel/internal/container_helpers.h>
20 #include <boost/random/uniform_real.hpp>
22 IMPCORE_BEGIN_NAMESPACE
25 #define IMP_CORE_HAS_MONTE_CARLO_MOVER 1
57 void set_return_best(
bool tf) { return_best_ = tf; }
65 void set_kt(
Float t) {
79 <<
" requires return best being on.");
88 IMPCORE_DEPRECATED_FUNCTION_DECL(2.2)
89 unsigned
int get_number_of_forward_steps()
const {
90 IMPCORE_DEPRECATED_FUNCTION_DEF(2.2,
91 "Use get_number_of_accepted_steps() instead.");
92 return get_number_of_accepted_steps();
97 return stat_downward_steps_taken_;
101 return stat_upward_steps_taken_;
105 return stat_downward_steps_taken_ + stat_upward_steps_taken_ +
110 return stat_downward_steps_taken_ + stat_upward_steps_taken_;
112 void reset_statistics() {
113 stat_downward_steps_taken_ = 0;
114 stat_upward_steps_taken_ = 0;
115 stat_num_failures_ = 0;
134 double get_maximum_difference()
const {
return max_difference_; }
142 IMP_LIST_ACTION(
public, Mover, Movers, mover, movers, MonteCarloMover *,
154 void set_incremental_scoring_function(IncrementalScoringFunction *isf);
155 bool get_use_incremental_scoring_function()
const {
return isf_; }
156 IncrementalScoringFunction *get_incremental_scoring_function()
const {
163 kernel::ParticleIndexes get_movable_particles()
const;
168 bool do_accept_or_reject_move(
double score,
double last,
169 double proposal_ratio);
170 bool do_accept_or_reject_move(
double score,
double proposal_ratio) {
171 return do_accept_or_reject_move(score, get_last_accepted_energy(),
175 MonteCarloMoverResult do_move();
177 virtual void do_step();
194 isf_->set_moved_particles(moved);
196 if (get_maximum_difference() <
NO_MAX) {
198 false, last_energy_ + max_difference_);
208 double max_difference_;
209 unsigned int stat_downward_steps_taken_;
210 unsigned int stat_upward_steps_taken_;
211 unsigned int stat_num_failures_;
215 ::boost::uniform_real<> rand_;
223 unsigned int num_local_;
228 unsigned int get_number_of_steps()
const {
return num_local_; }
230 Optimizer *get_local_optimizer()
const {
return opt_; }
233 virtual void do_step() IMP_OVERRIDE;
251 virtual void do_step() IMP_OVERRIDE;
255 IMPCORE_END_NAMESPACE
const double NO_MAX
Use this value when you want to turn off maximum for restraint evaluation.
unsigned int get_number_of_upward_steps() const
Return how many times the optimizer has stepped to higher score.
Simple Monte Carlo optimizer.
double get_score_threshold() const
Get the score threshold.
double evaluate(bool derivatives)
Evaluate and return the score.
A smart pointer to a ref-counted Object that is a class memeber.
void set_maximum_difference(double d)
double get_last_accepted_energy() const
#define IMP_UNUSED(variable)
A smart pointer to a reference counted object.
IMP::base::Vector< IMP::base::Pointer< MonteCarloMover > > MonteCarloMovers
Import IMP/kernel/Optimizer.h in the namespace.
void set_score_threshold(double s)
Set the score threshold.
Import IMP/kernel/container_macros.h in the namespace.
#define IMP_USAGE_CHECK(expr, message)
A runtime test for incorrect usage of a class or method.
virtual void do_step()
a class that inherits from this should override this method
This variant of Monte Carlo uses basis hopping.
unsigned int get_number_of_downward_steps() const
Return how many times the optimizer has stepped to lower score.
#define IMP_INTERNAL_CHECK(expr, message)
An assertion to check for internal errors in IMP. An IMP::ErrorException will be thrown.
The base class for movers for MC optimization.
unsigned int get_number_of_accepted_steps() const
Get number of accepted moves.
This variant of Monte Carlo that relaxes after each move.
#define IMP_OBJECT_METHODS(Name)
Define the basic things needed by any Object.
double get_best_accepted_energy() const
Base class for all optimizers.
virtual double do_evaluate(const kernel::ParticleIndexes &moved) const
Get the current energy.
double Float
Basic floating-point value (could be float, double...)
ScoringFunction * get_scoring_function() const
Return the scoring function that is being used.
virtual double do_optimize(unsigned int ns)=0
override this function to do actual optimization
virtual void do_step()
a class that inherits from this should override this method
Import IMP/kernel/Configuration.h in the namespace.
Functions to generate vectors.
unsigned int get_number_of_proposed_steps() const
Get number of proposed moves.
Class for storing model, its restraints, constraints, and particles.