1 files changed, 41 insertions, 0 deletions
diff --git a/test/ec/perf_mult.py b/test/ec/perf_mult.py
new file mode 100755
index 0000000..2ec82b0
--- /dev/null
+++ b/test/ec/perf_mult.py
@@ -0,0 +1,41 @@
+#!/usr/bin/env python
+import click
+
+from pyecsca.ec.mod import has_gmp
+from pyecsca.ec.mult import LTRMultiplier
+from pyecsca.ec.params import get_params
+from pyecsca.misc.cfg import TemporaryConfig
+from utils import Profiler
+
+
+@click.command()
+@click.option("-p", "--profiler", type=click.Choice(("py", "c")), default="py")
+@click.option("-m", "--mod", type=click.Choice(("python", "gmp")), default="gmp" if has_gmp else "python")
+@click.option("-o", "--operations", type=click.INT, default=50)
+@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
+        p256 = get_params("secg", "secp256r1", "projective")
+        coords = p256.curve.coordinate_model
+        add = coords.formulas["add-2016-rcb"]
+        dbl = coords.formulas["dbl-2016-rcb"]
+        mult = LTRMultiplier(add, dbl)
+        click.echo(f"Profiling {operations} {p256.curve.prime.bit_length()}-bit scalar multiplication executions...")
+        one_point = p256.generator
+        with Profiler(profiler, directory, f"mult_ltr_rcb_p256_{operations}_{mod}"):
+            for _ in range(operations):
+                mult.init(p256, one_point)
+                one_point = mult.multiply(0x71a55e0c1abb3a0e069419e0f837bc195f1b9545e69fc51e53c4d48d7fea3b1a)
+        # ed25519 = get_params("other", "Ed25519", "extended")
+        # ecoords = ed25519.curve.coordinate_model
+        # dblg = ecoords.formulas["mdbl-2008-hwcd"]
+        # click.echo(f"Profiling {operations} {ed25519.curve.prime.bit_length()}-bit doubling formula executions (with assumption)...")
+        # eone_point = ed25519.generator
+        # with Profiler(profiler) as pr:
+        #     for _ in range(operations):
+        #         dblg(ed25519.curve.prime, eone_point, **ed25519.curve.parameters)
+
+
+if __name__ == "__main__":
+    main()