blob: 6302306f63d860ece3a4d91d4111afa77f2dbeb2 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
|
/* $Header$ */
#ifndef NOISE_H
#define NOISE_H
#include <stdlib.h>
#include "cctk.h"
/* constants */
#define TRUE 1
#define FALSE 0
#define MAXDIM 3
#
/* macros */
#define RAND_VAL random()*amplitude - offset
#define INDEX_3D(lsh, i, j, k) ((i) + (lsh)[0]*((j) + (lsh)[1]*(k)))
#define NOISY_BOUNDARY(doBC, \
iend, jend, kend, \
ii, jj, kk, \
rand) \
{ \
if (doBC) \
{ \
srandom(time(0)); \
for (k = 0; k < kend; k++) \
{ \
for (j = 0; j < jend; j++) \
{ \
for (i = 0; i < iend; i++) \
{ \
int _index = CCTK_GFINDEX3D(GH, ii, jj, kk); \
CCTK_REAL* v= (CCTK_REAL *) GH->data[var][timelvl] + _index; \
v[0] += RAND_VAL; \
} \
} \
} \
} \
}
#endif /* !define(NOISE_H) */
|