OpenCL 浮点计算

1.浮点计算。标准定义了三种表示实数的方法，分别表示为float,double,和half等数据类型。opencl标准只定义了float数据类型。如果目标设备还能支持另外两种，编程时也可以使用他们来计算处理。

2.float数据另外i行可以表示的数值共分为下面四种：

（1）正规化数--数字完全可以用float数据类型来表示；
（2）去正规化数--数字的绝对值比最小所能表示的正规化还小；
（3）无穷数--数字的绝对值比最大所能表示的正规化数还大；
（4）不是一个数（NaN）--非法操作（例如0/0，对负数求平方根、取对数）得到的数字。

openCL兼容设备支持第三类数（INF）和第四类数（NaN），但不支持第二类数，去正规化数运行处理所需的时间一般比正规化所需的时间长，因此可以考虑给编译函数clBuildProgram加入编译选项_cl_denorms_are_zero，来提高浮点运算的性能，如果在除法中使用去正规化数，结果是有意义的，而如果使用0，结果变为NaN.

3.如果对浮点数进行近似取舍，标准也分别定义了四种模式：

（1）向最近的偶数取舍--将浮点数向最近可表示的数值做取舍。如果浮点数是在两数之间，则向最后一位是0（也就是偶数）做取舍（opencl标准只要求第一种）
（2）向正无穷取舍--向更接近正无穷的数做取舍；
（3）向负无穷取舍--向更接近负无穷的数做取舍；
（4）向0取舍--向更接近0的值做取舍（舍弃小数部分）。

4.双转度数据类型：double数据类型是用64位宽来表示浮点数。
从0开始，有效数位（0-51）52bits，指数位（51-63）11bits,还有一位符号位（63-64）。

5.Opencl兼容设备并不需要支持double型，需要设定CL_DEVICE_EXTENSIONS参数，调用clGetDeviceInfo函数的方式来获取设备关于这一方面的信息，一次来确定cl_khr_fp64扩展是否被支持。

6.如果你的设备支持double数据类型，opencl将会更多的使用double型，而非float型来做数据处理，原因可分为以下几点。其一，支持double型的设备也必定支持正规化数，有些计算结果就不需要简单的划为0；其二，设备页必定支持之前所述的四种取舍模式，而不是简单的向最近的偶数做取舍。其三，支持double型的设备也必定支持所谓的浮点操作（floating point multiply-and-add,FMA），这个操作比乘法和加法的分立运算要快得多。