137a6a455466e5b197311771a777ab241e471ed8aEdward O'Callaghan/* ===-- divxc3.c - Implement __divxc3 -------------------------------------=== 237a6a455466e5b197311771a777ab241e471ed8aEdward O'Callaghan * 337a6a455466e5b197311771a777ab241e471ed8aEdward O'Callaghan * The LLVM Compiler Infrastructure 437a6a455466e5b197311771a777ab241e471ed8aEdward O'Callaghan * 59ad441ffec97db647fee3725b3424284fb913e14Howard Hinnant * This file is dual licensed under the MIT and the University of Illinois Open 69ad441ffec97db647fee3725b3424284fb913e14Howard Hinnant * Source Licenses. See LICENSE.TXT for details. 737a6a455466e5b197311771a777ab241e471ed8aEdward O'Callaghan * 837a6a455466e5b197311771a777ab241e471ed8aEdward O'Callaghan * ===----------------------------------------------------------------------=== 937a6a455466e5b197311771a777ab241e471ed8aEdward O'Callaghan * 1037a6a455466e5b197311771a777ab241e471ed8aEdward O'Callaghan * This file implements __divxc3 for the compiler_rt library. 1137a6a455466e5b197311771a777ab241e471ed8aEdward O'Callaghan * 1237a6a455466e5b197311771a777ab241e471ed8aEdward O'Callaghan */ 13b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar 14b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar#if !_ARCH_PPC 15b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar 16b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar#include "int_lib.h" 178603024a01cfa82ab8dce397f40a2d6f246076e3Daniel Dunbar#include "int_math.h" 18b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar 1937a6a455466e5b197311771a777ab241e471ed8aEdward O'Callaghan/* Returns: the quotient of (a + ib) / (c + id) */ 20b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar 21b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbarlong double _Complex 22b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar__divxc3(long double __a, long double __b, long double __c, long double __d) 23b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar{ 24b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar int __ilogbw = 0; 25c25c6d10b18efa071f0016f31f6f3dd8a8fa6676Daniel Dunbar long double __logbw = crt_logbl(crt_fmaxl(crt_fabsl(__c), crt_fabsl(__d))); 268603024a01cfa82ab8dce397f40a2d6f246076e3Daniel Dunbar if (crt_isfinite(__logbw)) 27b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar { 28b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar __ilogbw = (int)__logbw; 29c25c6d10b18efa071f0016f31f6f3dd8a8fa6676Daniel Dunbar __c = crt_scalbnl(__c, -__ilogbw); 30c25c6d10b18efa071f0016f31f6f3dd8a8fa6676Daniel Dunbar __d = crt_scalbnl(__d, -__ilogbw); 31b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar } 32b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar long double __denom = __c * __c + __d * __d; 33b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar long double _Complex z; 34c25c6d10b18efa071f0016f31f6f3dd8a8fa6676Daniel Dunbar __real__ z = crt_scalbnl((__a * __c + __b * __d) / __denom, -__ilogbw); 35c25c6d10b18efa071f0016f31f6f3dd8a8fa6676Daniel Dunbar __imag__ z = crt_scalbnl((__b * __c - __a * __d) / __denom, -__ilogbw); 368603024a01cfa82ab8dce397f40a2d6f246076e3Daniel Dunbar if (crt_isnan(__real__ z) && crt_isnan(__imag__ z)) 37b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar { 388603024a01cfa82ab8dce397f40a2d6f246076e3Daniel Dunbar if ((__denom == 0) && (!crt_isnan(__a) || !crt_isnan(__b))) 39b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar { 40c25c6d10b18efa071f0016f31f6f3dd8a8fa6676Daniel Dunbar __real__ z = crt_copysignl(CRT_INFINITY, __c) * __a; 41c25c6d10b18efa071f0016f31f6f3dd8a8fa6676Daniel Dunbar __imag__ z = crt_copysignl(CRT_INFINITY, __c) * __b; 42b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar } 438603024a01cfa82ab8dce397f40a2d6f246076e3Daniel Dunbar else if ((crt_isinf(__a) || crt_isinf(__b)) && 448603024a01cfa82ab8dce397f40a2d6f246076e3Daniel Dunbar crt_isfinite(__c) && crt_isfinite(__d)) 45b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar { 46c25c6d10b18efa071f0016f31f6f3dd8a8fa6676Daniel Dunbar __a = crt_copysignl(crt_isinf(__a) ? 1 : 0, __a); 47c25c6d10b18efa071f0016f31f6f3dd8a8fa6676Daniel Dunbar __b = crt_copysignl(crt_isinf(__b) ? 1 : 0, __b); 48c25c6d10b18efa071f0016f31f6f3dd8a8fa6676Daniel Dunbar __real__ z = CRT_INFINITY * (__a * __c + __b * __d); 49c25c6d10b18efa071f0016f31f6f3dd8a8fa6676Daniel Dunbar __imag__ z = CRT_INFINITY * (__b * __c - __a * __d); 50b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar } 518603024a01cfa82ab8dce397f40a2d6f246076e3Daniel Dunbar else if (crt_isinf(__logbw) && __logbw > 0 && 528603024a01cfa82ab8dce397f40a2d6f246076e3Daniel Dunbar crt_isfinite(__a) && crt_isfinite(__b)) 53b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar { 54c25c6d10b18efa071f0016f31f6f3dd8a8fa6676Daniel Dunbar __c = crt_copysignl(crt_isinf(__c) ? 1 : 0, __c); 55c25c6d10b18efa071f0016f31f6f3dd8a8fa6676Daniel Dunbar __d = crt_copysignl(crt_isinf(__d) ? 1 : 0, __d); 56b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar __real__ z = 0 * (__a * __c + __b * __d); 57b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar __imag__ z = 0 * (__b * __c - __a * __d); 58b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar } 59b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar } 60b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar return z; 61b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar} 62b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar 63b3a6901e66f55b35aa9e01bcb24134e6a65ea004Daniel Dunbar#endif 64