| [66cfc7] | 1 | /*
 | 
|---|
 | 2 |  * FunctionArgument.hpp
 | 
|---|
 | 3 |  *
 | 
|---|
 | 4 |  *  Created on: 02.10.2012
 | 
|---|
 | 5 |  *      Author: heber
 | 
|---|
 | 6 |  */
 | 
|---|
 | 7 | 
 | 
|---|
 | 8 | #ifndef FUNCTIONARGUMENT_HPP_
 | 
|---|
 | 9 | #define FUNCTIONARGUMENT_HPP_
 | 
|---|
 | 10 | 
 | 
|---|
 | 11 | // include config.h
 | 
|---|
 | 12 | #ifdef HAVE_CONFIG_H
 | 
|---|
 | 13 | #include <config.h>
 | 
|---|
 | 14 | #endif
 | 
|---|
 | 15 | 
 | 
|---|
 | 16 | #include <utility>
 | 
|---|
| [9e903b] | 17 | #include <iosfwd>
 | 
|---|
| [66cfc7] | 18 | 
 | 
|---|
 | 19 | /** This class encapsulates all information with respect to a single argument
 | 
|---|
 | 20 |  *  for a high-dimensional model function.
 | 
|---|
 | 21 |  *
 | 
|---|
 | 22 |  *  We restrict ourselves here to a function that dependent on a set of
 | 
|---|
 | 23 |  *  three-dimensional vectors, i.e. a set of positions in space. And for
 | 
|---|
 | 24 |  *  the moment to distances in between these sets.
 | 
|---|
 | 25 |  *
 | 
|---|
 | 26 |  */
 | 
|---|
 | 27 | struct argument_t
 | 
|---|
 | 28 | {
 | 
|---|
| [9e903b] | 29 |   //!> grant operator access to private parts
 | 
|---|
 | 30 |   friend std::ostream& operator<<(std::ostream &ost, const argument_t &arg);
 | 
|---|
 | 31 | 
 | 
|---|
 | 32 |   //!> typedef for the two indices of the argument
 | 
|---|
| [5b5724] | 33 |   typedef std::pair<size_t, size_t> indices_t;
 | 
|---|
| [691be4] | 34 |   //!> typedef for the underlying type of the particle
 | 
|---|
 | 35 |   typedef int ParticleType_t;
 | 
|---|
 | 36 |   //!> typedef for the two particle types of the argument
 | 
|---|
 | 37 |   typedef std::pair<ParticleType_t, ParticleType_t> types_t;
 | 
|---|
| [5b5724] | 38 | 
 | 
|---|
 | 39 |   /** Default constructor for class argument_t.
 | 
|---|
 | 40 |    *
 | 
|---|
 | 41 |    */
 | 
|---|
 | 42 |   argument_t() :
 | 
|---|
 | 43 |     indices( std::make_pair(0,1) ),
 | 
|---|
| [691be4] | 44 |     types( std::make_pair(0,0) ),
 | 
|---|
| [eb1efe] | 45 |     distance(0.),
 | 
|---|
 | 46 |     globalid(-1)
 | 
|---|
| [5b5724] | 47 |   {}
 | 
|---|
 | 48 | 
 | 
|---|
 | 49 |   /** Constructor for class argument_t.
 | 
|---|
 | 50 |    *
 | 
|---|
 | 51 |    * This constructors uses the index pair (0,1) as default.
 | 
|---|
 | 52 |    *
 | 
|---|
 | 53 |    * \param _distance distance argument
 | 
|---|
 | 54 |    */
 | 
|---|
 | 55 |   argument_t(const double &_distance) :
 | 
|---|
 | 56 |     indices( std::make_pair(0,1) ),
 | 
|---|
| [691be4] | 57 |     types( std::make_pair(0,0) ),
 | 
|---|
| [eb1efe] | 58 |     distance(_distance),
 | 
|---|
 | 59 |     globalid(-1)
 | 
|---|
| [5b5724] | 60 |   {}
 | 
|---|
 | 61 | 
 | 
|---|
 | 62 |   /** Constructor for class argument_t.
 | 
|---|
 | 63 |    *
 | 
|---|
 | 64 |    * \param _indices pair of indices associated with the \a _distance
 | 
|---|
 | 65 |    * \param _distance distance argument
 | 
|---|
 | 66 |    */
 | 
|---|
 | 67 |   argument_t(const indices_t &_indices, const double &_distance) :
 | 
|---|
 | 68 |     indices( _indices ),
 | 
|---|
| [691be4] | 69 |     types( std::make_pair(0,0) ),
 | 
|---|
 | 70 |     distance(_distance),
 | 
|---|
 | 71 |     globalid(-1)
 | 
|---|
 | 72 |   {}
 | 
|---|
 | 73 | 
 | 
|---|
 | 74 |   /** Constructor for class argument_t.
 | 
|---|
 | 75 |    *
 | 
|---|
 | 76 |    * \param _indices pair of indices associated with the \a _distance
 | 
|---|
 | 77 |    * \param _types pair of particle type
 | 
|---|
 | 78 |    * \param _distance distance argument
 | 
|---|
 | 79 |    */
 | 
|---|
 | 80 |   argument_t(const indices_t &_indices, const types_t &_types, const double &_distance) :
 | 
|---|
 | 81 |     indices( _indices ),
 | 
|---|
 | 82 |     types( _types ),
 | 
|---|
| [eb1efe] | 83 |     distance(_distance),
 | 
|---|
 | 84 |     globalid(-1)
 | 
|---|
| [5b5724] | 85 |   {}
 | 
|---|
 | 86 | 
 | 
|---|
| [355af8] | 87 |   /** Comparator with respect to the distance.
 | 
|---|
 | 88 |    *
 | 
|---|
 | 89 |    * \note We'll have this as static function to allow usage in e.g. STL's sort.
 | 
|---|
 | 90 |    *
 | 
|---|
 | 91 |    * \param one first argument
 | 
|---|
 | 92 |    * \param other other argument to compare to \a one to
 | 
|---|
 | 93 |    * \return true - first distance is less
 | 
|---|
 | 94 |    */
 | 
|---|
 | 95 |   static bool DistanceComparator(const argument_t &one, const argument_t &other)
 | 
|---|
 | 96 |   {
 | 
|---|
 | 97 |     return one.distance < other.distance;
 | 
|---|
 | 98 |   }
 | 
|---|
 | 99 | 
 | 
|---|
| [64bdfd] | 100 |   /** Comparator with respect to the pair of types.
 | 
|---|
 | 101 |    *
 | 
|---|
 | 102 |    * \note We'll have this as static function to allow usage in e.g. STL's sort.
 | 
|---|
 | 103 |    *
 | 
|---|
 | 104 |    * \param one first argument
 | 
|---|
 | 105 |    * \param other other argument to compare to \a one to
 | 
|---|
 | 106 |    * \return true - first type is less or if equal, second type is less, else
 | 
|---|
 | 107 |    */
 | 
|---|
 | 108 |   bool static TypeComparator(const argument_t &one, const argument_t &other)
 | 
|---|
 | 109 |   {
 | 
|---|
 | 110 |     if (one.types.first < other.types.first)
 | 
|---|
 | 111 |       return true;
 | 
|---|
 | 112 |     else if (one.types.first > other.types.first)
 | 
|---|
 | 113 |       return false;
 | 
|---|
 | 114 |     else
 | 
|---|
 | 115 |       return one.types.second < other.types.second;
 | 
|---|
 | 116 |   }
 | 
|---|
 | 117 | 
 | 
|---|
| [355af8] | 118 |   /** Comparator with respect to the pair of indices.
 | 
|---|
 | 119 |    *
 | 
|---|
 | 120 |    * \note We'll have this as static function to allow usage in e.g. STL's sort.
 | 
|---|
 | 121 |    *
 | 
|---|
 | 122 |    * \param one first argument
 | 
|---|
 | 123 |    * \param other other argument to compare to \a one to
 | 
|---|
 | 124 |    * \return true - first index is less or if equal, second index is less, else
 | 
|---|
 | 125 |    */
 | 
|---|
 | 126 |   bool static IndexComparator(const argument_t &one, const argument_t &other)
 | 
|---|
 | 127 |   {
 | 
|---|
 | 128 |     if (one.indices.first < other.indices.first)
 | 
|---|
 | 129 |       return true;
 | 
|---|
 | 130 |     else if (one.indices.first > other.indices.first)
 | 
|---|
 | 131 |       return false;
 | 
|---|
 | 132 |     else
 | 
|---|
 | 133 |       return one.indices.second < other.indices.second;
 | 
|---|
 | 134 |   }
 | 
|---|
 | 135 | 
 | 
|---|
| [64bdfd] | 136 |   /** Less comparator for FunctionArgument.
 | 
|---|
 | 137 |    *
 | 
|---|
 | 138 |    * @param other other argument to compare to
 | 
|---|
 | 139 |    * @return true - this argument is less than \a other, false - else
 | 
|---|
 | 140 |    */
 | 
|---|
 | 141 |   bool operator<(const argument_t &other) const
 | 
|---|
 | 142 |   {
 | 
|---|
 | 143 |     if (types.first < other.types.first)
 | 
|---|
 | 144 |       return true;
 | 
|---|
 | 145 |     else if (types.first > other.types.first)
 | 
|---|
 | 146 |       return false;
 | 
|---|
 | 147 |     else if (types.second < other.types.second)
 | 
|---|
 | 148 |       return true;
 | 
|---|
 | 149 |     else if (types.second > other.types.second)
 | 
|---|
 | 150 |       return false;
 | 
|---|
 | 151 |     else {
 | 
|---|
 | 152 |       if (indices.first < other.indices.first)
 | 
|---|
 | 153 |         return true;
 | 
|---|
 | 154 |       else if (indices.first > other.indices.first)
 | 
|---|
 | 155 |         return false;
 | 
|---|
 | 156 |       else
 | 
|---|
 | 157 |         return indices.second < other.indices.second;
 | 
|---|
 | 158 |     }
 | 
|---|
 | 159 |   }
 | 
|---|
 | 160 | 
 | 
|---|
 | 161 |   /** Equality operator for FunctionArgument.
 | 
|---|
 | 162 |    *
 | 
|---|
 | 163 |    * \note This compares only types and indices.
 | 
|---|
 | 164 |    *
 | 
|---|
 | 165 |    * @param other other argument to compare to
 | 
|---|
 | 166 |    * @return true - this argument is equal to \a other, false - else
 | 
|---|
 | 167 |    */
 | 
|---|
 | 168 |   bool operator==(const argument_t &other) const
 | 
|---|
 | 169 |   {
 | 
|---|
 | 170 |     if (types.first != other.types.first)
 | 
|---|
 | 171 |       return false;
 | 
|---|
 | 172 |     else if (types.second != other.types.second)
 | 
|---|
 | 173 |       return false;
 | 
|---|
 | 174 |     if (indices.first != other.indices.first)
 | 
|---|
 | 175 |       return false;
 | 
|---|
 | 176 |     else if (indices.second != other.indices.second)
 | 
|---|
 | 177 |       return false;
 | 
|---|
 | 178 |     else
 | 
|---|
 | 179 |       return true;
 | 
|---|
 | 180 |   }
 | 
|---|
 | 181 | 
 | 
|---|
| [66cfc7] | 182 |   //!> indices between which the distance is given
 | 
|---|
| [5b5724] | 183 |   indices_t indices;
 | 
|---|
| [691be4] | 184 |   //!> indices between which the distance is given
 | 
|---|
 | 185 |   types_t types;
 | 
|---|
| [66cfc7] | 186 |   //!> distance
 | 
|---|
 | 187 |   double distance;
 | 
|---|
| [eb1efe] | 188 |   //!> global id refers to some global index, e.g. the configuration id in training set
 | 
|---|
 | 189 |   size_t globalid;
 | 
|---|
| [66cfc7] | 190 | };
 | 
|---|
 | 191 | 
 | 
|---|
| [9e903b] | 192 | /** Print given \a arg to stream \a ost.
 | 
|---|
 | 193 |  *
 | 
|---|
 | 194 |  * \param ost output stream to print to
 | 
|---|
 | 195 |  * \param arg argument to print
 | 
|---|
 | 196 |  * \return output stream for concatenation
 | 
|---|
 | 197 |  */
 | 
|---|
 | 198 | std::ostream& operator<<(std::ostream &ost, const argument_t &arg);
 | 
|---|
 | 199 | 
 | 
|---|
 | 200 | 
 | 
|---|
| [66cfc7] | 201 | #endif /* FUNCTIONARGUMENT_HPP_ */
 | 
|---|