diff options
| -rw-r--r-- | src/partial_sim.c | 20 | ||||
| -rw-r--r-- | src/post-refinement.c | 2 | ||||
| -rw-r--r-- | src/post-refinement.h | 10 |
3 files changed, 20 insertions, 12 deletions
diff --git a/src/partial_sim.c b/src/partial_sim.c index 2a412ae4..21068399 100644 --- a/src/partial_sim.c +++ b/src/partial_sim.c @@ -38,19 +38,25 @@ static void mess_up_cell(UnitCell *cell) double cx, cy, cz; /* Cell noise in percent */ - const double cnoise = 2.0; + const double cnoise = 0.5; + + STATUS("Real:\n"); + cell_print(cell); cell_get_reciprocal(cell, &ax, &ay, &az, &bx, &by, &bz, &cx, &cy, &cz); 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); +// 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); + + STATUS("Changed:\n"); + cell_print(cell); } diff --git a/src/post-refinement.c b/src/post-refinement.c index 7301b57b..64704598 100644 --- a/src/post-refinement.c +++ b/src/post-refinement.c @@ -532,4 +532,6 @@ void pr_refine(struct image *image, const RefList *full, const char *sym) i++; } while ( (max_shift > 0.01) && (i < MAX_CYCLES) ); + + cell_print(image->indexed_cell); } diff --git a/src/post-refinement.h b/src/post-refinement.h index 5b71cbd5..16db022c 100644 --- a/src/post-refinement.h +++ b/src/post-refinement.h @@ -27,15 +27,15 @@ /* Refineable parameters */ enum { REF_ASX, - REF_BSX, - REF_CSX, REF_ASY, - REF_BSY, - REF_CSY, REF_ASZ, + NUM_PARAMS, + REF_BSX, + REF_BSY, REF_BSZ, + REF_CSX, + REF_CSY, REF_CSZ, - NUM_PARAMS, REF_DIV, REF_R, }; |