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
110
111
112
113
114
115
116
117
118
119
120
121
122
123
|
#!/usr/bin/env python
import click
from pyecsca.ec.mod.base import Mod, mod as make_mod
from pyecsca.ec.mod.flint import has_flint
from pyecsca.ec.mod.gmp import has_gmp
from pyecsca.misc.cfg import TemporaryConfig
from test.utils import Profiler
@click.command()
@click.option(
"-p",
"--profiler",
type=click.Choice(("py", "c", "raw")),
default="py",
envvar="PROF",
)
@click.option(
"-m",
"--mod",
type=click.Choice(("python", "gmp", "flint")),
default="flint" if has_flint else "gmp" if has_gmp else "python",
envvar="MOD",
)
@click.option("-o", "--operations", type=click.INT, default=100000)
@click.option(
"-d",
"--directory",
type=click.Path(file_okay=False, dir_okay=True),
default=None,
envvar="DIR",
)
def main(profiler, mod, operations, directory):
with TemporaryConfig() as cfg:
cfg.ec.mod_implementation = mod
n = 0xFFFFFFFF00000001000000000000000000000000FFFFFFFFFFFFFFFFFFFFFFFF
a = make_mod(0x11111111111111111111111111111111, n)
b = make_mod(0xBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB, n)
click.echo(f"Profiling {operations} {n.bit_length()}-bit modular inverse...")
with Profiler(
profiler, directory, f"mod_256b_inverse_{operations}_{mod}", operations
):
for _ in range(operations):
a.inverse()
click.echo(
f"Profiling {operations} {n.bit_length()}-bit modular square root..."
)
with Profiler(
profiler, directory, f"mod_256b_sqrt_{operations}_{mod}", operations
):
for _ in range(operations):
a.sqrt()
click.echo(f"Profiling {operations} {n.bit_length()}-bit modular multiply...")
c = a
with Profiler(
profiler, directory, f"mod_256b_multiply_{operations}_{mod}", operations
):
for _ in range(operations):
c = c * b
click.echo(
f"Profiling {operations} {n.bit_length()}-bit constant modular multiply..."
)
c = a
with Profiler(
profiler,
directory,
f"mod_256b_constmultiply_{operations}_{mod}",
operations,
):
for _ in range(operations):
c = c * 48006
click.echo(f"Profiling {operations} {n.bit_length()}-bit modular square...")
c = a
with Profiler(
profiler, directory, f"mod_256b_square_{operations}_{mod}", operations
):
for _ in range(operations):
c = c**2
click.echo(f"Profiling {operations} {n.bit_length()}-bit modular add...")
c = a
with Profiler(
profiler, directory, f"mod_256b_add_{operations}_{mod}", operations
):
for _ in range(operations):
c = c + b
click.echo(f"Profiling {operations} {n.bit_length()}-bit modular subtract...")
c = a
with Profiler(
profiler, directory, f"mod_256b_subtract_{operations}_{mod}", operations
):
for _ in range(operations):
c = c - b
click.echo(
f"Profiling {operations} {n.bit_length()}-bit modular quadratic residue checks..."
)
with Profiler(
profiler, directory, f"mod_256b_isresidue_{operations}_{mod}", operations
):
for _ in range(operations):
a.is_residue()
click.echo(f"Profiling {operations} {n.bit_length()}-bit modular random...")
with Profiler(
profiler, directory, f"mod_256b_random_{operations}_{mod}", operations
):
for _ in range(operations):
Mod.random(n)
click.echo(f"Profiling {operations} {n.bit_length()}-bit cubic residue checks...")
with Profiler(
profiler, directory, f"mod_256b_iscubicresidue_{operations}_{mod}", operations
):
for _ in range(operations):
a.is_cubic_residue()
click.echo(f"Profiling {operations} {n.bit_length()}-bit cubic roots...")
with Profiler(
profiler, directory, f"mod_256b_cuberoot_{operations}_{mod}", operations
):
for _ in range(operations):
b.cube_root()
if __name__ == "__main__":
main()
|
