1// -*- C++ -*-
2//===---------------------------- ratio -----------------------------------===//
3//
4// Part of the LLVM Project, under the Apache License v2.0 with LLVM Exceptions.
5// See https://llvm.org/LICENSE.txt for license information.
6// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
7//
8//===----------------------------------------------------------------------===//
9
10#ifndef _LIBCPP_RATIO
11#define _LIBCPP_RATIO
12
13/*
14 ratio synopsis
15
16namespace std
17{
18
19template <intmax_t N, intmax_t D = 1>
20class ratio
21{
22public:
23 static constexpr intmax_t num;
24 static constexpr intmax_t den;
25 typedef ratio<num, den> type;
26};
27
28// ratio arithmetic
29template <class R1, class R2> using ratio_add = ...;
30template <class R1, class R2> using ratio_subtract = ...;
31template <class R1, class R2> using ratio_multiply = ...;
32template <class R1, class R2> using ratio_divide = ...;
33
34// ratio comparison
35template <class R1, class R2> struct ratio_equal;
36template <class R1, class R2> struct ratio_not_equal;
37template <class R1, class R2> struct ratio_less;
38template <class R1, class R2> struct ratio_less_equal;
39template <class R1, class R2> struct ratio_greater;
40template <class R1, class R2> struct ratio_greater_equal;
41
42// convenience SI typedefs
43typedef ratio<1, 1000000000000000000000000> yocto; // not supported
44typedef ratio<1, 1000000000000000000000> zepto; // not supported
45typedef ratio<1, 1000000000000000000> atto;
46typedef ratio<1, 1000000000000000> femto;
47typedef ratio<1, 1000000000000> pico;
48typedef ratio<1, 1000000000> nano;
49typedef ratio<1, 1000000> micro;
50typedef ratio<1, 1000> milli;
51typedef ratio<1, 100> centi;
52typedef ratio<1, 10> deci;
53typedef ratio< 10, 1> deca;
54typedef ratio< 100, 1> hecto;
55typedef ratio< 1000, 1> kilo;
56typedef ratio< 1000000, 1> mega;
57typedef ratio< 1000000000, 1> giga;
58typedef ratio< 1000000000000, 1> tera;
59typedef ratio< 1000000000000000, 1> peta;
60typedef ratio< 1000000000000000000, 1> exa;
61typedef ratio< 1000000000000000000000, 1> zetta; // not supported
62typedef ratio<1000000000000000000000000, 1> yotta; // not supported
63
64 // 20.11.5, ratio comparison
65 template <class R1, class R2> inline constexpr bool ratio_equal_v
66 = ratio_equal<R1, R2>::value; // C++17
67 template <class R1, class R2> inline constexpr bool ratio_not_equal_v
68 = ratio_not_equal<R1, R2>::value; // C++17
69 template <class R1, class R2> inline constexpr bool ratio_less_v
70 = ratio_less<R1, R2>::value; // C++17
71 template <class R1, class R2> inline constexpr bool ratio_less_equal_v
72 = ratio_less_equal<R1, R2>::value; // C++17
73 template <class R1, class R2> inline constexpr bool ratio_greater_v
74 = ratio_greater<R1, R2>::value; // C++17
75 template <class R1, class R2> inline constexpr bool ratio_greater_equal_v
76 = ratio_greater_equal<R1, R2>::value; // C++17
77}
78*/
79
80#include <__config>
81#include <cstdint>
82#include <climits>
83#include <type_traits>
84
85#if !defined(_LIBCPP_HAS_NO_PRAGMA_SYSTEM_HEADER)
86#pragma GCC system_header
87#endif
88
89_LIBCPP_PUSH_MACROS
90#include <__undef_macros>
91
92
93_LIBCPP_BEGIN_NAMESPACE_STD
94
95// __static_gcd
96
97template <intmax_t _Xp, intmax_t _Yp>
98struct __static_gcd
99{
100 static const intmax_t value = __static_gcd<_Yp, _Xp % _Yp>::value;
101};
102
103template <intmax_t _Xp>
104struct __static_gcd<_Xp, 0>
105{
106 static const intmax_t value = _Xp;
107};
108
109template <>
110struct __static_gcd<0, 0>
111{
112 static const intmax_t value = 1;
113};
114
115// __static_lcm
116
117template <intmax_t _Xp, intmax_t _Yp>
118struct __static_lcm
119{
120 static const intmax_t value = _Xp / __static_gcd<_Xp, _Yp>::value * _Yp;
121};
122
123template <intmax_t _Xp>
124struct __static_abs
125{
126 static const intmax_t value = _Xp < 0 ? -_Xp : _Xp;
127};
128
129template <intmax_t _Xp>
130struct __static_sign
131{
132 static const intmax_t value = _Xp == 0 ? 0 : (_Xp < 0 ? -1 : 1);
133};
134
135template <intmax_t _Xp, intmax_t _Yp, intmax_t = __static_sign<_Yp>::value>
136class __ll_add;
137
138template <intmax_t _Xp, intmax_t _Yp>
139class __ll_add<_Xp, _Yp, 1>
140{
141 static const intmax_t min = (1LL << (sizeof(intmax_t) * CHAR_BIT - 1)) + 1;
142 static const intmax_t max = -min;
143
144 static_assert(_Xp <= max - _Yp, "overflow in __ll_add");
145public:
146 static const intmax_t value = _Xp + _Yp;
147};
148
149template <intmax_t _Xp, intmax_t _Yp>
150class __ll_add<_Xp, _Yp, 0>
151{
152public:
153 static const intmax_t value = _Xp;
154};
155
156template <intmax_t _Xp, intmax_t _Yp>
157class __ll_add<_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_add");
163public:
164 static const intmax_t value = _Xp + _Yp;
165};
166
167template <intmax_t _Xp, intmax_t _Yp, intmax_t = __static_sign<_Yp>::value>
168class __ll_sub;
169
170template <intmax_t _Xp, intmax_t _Yp>
171class __ll_sub<_Xp, _Yp, 1>
172{
173 static const intmax_t min = (1LL << (sizeof(intmax_t) * CHAR_BIT - 1)) + 1;
174 static const intmax_t max = -min;
175
176 static_assert(min + _Yp <= _Xp, "overflow in __ll_sub");
177public:
178 static const intmax_t value = _Xp - _Yp;
179};
180
181template <intmax_t _Xp, intmax_t _Yp>
182class __ll_sub<_Xp, _Yp, 0>
183{
184public:
185 static const intmax_t value = _Xp;
186};
187
188template <intmax_t _Xp, intmax_t _Yp>
189class __ll_sub<_Xp, _Yp, -1>
190{
191 static const intmax_t min = (1LL << (sizeof(intmax_t) * CHAR_BIT - 1)) + 1;
192 static const intmax_t max = -min;
193
194 static_assert(_Xp <= max + _Yp, "overflow in __ll_sub");
195public:
196 static const intmax_t value = _Xp - _Yp;
197};
198
199template <intmax_t _Xp, intmax_t _Yp>
200class __ll_mul
201{
202 static const intmax_t nan = (1LL << (sizeof(intmax_t) * CHAR_BIT - 1));
203 static const intmax_t min = nan + 1;
204 static const intmax_t max = -min;
205 static const intmax_t __a_x = __static_abs<_Xp>::value;
206 static const intmax_t __a_y = __static_abs<_Yp>::value;
207
208 static_assert(_Xp != nan && _Yp != nan && __a_x <= max / __a_y, "overflow in __ll_mul");
209public:
210 static const intmax_t value = _Xp * _Yp;
211};
212
213template <intmax_t _Yp>
214class __ll_mul<0, _Yp>
215{
216public:
217 static const intmax_t value = 0;
218};
219
220template <intmax_t _Xp>
221class __ll_mul<_Xp, 0>
222{
223public:
224 static const intmax_t value = 0;
225};
226
227template <>
228class __ll_mul<0, 0>
229{
230public:
231 static const intmax_t value = 0;
232};
233
234// Not actually used but left here in case needed in future maintenance
235template <intmax_t _Xp, intmax_t _Yp>
236class __ll_div
237{
238 static const intmax_t nan = (1LL << (sizeof(intmax_t) * CHAR_BIT - 1));
239 static const intmax_t min = nan + 1;
240 static const intmax_t max = -min;
241
242 static_assert(_Xp != nan && _Yp != nan && _Yp != 0, "overflow in __ll_div");
243public:
244 static const intmax_t value = _Xp / _Yp;
245};
246
247template <intmax_t _Num, intmax_t _Den = 1>
248class _LIBCPP_TEMPLATE_VIS ratio
249{
250 static_assert(__static_abs<_Num>::value >= 0, "ratio numerator is out of range");
251 static_assert(_Den != 0, "ratio divide by 0");
252 static_assert(__static_abs<_Den>::value > 0, "ratio denominator is out of range");
253 static _LIBCPP_CONSTEXPR const intmax_t __na = __static_abs<_Num>::value;
254 static _LIBCPP_CONSTEXPR const intmax_t __da = __static_abs<_Den>::value;
255 static _LIBCPP_CONSTEXPR const intmax_t __s = __static_sign<_Num>::value * __static_sign<_Den>::value;
256 static _LIBCPP_CONSTEXPR const intmax_t __gcd = __static_gcd<__na, __da>::value;
257public:
258 static _LIBCPP_CONSTEXPR const intmax_t num = __s * __na / __gcd;
259 static _LIBCPP_CONSTEXPR const intmax_t den = __da / __gcd;
260
261 typedef ratio<num, den> type;
262};
263
264template <intmax_t _Num, intmax_t _Den>
265_LIBCPP_CONSTEXPR const intmax_t ratio<_Num, _Den>::num;
266
267template <intmax_t _Num, intmax_t _Den>
268_LIBCPP_CONSTEXPR const intmax_t ratio<_Num, _Den>::den;
269
270template <class _Tp> struct __is_ratio : false_type {};
271template <intmax_t _Num, intmax_t _Den> struct __is_ratio<ratio<_Num, _Den> > : true_type {};
272
273typedef ratio<1LL, 1000000000000000000LL> atto;
274typedef ratio<1LL, 1000000000000000LL> femto;
275typedef ratio<1LL, 1000000000000LL> pico;
276typedef ratio<1LL, 1000000000LL> nano;
277typedef ratio<1LL, 1000000LL> micro;
278typedef ratio<1LL, 1000LL> milli;
279typedef ratio<1LL, 100LL> centi;
280typedef ratio<1LL, 10LL> deci;
281typedef ratio< 10LL, 1LL> deca;
282typedef ratio< 100LL, 1LL> hecto;
283typedef ratio< 1000LL, 1LL> kilo;
284typedef ratio< 1000000LL, 1LL> mega;
285typedef ratio< 1000000000LL, 1LL> giga;
286typedef ratio< 1000000000000LL, 1LL> tera;
287typedef ratio< 1000000000000000LL, 1LL> peta;
288typedef ratio<1000000000000000000LL, 1LL> exa;
289
290template <class _R1, class _R2>
291struct __ratio_multiply
292{
293private:
294 static const intmax_t __gcd_n1_d2 = __static_gcd<_R1::num, _R2::den>::value;
295 static const intmax_t __gcd_d1_n2 = __static_gcd<_R1::den, _R2::num>::value;
296public:
297 typedef typename ratio
298 <
299 __ll_mul<_R1::num / __gcd_n1_d2, _R2::num / __gcd_d1_n2>::value,
300 __ll_mul<_R2::den / __gcd_n1_d2, _R1::den / __gcd_d1_n2>::value
301 >::type type;
302};
303
304#ifndef _LIBCPP_CXX03_LANG
305
306template <class _R1, class _R2> using ratio_multiply
307 = typename __ratio_multiply<_R1, _R2>::type;
308
309#else // _LIBCPP_CXX03_LANG
310
311template <class _R1, class _R2>
312struct _LIBCPP_TEMPLATE_VIS ratio_multiply
313 : public __ratio_multiply<_R1, _R2>::type {};
314
315#endif // _LIBCPP_CXX03_LANG
316
317template <class _R1, class _R2>
318struct __ratio_divide
319{
320private:
321 static const intmax_t __gcd_n1_n2 = __static_gcd<_R1::num, _R2::num>::value;
322 static const intmax_t __gcd_d1_d2 = __static_gcd<_R1::den, _R2::den>::value;
323public:
324 typedef typename ratio
325 <
326 __ll_mul<_R1::num / __gcd_n1_n2, _R2::den / __gcd_d1_d2>::value,
327 __ll_mul<_R2::num / __gcd_n1_n2, _R1::den / __gcd_d1_d2>::value
328 >::type type;
329};
330
331#ifndef _LIBCPP_CXX03_LANG
332
333template <class _R1, class _R2> using ratio_divide
334 = typename __ratio_divide<_R1, _R2>::type;
335
336#else // _LIBCPP_CXX03_LANG
337
338template <class _R1, class _R2>
339struct _LIBCPP_TEMPLATE_VIS ratio_divide
340 : public __ratio_divide<_R1, _R2>::type {};
341
342#endif // _LIBCPP_CXX03_LANG
343
344template <class _R1, class _R2>
345struct __ratio_add
346{
347private:
348 static const intmax_t __gcd_n1_n2 = __static_gcd<_R1::num, _R2::num>::value;
349 static const intmax_t __gcd_d1_d2 = __static_gcd<_R1::den, _R2::den>::value;
350public:
351 typedef typename ratio_multiply
352 <
353 ratio<__gcd_n1_n2, _R1::den / __gcd_d1_d2>,
354 ratio
355 <
356 __ll_add
357 <
358 __ll_mul<_R1::num / __gcd_n1_n2, _R2::den / __gcd_d1_d2>::value,
359 __ll_mul<_R2::num / __gcd_n1_n2, _R1::den / __gcd_d1_d2>::value
360 >::value,
361 _R2::den
362 >
363 >::type type;
364};
365
366#ifndef _LIBCPP_CXX03_LANG
367
368template <class _R1, class _R2> using ratio_add
369 = typename __ratio_add<_R1, _R2>::type;
370
371#else // _LIBCPP_CXX03_LANG
372
373template <class _R1, class _R2>
374struct _LIBCPP_TEMPLATE_VIS ratio_add
375 : public __ratio_add<_R1, _R2>::type {};
376
377#endif // _LIBCPP_CXX03_LANG
378
379template <class _R1, class _R2>
380struct __ratio_subtract
381{
382private:
383 static const intmax_t __gcd_n1_n2 = __static_gcd<_R1::num, _R2::num>::value;
384 static const intmax_t __gcd_d1_d2 = __static_gcd<_R1::den, _R2::den>::value;
385public:
386 typedef typename ratio_multiply
387 <
388 ratio<__gcd_n1_n2, _R1::den / __gcd_d1_d2>,
389 ratio
390 <
391 __ll_sub
392 <
393 __ll_mul<_R1::num / __gcd_n1_n2, _R2::den / __gcd_d1_d2>::value,
394 __ll_mul<_R2::num / __gcd_n1_n2, _R1::den / __gcd_d1_d2>::value
395 >::value,
396 _R2::den
397 >
398 >::type type;
399};
400
401#ifndef _LIBCPP_CXX03_LANG
402
403template <class _R1, class _R2> using ratio_subtract
404 = typename __ratio_subtract<_R1, _R2>::type;
405
406#else // _LIBCPP_CXX03_LANG
407
408template <class _R1, class _R2>
409struct _LIBCPP_TEMPLATE_VIS ratio_subtract
410 : public __ratio_subtract<_R1, _R2>::type {};
411
412#endif // _LIBCPP_CXX03_LANG
413
414// ratio_equal
415
416template <class _R1, class _R2>
417struct _LIBCPP_TEMPLATE_VIS ratio_equal
418 : public _LIBCPP_BOOL_CONSTANT((_R1::num == _R2::num && _R1::den == _R2::den)) {};
419
420template <class _R1, class _R2>
421struct _LIBCPP_TEMPLATE_VIS ratio_not_equal
422 : public _LIBCPP_BOOL_CONSTANT((!ratio_equal<_R1, _R2>::value)) {};
423
424// ratio_less
425
426template <class _R1, class _R2, bool _Odd = false,
427 intmax_t _Q1 = _R1::num / _R1::den, intmax_t _M1 = _R1::num % _R1::den,
428 intmax_t _Q2 = _R2::num / _R2::den, intmax_t _M2 = _R2::num % _R2::den>
429struct __ratio_less1
430{
431 static const bool value = _Odd ? _Q2 < _Q1 : _Q1 < _Q2;
432};
433
434template <class _R1, class _R2, bool _Odd, intmax_t _Qp>
435struct __ratio_less1<_R1, _R2, _Odd, _Qp, 0, _Qp, 0>
436{
437 static const bool value = false;
438};
439
440template <class _R1, class _R2, bool _Odd, intmax_t _Qp, intmax_t _M2>
441struct __ratio_less1<_R1, _R2, _Odd, _Qp, 0, _Qp, _M2>
442{
443 static const bool value = !_Odd;
444};
445
446template <class _R1, class _R2, bool _Odd, intmax_t _Qp, intmax_t _M1>
447struct __ratio_less1<_R1, _R2, _Odd, _Qp, _M1, _Qp, 0>
448{
449 static const bool value = _Odd;
450};
451
452template <class _R1, class _R2, bool _Odd, intmax_t _Qp, intmax_t _M1,
453 intmax_t _M2>
454struct __ratio_less1<_R1, _R2, _Odd, _Qp, _M1, _Qp, _M2>
455{
456 static const bool value = __ratio_less1<ratio<_R1::den, _M1>,
457 ratio<_R2::den, _M2>, !_Odd>::value;
458};
459
460template <class _R1, class _R2, intmax_t _S1 = __static_sign<_R1::num>::value,
461 intmax_t _S2 = __static_sign<_R2::num>::value>
462struct __ratio_less
463{
464 static const bool value = _S1 < _S2;
465};
466
467template <class _R1, class _R2>
468struct __ratio_less<_R1, _R2, 1LL, 1LL>
469{
470 static const bool value = __ratio_less1<_R1, _R2>::value;
471};
472
473template <class _R1, class _R2>
474struct __ratio_less<_R1, _R2, -1LL, -1LL>
475{
476 static const bool value = __ratio_less1<ratio<-_R2::num, _R2::den>, ratio<-_R1::num, _R1::den> >::value;
477};
478
479template <class _R1, class _R2>
480struct _LIBCPP_TEMPLATE_VIS ratio_less
481 : public _LIBCPP_BOOL_CONSTANT((__ratio_less<_R1, _R2>::value)) {};
482
483template <class _R1, class _R2>
484struct _LIBCPP_TEMPLATE_VIS ratio_less_equal
485 : public _LIBCPP_BOOL_CONSTANT((!ratio_less<_R2, _R1>::value)) {};
486
487template <class _R1, class _R2>
488struct _LIBCPP_TEMPLATE_VIS ratio_greater
489 : public _LIBCPP_BOOL_CONSTANT((ratio_less<_R2, _R1>::value)) {};
490
491template <class _R1, class _R2>
492struct _LIBCPP_TEMPLATE_VIS ratio_greater_equal
493 : public _LIBCPP_BOOL_CONSTANT((!ratio_less<_R1, _R2>::value)) {};
494
495template <class _R1, class _R2>
496struct __ratio_gcd
497{
498 typedef ratio<__static_gcd<_R1::num, _R2::num>::value,
499 __static_lcm<_R1::den, _R2::den>::value> type;
500};
501
502#if _LIBCPP_STD_VER > 14 && !defined(_LIBCPP_HAS_NO_VARIABLE_TEMPLATES)
503template <class _R1, class _R2>
504_LIBCPP_INLINE_VAR _LIBCPP_CONSTEXPR bool ratio_equal_v
505 = ratio_equal<_R1, _R2>::value;
506
507template <class _R1, class _R2>
508_LIBCPP_INLINE_VAR _LIBCPP_CONSTEXPR bool ratio_not_equal_v
509 = ratio_not_equal<_R1, _R2>::value;
510
511template <class _R1, class _R2>
512_LIBCPP_INLINE_VAR _LIBCPP_CONSTEXPR bool ratio_less_v
513 = ratio_less<_R1, _R2>::value;
514
515template <class _R1, class _R2>
516_LIBCPP_INLINE_VAR _LIBCPP_CONSTEXPR bool ratio_less_equal_v
517 = ratio_less_equal<_R1, _R2>::value;
518
519template <class _R1, class _R2>
520_LIBCPP_INLINE_VAR _LIBCPP_CONSTEXPR bool ratio_greater_v
521 = ratio_greater<_R1, _R2>::value;
522
523template <class _R1, class _R2>
524_LIBCPP_INLINE_VAR _LIBCPP_CONSTEXPR bool ratio_greater_equal_v
525 = ratio_greater_equal<_R1, _R2>::value;
526#endif
527
528_LIBCPP_END_NAMESPACE_STD
529
530_LIBCPP_POP_MACROS
531
532#endif // _LIBCPP_RATIO
533