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
|
//------------------------------------------------------------------------------
// GB_math_macros.h: simple math macros
//------------------------------------------------------------------------------
// SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2025, All Rights Reserved.
// SPDX-License-Identifier: Apache-2.0
//------------------------------------------------------------------------------
#ifndef GB_MATH_MACROS_H
#define GB_MATH_MACROS_H
//------------------------------------------------------------------------------
// min, max, and NaN handling
//------------------------------------------------------------------------------
// For floating-point computations, SuiteSparse:GraphBLAS relies on the IEEE
// 754 standard for the basic operations (+ - / *). Comparator also
// work as they should; any compare with NaN is always false, even
// eq(NaN,NaN) is false. This follows the IEEE 754 standard.
// For integer MIN and MAX, SuiteSparse:GraphBLAS relies on one comparator:
// z = min(x,y) = (x < y) ? x : y
// z = max(x,y) = (x > y) ? x : y
// However, this is not suitable for floating-point x and y. Compares with
// NaN always return false, so if either x or y are NaN, then z = y, for both
// min(x,y) and max(x,y).
// The C11 fmin, fminf, fmax, and fmaxf functions have the 'omitnan'
// behavior. These are used in SuiteSparse:GraphBLAS v2.3.0 and later.
// for integers only:
#define GB_IABS(x) (((x) >= 0) ? (x) : (-(x)))
// suitable for integers, and non-NaN floating point:
#define GB_IMAX(x,y) (((x) > (y)) ? (x) : (y))
#define GB_IMIN(x,y) (((x) < (y)) ? (x) : (y))
//------------------------------------------------------------------------------
// left and right shift
//------------------------------------------------------------------------------
#define GB_LSHIFT(x,k) (((uint64_t) (x)) << k)
#define GB_RSHIFT(x,k,b) (((x) >> (k)) & ((((uint64_t) 0x1) << (b)) - 1))
#endif
|
