rmf_links.h

Go to the documentation of this file.

/**
 *  \file npctransport/rmf_links.h
 *  \brief description
 *
 *  Copyright 2007-2021 IMP Inventors. All rights reserved.
 */

#ifndef IMPNPCTRANSPORT_RMF_LINKS_H
#define IMPNPCTRANSPORT_RMF_LINKS_H

#include "npctransport_config.h"
#include "SimulationData.h"
#include <IMP/rmf/atom_links.h>
#include <IMP/rmf/link_macros.h>
#include <boost/unordered_map.hpp>
#include <IMP/WeakPointer.h>

IMPNPCTRANSPORT_BEGIN_NAMESPACE

//! Load sites data from an RMF file
class IMPNPCTRANSPORTEXPORT HierarchyWithSitesLoadLink
    : public rmf::HierarchyLoadLink {
  WeakPointer<SimulationData> sd_;
  RMF::decorator::BallConstFactory bf_;
  // Transporting decorator related keys:
  RMF::IntKey is_last_entry_from_top_key_;
  RMF::IntKey n_entries_bottom_key_;
  RMF::IntKey n_entries_top_key_;
  // optimization data:
  RMF::IntKey coordinates_are_optimized_key_;
  // Slab decorator related keys:
  RMF::FloatKey pore_radius_key_;
  RMF::IntKey pore_radius_is_optimized_key_;
  // Spring decorator keys
  RMF::FloatKey rest_length_key_;

  boost::unordered_map<ParticleIndex, ParticleIndexes> particles_;

 protected:
  /** Called by rmf::HierarchyLoadLink::add_link_recursive()

      link the specified rmf node to particle cur. In addition, load
      the list of sites associated with the rmf node together with
      their coordinates (children of type RMF::Geometry decorated with
      RMF::BallConst) and associate them with the particle type of cur
      in the SimulationData of cur.
  */
  virtual void do_link_particle(Model *m, ParticleIndex root,
                                ParticleIndex cur,
                                RMF::NodeConstHandle node) IMP_OVERRIDE;

  virtual void do_setup_particle(Model* m, ParticleIndex root,
                                 ParticleIndex cur,
                                 RMF::NodeConstHandle node) IMP_OVERRIDE {
    IMP_NOT_IMPLEMENTED;
    IMP_UNUSED(m);     IMP_UNUSED(root);
    IMP_UNUSED(cur);    IMP_UNUSED(node);
  }

  /** load the values of a hierarchy. also
      loads dynamic Transporting decorator transport directionality
      information if needed (which is used in
      ParticleTransportStatisticsOptimizerState)
  */
  virtual void do_load_hierarchy(RMF::NodeConstHandle root_node,
                                 Model *m,
                                 ParticleIndex pi) IMP_OVERRIDE;

 public:
  HierarchyWithSitesLoadLink(RMF::FileConstHandle fh);
  static const char *get_name() {return "npctransport load";}
};

//! Save sites data to an RMF file
class IMPNPCTRANSPORTEXPORT HierarchyWithSitesSaveLink
    : public rmf::HierarchySaveLink {
  WeakPointer<SimulationData> sd_;
  RMF::decorator::BallFactory bf_;
  RMF::decorator::ColoredFactory cf_;
  // Transporting decorator related keys:
  RMF::IntKey is_last_entry_from_top_key_;
  RMF::IntKey n_entries_bottom_key_;
  RMF::IntKey n_entries_top_key_;
  // optimization data:
  RMF::IntKey coordinates_are_optimized_key_;
  // Slab decorator related keys:
  RMF::FloatKey pore_radius_key_;
  RMF::IntKey pore_radius_is_optimized_key_;
  // Spring decorator keys
  RMF::FloatKey rest_length_key_;

  boost::unordered_map<ParticleIndex, ParticleIndexes> particles_;

  // for testing without sd - site coords & radii per each particle type
  boost::unordered_map
    <core::ParticleType, algebra::Sphere3Ds>  test_sites_;

 private:
  // add sites associated with particle type t to node cur_node, using
  // sd_->get_sites() if sd_ is initialized, or using internal
  // test_sites_ table (used for testing only)
  void add_sites_to_node(RMF::NodeHandle cur_node, core::ParticleType t) const;

 protected:
  virtual void do_setup_node(Model *m, ParticleIndex root,
                           ParticleIndex cur,
                           RMF::NodeHandle cur_node) IMP_OVERRIDE;

  /** save the values of a hierarchy
  */
  virtual void do_save_hierarchy(Model *m, ParticleIndex root,
                                 RMF::NodeHandle root_node) IMP_OVERRIDE;

 public:
  HierarchyWithSitesSaveLink(RMF::FileHandle fh);
#ifndef SWIG
  // for testing
  void add_test_sites(core::ParticleType t,
                      algebra::Sphere3Ds sites) {
    test_sites_[t] = sites;
  }
#endif
  static const char *get_name() {return "npctransport save";}
};

//! for testing - adds the list of sites with specified radius, to be
//! associated with particle type t. The file handle fh relies on this list
//! only if it doesn't have particles with simulation data keys
IMPNPCTRANSPORTEXPORT void add_test_sites(RMF::FileHandle fh, core::ParticleType t,
                                     double display_radius, algebra::Vector3Ds sites);

//! for testing - adds the list of sites with specified display radius, to be
//! associated with particle type t. The file handle fh relies on this list
//! only if it doesn't have particles with simulation data keys
IMPNPCTRANSPORTEXPORT void add_test_sites(RMF::FileHandle fh,
                                          core::ParticleType t,
                                          algebra::Sphere3Ds sites);

// note that the corresponding define macro in the .cpp file implicitly
// uses the HierarchyWithSitesLoadLink and HierarchyWithSitesSaveLink
// classes
/**
   Functions for adding a hierarchy to an RMF file or linking an RMF
   file to an existing hierarchy, including support for particles with sites.

   These functions are practically identical to the add / link
   methods in modules/rmf/include/hierarchy_io.h, such as
   IMP::rmf::link_hierarchies(), except here NPC particle sites
   support is included.
 */
IMP_DECLARE_LINKERS(HierarchyWithSites,
                    hierarchy_with_sites, hierarchies_with_sites,
                    atom::Hierarchy, atom::Hierarchies,
                    (RMF::FileConstHandle fh, Model *m),
                    See IMP::rmf::link_hierarchies() for more details.
                    The only difference is the addition of particle sites
                    support.
                    );

IMPNPCTRANSPORT_END_NAMESPACE

#endif /* IMPNPCTRANSPORT_RMF_LINKS_H */