0.25的倍数可以精确地表示为双精度吗？

问题描述

我有以下用于查找四分位数的代码：

#include <stdio.h>
#include <stdlib.h>

typedef struct {
    double qrt[3];
    double *value;
    int count;
} t_data;

static void set_qrt(t_data *data, int qrt)
{
    int n, e;
    double d;

    d = qrt * 0.25 * data->count + 0.5;
    n = (int)d;
    e = n != d;
    data->qrt[qrt - 1] = data->value[n - 1];
    if (e) {
        data->qrt[qrt - 1] += data->value[n];
        data->qrt[qrt - 1] *= 0.5;
    }
}

static void set_qrts(t_data *data)
{
    set_qrt(data, 2);
    if (data->count > 1) {
        set_qrt(data, 1);
        set_qrt(data, 3);
    } else {
        data->qrt[0] = 0.0;
        data->qrt[2] = 0.0;
    }
}

static int comp(const void *pa, const void *pb)
{
    const double a = *(const double *)pa;
    const double b = *(const double *)pb;

    return (a > b) ? 1 : (a < b) ? -1 : 0;
}

int main(void)
{
    double values[] = {3.7, 8.9, 7.1, 5.4, 1.2, 6.8, 4.3, 2.7};
    t_data data;

    data.value = values;
    data.count = (int)(sizeof(values) / sizeof(double));
    qsort(data.value, data.count, sizeof(double), comp);
    set_qrts(&data);
    printf("Q1 = %.1f
Q2 = %.1f
Q3 = %.1f
", data.qrt[0], data.qrt[1], data.qrt[2]);
}

是

d = qrt * 0.25 * data->count + 0.5;
n = (int)d;
e = n != d;

保证按预期工作？(e==isinteger(D))

推荐答案

数字0.5、0.25、0.125等表示2的负幂，因此在IEEE 754 types中完全可以表示。使用这些数字不会导致表示错误。

这篇关于0.25的倍数可以精确地表示为双精度吗？的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持！