Update copyright to LGPL on all files
[dyninst.git] / common / src / irixKludges.C
1 /*
2  * Copyright (c) 1996-2009 Barton P. Miller
3  * 
4  * We provide the Paradyn Parallel Performance Tools (below
5  * described as "Paradyn") on an AS IS basis, and do not warrant its
6  * validity or performance.  We reserve the right to update, modify,
7  * or discontinue this software at any time.  We shall have no
8  * obligation to supply such updates or modifications or any other
9  * form of support to you.
10  * 
11  * By your use of Paradyn, you understand and agree that we (or any
12  * other person or entity with proprietary rights in Paradyn) are
13  * under no obligation to provide either maintenance services,
14  * update services, notices of latent defects, or correction of
15  * defects for Paradyn.
16  * 
17  * This library is free software; you can redistribute it and/or
18  * modify it under the terms of the GNU Lesser General Public
19  * License as published by the Free Software Foundation; either
20  * version 2.1 of the License, or (at your option) any later version.
21  * 
22  * This library is distributed in the hope that it will be useful,
23  * but WITHOUT ANY WARRANTY; without even the implied warranty of
24  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
25  * Lesser General Public License for more details.
26  * 
27  * You should have received a copy of the GNU Lesser General Public
28  * License along with this library; if not, write to the Free Software
29  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
30  */
31
32 #include "common/h/headers.h"
33 #include "common/h/irixKludges.h"
34
35 // copied from solarisKludges.C
36 unsigned long long PDYN_div1000(unsigned long long in) {
37    /* Divides by 1000 without an integer division instruction or library call, both of
38     * which are slow.
39     * We do only shifts, adds, and subtracts.
40     *
41     * We divide by 1000 in this way:
42     * multiply by 1/1000, or multiply by (1/1000)*2^30 and then right-shift by 30.
43     * So what is 1/1000 * 2^30?
44     * It is 1,073,742.   (actually this is rounded)
45     * So we can multiply by 1,073,742 and then right-shift by 30 (neat, eh?)
46     *
47     * Now for multiplying by 1,073,742...
48     * 1,073,742 = (1,048,576 + 16384 + 8192 + 512 + 64 + 8 + 4 + 2)
49     * or, slightly optimized:
50     * = (1,048,576 + 16384 + 8192 + 512 + 64 + 16 - 2)
51     * for a total of 8 shifts and 6 add/subs, or 14 operations.
52     *
53     */
54
55    unsigned long long temp = in << 20; // multiply by 1,048,576
56       // beware of overflow; left shift by 20 is quite a lot.
57       // If you know that the input fits in 32 bits (4 billion) then
58       // no problem.  But if it's much bigger then start worrying...
59
60    temp += in << 14; // 16384
61    temp += in << 13; // 8192
62    temp += in << 9;  // 512
63    temp += in << 6;  // 64
64    temp += in << 4;  // 16
65    temp -= in >> 2;  // 2
66
67    return (temp >> 30); // divide by 2^30
68 }
69
70 // copied from solarisKludges.C
71 unsigned long long PDYN_divMillion(unsigned long long in) {
72    /* Divides by 1,000,000 without an integer division instruction or library call,
73     * both of which are slow.
74     * We do only shifts, adds, and subtracts.
75     *
76     * We divide by 1,000,000 in this way:
77     * multiply by 1/1,000,000, or multiply by (1/1,000,000)*2^30 and then right-shift
78     * by 30.  So what is 1/1,000,000 * 2^30?
79     * It is 1,074.   (actually this is rounded)
80     * So we can multiply by 1,074 and then right-shift by 30 (neat, eh?)
81     *
82     * Now for multiplying by 1,074
83     * 1,074 = (1024 + 32 + 16 + 2)
84     * for a total of 4 shifts and 4 add/subs, or 8 operations.
85     *
86     * Note: compare with div1000 -- it's cheaper to divide by a million than
87     *       by a thousand (!)
88     *
89     */
90
91    unsigned long long temp = in << 10; // multiply by 1024
92       // beware of overflow...if the input arg uses more than 52 bits
93       // than start worrying about whether (in << 10) plus the smaller additions
94       // we're gonna do next will fit in 64...
95
96    temp += in << 5; // 32
97    temp += in << 4; // 16
98    temp += in << 1; // 2
99
100    return (temp >> 30); // divide by 2^30
101 }
102
103 // copied from solarisKludges.C
104 unsigned long long PDYN_mulMillion(unsigned long long in) {
105    unsigned long long result = in;
106
107    /* multiply by 125 by multiplying by 128 and subtracting 3x */
108    result = (result << 7) - result - result - result;
109
110    /* multiply by 125 again, for a total of 15625x */
111    result = (result << 7) - result - result - result;
112
113    /* multiply by 64, for a total of 1,000,000x */
114    result <<= 6;
115
116    /* cost was: 3 shifts and 6 subtracts
117     * cost of calling mul1000(mul1000()) would be: 6 shifts and 4 subtracts
118     *
119     * Another algorithm is to multiply by 2^6 and then 5^6.
120     * The former is super-cheap (one shift); the latter is more expensive.
121     * 5^6 = 15625 = 16384 - 512 - 256 + 8 + 1
122     * so multiplying by 5^6 means 4 shift operations and 4 add/sub ops
123     * so multiplying by 1000000 means 5 shift operations and 4 add/sub ops.
124     * That may or may not be cheaper than what we're doing (3 shifts; 6 subtracts);
125     * I'm not sure.  --ari
126     */
127
128    return result;
129 }