diff options
author | Thomas White <taw@physics.org> | 2013-04-11 18:15:34 +0200 |
---|---|---|
committer | Thomas White <taw@physics.org> | 2013-04-11 18:16:08 +0200 |
commit | f7314aa835d5379e3fdd76845f401c26790f171c (patch) | |
tree | 78cba94de8e870bbedf8c31890101c58a58f7e5a /src | |
parent | 2daf2a8fa448b22706ac67718084b88dcb6f613b (diff) |
check_hkl: Fix iteration limits for possible reflection calculation
Diffstat (limited to 'src')
-rw-r--r-- | src/check_hkl.c | 24 |
1 files changed, 12 insertions, 12 deletions
diff --git a/src/check_hkl.c b/src/check_hkl.c index a0acadb6..0f485654 100644 --- a/src/check_hkl.c +++ b/src/check_hkl.c @@ -89,9 +89,9 @@ static void plot_shells(RefList *list, UnitCell *cell, const SymOpList *sym, RefListIterator *iter; RefList *counted; int hmax, kmax, lmax; - double asx, asy, asz; - double bsx, bsy, bsz; - double csx, csy, csz; + double ax, ay, az; + double bx, by, bz; + double cx, cy, cz; possible = malloc(nshells*sizeof(int)); measurements = malloc(nshells*sizeof(unsigned int)); @@ -182,15 +182,15 @@ static void plot_shells(RefList *list, UnitCell *cell, const SymOpList *sym, /* Count the number of reflections possible in each shell */ counted = reflist_new(); - cell_get_reciprocal(cell, &asx, &asy, &asz, - &bsx, &bsy, &bsz, - &csx, &csy, &csz); - hmax = rmax / modulus(asx, asy, asz); - kmax = rmax / modulus(bsx, bsy, bsz); - lmax = rmax / modulus(csx, csy, csz); - for ( h=-hmax; h<hmax; h++ ) { - for ( k=-kmax; k<kmax; k++ ) { - for ( l=-lmax; l<lmax; l++ ) { + cell_get_cartesian(cell, &ax, &ay, &az, + &bx, &by, &bz, + &cx, &cy, &cz); + hmax = rmax * modulus(ax, ay, az); + kmax = rmax * modulus(bx, by, bz); + lmax = rmax * modulus(cx, cy, cz); + for ( h=-hmax; h<=hmax; h++ ) { + for ( k=-kmax; k<=kmax; k++ ) { + for ( l=-lmax; l<=lmax; l++ ) { double d; signed int hs, ks, ls; |