aboutsummaryrefslogtreecommitdiff
path: root/src
diff options
context:
space:
mode:
authorThomas White <taw@physics.org>2011-05-11 18:03:15 +0200
committerThomas White <taw@physics.org>2012-02-22 15:27:26 +0100
commita56aafbeae1b6438c84cb433bfb17023581d5527 (patch)
tree874163e8abef9bc60594c0b5fb0b2b37c5458ea0 /src
parent0910580f88ce9eeaab33f0b694a5fbac5884c642 (diff)
Formatting
Diffstat (limited to 'src')
-rw-r--r--src/partial_sim.c21
1 files changed, 12 insertions, 9 deletions
diff --git a/src/partial_sim.c b/src/partial_sim.c
index c0c44212..d9d7fc52 100644
--- a/src/partial_sim.c
+++ b/src/partial_sim.c
@@ -37,16 +37,19 @@ static void mess_up_cell(UnitCell *cell)
double bx, by, bz;
double cx, cy, cz;
+ /* Cell noise in percent */
+ const double cnoise = 2.0;
+
cell_get_reciprocal(cell, &ax, &ay, &az, &bx, &by, &bz, &cx, &cy, &cz);
- ax = gaussian_noise(ax, fabs(ax)/50.0);
- ay = gaussian_noise(ay, fabs(ay)/50.0);
- az = gaussian_noise(az, fabs(az)/50.0);
- bx = gaussian_noise(bx, fabs(bx)/50.0);
- by = gaussian_noise(by, fabs(by)/50.0);
- bz = gaussian_noise(bz, fabs(bz)/50.0);
- cx = gaussian_noise(cx, fabs(cx)/50.0);
- cy = gaussian_noise(cy, fabs(cy)/50.0);
- cz = gaussian_noise(cz, fabs(cz)/50.0);
+ ax = gaussian_noise(ax, cnoise*fabs(ax)/100.0);
+ ay = gaussian_noise(ay, cnoise*fabs(ay)/100.0);
+ az = gaussian_noise(az, cnoise*fabs(az)/100.0);
+ bx = gaussian_noise(bx, cnoise*fabs(bx)/100.0);
+ by = gaussian_noise(by, cnoise*fabs(by)/100.0);
+ bz = gaussian_noise(bz, cnoise*fabs(bz)/100.0);
+ cx = gaussian_noise(cx, cnoise*fabs(cx)/100.0);
+ cy = gaussian_noise(cy, cnoise*fabs(cy)/100.0);
+ cz = gaussian_noise(cz, cnoise*fabs(cz)/100.0);
cell_set_reciprocal(cell, ax, ay, az, bx, by, bz, cx, cy, cz);
}