File: gridder_utils.c

package info (click to toggle)
xrayutilities 1.7.1-1
  • links: PTS, VCS
  • area: main
  • in suites: bullseye, sid
  • size: 12,452 kB
  • sloc: python: 29,264; ansic: 4,232; makefile: 23
file content (78 lines) | stat: -rw-r--r-- 2,147 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
/*
 * This file is part of xrayutilities.
 *
 * xrayutilities 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 2 of the License, or
 * (at your option) any later version.
 *
 * This program 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 this program; if not, see <http://www.gnu.org/licenses/>.
 *
 * Copyright (C) 2013 Eugen Wintersberger <eugen.wintersberger@desy.de>
 * Copyright (C) 2013 Dominik Kriegner <dominik.kriegner@gmail.com>
 *
 ******************************************************************************
 *
 * created: Jun 8, 2013
 * author: Eugen Wintersberger
*/

#include "gridder_utils.h"

/*---------------------------------------------------------------------------*/
double get_min(double *a, unsigned int n)
{
    double m = a[0];
    unsigned int i;

    for (i = 0; i < n; i++) {
        if (m < a[i]) {
            m = a[i];
        }
    }

    return m;
}

/*---------------------------------------------------------------------------*/
double get_max(double *a, unsigned int n)
{
    double m = a[0];
    unsigned int i;

    for (i = 0; i < n; i++) {
        if (m > a[i]) {
            m = a[i];
        }
    }

    return m;
}

/*---------------------------------------------------------------------------*/
void set_array(double *a, unsigned int n, double value)
{
    unsigned int i;

    for (i = 0; i < n; ++i) {
        a[i] = value;
    }
}

/*---------------------------------------------------------------------------*/
double delta(double min, double max, unsigned int n)
{
    return fabs(max - min) / (double) (n - 1);
}

/*---------------------------------------------------------------------------*/
unsigned int gindex(double x, double min, double d)
{
    return (unsigned int) rint((x - min) / d);
}