mdtools/md_code/randomSPRNG.cpp

#include <iostream>
#include <cmath>

#include "randomSPRNG.hpp"
#include "simError.h"
#ifdef IS_MPI
#include "mpiSimulation.hpp"
#endif

/* randomStreamSPRNF creates a new SPRNG stream for random numbers
 */

int randomSPRNG::nStreamsInitialized;

randomSPRNG::randomSPRNG(iseed){
#ifdef IS_MPI
  nStreamsInitialized++;
  //Fix this
#error "Fix instance of mpiSimulation in randomSPRNG"
  nSPRNGStreams = getnumberProcessors();
  myStreamNumber = nStreamsInitialized*(getMyNode() + nSPRNGStreams);


#else
  nStreamsInitialized++;
  nSPRNGStreams = 1;
  myStreamNumber = nSPRNGStreams;

#endif

  thisStream = init_sprng(GTYPE,myStreamNumber,nSPRNGStreams,
                          iseed,SPRNG_DEFAULT);
}

randomSPRNG::~randomSPRNG(){

  free_sprng(int *thisStream);
  nStreamsInitialized--;

}


double randomSPRNG::getRandom(){
  double ranNum;
  ranNum = sprng(thisStream);
  return ranNum;
}


// Gaussian SPRNG class...

gaussianSPRNG::gaussianSPRNG() : randomSPRNG()
{
}

gaussianSPRNG::~gaussianSPRNG() : ~randomSPRNG()
{
}

double gaussianSPRNG::getGaussian(){
  double ranNum1;
  double ranNum2;
  double gaussianNumber;
  double R1;
  double R2;

  ranNum1 = getRandom();
  ranNum2 = getRandom();

  gaussianNumber = sqrt(-2.0 * log(ranNum1)) * cos(2 * M_PI * ranNum2);

  return gaussianNumber;
}
Revision:	221
Committed:	Thu Jan 2 20:14:08 2003 UTC (21 years, 6 months ago) by chuckv
File size:	1290 byte(s)
Log Message:	Thermo now can use SPRNG or rand48 (if not MPI). Finished up work on randomSPRNG. Still need to fix access to MPISimulation.
#	Content
1	#include <iostream>
2	#include <cmath>
3
4	#include "randomSPRNG.hpp"
5	#include "simError.h"
6	#ifdef IS_MPI
7	#include "mpiSimulation.hpp"
8	#endif
9
10	/* randomStreamSPRNF creates a new SPRNG stream for random numbers
11	*/
12
13	int randomSPRNG::nStreamsInitialized;
14
15	randomSPRNG::randomSPRNG(iseed){
16	#ifdef IS_MPI
17	nStreamsInitialized++;
18	//Fix this
19	#error "Fix instance of mpiSimulation in randomSPRNG"
20	nSPRNGStreams = getnumberProcessors();
21	myStreamNumber = nStreamsInitialized*(getMyNode() + nSPRNGStreams);
22
23
24	#else
25	nStreamsInitialized++;
26	nSPRNGStreams = 1;
27	myStreamNumber = nSPRNGStreams;
28
29	#endif
30
31	thisStream = init_sprng(GTYPE,myStreamNumber,nSPRNGStreams,
32	iseed,SPRNG_DEFAULT);
33	}
34
35	randomSPRNG::~randomSPRNG(){
36
37	free_sprng(int *thisStream);
38	nStreamsInitialized--;
39
40	}
41
42
43	double randomSPRNG::getRandom(){
44	double ranNum;
45	ranNum = sprng(thisStream);
46	return ranNum;
47	}
48
49
50	// Gaussian SPRNG class...
51
52	gaussianSPRNG::gaussianSPRNG() : randomSPRNG()
53	{
54	}
55
56	gaussianSPRNG::~gaussianSPRNG() : ~randomSPRNG()
57	{
58	}
59
60	double gaussianSPRNG::getGaussian(){
61	double ranNum1;
62	double ranNum2;
63	double gaussianNumber;
64	double R1;
65	double R2;
66
67	ranNum1 = getRandom();
68	ranNum2 = getRandom();
69
70	gaussianNumber = sqrt(-2.0 * log(ranNum1)) * cos(2 * M_PI * ranNum2);
71
72	return gaussianNumber;
73	}