blob: 749b41157ff1f6e55a09011ad1fdd8bb8dca9978 (
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
|
/* $Header$ */
#ifndef NOISE_H
#define NOISE_H
#include <assert.h>
#include <stdlib.h>
#include <string.h>
#include <time.h>
#include "cctk.h"
/* constants */
#define MAXDIM 3
/* macros */
#define RAND_VAL ((random()*(1.0/RAND_MAX)-0.5)*amplitude)
#define BOUNDARY_NOISE(doBC, \
iend, jend, kend, \
ii, jj, kk) \
{ \
if (doBC) \
{ \
CCTK_REAL* v= CCTK_VarDataPtrI(GH, timelvl, var); \
for (k = 0; k < kend; k++) \
{ \
for (j = 0; j < jend; j++) \
{ \
for (i = 0; i < iend; i++) \
{ \
const int _index = CCTK_GFINDEX3D(GH, (ii), (jj), (kk)); \
v[_index] += RAND_VAL; \
} \
} \
} \
} \
}
#endif /* !define(NOISE_H) */
|