big_pie 0.1.1 → 0.1.2
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- checksums.yaml +4 -4
- data/ext/bigpie.c +10 -31
- data/ext/extconf.rb +1 -1
- data/lib/big_pie.rb +4 -1
- data/lib/big_pie/version.rb +1 -1
- metadata +2 -2
checksums.yaml
CHANGED
|
@@ -1,7 +1,7 @@
|
|
|
1
1
|
---
|
|
2
2
|
SHA256:
|
|
3
|
-
metadata.gz:
|
|
4
|
-
data.tar.gz:
|
|
3
|
+
metadata.gz: 45a9f7bdb31f20d220aba497b88d80cb0bd74b5395db56818b36c133619bb222
|
|
4
|
+
data.tar.gz: 583aa41083110f3860cafe0364638baac775b3f472c911c663c86ebb68b36050
|
|
5
5
|
SHA512:
|
|
6
|
-
metadata.gz:
|
|
7
|
-
data.tar.gz:
|
|
6
|
+
metadata.gz: 817f9e39c9a346ab5e21d5dee883a6e1c62c2a9f0bce5ea046981ad4ffc126175339a2c5e44f035a654a635b11902535df5226735e0917c537b1bae65c4cd010
|
|
7
|
+
data.tar.gz: bbec218860433997d5bea3037b2cf103ab09a8599c7f40974a21a2d541059a036e955aecb097e7597236de7a8efa561773e6b504a1f8ad17d6a3d922e5ea2564
|
data/ext/bigpie.c
CHANGED
|
@@ -12,7 +12,7 @@
|
|
|
12
12
|
#pragma intel optimization_level 3
|
|
13
13
|
#endif
|
|
14
14
|
|
|
15
|
-
VALUE calculatePi(volatile VALUE obj, VALUE number) {
|
|
15
|
+
VALUE calculatePi(volatile VALUE obj, volatile VALUE number) {
|
|
16
16
|
VALUE ary = rb_ary_new() ;
|
|
17
17
|
VALUE shove = rb_intern("<<") ;
|
|
18
18
|
|
|
@@ -21,10 +21,9 @@ VALUE calculatePi(volatile VALUE obj, VALUE number) {
|
|
|
21
21
|
|
|
22
22
|
register uint64_t num = NUM2ULL(number) ;
|
|
23
23
|
register uint64_t index = 0 ;
|
|
24
|
-
register
|
|
24
|
+
register int8_t comp ;
|
|
25
25
|
|
|
26
26
|
mpz_t q, t, k, m, x, r ;
|
|
27
|
-
mpz_t m_t ;
|
|
28
27
|
mpz_t temp1 ;
|
|
29
28
|
|
|
30
29
|
mpz_init(q) ;
|
|
@@ -33,7 +32,6 @@ VALUE calculatePi(volatile VALUE obj, VALUE number) {
|
|
|
33
32
|
mpz_init(m) ;
|
|
34
33
|
mpz_init(x) ;
|
|
35
34
|
mpz_init(r) ;
|
|
36
|
-
mpz_init(m_t) ;
|
|
37
35
|
mpz_init(temp1) ;
|
|
38
36
|
|
|
39
37
|
mpz_set_ui(q, 1) ;
|
|
@@ -44,51 +42,33 @@ VALUE calculatePi(volatile VALUE obj, VALUE number) {
|
|
|
44
42
|
mpz_set_ui(r, 0) ;
|
|
45
43
|
|
|
46
44
|
while(index < num) {
|
|
47
|
-
|
|
45
|
+
mpz_sub(temp1, r, t) ;
|
|
48
46
|
mpz_addmul_ui(temp1, q, 4) ;
|
|
49
|
-
|
|
50
|
-
mpz_mul(m_t, m, t) ;
|
|
47
|
+
mpz_submul(temp1, m, t) ;
|
|
51
48
|
|
|
52
|
-
comp =
|
|
53
|
-
|
|
49
|
+
comp = mpz_cmp_ui(temp1, 0) ;
|
|
50
|
+
|
|
51
|
+
if(comp < 0) {
|
|
54
52
|
++index ;
|
|
55
53
|
VALUE m_to_ui = INT2FIX(mpz_get_ui(m)) ;
|
|
56
54
|
rb_funcallv_public(ary, shove, 1, &m_to_ui) ;
|
|
57
55
|
|
|
58
|
-
//m
|
|
59
|
-
mpz_mul_ui(temp1, q, 3) ;
|
|
60
|
-
mpz_addmul_ui(temp1, r, 10) ;
|
|
61
|
-
mpz_tdiv_q(temp1, temp1, t) ;
|
|
62
|
-
mpz_submul_ui(temp1, m, 10) ;
|
|
63
|
-
mpz_set(m, temp1) ;
|
|
64
|
-
|
|
65
56
|
// r
|
|
66
|
-
|
|
57
|
+
mpz_submul(r, m, t) ;
|
|
67
58
|
mpz_mul_ui(r, r, 10) ;
|
|
68
59
|
|
|
69
|
-
|
|
60
|
+
mpz_submul_ui(m, m, 10) ;
|
|
70
61
|
mpz_mul_ui(q, q, 10) ;
|
|
71
62
|
} else {
|
|
72
|
-
// t
|
|
73
63
|
mpz_mul(t, t, x) ;
|
|
74
64
|
|
|
75
|
-
// m
|
|
76
|
-
mpz_mul_ui(temp1, k, 7) ;
|
|
77
|
-
mpz_addmul_ui(temp1, q, 2) ;
|
|
78
|
-
mpz_addmul(temp1, r, x) ;
|
|
79
|
-
mpz_tdiv_q(m, temp1, t) ;
|
|
80
|
-
|
|
81
65
|
// r
|
|
82
66
|
mpz_addmul_ui(r, q, 2) ;
|
|
83
67
|
mpz_mul(r, r, x) ;
|
|
84
68
|
|
|
85
|
-
|
|
69
|
+
mpz_tdiv_q(m, r, t) ;
|
|
86
70
|
mpz_mul(q, q, k) ;
|
|
87
|
-
|
|
88
|
-
// k
|
|
89
71
|
mpz_add_ui(k, k, 1) ;
|
|
90
|
-
|
|
91
|
-
// x
|
|
92
72
|
mpz_add_ui(x, x, 2) ;
|
|
93
73
|
}
|
|
94
74
|
}
|
|
@@ -99,7 +79,6 @@ VALUE calculatePi(volatile VALUE obj, VALUE number) {
|
|
|
99
79
|
mpz_clear(m) ;
|
|
100
80
|
mpz_clear(x) ;
|
|
101
81
|
mpz_clear(r) ;
|
|
102
|
-
mpz_clear(m_t) ;
|
|
103
82
|
mpz_clear(temp1) ;
|
|
104
83
|
|
|
105
84
|
return ary ;
|
data/ext/extconf.rb
CHANGED
|
@@ -2,7 +2,7 @@ require 'mkmf'
|
|
|
2
2
|
|
|
3
3
|
abort "\e[38;2;255;255;0m\u2B21 You don't have gmp installed! Please follow the installation part on the documentation...\e[0m" unless have_header("gmp.h")
|
|
4
4
|
|
|
5
|
-
$CFLAGS = "-
|
|
5
|
+
$CFLAGS = "-Wall -Ofast -march=native -mtune=native -lgmp"
|
|
6
6
|
puts "\e[1;38;2;255;80;80m\u2B22 CFLAGS: #{$CFLAGS}\e[0m"
|
|
7
7
|
puts "\e[1;38;2;80;80;255m\u2B22 LDFLAGS: #{$LDFLAGS}\e[0m"
|
|
8
8
|
|
data/lib/big_pie.rb
CHANGED
data/lib/big_pie/version.rb
CHANGED
metadata
CHANGED
|
@@ -1,14 +1,14 @@
|
|
|
1
1
|
--- !ruby/object:Gem::Specification
|
|
2
2
|
name: big_pie
|
|
3
3
|
version: !ruby/object:Gem::Version
|
|
4
|
-
version: 0.1.
|
|
4
|
+
version: 0.1.2
|
|
5
5
|
platform: ruby
|
|
6
6
|
authors:
|
|
7
7
|
- Sourav Goswami
|
|
8
8
|
autorequire:
|
|
9
9
|
bindir: bin
|
|
10
10
|
cert_chain: []
|
|
11
|
-
date: 2021-01-
|
|
11
|
+
date: 2021-01-03 00:00:00.000000000 Z
|
|
12
12
|
dependencies: []
|
|
13
13
|
description: Calculate N Digits of Pi, argument upto unsigned long long (generally
|
|
14
14
|
2 ** 64 - 1)
|