17

Given the 32 bits that represent an IEEE 754 floating-point number, how can the number be converted to an integer, using integer or bit operations on the representation (rather than using a machine instruction or compiler operation to convert)?

I have the following function but it fails in some cases:

Input: int x (contains 32 bit single precision number in IEEE 754 format)

  if(x == 0) return x;

  unsigned int signBit = 0;
  unsigned int absX = (unsigned int)x;
  if (x < 0)
  {
      signBit = 0x80000000u;
      absX = (unsigned int)-x;
  }

  unsigned int exponent = 158;
  while ((absX & 0x80000000) == 0)
  {
      exponent--;
      absX <<= 1;
  }

  unsigned int mantissa = absX >> 8;

  unsigned int result = signBit | (exponent << 23) | (mantissa & 0x7fffff);
  printf("\nfor x: %x, result: %x",x,result);
  return result;
DaveyDaveDave
  • 8,579
  • 11
  • 62
  • 73
Anonymous
  • 179
  • 1
  • 1
  • 5
  • This don't cast a float into an int. It just copy bitwise their machine representation, without e.g. converting `2.03e1` to `20` [by rounding] as the `(int)2.03e1` cast will. – Basile Starynkevitch Sep 09 '12 at 20:59
  • 2
    You *want* do do it bitwise? Well, that's how you do it bitwise - it just reinterprets the bytes. No steps, really. – Ry- Sep 09 '12 at 21:00
  • But 0x7eff8965 = 1325268755 (after casting). If you use the HEX in IEEE 754 Calc, you get 1.6983327e+38 and HEX to decimal gives: 2130676069 - none of them give the correct result of 1325268755. – Anonymous Sep 09 '12 at 21:04
  • 2
    This code has undefined behavior in C. See section 6.5 in the standard. – Paul Hankin Sep 09 '12 at 21:30
  • Is your question this: Given the 32 bits that represent a float x, how can the conversion `(int) x` be implemented, using integer/bit operations on the representation (rather than using a machine instruction to convert floating-point to integer)? – Eric Postpischil Sep 09 '12 at 21:40
  • @EricPostpischil - Yes! exactly. – Anonymous Sep 09 '12 at 21:55
  • There's another related question by Anon: [Negate Floating Number in C](http://stackoverflow.com/questions/12336675/negate-floating-number-in-c-fails-in-some-cases), also about bitwise manipulation of IEEE 754 values. There was a second related question in the last 24 hours or so: [How to manually (bitwise) perform `(float)x`](http://stackoverflow.com/questions/12336314/how-to-manually-bitwise-perform-floatx). – Jonathan Leffler Sep 09 '12 at 22:19
  • Indeed, http://stackoverflow.com/questions/12336314/how-to-manually-bitwise-perform-floatx Having the same problem, it doesn't want to round correctly... very frustrating –  Sep 10 '12 at 01:08
  • @Silver - yes! I still have to work on float_times_four. That is time consuming too! – Anonymous Sep 10 '12 at 01:10
  • For float_times_four, you want to separate it into a bunch of cases (is NaN, is zero, is infinity, is normal, is denormalized (that last one was the part that took me a while)) –  Sep 10 '12 at 01:24
  • I have already done that. If less than 0x0071FFFF then just return uf*4, else just add to mantissa. But I am not sure what to do when both exponent and mantissa have to be changed. Also, are you converting the number and doing multiplication, or just manipulating bits in the IEEE form? – Anonymous Sep 10 '12 at 01:27
  • I asked whether you were trying to convert a float (given its representation) to an int, and you answered yes. But your code looks like you are trying to convert an int to a float. Which is it? (The latter is addressed [here](http://stackoverflow.com/questions/12336314/how-to-manually-bitwise-perform-floatx).) – Eric Postpischil Sep 10 '12 at 02:06
  • 1
    BTW, the code you posted converts a 32-bit signed int to its 32-bit IEEE 754 single-precision with rounding toward zero. I know because I wrote it yesterday. – Gabe Sep 10 '12 at 04:45
  • duplicates: [How to manually (bitwise) perform (float)x?](https://stackoverflow.com/q/12336314/995714), [Converting Int to Float or Float to Int using Bitwise operations (software floating point)](https://stackoverflow.com/q/20302904/995714), [How to convert an unsigned int to a float?](https://stackoverflow.com/q/19529356/995714) – phuclv May 15 '19 at 04:14

7 Answers7

23

C has the "union" to handle this type of view of data:

typedef union {
  int i;
  float f;
 } u;
 u u1;
 u1.f = 45.6789;
 /* now u1.i refers to the int version of the float */
 printf("%d",u1.i);
Carl
  • 239
  • 2
  • 2
18

&x gives the address of x so has float* type.

(int*)&x cast that pointer to a pointer to int ie to a int* thing.

*(int*)&x dereference that pointer into an int value. It won't do what you believe on machines where int and float have different sizes.

And there could be endianness issues.

This solution was used in the fast inverse square root algorithm.

mihaib
  • 15
  • 4
Basile Starynkevitch
  • 216,767
  • 17
  • 275
  • 509
  • So you are saying that the code just gets the location of x and prints it out? In that case, the value would change on each run. – Anonymous Sep 09 '12 at 21:09
  • 1
    No it gives the integer contained at the location of the float, so, when `sizeof(int) == sizeof[float]` it gives the `int` of the same machine bit representation as your `x` ; nothing is printed unless you call a printing routine like `printf` (which is not in your question) – Basile Starynkevitch Sep 09 '12 at 21:18
  • Ok, so it gives the value stored at the location in memory and casts it to an int type. How can I do this without casting? – Anonymous Sep 09 '12 at 21:27
  • @BasileStarynkevitch: what would the problem with endianness be? If you are just looking to pick out the bits of a float, I don't think it would matter if ints are stored big- or little-endian. – Björn Lindqvist Jan 21 '18 at 16:00
  • Endianness would be a problem if you were converting a float to an unsigned int, where you are using the bits as flags and the sending function/program/device can only send floats. – Mark Walsh Jan 15 '19 at 20:37
7

(Somebody should double-check this answer, especially border cases and the rounding of negative values. Also, I wrote it for round-to-nearest. To reproduce C’s conversion, this should be changed to round-toward-zero.)

Essentially, the process is:

Separate the 32 bits into one sign bit (s), eight exponent bits (e), and 23 significand bits (f). We will treat these as twos-complement integers.

If e is 255, the floating-point object is either infinity (if f is zero) or a NaN (otherwise). In this case, the conversion cannot be performed, and an error should be reported.

Otherwise, if e is not zero, add 224 to f. (If e is not zero, the significand implicitly has a 1 bit at its front. Adding 224 makes that bit explicit in f.)

Subtract 127 from e. (This converts the exponent from its biased/encoded form to the actual exponent. If we were doing a general conversion to any value, we would have to handle the special case when e is zero: Subtract 126 instead of 127. But, since we are only converting to an integer result, we can neglect this case, as long as the integer result is zero for these tiny input numbers.)

If s is 0 (the sign is positive) and e is 31 or more, then the value overflows a signed 32-bit integer (it is 231 or larger). The conversion cannot be performed, and an error should be reported.

If s is 1 (the sign is negative) and e is more than 31, then the value overflows a signed 32-bit integer (it is less than or equal to -232). If s is one, e is 32, and f is greater than 224 (any of the original significand bits were set), then the value overflows a signed 32-bit integer (it is less than -231; if the original f were zero, it would be exactly -231, which does not overflow). In any of these cases, the conversion cannot be performed, and an error should be reported.

Now we have an s, an e, and an f for a value which does not overflow, so we can prepare the final value.

If s is 1, set f to -f.

The exponent value is for a significand between 1 (inclusive) and 2 (exclusive), but our significand starts with a bit at 224. So we have to adjust for that. If e is 24, our significand is correct, and we are done, so return f as the result. If e is greater than 24 or less than 24, we have to shift the significand appropriately. Also, if we are going to shift f right, we may have to round it, to get a result rounded to the nearest integer.

If e is greater than 24, shift f left e-24 bits. Return f as the result.

If e is less than -1, the floating-point number is between -½ and ½, exclusive. Return 0 as the result.

Otherwise, we will shift f right 24-e bits. However, we will first save the bits we need for rounding. Set r to the result of casting f to an unsigned 32-bit integer and shifting it left by 32-(24-e) bits (equivalently, left by 8+e bits). This takes the bits that will be shifted out of f (below) and “left adjusts” them in the 32 bits, so we have a fixed position where they start.

Shift f right 24-e bits.

If r is less than 231, do nothing (this is rounding down; the shift truncated bits). If r is greater than 231, add one to f (this is rounding up). If r equals 231, add the low bit of f to f. (If f is odd, add one to f. Of the two equally near values, this rounds to the even value.) Return f.

Eric Postpischil
  • 168,892
  • 12
  • 149
  • 276
  • Thanks for the explaination. I wrote the function but it fails in some cases. – Anonymous Sep 10 '12 at 00:59
  • "Subtract 127 from e." happens when `e > 0`. Else "Subtract 126 from 0." – chux - Reinstate Monica Dec 02 '13 at 01:07
  • @chux: Yes, one would need to adjust when converting a floating-point encoding to a number in general. This question asks about the special case of converting a floating-point encoding to an integer. In that case, we can neglect proper handling of tiny values, since they will produce zero in the end. – Eric Postpischil Dec 02 '13 at 14:18
6
// With the proviso that your compiler implementation uses
// the same number of bytes for an int as for a float:
// example float
float f = 1.234f;
// get address of float, cast as pointer to int, reference
int i = *((int *)&f);
// get address of int, cast as pointer to float, reference
float g = *((float *)&i);
printf("%f %f %08x\n",f,g,i);
Dino Dini
  • 207
  • 2
  • 4
5
float x = 43.133;
int y;

assert (sizeof x == sizeof y);
memcpy (&y, &x, sizeof x);
...
wildplasser
  • 41,380
  • 7
  • 58
  • 102
  • memcpy did not work. int x (contains 32 bit float) is the input, then int result; memcpy(&result, &x, 4) does not work. (4 is ok as it will only run on 32bit machines) – Anonymous Sep 09 '12 at 21:37
  • Maybe your assert (or your sizeof) is broke? BTW:Oops, I should have used x instead of f. BRB. – wildplasser Sep 09 '12 at 21:53
1

You can cast the float using a reference. A cast like this should never generate any code.

C++

float f = 1.0f;
int i = (int &)f;
printf("Float %f is 0x%08x\n", f, i);

Output:

Float 1.000000 is 0x3f800000

If you want c++ style cast use a reinterpret_cast, like this.

int i = reinterpret_cast<int &>(f);

It does not work with expressions, you have to store it in a variable.

    int i_times_two;
    float f_times_two = f * 2.0f;
    i_times_two = (int &)f_times_two;

    i_times_two = (int &)(f * 2.0f);
main.cpp:25:13: error: C-style cast from rvalue to reference type 'int &'
  • Please add some informtion about how your code works – Koopakiller Aug 11 '15 at 21:01
  • Your second example will work if you use a rvalue reference, replace ``(int&)`` with ``(int&&)``. This is required as the expression returns an rvalue reference which lvalue references cannot bind to. I assume you could also use ``(const int &)`` to bind to both. – nitronoid Aug 19 '17 at 10:34
0

You cannot (meaningfully) convert a floating point number into an 'integer' (signed int or int) in this way.

It may end up having the integer type, but it's actually just an index into the encoding space of IEEE754, not a meaningful value in itself.

You might argue that an unsigned int serves dual purpose as a bit pattern and an integer value, but int does not.


Also there are platform issues with bit manipulation of signed ints.

Alex Brown
  • 40,336
  • 10
  • 90
  • 107
  • 5
    Meaningful use: receiving two `int16_t`s on a bus, that actually represent a `float32`. Reinterpret the two `int16_t` as a float. – Gauthier Mar 30 '15 at 10:46