File: gasdev.c

package info (click to toggle)
grass 8.4.2-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid
  • size: 277,040 kB
  • sloc: ansic: 460,798; python: 227,732; cpp: 42,026; sh: 11,262; makefile: 7,007; xml: 3,637; sql: 968; lex: 520; javascript: 484; yacc: 450; asm: 387; perl: 157; sed: 25; objc: 6; ruby: 4
file content (22 lines) | stat: -rw-r--r-- 579 bytes parent folder | download | duplicates (2) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
 
/* gasdev.c                                                             */
#include <math.h>
#include "ransurf.h"
#include "local_proto.h"

/* GasDev() returns a random double with a mean of 0.0 and a standard   */
/*      deviation of 1.0.                                               */
double GasDev(void)
{
    double fac, r, v1, v2;

    do {
        v1 = 2.0 * ran1() - 1.0;
        v2 = 2.0 * ran1() - 1.0;
        r = v1 * v1 + v2 * v2;
    } while (r >= 1.0);

    fac = sqrt(-2.0 * log(r) / r);
    G_debug(3, "(v2 * fac):%.12lf", v2 * fac);

    return (v2 * fac);
}