| 1 | /*
 | 
|---|
| 2 |  * Project: MoleCuilder
 | 
|---|
| 3 |  * Description: creates and alters molecular systems
 | 
|---|
| 4 |  * Copyright (C)  2010-2012 University of Bonn. All rights reserved.
 | 
|---|
| 5 |  * Please see the LICENSE file or "Copyright notice" in builder.cpp for details.
 | 
|---|
| 6 |  */
 | 
|---|
| 7 | 
 | 
|---|
| 8 | /*
 | 
|---|
| 9 |  * \file RandomNumberDistribution_Encapsulation_impl.hpp
 | 
|---|
| 10 |  *
 | 
|---|
| 11 |  * Herein we place specializations of specific functions for specific distributions.
 | 
|---|
| 12 |  *
 | 
|---|
| 13 |  *  Created on: Jan 5, 2011
 | 
|---|
| 14 |  *      Author: heber
 | 
|---|
| 15 |  */
 | 
|---|
| 16 | 
 | 
|---|
| 17 | // include config.h
 | 
|---|
| 18 | #ifdef HAVE_CONFIG_H
 | 
|---|
| 19 | #include <config.h>
 | 
|---|
| 20 | #endif
 | 
|---|
| 21 | 
 | 
|---|
| 22 | #include "CodePatterns/MemDebug.hpp"
 | 
|---|
| 23 | 
 | 
|---|
| 24 | #include <limits>
 | 
|---|
| 25 | 
 | 
|---|
| 26 | #include "RandomNumberDistribution_Encapsulation.hpp"
 | 
|---|
| 27 | 
 | 
|---|
| 28 | /* =========  manipulatedclone() ================ */
 | 
|---|
| 29 | 
 | 
|---|
| 30 | template <>
 | 
|---|
| 31 | RandomNumberDistribution_Encapsulation< boost::uniform_smallint<> >::
 | 
|---|
| 32 | RandomNumberDistribution_Encapsulation (const RandomNumberDistribution_Parameters&_params) :
 | 
|---|
| 33 |   distribution_type(_params.min, _params.max)
 | 
|---|
| 34 | {}
 | 
|---|
| 35 | 
 | 
|---|
| 36 | template <>
 | 
|---|
| 37 | RandomNumberDistribution_Encapsulation< boost::uniform_int<> >::
 | 
|---|
| 38 | RandomNumberDistribution_Encapsulation(const RandomNumberDistribution_Parameters&_params) :
 | 
|---|
| 39 |   distribution_type(_params.min, _params.max)
 | 
|---|
| 40 | {}
 | 
|---|
| 41 | 
 | 
|---|
| 42 | template <>
 | 
|---|
| 43 | RandomNumberDistribution_Encapsulation< boost::uniform_real<> >::
 | 
|---|
| 44 | RandomNumberDistribution_Encapsulation(const RandomNumberDistribution_Parameters&_params) :
 | 
|---|
| 45 |   distribution_type(_params.min, _params.max)
 | 
|---|
| 46 | {}
 | 
|---|
| 47 | 
 | 
|---|
| 48 | template <>
 | 
|---|
| 49 | RandomNumberDistribution_Encapsulation< boost::bernoulli_distribution<> >::
 | 
|---|
| 50 | RandomNumberDistribution_Encapsulation(const RandomNumberDistribution_Parameters&_params) :
 | 
|---|
| 51 |   distribution_type(_params.p)
 | 
|---|
| 52 | {}
 | 
|---|
| 53 | 
 | 
|---|
| 54 | template <>
 | 
|---|
| 55 | RandomNumberDistribution_Encapsulation< boost::binomial_distribution<> >::
 | 
|---|
| 56 | RandomNumberDistribution_Encapsulation(const RandomNumberDistribution_Parameters&_params) :
 | 
|---|
| 57 |   distribution_type((int)_params.t, _params.p)
 | 
|---|
| 58 | {}
 | 
|---|
| 59 | 
 | 
|---|
| 60 | template <>
 | 
|---|
| 61 | RandomNumberDistribution_Encapsulation< boost::cauchy_distribution<> >::
 | 
|---|
| 62 | RandomNumberDistribution_Encapsulation(const RandomNumberDistribution_Parameters&_params) :
 | 
|---|
| 63 |   distribution_type(_params.median, _params.sigma)
 | 
|---|
| 64 | {}
 | 
|---|
| 65 | 
 | 
|---|
| 66 | template <>
 | 
|---|
| 67 | RandomNumberDistribution_Encapsulation< boost::gamma_distribution<> >::
 | 
|---|
| 68 | RandomNumberDistribution_Encapsulation(const RandomNumberDistribution_Parameters&_params) :
 | 
|---|
| 69 |   distribution_type(_params.alpha)
 | 
|---|
| 70 | {}
 | 
|---|
| 71 | 
 | 
|---|
| 72 | template <>
 | 
|---|
| 73 | RandomNumberDistribution_Encapsulation< boost::poisson_distribution<> >::
 | 
|---|
| 74 | RandomNumberDistribution_Encapsulation(const RandomNumberDistribution_Parameters&_params) :
 | 
|---|
| 75 |   distribution_type(_params.mean)
 | 
|---|
| 76 | {}
 | 
|---|
| 77 | 
 | 
|---|
| 78 | template <>
 | 
|---|
| 79 | RandomNumberDistribution_Encapsulation< boost::geometric_distribution<> >::
 | 
|---|
| 80 | RandomNumberDistribution_Encapsulation(const RandomNumberDistribution_Parameters&_params) :
 | 
|---|
| 81 |   distribution_type(_params.p)
 | 
|---|
| 82 | {}
 | 
|---|
| 83 | 
 | 
|---|
| 84 | template <>
 | 
|---|
| 85 | RandomNumberDistribution_Encapsulation< boost::triangle_distribution<> >::
 | 
|---|
| 86 | RandomNumberDistribution_Encapsulation(const RandomNumberDistribution_Parameters&_params) :
 | 
|---|
| 87 |   distribution_type(_params.a, _params.b, _params.c)
 | 
|---|
| 88 | {}
 | 
|---|
| 89 | 
 | 
|---|
| 90 | template <>
 | 
|---|
| 91 | RandomNumberDistribution_Encapsulation< boost::exponential_distribution<> >::
 | 
|---|
| 92 | RandomNumberDistribution_Encapsulation(const RandomNumberDistribution_Parameters&_params) :
 | 
|---|
| 93 |   distribution_type(_params.lambda)
 | 
|---|
| 94 | {}
 | 
|---|
| 95 | 
 | 
|---|
| 96 | template <>
 | 
|---|
| 97 | RandomNumberDistribution_Encapsulation< boost::normal_distribution<> >::
 | 
|---|
| 98 | RandomNumberDistribution_Encapsulation(const RandomNumberDistribution_Parameters&_params) :
 | 
|---|
| 99 |   distribution_type(_params.mean, _params.sigma)
 | 
|---|
| 100 | {}
 | 
|---|
| 101 | 
 | 
|---|
| 102 | template <>
 | 
|---|
| 103 | RandomNumberDistribution_Encapsulation< boost::lognormal_distribution<> >::
 | 
|---|
| 104 | RandomNumberDistribution_Encapsulation(const RandomNumberDistribution_Parameters&_params) :
 | 
|---|
| 105 |   distribution_type(_params.mean, _params.sigma)
 | 
|---|
| 106 | {}
 | 
|---|
| 107 | 
 | 
|---|
| 108 | /* ===============  min() ======================= */
 | 
|---|
| 109 | 
 | 
|---|
| 110 | template <>
 | 
|---|
| 111 | double RandomNumberDistribution_Encapsulation< boost::uniform_smallint<> >::min() const
 | 
|---|
| 112 | {
 | 
|---|
| 113 |         return distribution_type.min();
 | 
|---|
| 114 | }
 | 
|---|
| 115 | 
 | 
|---|
| 116 | template <>
 | 
|---|
| 117 | double RandomNumberDistribution_Encapsulation< boost::uniform_int<> >::min() const
 | 
|---|
| 118 | {
 | 
|---|
| 119 |         return distribution_type.min();
 | 
|---|
| 120 | }
 | 
|---|
| 121 | 
 | 
|---|
| 122 | template <>
 | 
|---|
| 123 | double RandomNumberDistribution_Encapsulation< boost::uniform_01<> >::min() const
 | 
|---|
| 124 | {
 | 
|---|
| 125 |         return distribution_type.min();
 | 
|---|
| 126 | }
 | 
|---|
| 127 | 
 | 
|---|
| 128 | template <>
 | 
|---|
| 129 | double RandomNumberDistribution_Encapsulation< boost::uniform_real<> >::min() const
 | 
|---|
| 130 | {
 | 
|---|
| 131 |         return distribution_type.min();
 | 
|---|
| 132 | }
 | 
|---|
| 133 | 
 | 
|---|
| 134 | template <>
 | 
|---|
| 135 | double RandomNumberDistribution_Encapsulation< boost::bernoulli_distribution<> >::min() const
 | 
|---|
| 136 | {
 | 
|---|
| 137 |   return 0.;
 | 
|---|
| 138 | }
 | 
|---|
| 139 | 
 | 
|---|
| 140 | template <>
 | 
|---|
| 141 | double RandomNumberDistribution_Encapsulation< boost::binomial_distribution<> >::min() const
 | 
|---|
| 142 | {
 | 
|---|
| 143 |   return 0.;
 | 
|---|
| 144 | }
 | 
|---|
| 145 | 
 | 
|---|
| 146 | template <>
 | 
|---|
| 147 | double RandomNumberDistribution_Encapsulation< boost::cauchy_distribution<> >::min() const
 | 
|---|
| 148 | {
 | 
|---|
| 149 |   return std::numeric_limits<double>::min();
 | 
|---|
| 150 | }
 | 
|---|
| 151 | 
 | 
|---|
| 152 | template <>
 | 
|---|
| 153 | double RandomNumberDistribution_Encapsulation< boost::gamma_distribution<> >::min() const
 | 
|---|
| 154 | {
 | 
|---|
| 155 |   return 0.;
 | 
|---|
| 156 | }
 | 
|---|
| 157 | 
 | 
|---|
| 158 | template <>
 | 
|---|
| 159 | double RandomNumberDistribution_Encapsulation< boost::poisson_distribution<> >::min() const
 | 
|---|
| 160 | {
 | 
|---|
| 161 |   return 0.;
 | 
|---|
| 162 | }
 | 
|---|
| 163 | 
 | 
|---|
| 164 | template <>
 | 
|---|
| 165 | double RandomNumberDistribution_Encapsulation< boost::geometric_distribution<> >::min() const
 | 
|---|
| 166 | {
 | 
|---|
| 167 |   return 0.;
 | 
|---|
| 168 | }
 | 
|---|
| 169 | 
 | 
|---|
| 170 | template <>
 | 
|---|
| 171 | double RandomNumberDistribution_Encapsulation< boost::triangle_distribution<> >::min() const
 | 
|---|
| 172 | {
 | 
|---|
| 173 |   return distribution_type.a();
 | 
|---|
| 174 | }
 | 
|---|
| 175 | 
 | 
|---|
| 176 | template <>
 | 
|---|
| 177 | double RandomNumberDistribution_Encapsulation< boost::exponential_distribution<> >::min() const
 | 
|---|
| 178 | {
 | 
|---|
| 179 |   return 0.;
 | 
|---|
| 180 | }
 | 
|---|
| 181 | 
 | 
|---|
| 182 | template <>
 | 
|---|
| 183 | double RandomNumberDistribution_Encapsulation< boost::normal_distribution<> >::min() const
 | 
|---|
| 184 | {
 | 
|---|
| 185 |   return std::numeric_limits<double>::min();
 | 
|---|
| 186 | }
 | 
|---|
| 187 | 
 | 
|---|
| 188 | template <>
 | 
|---|
| 189 | double RandomNumberDistribution_Encapsulation< boost::lognormal_distribution<> >::min() const
 | 
|---|
| 190 | {
 | 
|---|
| 191 |   return 0.;
 | 
|---|
| 192 | }
 | 
|---|
| 193 | 
 | 
|---|
| 194 | /* ===============  max() ======================= */
 | 
|---|
| 195 | 
 | 
|---|
| 196 | template <>
 | 
|---|
| 197 | double RandomNumberDistribution_Encapsulation< boost::uniform_smallint<> >::max() const
 | 
|---|
| 198 | {
 | 
|---|
| 199 |   return distribution_type.max();
 | 
|---|
| 200 | }
 | 
|---|
| 201 | 
 | 
|---|
| 202 | template <>
 | 
|---|
| 203 | double RandomNumberDistribution_Encapsulation< boost::uniform_int<> >::max() const
 | 
|---|
| 204 | {
 | 
|---|
| 205 |   return distribution_type.max();
 | 
|---|
| 206 | }
 | 
|---|
| 207 | 
 | 
|---|
| 208 | template <>
 | 
|---|
| 209 | double RandomNumberDistribution_Encapsulation< boost::uniform_01<> >::max() const
 | 
|---|
| 210 | {
 | 
|---|
| 211 |   return distribution_type.max();
 | 
|---|
| 212 | }
 | 
|---|
| 213 | 
 | 
|---|
| 214 | template <>
 | 
|---|
| 215 | double RandomNumberDistribution_Encapsulation< boost::uniform_real<> >::max() const
 | 
|---|
| 216 | {
 | 
|---|
| 217 |   return distribution_type.max();
 | 
|---|
| 218 | }
 | 
|---|
| 219 | 
 | 
|---|
| 220 | template <>
 | 
|---|
| 221 | double RandomNumberDistribution_Encapsulation< boost::bernoulli_distribution<> >::max() const
 | 
|---|
| 222 | {
 | 
|---|
| 223 |   return 1.;
 | 
|---|
| 224 | }
 | 
|---|
| 225 | 
 | 
|---|
| 226 | template <>
 | 
|---|
| 227 | double RandomNumberDistribution_Encapsulation< boost::binomial_distribution<> >::max() const
 | 
|---|
| 228 | {
 | 
|---|
| 229 |   return distribution_type.t();
 | 
|---|
| 230 | }
 | 
|---|
| 231 | 
 | 
|---|
| 232 | template <>
 | 
|---|
| 233 | double RandomNumberDistribution_Encapsulation< boost::cauchy_distribution<> >::max() const
 | 
|---|
| 234 | {
 | 
|---|
| 235 |   return std::numeric_limits<double>::max();
 | 
|---|
| 236 | }
 | 
|---|
| 237 | 
 | 
|---|
| 238 | template <>
 | 
|---|
| 239 | double RandomNumberDistribution_Encapsulation< boost::gamma_distribution<> >::max() const
 | 
|---|
| 240 | {
 | 
|---|
| 241 |   return std::numeric_limits<double>::max();
 | 
|---|
| 242 | }
 | 
|---|
| 243 | 
 | 
|---|
| 244 | template <>
 | 
|---|
| 245 | double RandomNumberDistribution_Encapsulation< boost::poisson_distribution<> >::max() const
 | 
|---|
| 246 | {
 | 
|---|
| 247 |   return std::numeric_limits<double>::max();
 | 
|---|
| 248 | }
 | 
|---|
| 249 | 
 | 
|---|
| 250 | template <>
 | 
|---|
| 251 | double RandomNumberDistribution_Encapsulation< boost::geometric_distribution<> >::max() const
 | 
|---|
| 252 | {
 | 
|---|
| 253 |   return std::numeric_limits<double>::max();
 | 
|---|
| 254 | }
 | 
|---|
| 255 | 
 | 
|---|
| 256 | template <>
 | 
|---|
| 257 | double RandomNumberDistribution_Encapsulation< boost::triangle_distribution<> >::max() const
 | 
|---|
| 258 | {
 | 
|---|
| 259 |   return distribution_type.c();
 | 
|---|
| 260 | }
 | 
|---|
| 261 | 
 | 
|---|
| 262 | template <>
 | 
|---|
| 263 | double RandomNumberDistribution_Encapsulation< boost::exponential_distribution<> >::max() const
 | 
|---|
| 264 | {
 | 
|---|
| 265 |   return std::numeric_limits<double>::max();
 | 
|---|
| 266 | }
 | 
|---|
| 267 | 
 | 
|---|
| 268 | template <>
 | 
|---|
| 269 | double RandomNumberDistribution_Encapsulation< boost::normal_distribution<> >::max() const
 | 
|---|
| 270 | {
 | 
|---|
| 271 |   return std::numeric_limits<double>::max();
 | 
|---|
| 272 | }
 | 
|---|
| 273 | 
 | 
|---|
| 274 | template <>
 | 
|---|
| 275 | double RandomNumberDistribution_Encapsulation< boost::lognormal_distribution<> >::max() const
 | 
|---|
| 276 | {
 | 
|---|
| 277 |   return std::numeric_limits<double>::max();
 | 
|---|
| 278 | }
 | 
|---|
| 279 | 
 | 
|---|
| 280 | /* ===============  p() ======================= */
 | 
|---|
| 281 | 
 | 
|---|
| 282 | template <>
 | 
|---|
| 283 | double RandomNumberDistribution_Encapsulation< boost::bernoulli_distribution<> >::p() const
 | 
|---|
| 284 | {
 | 
|---|
| 285 |   return distribution_type.p();
 | 
|---|
| 286 | }
 | 
|---|
| 287 | 
 | 
|---|
| 288 | template <>
 | 
|---|
| 289 | double RandomNumberDistribution_Encapsulation< boost::binomial_distribution<> >::p() const
 | 
|---|
| 290 | {
 | 
|---|
| 291 |   return distribution_type.p();
 | 
|---|
| 292 | }
 | 
|---|
| 293 | 
 | 
|---|
| 294 | template <>
 | 
|---|
| 295 | double RandomNumberDistribution_Encapsulation< boost::geometric_distribution<> >::p() const
 | 
|---|
| 296 | {
 | 
|---|
| 297 |   return distribution_type.p();
 | 
|---|
| 298 | }
 | 
|---|
| 299 | 
 | 
|---|
| 300 | 
 | 
|---|
| 301 | /* ===============  t() ======================= */
 | 
|---|
| 302 | 
 | 
|---|
| 303 | template <>
 | 
|---|
| 304 | double RandomNumberDistribution_Encapsulation< boost::binomial_distribution<> >::t() const
 | 
|---|
| 305 | {
 | 
|---|
| 306 |   return distribution_type.t();
 | 
|---|
| 307 | }
 | 
|---|
| 308 | 
 | 
|---|
| 309 | /* ===============  median() ======================= */
 | 
|---|
| 310 | 
 | 
|---|
| 311 | template <>
 | 
|---|
| 312 | double RandomNumberDistribution_Encapsulation< boost::cauchy_distribution<> >::median() const
 | 
|---|
| 313 | {
 | 
|---|
| 314 |   return distribution_type.median();
 | 
|---|
| 315 | }
 | 
|---|
| 316 | 
 | 
|---|
| 317 | /* ===============  sigma() ======================= */
 | 
|---|
| 318 | 
 | 
|---|
| 319 | template <>
 | 
|---|
| 320 | double RandomNumberDistribution_Encapsulation< boost::cauchy_distribution<> >::sigma() const
 | 
|---|
| 321 | {
 | 
|---|
| 322 |   return distribution_type.sigma();
 | 
|---|
| 323 | }
 | 
|---|
| 324 | 
 | 
|---|
| 325 | template <>
 | 
|---|
| 326 | double RandomNumberDistribution_Encapsulation< boost::normal_distribution<> >::sigma() const
 | 
|---|
| 327 | {
 | 
|---|
| 328 |   return distribution_type.sigma();
 | 
|---|
| 329 | }
 | 
|---|
| 330 | 
 | 
|---|
| 331 | template <>
 | 
|---|
| 332 | double RandomNumberDistribution_Encapsulation< boost::lognormal_distribution<> >::sigma() const
 | 
|---|
| 333 | {
 | 
|---|
| 334 |   return distribution_type.sigma();
 | 
|---|
| 335 | }
 | 
|---|
| 336 | 
 | 
|---|
| 337 | /* ===============  alpha() ======================= */
 | 
|---|
| 338 | 
 | 
|---|
| 339 | template <>
 | 
|---|
| 340 | double RandomNumberDistribution_Encapsulation< boost::gamma_distribution<> >::alpha() const
 | 
|---|
| 341 | {
 | 
|---|
| 342 |   return distribution_type.alpha();
 | 
|---|
| 343 | }
 | 
|---|
| 344 | 
 | 
|---|
| 345 | 
 | 
|---|
| 346 | /* ===============  mean() ======================= */
 | 
|---|
| 347 | 
 | 
|---|
| 348 | template <>
 | 
|---|
| 349 | double RandomNumberDistribution_Encapsulation< boost::poisson_distribution<> >::mean() const
 | 
|---|
| 350 | {
 | 
|---|
| 351 |   return distribution_type.mean();
 | 
|---|
| 352 | }
 | 
|---|
| 353 | 
 | 
|---|
| 354 | template <>
 | 
|---|
| 355 | double RandomNumberDistribution_Encapsulation< boost::normal_distribution<> >::mean() const
 | 
|---|
| 356 | {
 | 
|---|
| 357 |   return distribution_type.mean();
 | 
|---|
| 358 | }
 | 
|---|
| 359 | 
 | 
|---|
| 360 | template <>
 | 
|---|
| 361 | double RandomNumberDistribution_Encapsulation< boost::lognormal_distribution<> >::mean() const
 | 
|---|
| 362 | {
 | 
|---|
| 363 |   return distribution_type.mean();
 | 
|---|
| 364 | }
 | 
|---|
| 365 | 
 | 
|---|
| 366 | /* ===============  a() ======================= */
 | 
|---|
| 367 | 
 | 
|---|
| 368 | template <>
 | 
|---|
| 369 | double RandomNumberDistribution_Encapsulation< boost::triangle_distribution<> >::a() const
 | 
|---|
| 370 | {
 | 
|---|
| 371 |   return distribution_type.a();
 | 
|---|
| 372 | }
 | 
|---|
| 373 | 
 | 
|---|
| 374 | 
 | 
|---|
| 375 | /* ===============  b() ======================= */
 | 
|---|
| 376 | 
 | 
|---|
| 377 | template <>
 | 
|---|
| 378 | double RandomNumberDistribution_Encapsulation< boost::triangle_distribution<> >::b() const
 | 
|---|
| 379 | {
 | 
|---|
| 380 |   return distribution_type.b();
 | 
|---|
| 381 | }
 | 
|---|
| 382 | 
 | 
|---|
| 383 | 
 | 
|---|
| 384 | /* ===============  c() ======================= */
 | 
|---|
| 385 | 
 | 
|---|
| 386 | template <>
 | 
|---|
| 387 | double RandomNumberDistribution_Encapsulation< boost::triangle_distribution<> >::c() const
 | 
|---|
| 388 | {
 | 
|---|
| 389 |   return distribution_type.c();
 | 
|---|
| 390 | }
 | 
|---|
| 391 | 
 | 
|---|
| 392 | 
 | 
|---|
| 393 | /* ===============  lambda() ======================= */
 | 
|---|
| 394 | 
 | 
|---|
| 395 | template <>
 | 
|---|
| 396 | double RandomNumberDistribution_Encapsulation< boost::exponential_distribution<> >::lambda() const
 | 
|---|
| 397 | {
 | 
|---|
| 398 |   return distribution_type.lambda();
 | 
|---|
| 399 | }
 | 
|---|
| 400 | 
 | 
|---|
| 401 | 
 | 
|---|