| [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_ */ | 
|---|