| [081e5d] | 1 | /*
 | 
|---|
 | 2 |  * Project: MoleCuilder
 | 
|---|
 | 3 |  * Description: creates and alters molecular systems
 | 
|---|
 | 4 |  * Copyright (C)  2010 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 | 
 | 
|---|
| [ecb6c5] | 24 | #include <limits>
 | 
|---|
 | 25 | 
 | 
|---|
| [081e5d] | 26 | #include "RandomNumberDistribution_Encapsulation.hpp"
 | 
|---|
 | 27 | 
 | 
|---|
| [0275ad] | 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 | {}
 | 
|---|
| [081e5d] | 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 | 
 | 
|---|
| [ecb6c5] | 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 | }
 | 
|---|
| [081e5d] | 187 | 
 | 
|---|
| [ecb6c5] | 188 | template <>
 | 
|---|
 | 189 | double RandomNumberDistribution_Encapsulation< boost::lognormal_distribution<> >::min() const
 | 
|---|
 | 190 | {
 | 
|---|
 | 191 |   return 0.;
 | 
|---|
 | 192 | }
 | 
|---|
| [081e5d] | 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 | 
 | 
|---|
| [ecb6c5] | 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 | }
 | 
|---|
| [081e5d] | 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 | 
 | 
|---|