AskOverflow.Dev

AskOverflow.Dev Logo AskOverflow.Dev Logo

AskOverflow.Dev Navigation

  • 主页
  • 系统&网络
  • Ubuntu
  • Unix
  • DBA
  • Computer
  • Coding
  • LangChain

Mobile menu

Close
  • 主页
  • 系统&网络
    • 最新
    • 热门
    • 标签
  • Ubuntu
    • 最新
    • 热门
    • 标签
  • Unix
    • 最新
    • 标签
  • DBA
    • 最新
    • 标签
  • Computer
    • 最新
    • 标签
  • Coding
    • 最新
    • 标签
主页 / unix / 问题 / 564397
Accepted
HerpDerpington
HerpDerpington
Asked: 2020-01-28 07:54:10 +0800 CST2020-01-28 07:54:10 +0800 CST 2020-01-28 07:54:10 +0800 CST

为什么 gawk(有时?)认为 2.0e-318 > 2.0?

  • 772

我正在尝试使用 gawk 找到一列数据的最大值:

gawk 'BEGIN{max=0} {if($1>0+max) max=$1} END {print max}' dataset.dat

dataset.dat看起来像这样:

2.0
2.0e-318

该命令的输出是

2.0e-318

明显小于 2。

我的错误在哪里?

编辑

有趣的是,如果交换输入文件的行,输出变为

2.0

编辑 2

我的 gawk 版本是 GNU Awk 4.2.1,API:2.0(GNU MPFR 4.0.2,GNU MP 6.1.2)。

ubuntu gawk
  • 2 2 个回答
  • 93 Views

2 个回答

  • Voted
  1. user232326
    2020-01-28T13:39:48+08:002020-01-28T13:39:48+08:00

    2e-318在 awk 中处理如此小的数字 ( ) 有几个问题。

    • 首先,输入需要在使用前转换为数字。这通常是通过添加 0 来完成的。因此,您需要以下内容:

      val=0+$1
      
    • 其次,普通双精度浮点数(53 位尾数和 11 位指数)11 位宽度的指数允许表示 10e-308 和 10e308 之间的数字,因此,普通浮点数将无法表示此类数字。

      $ echo '1e-307 1e-308' | awk '{print $1,$1+0,$2,$2+0}'
      1e-307 1e-307 1e-308 0
      

      默认 GNU awk 将不接受以下(正常)值1e-308。

    • 第三,awk(CNVFMT 和 OFMT)的默认转换格式设置为"%.6g". 超过 6 位有效数字的数字将被截断。要获得更重要的数字:要求他们。就像%.15g15 一样(对于 53 位尾数,不要要求超过 17,它可能会撒谎)。

    • 第四,最好将第一个值设置为第max一个输入。如果输入的最大值为负,将最大值设置为 0 将失败。


    如果您使用的是 GNU awk 并且它已经以任意精度编译,您可以使用:

    $ printf '%s\n' 2e-318 2e-317 2e-307 2e-308 2e-319 | 
        awk -M -v PREC=100     'BEGIN{OFMT="%.15g"};
            {val=0+$1};
            NR==1{max=val};
            {print($1,val,max)};
            val>max{max=val}
            END{print max}'
    
    2e-318 2e-318 2e-318
    2e-317 2e-317 2e-318
    2e-307 2e-307 2e-317
    2e-308 2e-308 2e-307
    2e-319 2e-319 2e-307
    2e-307 
    

    或简化为您的用例:

    awk -M -v PREC=100 '
        BEGIN{OFMT="%.15g"};    # allow more than 6 figures
        {val=0+$1};             # convert input to a (float) number.
        NR==1{max=val};         # On the first line, set the max value.
        val>max{max=val}        # On every entry keep track of the max.
        END{print max}          # At the end, print the max.
        '  file                 # file with input (one per line).
    
    • 3
  2. Best Answer
    Paul_Pedant
    2020-01-28T08:09:20+08:002020-01-28T08:09:20+08:00

    0+ 需要作为每个 $1 的前缀以强制进行数字转换。max 不需要 0+ - 它在存储时已经转换为数字。

    Paul--) AWK='
    > BEGIN { max = 0; }
    > 0+$1 > max { max = 0 + $1; }
    > END { print max; }
    > '
    Paul--) awk "${AWK}" <<[][]
    > 2.0
    > 2.0e-318
    > [][]
    2
    Paul--) awk "${AWK}" <<[][]
    > 2.0e-318
    > 2.0
    > [][]
    2
    
    • 1

相关问题

  • 使 Apple 键盘设置在 ../hid_apple/parameters/ 中持久化

  • 用文件中的点绘制二维图形

  • 如何验证我使用 9.9.9.9 作为 DNS?

  • 并行安装多个 linux 发行版

  • 自编译 OpenVPN 不会从 systemd 启动

Sidebar

Stats

  • 问题 205573
  • 回答 270741
  • 最佳答案 135370
  • 用户 68524
  • 热门
  • 回答
  • Marko Smith

    模块 i915 可能缺少固件 /lib/firmware/i915/*

    • 3 个回答
  • Marko Smith

    无法获取 jessie backports 存储库

    • 4 个回答
  • Marko Smith

    如何将 GPG 私钥和公钥导出到文件

    • 4 个回答
  • Marko Smith

    我们如何运行存储在变量中的命令?

    • 5 个回答
  • Marko Smith

    如何配置 systemd-resolved 和 systemd-networkd 以使用本地 DNS 服务器来解析本地域和远程 DNS 服务器来解析远程域?

    • 3 个回答
  • Marko Smith

    dist-upgrade 后 Kali Linux 中的 apt-get update 错误 [重复]

    • 2 个回答
  • Marko Smith

    如何从 systemctl 服务日志中查看最新的 x 行

    • 5 个回答
  • Marko Smith

    Nano - 跳转到文件末尾

    • 8 个回答
  • Marko Smith

    grub 错误:你需要先加载内核

    • 4 个回答
  • Marko Smith

    如何下载软件包而不是使用 apt-get 命令安装它?

    • 7 个回答
  • Martin Hope
    user12345 无法获取 jessie backports 存储库 2019-03-27 04:39:28 +0800 CST
  • Martin Hope
    Carl 为什么大多数 systemd 示例都包含 WantedBy=multi-user.target? 2019-03-15 11:49:25 +0800 CST
  • Martin Hope
    rocky 如何将 GPG 私钥和公钥导出到文件 2018-11-16 05:36:15 +0800 CST
  • Martin Hope
    Evan Carroll systemctl 状态显示:“状态:降级” 2018-06-03 18:48:17 +0800 CST
  • Martin Hope
    Tim 我们如何运行存储在变量中的命令? 2018-05-21 04:46:29 +0800 CST
  • Martin Hope
    Ankur S 为什么 /dev/null 是一个文件?为什么它的功能不作为一个简单的程序来实现? 2018-04-17 07:28:04 +0800 CST
  • Martin Hope
    user3191334 如何从 systemctl 服务日志中查看最新的 x 行 2018-02-07 00:14:16 +0800 CST
  • Martin Hope
    Marko Pacak Nano - 跳转到文件末尾 2018-02-01 01:53:03 +0800 CST
  • Martin Hope
    Kidburla 为什么真假这么大? 2018-01-26 12:14:47 +0800 CST
  • Martin Hope
    Christos Baziotis 在一个巨大的(70GB)、一行、文本文件中替换字符串 2017-12-30 06:58:33 +0800 CST

热门标签

linux bash debian shell-script text-processing ubuntu centos shell awk ssh

Explore

  • 主页
  • 问题
    • 最新
    • 热门
  • 标签
  • 帮助

Footer

AskOverflow.Dev

关于我们

  • 关于我们
  • 联系我们

Legal Stuff

  • Privacy Policy

Language

  • Pt
  • Server
  • Unix

© 2023 AskOverflow.DEV All Rights Reserve