diff options
author | Thomas White <taw@physics.org> | 2021-10-28 15:00:56 +0200 |
---|---|---|
committer | Thomas White <taw@physics.org> | 2021-10-28 15:00:56 +0200 |
commit | 3dfb24e9cacf8bc62f39b0de37c5b0ddceba9e71 (patch) | |
tree | ef709dab75a8527675751424d8e55b06070aa098 /tests/symop_parse.c | |
parent | 86633959138c6ec528c8cf635c6f0d39d79b61cc (diff) |
Make symmetry operation parser re-entrant
Fixes: https://gitlab.desy.de/thomas.white/crystfel/-/issues/1
Diffstat (limited to 'tests/symop_parse.c')
-rw-r--r-- | tests/symop_parse.c | 68 |
1 files changed, 68 insertions, 0 deletions
diff --git a/tests/symop_parse.c b/tests/symop_parse.c new file mode 100644 index 00000000..0865ff9b --- /dev/null +++ b/tests/symop_parse.c @@ -0,0 +1,68 @@ +/* + * symop_parse.c + * + * Check that symmetry operation parsing works + * + * Copyright © 2021 Deutsches Elektronen-Synchrotron DESY, + * a research centre of the Helmholtz Association. + * + * Authors: + * 2021 Thomas White <taw@physics.org> + * + * This file is part of CrystFEL. + * + * CrystFEL is free software: you can redistribute it and/or modify + * it under the terms of the GNU General Public License as published by + * the Free Software Foundation, either version 3 of the License, or + * (at your option) any later version. + * + * CrystFEL is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the + * GNU General Public License for more details. + * + * You should have received a copy of the GNU General Public License + * along with CrystFEL. If not, see <http://www.gnu.org/licenses/>. + * + */ + +#include <stdio.h> +#include <string.h> +#include <stdarg.h> + +#include <symmetry.h> +#include <rational.h> + +int main(int argc, char *argv[]) +{ + int r = 0; + RationalMatrix *mtx; + SymOpList *sym; + + mtx = parse_symmetry_operation("h,k,l"); + if ( !rtnl_mtx_is_identity(mtx) ) { + printf("h,k,l not an identity:\n"); + rtnl_mtx_print(mtx); + r = 1; + } + + mtx = parse_symmetry_operation("k,h,-l"); + + mtx = parse_symmetry_operation("h,k,l"); + if ( !rtnl_mtx_is_identity(mtx) ) { + printf("h,k,l not an identity on second attempt:\n"); + rtnl_mtx_print(mtx); + r = 1; + } + + sym = parse_symmetry_operations("h,k,l;k,h,-l;-h,-k,l"); + if ( sym == NULL ) r = 1; + + mtx = parse_symmetry_operation("h,k,fail"); + if ( mtx != NULL ) r = 1; + + sym = parse_symmetry_operations("k,h,-l;h,k,fail;h,k,l"); + if ( sym != NULL ) r = 1; + + return r; +} |