tests/hkl-source-t.c

   1 /* This file is part of the hkl library.
   2  *
   3  * The hkl library is free software: you can redistribute it and/or modify
   4  * it under the terms of the GNU General Public License as published by
   5  * the Free Software Foundation, either version 3 of the License, or
   6  * (at your option) any later version.
   7  *
   8  * The hkl library is distributed in the hope that it will be useful,
   9  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  10  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  11  * GNU General Public License for more details.
  12  *
  13  * You should have received a copy of the GNU General Public License
  14  * along with the hkl library.  If not, see <http://www.gnu.org/licenses/>.
  15  *
  16  * Copyright (C) 2003-2019 Synchrotron SOLEIL
  17  *                         L'Orme des Merisiers Saint-Aubin
  18  *                         BP 48 91192 GIF-sur-YVETTE CEDEX
  19  *
  20  * Authors: Picca Frédéric-Emmanuel <picca@synchrotron-soleil.fr>
  21  */
  22 #include "hkl.h"
  23 #include <tap/basic.h>
  24 #include <tap/float.h>
  25
  26 #include "hkl-source-private.h"
  27
  28 static void new_copy(void)
  29 {
  30         HklSource s, c;
  31
  32         hkl_source_init(&s, 1.54, 1, 0, 0);
  33         c = s;
  34
  35         is_double(c.wave_length, s.wave_length, HKL_EPSILON, __func__);
  36         ok(FALSE == hkl_vector_cmp(&c.direction, &s.direction), __func__);
  37 }
  38
  39 static void init(void)
  40 {
  41         HklSource s;
  42
  43         hkl_source_init(&s, 1, 1, 0, 0);
  44
  45         is_double(1., s.wave_length, HKL_EPSILON, __func__);
  46         is_double(1., s.direction.data[0], HKL_EPSILON, __func__);
  47         is_double(0., s.direction.data[1], HKL_EPSILON, __func__);
  48         is_double(0., s.direction.data[2], HKL_EPSILON, __func__);
  49 }
  50
  51 static void cmp(void)
  52 {
  53         HklSource ref, s1, s2;
  54
  55         hkl_source_init(&ref, 1.54, 1, 0, 0);
  56         hkl_source_init(&s1, 1.54, 1, 0, 0);
  57         hkl_source_init(&s2, 1, 1, 0, 0);
  58
  59         ok(TRUE == hkl_source_cmp(&ref, &s1), __func__);
  60         ok(FALSE == hkl_source_cmp(&ref, &s2), __func__);
  61 }
  62
  63 static void compute_ki(void)
  64 {
  65         HklSource s;
  66         HklVector ki_ref = {{HKL_TAU / 1.54, 0, 0}};
  67         HklVector ki;
  68
  69         hkl_source_init(&s, 1.54, 1, 0, 0);
  70
  71         hkl_source_compute_ki(&s, &ki);
  72         ok(0 == hkl_vector_cmp(&ki_ref, &ki), __func__);
  73 }
  74
  75 static void get_wavelength(void)
  76 {
  77         HklSource s;
  78
  79         hkl_source_init(&s, 1, 1, 0, 0);
  80
  81         is_double(1., hkl_source_get_wavelength(&s), HKL_EPSILON, __func__);
  82 }
  83
  84 int main(void)
  85 {
  86         plan(10);
  87
  88         new_copy();
  89         init();
  90         cmp();
  91         compute_ki();
  92         get_wavelength();
  93
  94         return 0;
  95 }