10b57cec5SDimitry Andric //===-- powitf2.cpp - Implement __powitf2 ---------------------------------===// 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 __powitf2 for the compiler_rt library. 100b57cec5SDimitry Andric // 110b57cec5SDimitry Andric //===----------------------------------------------------------------------===// 120b57cec5SDimitry Andric 135ffd83dbSDimitry Andric #define QUAD_PRECISION 145ffd83dbSDimitry Andric #include "fp_lib.h" 150b57cec5SDimitry Andric 1606c3fb27SDimitry Andric #if defined(CRT_HAS_TF_MODE) 170b57cec5SDimitry Andric 180b57cec5SDimitry Andric // Returns: a ^ b 190b57cec5SDimitry Andric 20*5f757f3fSDimitry Andric COMPILER_RT_ABI fp_t __powitf2(fp_t a, int b) { 210b57cec5SDimitry Andric const int recip = b < 0; 22*5f757f3fSDimitry Andric fp_t r = 1; 230b57cec5SDimitry Andric while (1) { 240b57cec5SDimitry Andric if (b & 1) 250b57cec5SDimitry Andric r *= a; 260b57cec5SDimitry Andric b /= 2; 270b57cec5SDimitry Andric if (b == 0) 280b57cec5SDimitry Andric break; 290b57cec5SDimitry Andric a *= a; 300b57cec5SDimitry Andric } 310b57cec5SDimitry Andric return recip ? 1 / r : r; 320b57cec5SDimitry Andric } 330b57cec5SDimitry Andric 340b57cec5SDimitry Andric #endif 35