C如何处理混合精度运算

Question

float  FLOAT  = 1.0f;
double DOUBLE = 2.0;

float  a = FLOAT / DOUBLE;
double b = FLOAT / DOUBLE;

a 和 b 的计算方式是否相同？

FLOAT 和 DOUBLE 在编译时如何转换？

看来默认转换是向上转换，以防丢失。

其实我是在 GPU 上做一些对精度敏感的计算，代码应该是这样的：

float a = FLOAT / 2.0 + 1.0/3.0;

其中代码包含很长的表达式，包含许多数字和变量（实际上是从 Matlab 代码生成的）。

那我该如何控制这种转换行为呢？ 除了在2.0f 中写下所有数字（表达式中的数以千计）。

Answer 1

根据 C 标准（6.3.1.8 常用算术转换）

否则，如果任一操作数对应的实数类型为双精度， 另一个操作数在不改变类型域的情况下转换为 对应的实数类型为double的类型。

所以在声明的两个初始化表达式中

float  a = FLOAT / DOUBLE;
double b = FLOAT / DOUBLE;

float 类型的操作数被转换为double 类型。因此，表达式的结果具有double 类型。在第一个声明中，结果进一步转换为 float，因为变量 a 的类型为 float。