1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
|
/* Copyright (c) <2003-2011> <Julio Jerez, Newton Game Dynamics>
*
* This software is provided 'as-is', without any express or implied
* warranty. In no event will the authors be held liable for any damages
* arising from the use of this software.
*
* Permission is granted to anyone to use this software for any purpose,
* including commercial applications, and to alter it and redistribute it
* freely, subject to the following restrictions:
*
* 1. The origin of this software must not be misrepresented; you must not
* claim that you wrote the original software. If you use this software
* in a product, an acknowledgment in the product documentation would be
* appreciated but is not required.
*
* 2. Altered source versions must be plainly marked as such, and must not be
* misrepresented as being the original software.
*
* 3. This notice may not be removed or altered from any source distribution.
*/
#ifndef __dgGoogol__
#define __dgGoogol__
#include "dgArray.h"
#include "dgVector.h"
//#define DG_GOOGOL_SIZE 16
#define DG_GOOGOL_SIZE 4
class dgGoogol
{
public:
dgGoogol(void);
dgGoogol(double value);
~dgGoogol(void);
double GetAproximateValue() const;
void InitFloatFloat (double value);
dgGoogol operator+ (const dgGoogol &A) const;
dgGoogol operator- (const dgGoogol &A) const;
dgGoogol operator* (const dgGoogol &A) const;
dgGoogol operator/ (const dgGoogol &A) const;
dgGoogol operator+= (const dgGoogol &A);
dgGoogol operator-= (const dgGoogol &A);
dgGoogol Floor () const;
#ifdef _DEBUG
void ToString (char* const string) const;
#endif
private:
void NegateMantissa (uint64_t* const mantissa) const;
void CopySignedMantissa (uint64_t* const mantissa) const;
int32_t NormalizeMantissa (uint64_t* const mantissa) const;
uint64_t CheckCarrier (uint64_t a, uint64_t b) const;
void ShiftRightMantissa (uint64_t* const mantissa, int32_t bits) const;
int32_t LeadinZeros (uint64_t a) const;
void ExtendeMultiply (uint64_t a, uint64_t b, uint64_t& high, uint64_t& low) const;
void ScaleMantissa (uint64_t* const out, uint64_t scale) const;
int8_t m_sign;
int16_t m_exponent;
uint64_t m_mantissa[DG_GOOGOL_SIZE];
};
class dgHugeVector: public dgTemplateVector<dgGoogol>
{
public:
dgHugeVector ()
:dgTemplateVector<dgGoogol>()
{
}
dgHugeVector (const dgBigVector& a)
:dgTemplateVector<dgGoogol>(dgGoogol (a.m_x), dgGoogol (a.m_y), dgGoogol (a.m_z), dgGoogol (a.m_w))
{
}
dgHugeVector (const dgTemplateVector<dgGoogol>& a)
:dgTemplateVector<dgGoogol>(a)
{
}
dgHugeVector (double x, double y, double z, double w)
:dgTemplateVector<dgGoogol>(x, y, z, w)
{
}
dgGoogol EvaluePlane (const dgHugeVector& point) const
{
return (point % (*this)) + m_w;
}
};
#endif
|