为何std::stof()对合法非规格化float抛出std::out_of_range异常？

问题现象 #

使用g++ 12.2（运行于amd64 Debian 12）编译执行以下代码时，std::stof()会抛出std::out_of_range异常：

#include <iostream>

int main()
{
  const std::string s("5.87747175e-39");
  float f = std::stof(s);
  std::cout << s << " - " << f << std::endl;
  return 0;
}

尽管5.87747175e-39是合法的IEEE-754单精度非规格化浮点数，且Python可正常转换该值，但所有非规格化float都会触发此问题，而FLT_MIN对应的最小规格化float则能正常转换。

原因 #

这是由于g的标准库实现中，std::stof()的范围检查逻辑将非规格化数判定为超出float的有效范围。根据C标准，std::stof()在输入值小于FLT_MIN（最小规格化正数）时，允许实现选择抛出std::out_of_range异常，而非强制要求支持非规格化数的转换。

解决方法 #

方法1：通过long double中转转换 #

利用std::stold()先将字符串转换为范围更大的long double，再强制转换为float，这样可以正确保留非规格化数的表示：

#include <iostream>
#include <string>

int main()
{
  const std::string s("5.87747175e-39");
  long double ld = std::stold(s);
  float f = static_cast<float>(ld);
  std::cout << s << " - " << f << std::endl;
  return 0;
}

方法2：手动构造非规格化float #

直接解析字符串的尾数和指数，通过操作float的二进制内存布局来构造对应的非规格化值。这种方式需要熟悉IEEE-754单精度浮点数的格式，示例简化实现：

#include <iostream>
#include <string>
#include <cstdint>
#include <cmath>

float string_to_float_non_normal(const std::string& s) {
    size_t exp_pos = s.find('e');
    std::string mantissa_str = s.substr(0, exp_pos);
    int exponent = std::stoi(s.substr(exp_pos + 1));

    double mantissa = std::stod(mantissa_str);
    // 转换为IEEE-754单精度非规格化数的二进制表示
    float result;
    uint32_t* ptr = reinterpret_cast<uint32_t*>(&result);
    
    // 非规格化数指数域全为0，尾数直接存储有效数字（需补充完整转换逻辑）
    *ptr = static_cast<uint32_t>(std::ldexp(mantissa, exponent + 126)) & 0x7FFFFF;
    return result;
}

int main()
{
  const std::string s("5.87747175e-39");
  float f = string_to_float_non_normal(s);
  std::cout << s << " - " << f << std::endl;
  return 0;
}

注意：手动构造需严格遵循IEEE-754规范，处理各类边界场景，适合对浮点数格式有深入了解的开发场景。

内容的提问来源于stack exchange，提问作者Paul Grinberg