10b57cec5SDimitry Andric //===-- powixf2.cpp - Implement __powixf2 ---------------------------------===// 20b57cec5SDimitry Andric // 30b57cec5SDimitry Andric // Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions. 40b57cec5SDimitry Andric // See https://llvm.org/LICENSE.txt for license information. 50b57cec5SDimitry Andric // SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception 60b57cec5SDimitry Andric // 70b57cec5SDimitry Andric //===----------------------------------------------------------------------===// 80b57cec5SDimitry Andric // 90b57cec5SDimitry Andric // This file implements __powixf2 for the compiler_rt library. 100b57cec5SDimitry Andric // 110b57cec5SDimitry Andric //===----------------------------------------------------------------------===// 120b57cec5SDimitry Andric 130b57cec5SDimitry Andric #if !_ARCH_PPC 140b57cec5SDimitry Andric 150b57cec5SDimitry Andric #include "int_lib.h" 160b57cec5SDimitry Andric 170b57cec5SDimitry Andric // Returns: a ^ b 180b57cec5SDimitry Andric 19*5f757f3fSDimitry Andric COMPILER_RT_ABI xf_float __powixf2(xf_float a, int b) { 200b57cec5SDimitry Andric const int recip = b < 0; 21*5f757f3fSDimitry Andric xf_float r = 1; 220b57cec5SDimitry Andric while (1) { 230b57cec5SDimitry Andric if (b & 1) 240b57cec5SDimitry Andric r *= a; 250b57cec5SDimitry Andric b /= 2; 260b57cec5SDimitry Andric if (b == 0) 270b57cec5SDimitry Andric break; 280b57cec5SDimitry Andric a *= a; 290b57cec5SDimitry Andric } 300b57cec5SDimitry Andric return recip ? 1 / r : r; 310b57cec5SDimitry Andric } 320b57cec5SDimitry Andric 330b57cec5SDimitry Andric #endif 34