From a56aafbeae1b6438c84cb433bfb17023581d5527 Mon Sep 17 00:00:00 2001 From: Thomas White Date: Wed, 11 May 2011 18:03:15 +0200 Subject: Formatting --- src/partial_sim.c | 21 ++++++++++++--------- 1 file changed, 12 insertions(+), 9 deletions(-) (limited to 'src/partial_sim.c') 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); } -- cgit v1.2.3