Project

General

Profile

Statistics
| Revision:

root / lab4 / .minix-src / include / c++ / ratio @ 13

History | View | Annotate | Download (14.4 KB)

1
// -*- C++ -*-
2
//===---------------------------- ratio -----------------------------------===//
3
//
4
//                     The LLVM Compiler Infrastructure
5
//
6
// This file is dual licensed under the MIT and the University of Illinois Open
7
// Source Licenses. See LICENSE.TXT for details.
8
//
9
//===----------------------------------------------------------------------===//
10

    
11
#ifndef _LIBCPP_RATIO
12
#define _LIBCPP_RATIO
13

    
14
/*
15
    ratio synopsis
16

    
17
namespace std
18
{
19

    
20
template <intmax_t N, intmax_t D = 1>
21
class ratio
22
{
23
public:
24
    static constexpr intmax_t num;
25
    static constexpr intmax_t den;
26
    typedef ratio<num, den> type;
27
};
28

    
29
// ratio arithmetic
30
template <class R1, class R2> using ratio_add = ...;
31
template <class R1, class R2> using ratio_subtract = ...;
32
template <class R1, class R2> using ratio_multiply = ...;
33
template <class R1, class R2> using ratio_divide = ...;
34

    
35
// ratio comparison
36
template <class R1, class R2> struct ratio_equal;
37
template <class R1, class R2> struct ratio_not_equal;
38
template <class R1, class R2> struct ratio_less;
39
template <class R1, class R2> struct ratio_less_equal;
40
template <class R1, class R2> struct ratio_greater;
41
template <class R1, class R2> struct ratio_greater_equal;
42

    
43
// convenience SI typedefs
44
typedef ratio<1, 1000000000000000000000000> yocto;  // not supported
45
typedef ratio<1,    1000000000000000000000> zepto;  // not supported
46
typedef ratio<1,       1000000000000000000> atto;
47
typedef ratio<1,          1000000000000000> femto;
48
typedef ratio<1,             1000000000000> pico;
49
typedef ratio<1,                1000000000> nano;
50
typedef ratio<1,                   1000000> micro;
51
typedef ratio<1,                      1000> milli;
52
typedef ratio<1,                       100> centi;
53
typedef ratio<1,                        10> deci;
54
typedef ratio<                       10, 1> deca;
55
typedef ratio<                      100, 1> hecto;
56
typedef ratio<                     1000, 1> kilo;
57
typedef ratio<                  1000000, 1> mega;
58
typedef ratio<               1000000000, 1> giga;
59
typedef ratio<            1000000000000, 1> tera;
60
typedef ratio<         1000000000000000, 1> peta;
61
typedef ratio<      1000000000000000000, 1> exa;
62
typedef ratio<   1000000000000000000000, 1> zetta;  // not supported
63
typedef ratio<1000000000000000000000000, 1> yotta;  // not supported
64

    
65
}
66
*/
67

    
68
#include <__config>
69
#include <cstdint>
70
#include <climits>
71
#include <type_traits>
72

    
73
#include <__undef_min_max>
74

    
75
#if !defined(_LIBCPP_HAS_NO_PRAGMA_SYSTEM_HEADER)
76
#pragma GCC system_header
77
#endif
78

    
79
_LIBCPP_BEGIN_NAMESPACE_STD
80

    
81
// __static_gcd
82

    
83
template <intmax_t _Xp, intmax_t _Yp>
84
struct __static_gcd
85
{
86
    static const intmax_t value = __static_gcd<_Yp, _Xp % _Yp>::value;
87
};
88

    
89
template <intmax_t _Xp>
90
struct __static_gcd<_Xp, 0>
91
{
92
    static const intmax_t value = _Xp;
93
};
94

    
95
template <>
96
struct __static_gcd<0, 0>
97
{
98
    static const intmax_t value = 1;
99
};
100

    
101
// __static_lcm
102

    
103
template <intmax_t _Xp, intmax_t _Yp>
104
struct __static_lcm
105
{
106
    static const intmax_t value = _Xp / __static_gcd<_Xp, _Yp>::value * _Yp;
107
};
108

    
109
template <intmax_t _Xp>
110
struct __static_abs
111
{
112
    static const intmax_t value = _Xp < 0 ? -_Xp : _Xp;
113
};
114

    
115
template <intmax_t _Xp>
116
struct __static_sign
117
{
118
    static const intmax_t value = _Xp == 0 ? 0 : (_Xp < 0 ? -1 : 1);
119
};
120

    
121
template <intmax_t _Xp, intmax_t _Yp, intmax_t = __static_sign<_Yp>::value>
122
class __ll_add;
123

    
124
template <intmax_t _Xp, intmax_t _Yp>
125
class __ll_add<_Xp, _Yp, 1>
126
{
127
    static const intmax_t min = (1LL << (sizeof(intmax_t) * CHAR_BIT - 1)) + 1;
128
    static const intmax_t max = -min;
129

    
130
    static_assert(_Xp <= max - _Yp, "overflow in __ll_add");
131
public:
132
    static const intmax_t value = _Xp + _Yp;
133
};
134

    
135
template <intmax_t _Xp, intmax_t _Yp>
136
class __ll_add<_Xp, _Yp, 0>
137
{
138
public:
139
    static const intmax_t value = _Xp;
140
};
141

    
142
template <intmax_t _Xp, intmax_t _Yp>
143
class __ll_add<_Xp, _Yp, -1>
144
{
145
    static const intmax_t min = (1LL << (sizeof(intmax_t) * CHAR_BIT - 1)) + 1;
146
    static const intmax_t max = -min;
147

    
148
    static_assert(min - _Yp <= _Xp, "overflow in __ll_add");
149
public:
150
    static const intmax_t value = _Xp + _Yp;
151
};
152

    
153
template <intmax_t _Xp, intmax_t _Yp, intmax_t = __static_sign<_Yp>::value>
154
class __ll_sub;
155

    
156
template <intmax_t _Xp, intmax_t _Yp>
157
class __ll_sub<_Xp, _Yp, 1>
158
{
159
    static const intmax_t min = (1LL << (sizeof(intmax_t) * CHAR_BIT - 1)) + 1;
160
    static const intmax_t max = -min;
161

    
162
    static_assert(min + _Yp <= _Xp, "overflow in __ll_sub");
163
public:
164
    static const intmax_t value = _Xp - _Yp;
165
};
166

    
167
template <intmax_t _Xp, intmax_t _Yp>
168
class __ll_sub<_Xp, _Yp, 0>
169
{
170
public:
171
    static const intmax_t value = _Xp;
172
};
173

    
174
template <intmax_t _Xp, intmax_t _Yp>
175
class __ll_sub<_Xp, _Yp, -1>
176
{
177
    static const intmax_t min = (1LL << (sizeof(intmax_t) * CHAR_BIT - 1)) + 1;
178
    static const intmax_t max = -min;
179

    
180
    static_assert(_Xp <= max + _Yp, "overflow in __ll_sub");
181
public:
182
    static const intmax_t value = _Xp - _Yp;
183
};
184

    
185
template <intmax_t _Xp, intmax_t _Yp>
186
class __ll_mul
187
{
188
    static const intmax_t nan = (1LL << (sizeof(intmax_t) * CHAR_BIT - 1));
189
    static const intmax_t min = nan + 1;
190
    static const intmax_t max = -min;
191
    static const intmax_t __a_x = __static_abs<_Xp>::value;
192
    static const intmax_t __a_y = __static_abs<_Yp>::value;
193

    
194
    static_assert(_Xp != nan && _Yp != nan && __a_x <= max / __a_y, "overflow in __ll_mul");
195
public:
196
    static const intmax_t value = _Xp * _Yp;
197
};
198

    
199
template <intmax_t _Yp>
200
class __ll_mul<0, _Yp>
201
{
202
public:
203
    static const intmax_t value = 0;
204
};
205

    
206
template <intmax_t _Xp>
207
class __ll_mul<_Xp, 0>
208
{
209
public:
210
    static const intmax_t value = 0;
211
};
212

    
213
template <>
214
class __ll_mul<0, 0>
215
{
216
public:
217
    static const intmax_t value = 0;
218
};
219

    
220
// Not actually used but left here in case needed in future maintenance
221
template <intmax_t _Xp, intmax_t _Yp>
222
class __ll_div
223
{
224
    static const intmax_t nan = (1LL << (sizeof(intmax_t) * CHAR_BIT - 1));
225
    static const intmax_t min = nan + 1;
226
    static const intmax_t max = -min;
227

    
228
    static_assert(_Xp != nan && _Yp != nan && _Yp != 0, "overflow in __ll_div");
229
public:
230
    static const intmax_t value = _Xp / _Yp;
231
};
232

    
233
template <intmax_t _Num, intmax_t _Den = 1>
234
class _LIBCPP_TYPE_VIS_ONLY ratio
235
{
236
    static_assert(__static_abs<_Num>::value >= 0, "ratio numerator is out of range");
237
    static_assert(_Den != 0, "ratio divide by 0");
238
    static_assert(__static_abs<_Den>::value >  0, "ratio denominator is out of range");
239
    static _LIBCPP_CONSTEXPR const intmax_t __na = __static_abs<_Num>::value;
240
    static _LIBCPP_CONSTEXPR const intmax_t __da = __static_abs<_Den>::value;
241
    static _LIBCPP_CONSTEXPR const intmax_t __s = __static_sign<_Num>::value * __static_sign<_Den>::value;
242
    static _LIBCPP_CONSTEXPR const intmax_t __gcd = __static_gcd<__na, __da>::value;
243
public:
244
    static _LIBCPP_CONSTEXPR const intmax_t num = __s * __na / __gcd;
245
    static _LIBCPP_CONSTEXPR const intmax_t den = __da / __gcd;
246

    
247
    typedef ratio<num, den> type;
248
};
249

    
250
template <intmax_t _Num, intmax_t _Den>
251
_LIBCPP_CONSTEXPR const intmax_t ratio<_Num, _Den>::num;
252

    
253
template <intmax_t _Num, intmax_t _Den>
254
_LIBCPP_CONSTEXPR const intmax_t ratio<_Num, _Den>::den;
255

    
256
template <class _Tp>                    struct __is_ratio                     : false_type {};
257
template <intmax_t _Num, intmax_t _Den> struct __is_ratio<ratio<_Num, _Den> > : true_type  {};
258

    
259
typedef ratio<1LL, 1000000000000000000LL> atto;
260
typedef ratio<1LL,    1000000000000000LL> femto;
261
typedef ratio<1LL,       1000000000000LL> pico;
262
typedef ratio<1LL,          1000000000LL> nano;
263
typedef ratio<1LL,             1000000LL> micro;
264
typedef ratio<1LL,                1000LL> milli;
265
typedef ratio<1LL,                 100LL> centi;
266
typedef ratio<1LL,                  10LL> deci;
267
typedef ratio<                 10LL, 1LL> deca;
268
typedef ratio<                100LL, 1LL> hecto;
269
typedef ratio<               1000LL, 1LL> kilo;
270
typedef ratio<            1000000LL, 1LL> mega;
271
typedef ratio<         1000000000LL, 1LL> giga;
272
typedef ratio<      1000000000000LL, 1LL> tera;
273
typedef ratio<   1000000000000000LL, 1LL> peta;
274
typedef ratio<1000000000000000000LL, 1LL> exa;
275

    
276
template <class _R1, class _R2>
277
struct __ratio_multiply
278
{
279
private:
280
    static const intmax_t __gcd_n1_d2 = __static_gcd<_R1::num, _R2::den>::value;
281
    static const intmax_t __gcd_d1_n2 = __static_gcd<_R1::den, _R2::num>::value;
282
public:
283
    typedef typename ratio
284
        <
285
            __ll_mul<_R1::num / __gcd_n1_d2, _R2::num / __gcd_d1_n2>::value,
286
            __ll_mul<_R2::den / __gcd_n1_d2, _R1::den / __gcd_d1_n2>::value
287
        >::type type;
288
};
289

    
290
#ifndef _LIBCPP_HAS_NO_TEMPLATE_ALIASES
291

    
292
template <class _R1, class _R2> using ratio_multiply
293
                                    = typename __ratio_multiply<_R1, _R2>::type;
294

    
295
#else  // _LIBCPP_HAS_NO_TEMPLATE_ALIASES
296

    
297
template <class _R1, class _R2>
298
struct _LIBCPP_TYPE_VIS_ONLY ratio_multiply
299
    : public __ratio_multiply<_R1, _R2>::type {};
300

    
301
#endif  // _LIBCPP_HAS_NO_TEMPLATE_ALIASES
302

    
303
template <class _R1, class _R2>
304
struct __ratio_divide
305
{
306
private:
307
    static const intmax_t __gcd_n1_n2 = __static_gcd<_R1::num, _R2::num>::value;
308
    static const intmax_t __gcd_d1_d2 = __static_gcd<_R1::den, _R2::den>::value;
309
public:
310
    typedef typename ratio
311
        <
312
            __ll_mul<_R1::num / __gcd_n1_n2, _R2::den / __gcd_d1_d2>::value,
313
            __ll_mul<_R2::num / __gcd_n1_n2, _R1::den / __gcd_d1_d2>::value
314
        >::type type;
315
};
316

    
317
#ifndef _LIBCPP_HAS_NO_TEMPLATE_ALIASES
318

    
319
template <class _R1, class _R2> using ratio_divide
320
                                      = typename __ratio_divide<_R1, _R2>::type;
321

    
322
#else  // _LIBCPP_HAS_NO_TEMPLATE_ALIASES
323

    
324
template <class _R1, class _R2>
325
struct _LIBCPP_TYPE_VIS_ONLY ratio_divide
326
    : public __ratio_divide<_R1, _R2>::type {};
327

    
328
#endif  // _LIBCPP_HAS_NO_TEMPLATE_ALIASES
329

    
330
template <class _R1, class _R2>
331
struct __ratio_add
332
{
333
private:
334
    static const intmax_t __gcd_n1_n2 = __static_gcd<_R1::num, _R2::num>::value;
335
    static const intmax_t __gcd_d1_d2 = __static_gcd<_R1::den, _R2::den>::value;
336
public:
337
    typedef typename ratio_multiply
338
        <
339
            ratio<__gcd_n1_n2, _R1::den / __gcd_d1_d2>,
340
            ratio
341
            <
342
                __ll_add
343
                <
344
                    __ll_mul<_R1::num / __gcd_n1_n2, _R2::den / __gcd_d1_d2>::value,
345
                    __ll_mul<_R2::num / __gcd_n1_n2, _R1::den / __gcd_d1_d2>::value
346
                >::value,
347
                _R2::den
348
            >
349
        >::type type;
350
};
351

    
352
#ifndef _LIBCPP_HAS_NO_TEMPLATE_ALIASES
353

    
354
template <class _R1, class _R2> using ratio_add
355
                                         = typename __ratio_add<_R1, _R2>::type;
356

    
357
#else  // _LIBCPP_HAS_NO_TEMPLATE_ALIASES
358

    
359
template <class _R1, class _R2>
360
struct _LIBCPP_TYPE_VIS_ONLY ratio_add
361
    : public __ratio_add<_R1, _R2>::type {};
362

    
363
#endif  // _LIBCPP_HAS_NO_TEMPLATE_ALIASES
364

    
365
template <class _R1, class _R2>
366
struct __ratio_subtract
367
{
368
private:
369
    static const intmax_t __gcd_n1_n2 = __static_gcd<_R1::num, _R2::num>::value;
370
    static const intmax_t __gcd_d1_d2 = __static_gcd<_R1::den, _R2::den>::value;
371
public:
372
    typedef typename ratio_multiply
373
        <
374
            ratio<__gcd_n1_n2, _R1::den / __gcd_d1_d2>,
375
            ratio
376
            <
377
                __ll_sub
378
                <
379
                    __ll_mul<_R1::num / __gcd_n1_n2, _R2::den / __gcd_d1_d2>::value,
380
                    __ll_mul<_R2::num / __gcd_n1_n2, _R1::den / __gcd_d1_d2>::value
381
                >::value,
382
                _R2::den
383
            >
384
        >::type type;
385
};
386

    
387
#ifndef _LIBCPP_HAS_NO_TEMPLATE_ALIASES
388

    
389
template <class _R1, class _R2> using ratio_subtract
390
                                    = typename __ratio_subtract<_R1, _R2>::type;
391

    
392
#else  // _LIBCPP_HAS_NO_TEMPLATE_ALIASES
393

    
394
template <class _R1, class _R2>
395
struct _LIBCPP_TYPE_VIS_ONLY ratio_subtract
396
    : public __ratio_subtract<_R1, _R2>::type {};
397

    
398
#endif  // _LIBCPP_HAS_NO_TEMPLATE_ALIASES
399

    
400
// ratio_equal
401

    
402
template <class _R1, class _R2>
403
struct _LIBCPP_TYPE_VIS_ONLY ratio_equal
404
    : public _LIBCPP_BOOL_CONSTANT((_R1::num == _R2::num && _R1::den == _R2::den)) {};
405

    
406
template <class _R1, class _R2>
407
struct _LIBCPP_TYPE_VIS_ONLY ratio_not_equal
408
    : public _LIBCPP_BOOL_CONSTANT((!ratio_equal<_R1, _R2>::value)) {};
409

    
410
// ratio_less
411

    
412
template <class _R1, class _R2, bool _Odd = false,
413
          intmax_t _Q1 = _R1::num / _R1::den, intmax_t _M1 = _R1::num % _R1::den,
414
          intmax_t _Q2 = _R2::num / _R2::den, intmax_t _M2 = _R2::num % _R2::den>
415
struct __ratio_less1
416
{
417
    static const bool value = _Odd ? _Q2 < _Q1 : _Q1 < _Q2;
418
};
419

    
420
template <class _R1, class _R2, bool _Odd, intmax_t _Qp>
421
struct __ratio_less1<_R1, _R2, _Odd, _Qp, 0, _Qp, 0>
422
{
423
    static const bool value = false;
424
};
425

    
426
template <class _R1, class _R2, bool _Odd, intmax_t _Qp, intmax_t _M2>
427
struct __ratio_less1<_R1, _R2, _Odd, _Qp, 0, _Qp, _M2>
428
{
429
    static const bool value = !_Odd;
430
};
431

    
432
template <class _R1, class _R2, bool _Odd, intmax_t _Qp, intmax_t _M1>
433
struct __ratio_less1<_R1, _R2, _Odd, _Qp, _M1, _Qp, 0>
434
{
435
    static const bool value = _Odd;
436
};
437

    
438
template <class _R1, class _R2, bool _Odd, intmax_t _Qp, intmax_t _M1,
439
                                                        intmax_t _M2>
440
struct __ratio_less1<_R1, _R2, _Odd, _Qp, _M1, _Qp, _M2>
441
{
442
    static const bool value = __ratio_less1<ratio<_R1::den, _M1>,
443
                                            ratio<_R2::den, _M2>, !_Odd>::value;
444
};
445

    
446
template <class _R1, class _R2, intmax_t _S1 = __static_sign<_R1::num>::value,
447
                                intmax_t _S2 = __static_sign<_R2::num>::value>
448
struct __ratio_less
449
{
450
    static const bool value = _S1 < _S2;
451
};
452

    
453
template <class _R1, class _R2>
454
struct __ratio_less<_R1, _R2, 1LL, 1LL>
455
{
456
    static const bool value = __ratio_less1<_R1, _R2>::value;
457
};
458

    
459
template <class _R1, class _R2>
460
struct __ratio_less<_R1, _R2, -1LL, -1LL>
461
{
462
    static const bool value = __ratio_less1<ratio<-_R2::num, _R2::den>, ratio<-_R1::num, _R1::den> >::value;
463
};
464

    
465
template <class _R1, class _R2>
466
struct _LIBCPP_TYPE_VIS_ONLY ratio_less
467
    : public _LIBCPP_BOOL_CONSTANT((__ratio_less<_R1, _R2>::value)) {};
468

    
469
template <class _R1, class _R2>
470
struct _LIBCPP_TYPE_VIS_ONLY ratio_less_equal
471
    : public _LIBCPP_BOOL_CONSTANT((!ratio_less<_R2, _R1>::value)) {};
472

    
473
template <class _R1, class _R2>
474
struct _LIBCPP_TYPE_VIS_ONLY ratio_greater
475
    : public _LIBCPP_BOOL_CONSTANT((ratio_less<_R2, _R1>::value)) {};
476

    
477
template <class _R1, class _R2>
478
struct _LIBCPP_TYPE_VIS_ONLY ratio_greater_equal
479
    : public _LIBCPP_BOOL_CONSTANT((!ratio_less<_R1, _R2>::value)) {};
480

    
481
template <class _R1, class _R2>
482
struct __ratio_gcd
483
{
484
    typedef ratio<__static_gcd<_R1::num, _R2::num>::value,
485
                  __static_lcm<_R1::den, _R2::den>::value> type;
486
};
487

    
488
_LIBCPP_END_NAMESPACE_STD
489

    
490
#endif  // _LIBCPP_RATIO