author | Tom Henderson <tomh@tomh.org> |
Wed, 18 Aug 2010 15:49:18 -0700 | |
changeset 6584 | def4153e27cd |
parent 6583 | b47b2f59d6a8 |
permissions | -rw-r--r-- |
168
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
1 |
/* -*- Mode:C++; c-file-style:"gnu"; indent-tabs-mode:nil; -*- */ |
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
2 |
/* |
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
3 |
* Copyright (c) 2006 INRIA |
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
4 |
* |
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
5 |
* This program is free software; you can redistribute it and/or modify |
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
6 |
* it under the terms of the GNU General Public License version 2 as |
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
7 |
* published by the Free Software Foundation; |
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
8 |
* |
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
9 |
* This program is distributed in the hope that it will be useful, |
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
10 |
* but WITHOUT ANY WARRANTY; without even the implied warranty of |
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
11 |
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
12 |
* GNU General Public License for more details. |
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
13 |
* |
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
14 |
* You should have received a copy of the GNU General Public License |
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
15 |
* along with this program; if not, write to the Free Software |
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
16 |
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA |
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
17 |
* |
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
18 |
* Author: Mathieu Lacage <mathieu.lacage@sophia.inria.fr> |
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
19 |
*/ |
6404
5e271cbf9d16
introduce new hp implementation for 64bit boxes
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6180
diff
changeset
|
20 |
#ifndef HIGH_PRECISION_CAIRO_H |
5e271cbf9d16
introduce new hp implementation for 64bit boxes
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6180
diff
changeset
|
21 |
#define HIGH_PRECISION_CAIRO_H |
168
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
22 |
|
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
23 |
#include <stdint.h> |
6419
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
24 |
#include <math.h> |
6569
08154c0542ae
forgot include
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6566
diff
changeset
|
25 |
#include <iostream> |
170
9f68cb694f40
implement 128 bit arithmetic with cairo's internal types
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
168
diff
changeset
|
26 |
#include "cairo-wideint-private.h" |
168
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
27 |
|
216
27c5e5be4607
add implementation comment
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
214
diff
changeset
|
28 |
/** |
27c5e5be4607
add implementation comment
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
214
diff
changeset
|
29 |
* This file contains an implementation of the HighPrecision class. |
27c5e5be4607
add implementation comment
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
214
diff
changeset
|
30 |
* Each instance of the Time class also contains an instance of this |
27c5e5be4607
add implementation comment
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
214
diff
changeset
|
31 |
* class which is used to perform all the arithmetic operations of |
27c5e5be4607
add implementation comment
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
214
diff
changeset
|
32 |
* the Time class. |
27c5e5be4607
add implementation comment
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
214
diff
changeset
|
33 |
* |
27c5e5be4607
add implementation comment
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
214
diff
changeset
|
34 |
* This code is a bit ugly with a lot of inline methods for speed: |
27c5e5be4607
add implementation comment
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
214
diff
changeset
|
35 |
* profiling this code on anything but the simplest scenarios shows |
27c5e5be4607
add implementation comment
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
214
diff
changeset
|
36 |
* that it is a big bottleneck if great care in its implementation |
6180
cd0d8ba00e6c
coding style
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6177
diff
changeset
|
37 |
* is not performed. My observations are that what dominates are |
216
27c5e5be4607
add implementation comment
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
214
diff
changeset
|
38 |
* Division operations (there are really really super costly) |
27c5e5be4607
add implementation comment
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
214
diff
changeset
|
39 |
* and Comparison operations (because there are typically a lot of |
27c5e5be4607
add implementation comment
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
214
diff
changeset
|
40 |
* these in any complex timekeeping code). |
27c5e5be4607
add implementation comment
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
214
diff
changeset
|
41 |
*/ |
27c5e5be4607
add implementation comment
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
214
diff
changeset
|
42 |
|
168
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
43 |
namespace ns3 { |
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
44 |
|
6180
cd0d8ba00e6c
coding style
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6177
diff
changeset
|
45 |
class HighPrecision |
168
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
46 |
{ |
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
47 |
public: |
214
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
48 |
inline HighPrecision (); |
6505
2d59c219852b
another missing explicit keyword
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6500
diff
changeset
|
49 |
explicit inline HighPrecision (int64_t value, bool dummy); |
6496
ed210e7279a7
rewrite time integer support for 32bit systems.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6419
diff
changeset
|
50 |
explicit inline HighPrecision (double value); |
214
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
51 |
|
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
52 |
inline int64_t GetInteger (void) const; |
6419
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
53 |
inline double GetDouble (void) const; |
6414
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
54 |
inline void Add (HighPrecision const &o); |
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
55 |
inline void Sub (HighPrecision const &o); |
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
56 |
void Mul (HighPrecision const &o); |
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
57 |
void Div (HighPrecision const &o); |
6496
ed210e7279a7
rewrite time integer support for 32bit systems.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6419
diff
changeset
|
58 |
void MulByInvert (const HighPrecision &o); |
ed210e7279a7
rewrite time integer support for 32bit systems.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6419
diff
changeset
|
59 |
static HighPrecision Invert (uint64_t v); |
168
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
60 |
|
214
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
61 |
inline int Compare (HighPrecision const &o) const; |
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
62 |
inline static HighPrecision Zero (void); |
6581
b3597c013780
Provide stubbed HighPrecision GetHigh(), GetLow() implementations to fix builds across platforms
Tom Henderson <tomh@tomh.org>
parents:
6573
diff
changeset
|
63 |
int64_t GetHigh (void) const; |
b3597c013780
Provide stubbed HighPrecision GetHigh(), GetLow() implementations to fix builds across platforms
Tom Henderson <tomh@tomh.org>
parents:
6573
diff
changeset
|
64 |
uint64_t GetLow (void) const; |
168
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
65 |
private: |
6496
ed210e7279a7
rewrite time integer support for 32bit systems.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6419
diff
changeset
|
66 |
static cairo_uint128_t Umul (cairo_uint128_t a, cairo_uint128_t b); |
ed210e7279a7
rewrite time integer support for 32bit systems.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6419
diff
changeset
|
67 |
static cairo_uint128_t Udiv (cairo_uint128_t a, cairo_uint128_t b); |
ed210e7279a7
rewrite time integer support for 32bit systems.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6419
diff
changeset
|
68 |
static cairo_uint128_t UmulByInvert (cairo_uint128_t a, cairo_uint128_t b); |
6414
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
69 |
inline bool IsNegative (void) const; |
214
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
70 |
|
6414
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
71 |
cairo_int128_t m_value; |
168
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
72 |
}; |
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
73 |
|
6566
eb0098cd929c
allow serialization/deserialization for debugging output
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6519
diff
changeset
|
74 |
std::ostream &operator << (std::ostream &os, const HighPrecision &hp); |
eb0098cd929c
allow serialization/deserialization for debugging output
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6519
diff
changeset
|
75 |
std::istream &operator >> (std::istream &is, HighPrecision &hp); |
eb0098cd929c
allow serialization/deserialization for debugging output
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6519
diff
changeset
|
76 |
|
eb0098cd929c
allow serialization/deserialization for debugging output
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6519
diff
changeset
|
77 |
|
6414
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
78 |
} // namespace ns3 |
168
037cd2b37c67
split high precision implementations in different files
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
diff
changeset
|
79 |
|
214
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
80 |
namespace ns3 { |
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
81 |
|
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
82 |
HighPrecision::HighPrecision () |
6180
cd0d8ba00e6c
coding style
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6177
diff
changeset
|
83 |
{ |
6414
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
84 |
m_value.hi = 0; |
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
85 |
m_value.lo = 0; |
6180
cd0d8ba00e6c
coding style
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6177
diff
changeset
|
86 |
} |
214
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
87 |
|
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
88 |
HighPrecision::HighPrecision (int64_t value, bool dummy) |
6180
cd0d8ba00e6c
coding style
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6177
diff
changeset
|
89 |
{ |
6414
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
90 |
m_value.hi = value; |
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
91 |
m_value.lo = 0; |
6180
cd0d8ba00e6c
coding style
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6177
diff
changeset
|
92 |
} |
214
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
93 |
|
6414
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
94 |
bool |
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
95 |
HighPrecision::IsNegative (void) const |
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
96 |
{ |
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
97 |
int64_t hi = m_value.hi; |
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
98 |
return hi < 0; |
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
99 |
} |
214
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
100 |
|
6180
cd0d8ba00e6c
coding style
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6177
diff
changeset
|
101 |
int64_t |
214
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
102 |
HighPrecision::GetInteger (void) const |
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
103 |
{ |
6414
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
104 |
return m_value.hi; |
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
105 |
} |
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
106 |
void |
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
107 |
HighPrecision::Add (HighPrecision const &o) |
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
108 |
{ |
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
109 |
m_value.hi += o.m_value.hi; |
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
110 |
m_value.lo += o.m_value.lo; |
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
111 |
if (m_value.lo < o.m_value.lo) |
214
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
112 |
{ |
6414
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
113 |
m_value.hi++; |
214
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
114 |
} |
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
115 |
} |
6414
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
116 |
void |
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
117 |
HighPrecision::Sub (HighPrecision const &o) |
214
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
118 |
{ |
6414
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
119 |
m_value.hi -= o.m_value.hi; |
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
120 |
m_value.lo -= o.m_value.lo; |
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
121 |
if (m_value.lo > o.m_value.lo) |
214
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
122 |
{ |
6414
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
123 |
m_value.hi--; |
214
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
124 |
} |
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
125 |
} |
6414
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
126 |
int |
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
127 |
HighPrecision::Compare (HighPrecision const &o) const |
214
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
128 |
{ |
6419
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
129 |
HighPrecision tmp = *this; |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
130 |
tmp.Sub (o); |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
131 |
return tmp.IsNegative ()?-1:(tmp.m_value.hi == 0 && tmp.m_value.lo == 0)?0:1; |
214
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
132 |
} |
6180
cd0d8ba00e6c
coding style
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6177
diff
changeset
|
133 |
HighPrecision |
214
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
134 |
HighPrecision::Zero (void) |
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
135 |
{ |
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
136 |
return HighPrecision (); |
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
137 |
} |
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
138 |
|
6419
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
139 |
|
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
140 |
#define HPCAIRO_MAX_64 18446744073709551615.0 |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
141 |
double HighPrecision::GetDouble (void) const |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
142 |
{ |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
143 |
bool is_negative = IsNegative (); |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
144 |
cairo_int128_t value = is_negative ? _cairo_int128_negate (m_value) : m_value; |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
145 |
double flo = value.lo; |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
146 |
flo /= HPCAIRO_MAX_64; |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
147 |
double retval = value.hi; |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
148 |
retval += flo; |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
149 |
retval = is_negative ? -retval: retval; |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
150 |
return retval; |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
151 |
} |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
152 |
#undef HPCAIRO_MAX_64 |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
153 |
|
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
154 |
#define HPCAIRO_MAX_64 18446744073709551615.0 |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
155 |
HighPrecision::HighPrecision (double value) |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
156 |
{ |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
157 |
double fhi = floor (value); |
6573
9d414b94c46d
add explicit cast for building on gcc-4.1.1
Tom Henderson <tomh@tomh.org>
parents:
6569
diff
changeset
|
158 |
int64_t hi = static_cast<int64_t> (fhi); |
6419
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
159 |
uint64_t lo = (uint64_t) ((value - fhi) * HPCAIRO_MAX_64); |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
160 |
m_value.hi = hi; |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
161 |
m_value.lo = lo; |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
162 |
} |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
163 |
#undef HPCAIRO_MAX_64 |
c08597bde2a0
inline a few more functions for optimization
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6414
diff
changeset
|
164 |
|
6414
e1ddd0c5cd22
bug 826: optimize hp implementation for 32bit systems
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6404
diff
changeset
|
165 |
} // namespace ns3 |
214
ac309a63f850
optimize Time arithmetic operations.
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
170
diff
changeset
|
166 |
|
6404
5e271cbf9d16
introduce new hp implementation for 64bit boxes
Mathieu Lacage <mathieu.lacage@sophia.inria.fr>
parents:
6180
diff
changeset
|
167 |
#endif /* HIGH_PRECISION_CAIRO_H */ |