2009-08-06 03:06:50 +08:00
|
|
|
/* ===-- powitf2.cpp - Implement __powitf2 ---------------------------------===
|
|
|
|
*
|
|
|
|
* The LLVM Compiler Infrastructure
|
|
|
|
*
|
2010-11-17 06:13:33 +08:00
|
|
|
* This file is dual licensed under the MIT and the University of Illinois Open
|
|
|
|
* Source Licenses. See LICENSE.TXT for details.
|
2009-08-06 03:06:50 +08:00
|
|
|
*
|
|
|
|
* ===----------------------------------------------------------------------===
|
|
|
|
*
|
|
|
|
* This file implements __powitf2 for the compiler_rt library.
|
|
|
|
*
|
|
|
|
* ===----------------------------------------------------------------------===
|
|
|
|
*/
|
2009-06-27 00:47:03 +08:00
|
|
|
|
|
|
|
#if _ARCH_PPC
|
|
|
|
|
|
|
|
#include "int_lib.h"
|
|
|
|
|
2009-08-06 03:06:50 +08:00
|
|
|
/* Returns: a ^ b */
|
2009-06-27 00:47:03 +08:00
|
|
|
|
|
|
|
long double
|
|
|
|
__powitf2(long double a, si_int b)
|
|
|
|
{
|
|
|
|
const int recip = b < 0;
|
|
|
|
long double r = 1;
|
|
|
|
while (1)
|
|
|
|
{
|
|
|
|
if (b & 1)
|
|
|
|
r *= a;
|
|
|
|
b /= 2;
|
|
|
|
if (b == 0)
|
|
|
|
break;
|
|
|
|
a *= a;
|
|
|
|
}
|
|
|
|
return recip ? 1/r : r;
|
|
|
|
}
|
|
|
|
|
|
|
|
#endif
|