500

I want a function that returns -1 for negative numbers and +1 for positive numbers. http://en.wikipedia.org/wiki/Sign_function It's easy enough to write my own, but it seems like something that ought to be in a standard library somewhere.

Edit: Specifically, I was looking for a function working on floats.

12
  • 26
    What should it return for 0? Dec 14, 2009 at 23:26
  • 84
    @Craig McQueen; that depends on if it is a positive zero or negative zero.
    – ysth
    Dec 15, 2009 at 5:53
  • 11
    @ysth @Craig McQueen, false for floats too, no? sgn(x)'s definition says to return 0 if x==0. According to IEEE 754, negative zero and positive zero should compare as equal.
    – RJFalconer
    Jun 4, 2014 at 11:28
  • 9
    @ysth "it depends on positive zero or negative zero". In fact, it does not.
    – RJFalconer
    Jun 6, 2014 at 8:10
  • 6
    Commenting late, but regarding signed zeros, another reasonable option is that sgn(x) returns x, when x is zero. In other words, you get 0 out, but it's a signed zero with the same sign as the input. @RJFalconer In the relatively few cases that signed zeros matter, you get a sensible answer, and in the other cases it makes no difference.
    – Brick
    Feb 15, 2018 at 19:14

21 Answers 21

606

The type-safe C++ version:

template <typename T> int sgn(T val) {
    return (T(0) < val) - (val < T(0));
}

Benefits:

  • Actually implements signum (-1, 0, or 1). Implementations here using copysign only return -1 or 1, which is not signum. Also, some implementations here are returning a float (or T) rather than an int, which seems wasteful.
  • Works for ints, floats, doubles, unsigned shorts, or any custom types constructible from integer 0 and orderable.
  • Fast! copysign is slow, especially if you need to promote and then narrow again. This is branchless and optimizes excellently
  • Standards-compliant! The bitshift hack is neat, but only works for some bit representations, and doesn't work when you have an unsigned type. It could be provided as a manual specialization when appropriate.
  • Accurate! Simple comparisons with zero can maintain the machine's internal high-precision representation (e.g. 80 bit on x87), and avoid a premature round to zero.

Caveats:

  • It's a template so it might take longer to compile in some circumstances.

  • Apparently some people think use of a new, somewhat esoteric, and very slow standard library function that doesn't even really implement signum is more understandable.

  • The < 0 part of the check triggers GCC's -Wtype-limits warning when instantiated for an unsigned type. You can avoid this by using some overloads:

     template <typename T> inline constexpr
     int signum(T x, std::false_type is_signed) {
         return T(0) < x;
     }
    
     template <typename T> inline constexpr
     int signum(T x, std::true_type is_signed) {
         return (T(0) < x) - (x < T(0));
     }
    
     template <typename T> inline constexpr
     int signum(T x) {
         return signum(x, std::is_signed<T>());
     }
    

    (Which is a good example of the first caveat.)

29
  • 20
    @GMan: GCC only just now (4.5) stopped having cost quadratic to the number of instantiations for template functions, and they are still drastically more expensive to parse and instantiate than manually written functions or the standard C preprocessor. The linker also has to do more work to remove duplicate instantiations. Templates also encourage #includes-in-#includes, which makes dependency calculation take longer and small (often implementation, not interface) changes to force more files to be recompiled.
    – user79758
    Jan 5, 2011 at 22:42
  • 16
    @Joe: Yes, and there's still no noticeable cost. C++ uses templates, that's just something we all have to understand, accept, and get over.
    – GManNickG
    Jan 5, 2011 at 22:54
  • 53
    Wait, what's this "copysign is slow" business...? Using current compilers (g++ 4.6+, clang++ 3.0), std::copysign seems to result in excellent code for me: 4 instructions (inlined), no branching, entirely using the FPU. The recipe given in this answer, by contrast, generates much worse code (many more instructions, including a multiply, moving back and forth between integer unit and FPU)... Jan 23, 2012 at 6:35
  • 14
    @snogglethorpe: If you're calling copysign on an int it promotes to float/double, and must narrow again on return. Your compiler may optimize that promotion out but I can't find anything suggesting that's guaranteed by the standard. Also to implement signum via copysign you need to manually handle the 0 case - please make sure you include that in any performance comparison.
    – user79758
    Jan 23, 2012 at 9:31
  • 60
    The first version is not branchless. Why do people think that a comparison used in an expression will not generate a branch? It will on most architectures. Only processors that have a cmove (or predication) will generate branchless code, but they'll do it also for ternaries or if/else if it is a win. Mar 12, 2012 at 15:41
321

I don't know of a standard function for it. Here's an interesting way to write it though:

(x > 0) - (x < 0)

Here's a more readable way to do it:

if (x > 0) return 1;
if (x < 0) return -1;
return 0;

If you like the ternary operator you can do this:

(x > 0) ? 1 : ((x < 0) ? -1 : 0)
23
  • 9
    Mark Ransom, your expressions give incorrect results for x==0.
    – avakar
    Dec 14, 2009 at 22:45
  • 5
    @Svante: "Each of the operators <, > ... shall yield 1 if the specified relation is true and 0 if it is false" Dec 14, 2009 at 23:15
  • 11
    @Svante: not exactly. A value of 0 is "false"; any other value is "true"; however, the relational and equality operators always return 0 or 1 (see Standard 6.5.8 and 6.5.9). -- the value of the expression a * (x == 42) is either 0 or a.
    – pmg
    Dec 14, 2009 at 23:21
  • 21
    High-Performance Mark, I'm amazed that you missed the C++ tag. This answer is very much valid and doesn't deserve a down-vote. Moreover, I wouldn't use copysign for integral x even if I had it available.
    – avakar
    Dec 15, 2009 at 8:39
  • 6
    Has anyone actually checked what code GCC/G++/any other compiler emits on a real platform? My guess is that the "branchless" version uses two branches instead of one. Bitshifting is probably a lot faster - and more portable in terms of performance. Sep 2, 2011 at 11:29
246

There is a C99 math library function called copysign(), which takes the sign from one argument and the absolute value from the other:

result = copysign(1.0, value) // double
result = copysignf(1.0, value) // float
result = copysignl(1.0, value) // long double

will give you a result of +/- 1.0, depending on the sign of value. Note that floating point zeroes are signed: (+0) will yield +1, and (-0) will yield -1.

9
  • 77
    Upvoted this one, downvoted most popular answer. Left reeling in amazement that SO community seems to prefer a hack to use of a standard library function. May the gods of programming condemn you all to trying to decipher hacks used by clever programmers unfamiliar with language standards. Yeah, I know this is going to cost me a ton of rep on SO, but I'd rather side with comingstorm than the rest of you ... Dec 15, 2009 at 7:42
  • 40
    This is close, but it gives the wrong answer for zero (according to the Wikipedia article in the question at least). Nice suggestion though. +1 anyway.
    – Mark Byers
    Dec 15, 2009 at 8:25
  • 11
    1) C99 is not fully supported everywhere (consider VC++); 2) this is also a C++ question. This is a good answer, but the upvoted one also works, and is more widely applicable. Dec 31, 2009 at 9:17
  • 10
    I wouldn't use copysign() on an AVR microcontroller, it adds an amazing 334 bytes to the program size compared to the "hacks" (if not already using anything else from math.h). Feb 29, 2016 at 12:53
  • 5
    I'm generally for using standard library functions, but this really does not do what was requested precisely because of the note at the end about signed floating-point 0. If your use case really wants sgn(0) to give +1 or -1, then this is ok, but I think that most people looking for a sgn function are going to want that to always give 0 as that is the usual mathematical convention and it matches other languages.
    – Brick
    Mar 14, 2016 at 11:57
100

It seems that most of the answers missed the original question.

Is there a standard sign function (signum, sgn) in C/C++?

Not in the standard library, however there is copysign which can be used almost the same way via copysign(1.0, arg) and there is a true sign function in boost, which might as well be part of the standard.

    #include <boost/math/special_functions/sign.hpp>

    //Returns 1 if x > 0, -1 if x < 0, and 0 if x is zero.
    template <class T>
    inline int sign (const T& z);
3
  • I have been wondering for the past few minutes why the standard library doesn't have sign function. It is just so common -- definitely more commonly used than gamma function which could be found in cmath header.
    – Taozi
    Feb 22, 2016 at 22:33
  • 5
    The explanation I often get for similar questions is "it's easy enough to implement yourself" Which IMO is not a good reason. It completely belies the problems of where standardization, unobvious edge cases, and where to put such a widely used tool.
    – Catskul
    Feb 22, 2016 at 22:35
  • I would not expect to see this marked as the answer because it says to use an external non-standard library. I don't use Boost and cannot use Boost so this is not helpful. Jan 19, 2021 at 2:34
93

Apparently, the answer to the original poster's question is no. There is no standard C++ sgn function.

0
48

Is there a standard sign function (signum, sgn) in C/C++?

Yes, depending on definition.

C99 and later has the signbit() macro in <math.h>

int signbit(real-floating x);
The signbit macro returns a nonzero value if and only if the sign of its argument value is negative. C11 §7.12.3.6


Yet OP wants something a little different.

I want a function that returns -1 for negative numbers and +1 for positive numbers. ... a function working on floats.

#define signbit_p1_or_n1(x)  ((signbit(x) ?  -1 : 1)

Deeper:

OP's question is not specific in the following cases: x = 0.0, -0.0, +NaN, -NaN.

A classic signum() returns +1 on x>0, -1 on x<0 and 0 on x==0.

Many answers have already covered that, but do not address x = -0.0, +NaN, -NaN. Many are geared for an integer point-of-view that usually lacks Not-a-Numbers (NaN) and -0.0.

Typical answers function like signnum_typical() On -0.0, +NaN, -NaN, they return 0.0, 0.0, 0.0.

int signnum_typical(double x) {
  if (x > 0.0) return 1;
  if (x < 0.0) return -1;
  return 0;
}

Instead, I propose this functionality: On -0.0, +NaN, -NaN, it returns -0.0, +NaN, -NaN.

double signnum_c(double x) {
  if (x > 0.0) return 1.0;
  if (x < 0.0) return -1.0;
  return x;
}
1
30

Faster than the above solutions, including the highest rated one:

(x < 0) ? -1 : (x > 0)
6
  • 1
    What type is x? Or are you using a #define?
    – Chance
    Feb 20, 2012 at 18:11
  • 3
    Your type is not faster. It will cause a cache miss quite often. Dec 30, 2012 at 3:19
  • 20
    Cache miss? I'm not sure how. Perhaps you meant branch misprediction?
    – Catskul
    Jun 1, 2013 at 5:10
  • 3
    It seems to me this will result in a warning of confusing integer and boolean types!
    – sergiol
    Sep 2, 2015 at 11:55
  • how this will be fast with the branch?
    – Nick
    Jul 15, 2017 at 19:55
17

There's a way to do it without branching, but it's not very pretty.

sign = -(int)((unsigned int)((int)v) >> (sizeof(int) * CHAR_BIT - 1));

http://graphics.stanford.edu/~seander/bithacks.html

Lots of other interesting, overly-clever stuff on that page, too...

2
  • 1
    If I read the link correctly that only returns -1 or 0. If you want -1, 0, or +1 then it's sign = (v != 0) | -(int)((unsigned int)((int)v) >> (sizeof(int) * CHAR_BIT - 1)); or sign = (v > 0) - (v < 0);.
    – Z boson
    Apr 21, 2015 at 7:58
  • 2
    this implies that v is an integer type not wider than int
    – phuclv
    Aug 10, 2019 at 15:36
13

If all you want is to test the sign, use signbit (returns true if its argument has a negative sign). Not sure why you would particularly want -1 or +1 returned; copysign is more convenient for that, but it sounds like it will return +1 for negative zero on some platforms with only partial support for negative zero, where signbit presumably would return true.

3
  • 9
    There's many mathematical applications in which the sign(x) is necessary. Otherwise I'd just do if (x < 0).
    – Chance
    Mar 9, 2012 at 16:42
  • It looks like this function takes a float as parameter... I wouldn't want to pay the conversion price just to get the sign bit... Nov 19, 2022 at 16:14
  • correct. but this question is about floats
    – ysth
    Nov 20, 2022 at 2:24
10

In general, there is no standard signum function in C/C++, and the lack of such a fundamental function tells you a lot about these languages.

Apart from that, I believe both majority viewpoints about the right approach to define such a function are in a way correct, and the "controversy" about it is actually a non-argument once you take into account two important caveats:

  • A signum function should always return the type of its operand, similarly to an abs() function, because signum is usually used for multiplication with an absolute value after the latter has been processed somehow. Therefore, the major use case of signum is not comparisons but arithmetic, and the latter shouldn't involve any expensive integer-to/from-floating-point conversions.

  • Floating point types do not feature a single exact zero value: +0.0 can be interpreted as "infinitesimally above zero", and -0.0 as "infinitesimally below zero". That's the reason why comparisons involving zero must internally check against both values, and an expression like x == 0.0 can be dangerous.

Regarding C, I think the best way forward with integral types is indeed to use the (x > 0) - (x < 0) expression, as it should be translated in a branch-free fashion, and requires only three basic operations. Best define inline functions that enforce a return type matching the argument type, and add a C11 define _Generic to map these functions to a common name.

With floating point values, I think inline functions based on C11 copysignf(1.0f, x), copysign(1.0, x), and copysignl(1.0l, x) are the way to go, simply because they're also highly likely to be branch-free, and additionally do not require casting the result from integer back into a floating point value. You should probably comment prominently that your floating point implementations of signum will not return zero because of the peculiarities of floating point zero values, processing time considerations, and also because it is often very useful in floating point arithmetic to receive the correct -1/+1 sign, even for zero values.

5

My copy of C in a Nutshell reveals the existence of a standard function called copysign which might be useful. It looks as if copysign(1.0, -2.0) would return -1.0 and copysign(1.0, 2.0) would return +1.0.

Pretty close huh?

3
  • Not standard, but may be widely available. Microsoft's starts with an underscore, which is the convention they use for non-standard extensions. Not the best choice when you're working with integers, though. Dec 14, 2009 at 23:33
  • 6
    copysign is both in the ISO C (C99) and POSIX standards. See opengroup.org/onlinepubs/000095399/functions/copysign.html
    – lhf
    Dec 15, 2009 at 1:19
  • 3
    What lhf said. Visual Studio is not a reference for the C standard. Dec 15, 2009 at 1:42
4

The question is old but there is now this kind of desired function. I added a wrapper with not, left shift and dec.

You can use a wrapper function based on signbit from C99 in order to get the exact desired behavior (see code further below).

Returns whether the sign of x is negative.
This can be also applied to infinites, NaNs and zeroes (if zero is unsigned, it is considered positive

#include <math.h>

int signValue(float a) {
    return ((!signbit(a)) << 1) - 1;
}

NB: I use operand not ("!") because the return value of signbit is not specified to be 1 (even though the examples let us think it would always be this way) but true for a negative number:

Return value
A non-zero value (true) if the sign of x is negative; and zero (false) otherwise.

Then I multiply by two with left shift (" << 1") which will give us 2 for a positive number and 0 for a negative one and finally decrement by 1 to obtain 1 and -1 for respectively positive and negative numbers as requested by OP.

2
4

The accepted answer with the overload below does indeed not trigger -Wtype-limits. But it does trigger unused argument warnings (on the is_signed variable). To avoid these the second argument should not be named like so:

template <typename T> inline constexpr
  int signum(T x, std::false_type) {
  return T(0) < x;
}

template <typename T> inline constexpr
  int signum(T x, std::true_type) {
  return (T(0) < x) - (x < T(0));
}

template <typename T> inline constexpr
  int signum(T x) {
  return signum(x, std::is_signed<T>());
}

For C++11 and higher an alternative could be.

template <typename T>
typename std::enable_if<std::is_unsigned<T>::value, int>::type
inline constexpr signum(T const x) {
    return T(0) < x;  
}

template <typename T>
typename std::enable_if<std::is_signed<T>::value, int>::type
inline constexpr signum(T const x) {
    return (T(0) < x) - (x < T(0));  
}

For me it does not trigger any warnings on GCC 5.3.1.

2
  • To avoid the -Wunused-parameter warning just use unnamed parameters. Jan 24, 2019 at 12:47
  • That is actually very true. I missed that. However, I like the C++11 alternative more either way. Jan 29, 2019 at 19:04
2

No, it doesn't exist in c++, like in matlab. I use a macro in my programs for this.

#define sign(a) ( ( (a) < 0 )  ?  -1   : ( (a) > 0 ) )
4
  • 6
    One should prefer templates over macros in C++.
    – Ruslan
    Jan 27, 2016 at 11:06
  • In C, there is no template ...... helloacm.com/how-to-implement-the-sgn-function-in-c
    – justyy
    Oct 20, 2016 at 10:55
  • I thought this was a good answer then I looked at my own code and found this: #define sign(x) (((x) > 0) - ((x) < 0)) which is good too. Feb 25, 2018 at 12:12
  • 1
    an inline function is better than a macro in C, and in C++ template is better
    – phuclv
    Apr 10, 2018 at 16:05
2

Bit off-topic, but I use this:

template<typename T>
constexpr int sgn(const T &a, const T &b) noexcept{
    return (a > b) - (a < b);
}

template<typename T>
constexpr int sgn(const T &a) noexcept{
    return sgn(a, T(0));
}

and I found first function - the one with two arguments, to be much more useful from "standard" sgn(), because it is most often used in code like this:

int comp(unsigned a, unsigned b){
   return sgn( int(a) - int(b) );
}

vs.

int comp(unsigned a, unsigned b){
   return sgn(a, b);
}

there is no cast for unsigned types and no additional minus.

in fact i have this piece of code using sgn()

template <class T>
int comp(const T &a, const T &b){
    log__("all");
    if (a < b)
        return -1;

    if (a > b)
        return +1;

    return 0;
}

inline int comp(int const a, int const b){
    log__("int");
    return a - b;
}

inline int comp(long int const a, long int const b){
    log__("long");
    return sgn(a, b);
}
2

You can use boost::math::sign() method from boost/math/special_functions/sign.hpp if boost is available.

2
1

Here's a branching-friendly implementation:

inline int signum(const double x) {
    if(x == 0) return 0;
    return (1 - (static_cast<int>((*reinterpret_cast<const uint64_t*>(&x)) >> 63) << 1));
}

Unless your data has zeros as half of the numbers, here the branch predictor will choose one of the branches as the most common. Both branches only involve simple operations.

Alternatively, on some compilers and CPU architectures a completely branchless version may be faster:

inline int signum(const double x) {
    return (x != 0) * 
        (1 - (static_cast<int>((*reinterpret_cast<const uint64_t*>(&x)) >> 63) << 1));
}

This works for IEEE 754 double-precision binary floating-point format: binary64 .

0

While the integer solution in the accepted answer is quite elegant it bothered me that it wouldn't be able to return NAN for double types, so I modified it slightly.

template <typename T> double sgn(T val) {
    return double((T(0) < val) - (val < T(0)))/(val == val);
}

Note that returning a floating point NAN as opposed to a hard coded NAN causes the sign bit to be set in some implementations, so the output for val = -NAN and val = NAN are going to be identical no matter what (if you prefer a "nan" output over a -nan you can put an abs(val) before the return...)

-1
int sign(float n)
{     
  union { float f; std::uint32_t i; } u { n };
  return 1 - ((u.i >> 31) << 1);
}

This function assumes:

  • binary32 representation of floating point numbers
  • a compiler that make an exception about the strict aliasing rule when using a named union
1
  • 3
    There are still some bad assumptions here. For example I don't believe the endianness of the float is guaranteed to be the endianness of the integer. Your check also fails on any architectures using ILP64. Really, you're just reimplementing copysign; if you're using static_assert you've got C++11, and might as well really use copysign.
    – user79758
    Mar 15, 2012 at 18:24
-3
double signof(double a) { return (a == 0) ? 0 : (a<0 ? -1 : 1); }
-3

Why use ternary operators and if-else when you can simply do this

#define sgn(x) x==0 ? 0 : x/abs(x)
4
  • 3
    Your definition uses a ternary operator as well.
    – Martin R
    Jun 30, 2018 at 14:26
  • Yes Definitely, but it just uses one ternary operator to separate zero and non-zero numbers. Others' versions include nested ternary ops to separate positive, negative and zero.
    – Jagreet
    Jun 30, 2018 at 16:15
  • Using an integer division is very inefficient and abs() is only for integers. Mar 19, 2019 at 12:18
  • Undefined behavior possible when x == INT_MIN. Feb 3, 2020 at 22:52

Not the answer you're looking for? Browse other questions tagged or ask your own question.