-
Notifications
You must be signed in to change notification settings - Fork 375
/
Copy pathgmt_common_math.h
107 lines (93 loc) · 3.58 KB
/
gmt_common_math.h
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
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
/*--------------------------------------------------------------------
*
* Copyright (c) 1991-2025 by the GMT Team (https://www.generic-mapping-tools.org/team.html)
* See LICENSE.TXT file for copying and redistribution conditions.
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU Lesser General Public License as published by
* the Free Software Foundation; version 3 or 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 Lesser General Public License for more details.
*
* Contact info: www.generic-mapping-tools.org
*--------------------------------------------------------------------*/
/*
* gmt_common_math.h declares shared math functions
*
* Author: Florian Wobbe
* Date: 10-MAR-2012
* Version: 5
*/
/*!
* \file gmt_common_math.h
* \brief declares shared math functions
*/
#pragma once
#ifndef GMT_COMMON_MATH_H
#define GMT_COMMON_MATH_H
#ifdef __cplusplus /* Basic C++ support */
extern "C" {
#endif
/* CMake definitions: This must be first! */
#include "gmt_config.h"
/* Declaration modifiers for DLL support (MSC et al) */
#include "declspec.h"
#include <stdint.h>
#include <stdbool.h>
/* int32_abs function that works with int32_t */
#if defined(SIZEOF_INT) && SIZEOF_INT == 4 && defined(HAVE_ABS)
# define int32_abs abs
#elif defined(SIZEOF_LONG) && SIZEOF_LONG == 4
# define int32_abs labs
#else
# define int32_abs(x) ((int32_t)(((x) >= 0) ? (x) : -(x)))
#endif
/* int64_abs function that works with int64_t */
#if defined(_WIN64)
# define int64_abs _abs64
#elif defined(SIZEOF_INT) && SIZEOF_INT == 8 && defined(HAVE_ABS)
# define int64_abs abs
#elif defined(SIZEOF_LONG) && SIZEOF_LONG == 8
# define int64_abs labs
#elif defined(SIZEOF_LONG_LONG) && SIZEOF_LONG_LONG == 8 && defined(HAVE_LLABS)
# define int64_abs llabs
#else
# define int64_abs(x) ((int64_t)(((x) >= 0) ? (x) : -(x)))
#endif
#ifdef __CYGWIN__ /* See http://gmt.soest.hawaii.edu/boards/1/topics/5428 */
#ifdef __x86_64
#define lrint(x) ((long int)(int)lrint(x))
#endif
#endif
/* Limit casting to one place (here) for dropping lrint output to signed or unsigned ints */
#define irint(x) ((int)lrint(x))
#define urint(x) ((unsigned int)lrint(x))
#define irintf(x) ((int)lrintf(x))
#define urintf(x) ((unsigned int)lrintf(x))
/* Safe rounding of float and double to signed and unsigned 64 bit ints */
#if defined(SIZEOF_INT) && SIZEOF_LONG == 8
# define irint64(x) lrint(x)
# define urint64(x) ((uint64_t)lrint(x))
# define irint64f(x) lrintf(x)
# define urint64f(x) ((uint64_t)lrintf(x))
#else /* SIZEOF_LONG_LONG == 8 by ISO C definition */
# define irint64(x) llrint(x)
# define urint64(x) ((uint64_t)llrint(x))
# define irint64f(x) llrintf(x)
# define urint64f(x) ((uint64_t)llrintf(x))
#endif
EXTERN_MSC bool floatAlmostEqualUlpsAndAbs(float A, float B, float maxDiff, int maxUlpsDiff);
EXTERN_MSC bool doubleAlmostEqualUlpsAndAbs(double A, double B, double maxDiff, int maxUlpsDiff);
EXTERN_MSC bool floatAlmostEqualUlps(float A, float B, int maxUlpsDiff);
EXTERN_MSC bool doubleAlmostEqualUlps(double A, double B, int maxUlpsDiff);
# define floatAlmostEqual(A, B) (floatAlmostEqualUlps(A, B, 5))
# define doubleAlmostEqual(A, B) (doubleAlmostEqualUlps(A, B, 5))
# define floatAlmostEqualZero(A, B) (floatAlmostEqualUlpsAndAbs(A, B, 5*FLT_EPSILON, 5))
# define doubleAlmostEqualZero(A, B) (doubleAlmostEqualUlpsAndAbs(A, B, 5*DBL_EPSILON, 5))
#ifdef __cplusplus
}
#endif
#endif /* !GMT_COMMON_MATH_H */