1*0b57cec5SDimitry Andric //===-- powidf2.cpp - Implement __powidf2 ---------------------------------===// 2*0b57cec5SDimitry Andric // 3*0b57cec5SDimitry Andric // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions. 4*0b57cec5SDimitry Andric // See https://llvm.org/LICENSE.txt for license information. 5*0b57cec5SDimitry Andric // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception 6*0b57cec5SDimitry Andric // 7*0b57cec5SDimitry Andric //===----------------------------------------------------------------------===// 8*0b57cec5SDimitry Andric // 9*0b57cec5SDimitry Andric // This file implements __powidf2 for the compiler_rt library. 10*0b57cec5SDimitry Andric // 11*0b57cec5SDimitry Andric //===----------------------------------------------------------------------===// 12*0b57cec5SDimitry Andric 13*0b57cec5SDimitry Andric #include "int_lib.h" 14*0b57cec5SDimitry Andric 15*0b57cec5SDimitry Andric // Returns: a ^ b 16*0b57cec5SDimitry Andric 17*0b57cec5SDimitry Andric COMPILER_RT_ABI double __powidf2(double a, si_int b) { 18*0b57cec5SDimitry Andric const int recip = b < 0; 19*0b57cec5SDimitry Andric double r = 1; 20*0b57cec5SDimitry Andric while (1) { 21*0b57cec5SDimitry Andric if (b & 1) 22*0b57cec5SDimitry Andric r *= a; 23*0b57cec5SDimitry Andric b /= 2; 24*0b57cec5SDimitry Andric if (b == 0) 25*0b57cec5SDimitry Andric break; 26*0b57cec5SDimitry Andric a *= a; 27*0b57cec5SDimitry Andric } 28*0b57cec5SDimitry Andric return recip ? 1 / r : r; 29*0b57cec5SDimitry Andric } 30