std::nextafter, std::nextafterf, std::nextafterl, std::nexttoward, std::nexttowardf, std::nexttowardl

#include <cfenv>
#include <cfloat>
#include <cmath>
#include <concepts>
#include <iomanip>
#include <iostream>
 
int main()
{
    float from1 = 0, to1 = std::nextafter(from1, 1.f);
    std::cout << "El siguiente float representable después de " << std::setprecision(20) << from1
              << " es " << to1
              << std::hexfloat << " (" << to1 << ")\n" << std::defaultfloat;
 
    float from2 = 1, to2 = std::nextafter(from2, 2.f);
    std::cout << "El siguiente float representable después de " << from2 << " es " << to2
              << std::hexfloat << " (" << to2 << ")\n" << std::defaultfloat;
 
    double from3 = std::nextafter(0.1, 0), to3 = 0.1;
    std::cout << "El número 0.1 se encuentra entre dos double válidos:\n"
              << std::setprecision(56) << "    " << from3
              << std::hexfloat << " (" << from3 << ')' << std::defaultfloat
              << "\ny   " << to3 << std::hexfloat << "  (" << to3 << ")\n"
              << std::defaultfloat << std::setprecision(20);
 
    std::cout << "\nDiferencia entre nextafter y nexttoward:\n";
    long double dir = std::nextafter(from1, 1.0L); // primer long double subnormal
    float x = std::nextafter(from1, dir); // primero convierte dir to float, dando 0
    std::cout << "Con nextafter, el próximo float después de " << from1 << " es " << x << '\n';
    x = std::nexttoward(from1, dir);
    std::cout << "Con nexttoward, el próximo float después de " << from1 << " es " << x << '\n';
 
    std::cout << "\nValores especiales:\n";
    {
        // #pragma STDC FENV_ACCESS ON
        std::feclearexcept(FE_ALL_EXCEPT);
        double from4 = DBL_MAX, to4 = std::nextafter(from4, INFINITY);
        std::cout << "El siguiente double representable después de " << std::setprecision(6)
                  << from4 << std::hexfloat << " (" << from4 << ')'
                  << std::defaultfloat << " es " << to4
                  << std::hexfloat << " (" << to4 << ")\n" << std::defaultfloat;
        if(std::fetestexcept(FE_OVERFLOW)) std::cout << "   generó FE_OVERFLOW\n";
        if(std::fetestexcept(FE_INEXACT)) std::cout << "   generó FE_INEXACT\n";
    } // fin del bloque FENV_ACCESS 
 
    float from5 = 0.0, to5 = std::nextafter(from5, -0.0);
    std::cout << "std::nextafter(+0.0, -0.0) da " << std::fixed << to5 << '\n';
 
    auto precision_loss_demo = []<std::floating_point Fp>(const auto rem, const Fp start) {
        std::cout << rem;
        for (Fp from = start, to, Δ;
            (Δ = (to = std::nextafter(from, +INFINITY)) - from) < Fp(10.0);
            from *= Fp(10.0))
            std::cout << "nextafter(" << std::scientific << std::setprecision(0) << from 
                      << ", INF) da " << std::fixed << std::setprecision(6) << to
                      << "; Δ = " << Δ << '\n';
    };
 
    precision_loss_demo("\nDemo de pérdida de precisión para float:\n", 10.0f);
    precision_loss_demo("\nDemo de pérdida de precisión para double:\n", 10.0e9);
    precision_loss_demo("\nDemo de pérdida de precisión para long double:\n", 10.0e17L);
}