IMP  2.4.0
The Integrative Modeling Platform
vonMisesKappaConjugateRestraint.h
Go to the documentation of this file.
1 /**
2  * \file IMP/isd/vonMisesKappaConjugateRestraint.h
3  * \brief Conjugate prior for \f$\kappa\f$ in the von Mises distribution.
4  *
5  * Copyright 2007-2015 IMP Inventors. All rights reserved.
6  *
7  */
8 
9 #ifndef IMPISD_VON_MISES_KAPPA_CONJUGATE_RESTRAINT_H
10 #define IMPISD_VON_MISES_KAPPA_CONJUGATE_RESTRAINT_H
11 
12 #include <IMP/isd/isd_config.h>
13 #include <IMP/SingletonScore.h>
14 #include <IMP/kernel/Restraint.h>
15 
16 IMPISD_BEGIN_NAMESPACE
17 
18 //! Conjugate prior for the concentration parameter of a von Mises distribution.
19 /** \f[ p(\kappa|c,R_0) = \frac{\exp(\kappa R_0)}{I_0(\kappa)^c} \f]
20  \f$0 < R_0 \le c\f$
21  where \f$kappa\f$ is the concentration parameter.
22  Default values: \f$R_0=1\f$ and \f$c=10\f$
23  */
26  double old_kappaval;
27  bool bessel_init_;
28  double I0_, I1_;
29  double c_, R0_;
30 
31  public:
32  //! Create the restraint.
34  double c = 10.0, double R0 = 0.0);
35 
36  /** This macro declares the basic needed methods: evaluate and show
37  */
38  virtual double unprotected_evaluate(IMP::kernel::DerivativeAccumulator *accum)
39  const IMP_OVERRIDE;
42 
43  virtual double get_probability() const;
44 
45  private:
46  void update_bessel(double kappaval); // update memoized bessel value
47  double get_kappa() const;
48 };
49 
50 IMPISD_END_NAMESPACE
51 
52 #endif /* IMPISD_VON_MISES_KAPPA_CONJUGATE_RESTRAINT_H */
Class for adding derivatives from restraints to the model.
#define IMP_OBJECT_METHODS(Name)
Define the basic things needed by any Object.
Definition: object_macros.h:25
A smart pointer to a reference counted object.
Definition: Pointer.h:87
Conjugate prior for the concentration parameter of a von Mises distribution.
Abstract base class for all restraints.
A restraint is a term in an IMP ScoringFunction.
Class to handle individual model particles.
Import IMP/kernel/SingletonScore.h in the namespace.
virtual ModelObjectsTemp do_get_inputs() const =0
#define IMP_OVERRIDE
Cause a compile error if this method does not override a parent method.
Class for storing model, its restraints, constraints, and particles.
Definition: kernel/Model.h:73