PR libgomp/64635
[official-gcc.git] / gcc / ada / s-expllu.adb
blob47192b9b0a8600dfbd6580c3df078df91d1f7dc3
1 ------------------------------------------------------------------------------
2 -- --
3 -- GNAT RUN-TIME COMPONENTS --
4 -- --
5 -- S Y S T E M . X P _ B M L --
6 -- --
7 -- B o d y --
8 -- --
9 -- Copyright (C) 1992-2014, Free Software Foundation, Inc. --
10 -- --
11 -- GNAT is free software; you can redistribute it and/or modify it under --
12 -- terms of the GNU General Public License as published by the Free Soft- --
13 -- ware Foundation; either version 3, or (at your option) any later ver- --
14 -- sion. GNAT is distributed in the hope that it will be useful, but WITH- --
15 -- OUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY --
16 -- or FITNESS FOR A PARTICULAR PURPOSE. --
17 -- --
18 -- As a special exception under Section 7 of GPL version 3, you are granted --
19 -- additional permissions described in the GCC Runtime Library Exception, --
20 -- version 3.1, as published by the Free Software Foundation. --
21 -- --
22 -- You should have received a copy of the GNU General Public License and --
23 -- a copy of the GCC Runtime Library Exception along with this program; --
24 -- see the files COPYING3 and COPYING.RUNTIME respectively. If not, see --
25 -- <http://www.gnu.org/licenses/>. --
26 -- --
27 -- GNAT was originally developed by the GNAT team at New York University. --
28 -- Extensive contributions were provided by Ada Core Technologies Inc. --
29 -- --
30 ------------------------------------------------------------------------------
32 with System.Unsigned_Types; use System.Unsigned_Types;
34 package body System.Exp_LLU is
36 ----------------------------
37 -- Exp_Long_Long_Unsigned --
38 ----------------------------
40 function Exp_Long_Long_Unsigned
41 (Left : Long_Long_Unsigned;
42 Right : Natural)
43 return Long_Long_Unsigned
45 Result : Long_Long_Unsigned := 1;
46 Factor : Long_Long_Unsigned := Left;
47 Exp : Natural := Right;
49 begin
50 -- We use the standard logarithmic approach, Exp gets shifted right
51 -- testing successive low order bits and Factor is the value of the
52 -- base raised to the next power of 2.
54 -- Note: it is not worth special casing the cases of base values -1,0,+1
55 -- since the expander does this when the base is a literal, and other
56 -- cases will be extremely rare.
58 if Exp /= 0 then
59 loop
60 if Exp rem 2 /= 0 then
61 Result := Result * Factor;
62 end if;
64 Exp := Exp / 2;
65 exit when Exp = 0;
66 Factor := Factor * Factor;
67 end loop;
68 end if;
70 return Result;
72 end Exp_Long_Long_Unsigned;
74 end System.Exp_LLU;