nexttoward(), nexttowardf(), nexttowardl()

Compute the next representable floating-point number


#include <math.h>

double nexttoward ( double x,
                    long double y);

float nexttowardf ( float x,
                    long double y );

long double nexttowardl ( long double x,
                          long double y);


The number that you want the next number after.
A number that specifies the direction you want to go; see below.



Use the -l m option to qcc to link against this library.


The nexttoward(), nexttowardf(), and nexttowardl() functions compute the next representable double-precision floating-point value following x in the direction of y. These functions are equivalent to nextafter(), nextafterf(), and nextafterl(), respectively, except that:

To check for error situations, use feclearexcept() and fetestexcept(). For example:


The next machine floating-point number of x in the direction towards y, converted to the type that the function returns.

If These functions return: Errors:
x equals y y
x is finite, and the correct value would overflow Inf FE_OVERFLOW
x or y is NaN NaN
x != y, and the correct value is subnormal, zero, or underflows The correct value (if representable), or 0.0 FE_UNDERFLOW

These functions raise FE_INEXACT if the FPU reports that the result can't be exactly represented as a floating-point number.


#include <stdio.h>
#include <inttypes.h>
#include <math.h>
#include <fenv.h>
#include <stdlib.h>

void dump_to_hex(double d) {
        printf("0x%08x %08x \n",
           (uint32_t)(*((uint64_t*)&d) >> 32),

int main(int argc, char** argv) 
    double a, b, c;
    int except_flags;

    a = 0;

    b = nexttoward(a, -1);
    except_flags = fetestexcept(FE_ALL_EXCEPT);
    if(except_flags) {
        /* An error occurred; handle it appropriately. */

    c = nexttoward(a, 1);
    except_flags = fetestexcept(FE_ALL_EXCEPT);
    if(except_flags) {
        /* An error occurred; handle it appropriately. */

    printf("Next possible value before %f is %f \n", a, b);
    printf("-->"); dump_to_hex(a);
    printf("-->"); dump_to_hex(b);
    printf("Next possible value after %f is %f \n", a, c);
    printf("-->"); dump_to_hex(a);
    printf("-->"); dump_to_hex(c);

    return EXIT_SUCCESS;

produces the output:

Next possible value before 0.000000 is 0.000000
-->0x00000000 00000000
-->0x80000000 00000001
Next possible value after 0.000000 is 0.000000
-->0x00000000 00000000


ANSI, POSIX 1003.1

Cancellation point No
Interrupt handler No
Signal handler No
Thread Yes