Random vector on a sphere: Difference between revisions
Jump to navigation
Jump to search
Carl McBride (talk | contribs) No edit summary |
Carl McBride (talk | contribs) m (Slight tidy) |
||
| Line 1: | Line 1: | ||
The ability to generate a randomly orientated vector is very useful in [[Monte Carlo]] simulations of anisotropic [[models]] | The ability to generate a randomly orientated vector is very useful in [[Monte Carlo]] simulations of anisotropic [[models]] | ||
or molecular systems. | or molecular systems. | ||
==Marsaglia algorithm== | ==Marsaglia algorithm== | ||
This is the algorithm proposed by George Marsaglia ( | This is the algorithm proposed by George Marsaglia <ref name="Marsaglia">[http://dx.doi.org/10.1214/aoms/1177692644 George Marsaglia "Choosing a Point from the Surface of a Sphere", The Annals of Mathematical Statistics '''43''' pp. 645-646 (1972)]</ref> | ||
*Independently generate V<sub>1</sub> and V<sub>2</sub>, taken from a uniform distribution on (-1,1) such that | *Independently generate V<sub>1</sub> and V<sub>2</sub>, taken from a uniform distribution on (-1,1) such that | ||
:<math>S=(V_1^2+V_2^2) < 1</math> | :<math>S=(V_1^2+V_2^2) < 1</math> | ||
*The random vector is then ( | *The random vector is then (Eq. 4 in <ref name="Marsaglia"> </ref> ): | ||
:<math>\left(2V_1 \sqrt{1-S},~ 2V_2 \sqrt{1-S},~ 1-2S\right)</math> | :<math>\left(2V_1 \sqrt{1-S},~ 2V_2 \sqrt{1-S},~ 1-2S\right)</math> | ||
==Fortran 90 implementation== | ==Fortran 90 implementation== | ||
This Fortran 90 implementation is adapted from Ref. | This Fortran 90 implementation is adapted from Ref. <ref>Daan Frenkel and Berend Smit "Understanding Molecular Simulation: From Algorithms to Applications" p. 410 Academic Press (1996)</ref> The function '''ran()''' calls a [[Random_numbers | randon number generator]]: | ||
<small><pre> | <small><pre> | ||
! The following is taken from Allen & Tildesley, p. 349 | ! The following is taken from Allen & Tildesley, p. 349 | ||
| Line 41: | Line 40: | ||
==References== | ==References== | ||
<references/> | |||
'''Related reading''' | |||
* M. P. Allen and D. J. Tildesley "Computer Simulation of Liquids", p. 349 Clarendon Press (1989) | |||
{{Source}} | |||
[[category: random numbers]] | [[category: random numbers]] | ||
[[category: computer simulation techniques]] | [[category: computer simulation techniques]] | ||
Revision as of 12:54, 21 October 2010
The ability to generate a randomly orientated vector is very useful in Monte Carlo simulations of anisotropic models or molecular systems.
Marsaglia algorithm
This is the algorithm proposed by George Marsaglia [1]
- Independently generate V1 and V2, taken from a uniform distribution on (-1,1) such that
- The random vector is then (Eq. 4 in [1] ):
Fortran 90 implementation
This Fortran 90 implementation is adapted from Ref. [2] The function ran() calls a randon number generator:
! The following is taken from Allen & Tildesley, p. 349
! Generate a random vector towards a point in the unit sphere
! Daniel Duque 2004
subroutine random_vector(vctr)
implicit none
real, dimension(3) :: vctr
real:: ran1,ran2,ransq,ranh
real:: ran
do
ran1=1.0-2.0*ran()
ran2=1.0-2.0*ran()
ransq=ran1**2+ran2**2
if(ransq.le.1.0) exit
enddo
ranh=2.0*sqrt(1.0-ransq)
vctr(1)=ran1*ranh
vctr(2)=ran2*ranh
vctr(3)=(1.0-2.0*ransq)
end subroutine random_vector
References
- ↑ 1.0 1.1 George Marsaglia "Choosing a Point from the Surface of a Sphere", The Annals of Mathematical Statistics 43 pp. 645-646 (1972)
- ↑ Daan Frenkel and Berend Smit "Understanding Molecular Simulation: From Algorithms to Applications" p. 410 Academic Press (1996)
Related reading
- M. P. Allen and D. J. Tildesley "Computer Simulation of Liquids", p. 349 Clarendon Press (1989)
|
This page contains computer source code. If you intend to compile and use this code you must check for yourself the validity of the code. Please read the SklogWiki disclaimer. |