From cad58fc5e7aba22dee7f9099e9fc46634b6a439c Mon Sep 17 00:00:00 2001 From: Elie Tournier Date: Wed, 9 Aug 2017 23:11:08 +0100 Subject: glsl: Add "built-in" functions to do fp32_to_fp64(fp32) Signed-off-by: Elie Tournier --- src/compiler/glsl/float64.glsl | 38 ++++++++++++++++++++++++++++++++++++++ 1 file changed, 38 insertions(+) (limited to 'src/compiler/glsl/float64.glsl') diff --git a/src/compiler/glsl/float64.glsl b/src/compiler/glsl/float64.glsl index 8cc0a278904..05e07d39ef6 100644 --- a/src/compiler/glsl/float64.glsl +++ b/src/compiler/glsl/float64.glsl @@ -1045,3 +1045,41 @@ __fp64_to_fp32(uint64_t __a) zFrac = mix(zFrac, zFrac | 0x40000000u, aExp != 0); return __roundAndPackFloat32(aSign, aExp - 0x381, zFrac); } + +/* Returns the result of converting the single-precision floating-point value + * `a' to the double-precision floating-point format. + */ +uint64_t +__fp32_to_fp64(float f) +{ + uint a = floatBitsToUint(f); + uint aFrac = a & 0x007FFFFFu; + int aExp = int((a>>23) & 0xFFu); + uint aSign = a>>31; + uint zFrac0 = 0u; + uint zFrac1 = 0u; + + if (aExp == 0xFF) { + if (aFrac != 0u) { + uint nanLo = 0u; + uint nanHi = a<<9; + __shift64Right(nanHi, nanLo, 12, nanHi, nanLo); + nanHi |= ((aSign<<31) | 0x7FF80000u); + return packUint2x32(uvec2(nanLo, nanHi)); + } + return __packFloat64(aSign, 0x7FF, 0u, 0u); + } + + if (aExp == 0) { + if (aFrac == 0u) + return __packFloat64(aSign, 0, 0u, 0u); + /* Normalize subnormal */ + int shiftCount = __countLeadingZeros32(aFrac) - 8; + aFrac <<= shiftCount; + aExp = 1 - shiftCount; + --aExp; + } + + __shift64Right(aFrac, 0u, 3, zFrac0, zFrac1); + return __packFloat64(aSign, aExp + 0x380, zFrac0, zFrac1); +} -- cgit v1.2.3