src/noise.h


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

/* $Header$ */

#ifndef NOISE_H
#define NOISE_H

#include <stdlib.h>
#include <time.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 BOUNDARY_NOISE(doBC,                                                 \
                       iend, jend, kend,                                     \
                       ii, jj, kk)            \
{                                                                             \
if (doBC)                                        \
  {                                              \
    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) */