hkl/hkl-pseudoaxis-zaxis-hkl.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-2010 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 <gsl/gsl_math.h>
  23 #include <gsl/gsl_vector.h>
  24
  25 #include <ccan/array_size/array_size.h>
  26
  27 #include "hkl-pseudoaxis-auto-private.h"
  28 #include "hkl-pseudoaxis-common-hkl-private.h"
  29
  30 /* #define DEBUG */
  31
  32 /***********************/
  33 /* numerical functions */
  34 /***********************/
  35
  36 static int _reflectivity_func(const gsl_vector *x, void *params, gsl_vector *f)
  37 {
  38         const double mu = x->data[0];
  39         const double gamma = x->data[3];
  40
  41         CHECK_NAN(x->data, x->size);
  42
  43         RUBh_minus_Q(x->data, params, f->data);
  44         f->data[3] = mu - gamma;
  45
  46         return  GSL_SUCCESS;
  47 }
  48
  49 static const HklFunction reflectivity_func = {
  50         .function = _reflectivity_func,
  51         .size = 4,
  52 };
  53
  54 /********/
  55 /* mode */
  56 /********/
  57
  58 static HklMode* zaxis()
  59 {
  60         static const char* axes[] = {"omega", "delta", "gamma"};
  61         static const HklFunction *functions[] = {&RUBh_minus_Q_func};
  62         static const HklModeAutoInfo info = {
  63                 INFO_AUTO(__func__, axes, functions),
  64         };
  65
  66         return hkl_mode_auto_new(&info,
  67                                                     &hkl_full_mode_operations);
  68 }
  69
  70 static HklMode* reflectivity()
  71 {
  72         static const char* axes[] = {"mu", "omega", "delta", "gamma"};
  73         static const HklFunction *functions[] = {&reflectivity_func};
  74         static const HklModeAutoInfo info = {
  75                 INFO_AUTO(__func__, axes, functions),
  76         };
  77
  78         return hkl_mode_auto_new(&info,
  79                                                     &hkl_full_mode_operations);
  80 }
  81
  82 /**********************/
  83 /* pseudo axis engine */
  84 /**********************/
  85
  86 HklEngine *hkl_engine_zaxis_hkl_new(void)
  87 {
  88         HklEngine *self;
  89         HklMode *default_mode;
  90
  91         self = hkl_engine_hkl_new();
  92
  93         default_mode = zaxis();
  94         hkl_engine_add_mode(self, default_mode);
  95         hkl_engine_select_mode(self, default_mode);
  96
  97         hkl_engine_add_mode(self, reflectivity());
  98
  99         return self;
 100 }