9 #ifndef IMPISD_GAUSSIAN_EM_RESTRAINT_H
10 #define IMPISD_GAUSSIAN_EM_RESTRAINT_H
12 #include "isd_config.h"
23 #include <IMP/algebra/eigen3/Eigen/Dense>
24 #include <boost/unordered_map.hpp>
26 IMPISD_BEGIN_NAMESPACE
28 #if !defined(SWIG) && !defined(IMP_DOXYGEN)
29 struct KahanAccumulation{
33 sum(0.0),correction(0.0)
36 struct KahanVectorAccumulation{
37 IMP_Eigen::Vector3d sum;
38 IMP_Eigen::Vector3d correction;
39 KahanVectorAccumulation():
40 sum(IMP_Eigen::Vector3d(0,0,0)),
41 correction(IMP_Eigen::Vector3d(0,0,0))
44 inline KahanAccumulation KahanSum(KahanAccumulation accumulation,
double value){
45 KahanAccumulation result;
46 double y = value - accumulation.correction;
47 double t = accumulation.sum + y;
48 result.correction = (t - accumulation.sum) - y;
52 inline KahanVectorAccumulation
53 KahanVectorSum(KahanVectorAccumulation accumulation, IMP_Eigen::Vector3d value){
54 KahanVectorAccumulation result;
55 IMP_Eigen::Vector3d y = value - accumulation.correction;
56 IMP_Eigen::Vector3d t = accumulation.sum + y;
57 result.correction = (t - accumulation.sum) - y;
99 bool update_model=
true,
bool backbone_slope=
false,
100 std::string name=
"GaussianEMRestraint%1%");
104 return exp(-unprotected_evaluate(NULL));
111 void compute_initial_scores();
122 void show(std::ostream &out)
const { out <<
"GEM restraint"; }
132 Float normalization_;
134 Float self_mm_score_;
Declare an efficient stl-compatible map.
A decorator for particles with mass.
Store a list of ParticleIndexes.
Return all pairs from a SingletonContainer.
#define IMP_OBJECT_METHODS(Name)
Define the basic things needed by any Object.
Float get_slope()
Get restraint slope.
void set_slope(Float s)
Set restraint slope.
Macros to define containers of objects.
Class for storing model, its restraints, constraints, and particles.
double get_probability() const
Returns exp(score)
Creates a restraint between two Gaussian Mixture Models, "model" and "density".
Return all pairs from a SingletonContainer.
Decorator to hold Gaussian3D.
std::ostream & show(Hierarchy h, std::ostream &out=std::cout)
Print the hierarchy using a given decorator to display each node.
double Float
Basic floating-point value (could be float, double...)
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 adding derivatives from restraints to the model.
A restraint is a term in an IMP ScoringFunction.