ruby_nacl 0.1.0

data/ext/ruby_nacl/NaCl/crypto_scalarmult/curve25519/athlon/mult.s

@@ -0,0 +1,410 @@
+.text
+.p2align 5
+.globl _crypto_scalarmult_curve25519_athlon_mult
+.globl crypto_scalarmult_curve25519_athlon_mult
+_crypto_scalarmult_curve25519_athlon_mult:
+crypto_scalarmult_curve25519_athlon_mult:
+mov %esp,%eax
+and $31,%eax
+add $32,%eax
+sub %eax,%esp
+movl %ebp,0(%esp)
+movl 4(%esp,%eax),%ecx
+movl 8(%esp,%eax),%edx
+movl 12(%esp,%eax),%ebp
+fldl 72(%edx)
+fmull 72(%ebp)
+fmull crypto_scalarmult_curve25519_athlon_scale
+fldl 0(%edx)
+fmull 64(%ebp)
+faddp %st(0),%st(1)
+fldl 8(%edx)
+fmull 56(%ebp)
+faddp %st(0),%st(1)
+fldl 0(%edx)
+fmull 72(%ebp)
+fldl 16(%edx)
+fmull 48(%ebp)
+faddp %st(0),%st(2)
+fldl 8(%edx)
+fmull 64(%ebp)
+faddp %st(0),%st(1)
+fldl 24(%edx)
+fmull 40(%ebp)
+faddp %st(0),%st(2)
+fldl 16(%edx)
+fmull 56(%ebp)
+faddp %st(0),%st(1)
+fldl 32(%edx)
+fmull 32(%ebp)
+faddp %st(0),%st(2)
+fldl 24(%edx)
+fmull 48(%ebp)
+faddp %st(0),%st(1)
+fldl 8(%edx)
+fmull 72(%ebp)
+fldl 40(%edx)
+fmull 24(%ebp)
+faddp %st(0),%st(3)
+fldl 32(%edx)
+fmull 40(%ebp)
+faddp %st(0),%st(2)
+fldl 16(%edx)
+fmull 64(%ebp)
+faddp %st(0),%st(1)
+fldl 48(%edx)
+fmull 16(%ebp)
+faddp %st(0),%st(3)
+fldl 40(%edx)
+fmull 32(%ebp)
+faddp %st(0),%st(2)
+fldl 24(%edx)
+fmull 56(%ebp)
+faddp %st(0),%st(1)
+fldl 56(%edx)
+fmull 8(%ebp)
+faddp %st(0),%st(3)
+fldl 48(%edx)
+fmull 24(%ebp)
+faddp %st(0),%st(2)
+fldl 32(%edx)
+fmull 48(%ebp)
+faddp %st(0),%st(1)
+fldl 64(%edx)
+fmull 0(%ebp)
+faddp %st(0),%st(3)
+fldl 56(%edx)
+fmull 16(%ebp)
+faddp %st(0),%st(2)
+fldl 40(%edx)
+fmull 40(%ebp)
+faddp %st(0),%st(1)
+fldl 16(%edx)
+fmull 72(%ebp)
+fldl 48(%edx)
+fmull 32(%ebp)
+faddp %st(0),%st(2)
+fldl 64(%edx)
+fmull 8(%ebp)
+faddp %st(0),%st(3)
+fldl crypto_scalarmult_curve25519_athlon_alpha230
+fadd %st(4),%st(0)
+fldl 24(%edx)
+fmull 64(%ebp)
+faddp %st(0),%st(2)
+fldl 56(%edx)
+fmull 24(%ebp)
+faddp %st(0),%st(3)
+fldl 72(%edx)
+fmull 0(%ebp)
+faddp %st(0),%st(4)
+fsubl crypto_scalarmult_curve25519_athlon_alpha230
+fldl 32(%edx)
+fmull 56(%ebp)
+faddp %st(0),%st(2)
+fldl 24(%edx)
+fmull 72(%ebp)
+fldl 64(%edx)
+fmull 16(%ebp)
+faddp %st(0),%st(4)
+fxch %st(1)
+fadd %st(0),%st(4)
+fldl 40(%edx)
+fmull 48(%ebp)
+faddp %st(0),%st(3)
+fldl 32(%edx)
+fmull 64(%ebp)
+faddp %st(0),%st(2)
+fsubrp %st(0),%st(5)
+fldl crypto_scalarmult_curve25519_athlon_alpha255
+fadd %st(4),%st(0)
+fldl 48(%edx)
+fmull 40(%ebp)
+faddp %st(0),%st(3)
+fldl 40(%edx)
+fmull 56(%ebp)
+faddp %st(0),%st(2)
+fldl 72(%edx)
+fmull 8(%ebp)
+faddp %st(0),%st(4)
+fsubl crypto_scalarmult_curve25519_athlon_alpha255
+fldl 56(%edx)
+fmull 32(%ebp)
+faddp %st(0),%st(3)
+fldl 32(%edx)
+fmull 72(%ebp)
+fldl 48(%edx)
+fmull 48(%ebp)
+faddp %st(0),%st(3)
+fldl 64(%edx)
+fmull 24(%ebp)
+faddp %st(0),%st(4)
+fxch %st(1)
+fadd %st(0),%st(4)
+fsubrp %st(0),%st(5)
+fxch %st(5)
+fstpl 64(%ecx)
+fldl 40(%edx)
+fmull 64(%ebp)
+faddp %st(0),%st(5)
+fldl 56(%edx)
+fmull 40(%ebp)
+faddp %st(0),%st(1)
+fldl 72(%edx)
+fmull 16(%ebp)
+faddp %st(0),%st(2)
+fxch %st(2)
+fmull crypto_scalarmult_curve25519_athlon_scale
+fldl 48(%edx)
+fmull 56(%ebp)
+faddp %st(0),%st(5)
+fldl 64(%edx)
+fmull 32(%ebp)
+faddp %st(0),%st(3)
+fxch %st(1)
+fmull crypto_scalarmult_curve25519_athlon_scale
+fldl 0(%edx)
+fmull 0(%ebp)
+faddp %st(0),%st(2)
+fxch %st(3)
+fstpl 72(%ecx)
+fldl 40(%edx)
+fmull 72(%ebp)
+fldl 56(%edx)
+fmull 48(%ebp)
+faddp %st(0),%st(5)
+fldl 72(%edx)
+fmull 24(%ebp)
+faddp %st(0),%st(3)
+fldl 0(%edx)
+fmull 8(%ebp)
+faddp %st(0),%st(4)
+fldl crypto_scalarmult_curve25519_athlon_alpha26
+fadd %st(2),%st(0)
+fldl 48(%edx)
+fmull 64(%ebp)
+faddp %st(0),%st(2)
+fldl 64(%edx)
+fmull 40(%ebp)
+faddp %st(0),%st(6)
+fxch %st(3)
+fmull crypto_scalarmult_curve25519_athlon_scale
+fldl 8(%edx)
+fmull 0(%ebp)
+faddp %st(0),%st(5)
+fxch %st(3)
+fsubl crypto_scalarmult_curve25519_athlon_alpha26
+fldl 56(%edx)
+fmull 56(%ebp)
+faddp %st(0),%st(2)
+fldl 72(%edx)
+fmull 32(%ebp)
+faddp %st(0),%st(6)
+fldl 0(%edx)
+fmull 16(%ebp)
+faddp %st(0),%st(4)
+fadd %st(0),%st(4)
+fsubrp %st(0),%st(2)
+fldl 48(%edx)
+fmull 72(%ebp)
+fldl 64(%edx)
+fmull 48(%ebp)
+faddp %st(0),%st(2)
+fxch %st(5)
+fmull crypto_scalarmult_curve25519_athlon_scale
+fldl 8(%edx)
+fmull 8(%ebp)
+faddp %st(0),%st(4)
+fldl crypto_scalarmult_curve25519_athlon_alpha51
+fadd %st(5),%st(0)
+fldl 56(%edx)
+fmull 64(%ebp)
+faddp %st(0),%st(7)
+fldl 72(%edx)
+fmull 40(%ebp)
+faddp %st(0),%st(3)
+fldl 0(%edx)
+fmull 24(%ebp)
+faddp %st(0),%st(2)
+fldl 16(%edx)
+fmull 0(%ebp)
+faddp %st(0),%st(5)
+fsubl crypto_scalarmult_curve25519_athlon_alpha51
+fxch %st(3)
+fstpl 0(%ecx)
+fldl 64(%edx)
+fmull 56(%ebp)
+faddp %st(0),%st(6)
+fxch %st(1)
+fmull crypto_scalarmult_curve25519_athlon_scale
+fldl 8(%edx)
+fmull 16(%ebp)
+faddp %st(0),%st(2)
+fxch %st(2)
+fadd %st(0),%st(3)
+fsubrp %st(0),%st(4)
+fldl 56(%edx)
+fmull 72(%ebp)
+fldl 72(%edx)
+fmull 48(%ebp)
+faddp %st(0),%st(6)
+fldl 0(%edx)
+fmull 32(%ebp)
+faddp %st(0),%st(3)
+fldl 16(%edx)
+fmull 8(%ebp)
+faddp %st(0),%st(2)
+fldl crypto_scalarmult_curve25519_athlon_alpha77
+fadd %st(4),%st(0)
+fldl 64(%edx)
+fmull 64(%ebp)
+faddp %st(0),%st(2)
+fxch %st(6)
+fmull crypto_scalarmult_curve25519_athlon_scale
+fldl 8(%edx)
+fmull 24(%ebp)
+faddp %st(0),%st(4)
+fldl 24(%edx)
+fmull 0(%ebp)
+faddp %st(0),%st(3)
+fxch %st(6)
+fsubl crypto_scalarmult_curve25519_athlon_alpha77
+fxch %st(5)
+fstpl 8(%ecx)
+fldl 72(%edx)
+fmull 56(%ebp)
+faddp %st(0),%st(1)
+fldl 0(%edx)
+fmull 40(%ebp)
+faddp %st(0),%st(6)
+fldl 16(%edx)
+fmull 16(%ebp)
+faddp %st(0),%st(3)
+fxch %st(4)
+fadd %st(0),%st(1)
+fsubrp %st(0),%st(3)
+fldl 64(%edx)
+fmull 72(%ebp)
+fxch %st(4)
+fmull crypto_scalarmult_curve25519_athlon_scale
+fldl 8(%edx)
+fmull 32(%ebp)
+faddp %st(0),%st(6)
+fldl 24(%edx)
+fmull 8(%ebp)
+faddp %st(0),%st(3)
+fldl crypto_scalarmult_curve25519_athlon_alpha102
+fadd %st(2),%st(0)
+fldl 72(%edx)
+fmull 64(%ebp)
+faddp %st(0),%st(6)
+fldl 0(%edx)
+fmull 48(%ebp)
+faddp %st(0),%st(2)
+fldl 16(%edx)
+fmull 24(%ebp)
+faddp %st(0),%st(7)
+fldl 32(%edx)
+fmull 0(%ebp)
+faddp %st(0),%st(4)
+fsubl crypto_scalarmult_curve25519_athlon_alpha102
+fxch %st(4)
+fstpl 16(%ecx)
+fxch %st(4)
+fmull crypto_scalarmult_curve25519_athlon_scale
+fldl 8(%edx)
+fmull 40(%ebp)
+faddp %st(0),%st(5)
+fldl 24(%edx)
+fmull 16(%ebp)
+faddp %st(0),%st(6)
+fxch %st(3)
+fadd %st(0),%st(2)
+fsubrp %st(0),%st(1)
+fldl 0(%edx)
+fmull 56(%ebp)
+faddp %st(0),%st(3)
+fldl 16(%edx)
+fmull 32(%ebp)
+faddp %st(0),%st(4)
+fldl 32(%edx)
+fmull 8(%ebp)
+faddp %st(0),%st(5)
+fldl crypto_scalarmult_curve25519_athlon_alpha128
+fadd %st(2),%st(0)
+fldl 8(%edx)
+fmull 48(%ebp)
+faddp %st(0),%st(4)
+fldl 24(%edx)
+fmull 24(%ebp)
+faddp %st(0),%st(5)
+fldl 40(%edx)
+fmull 0(%ebp)
+faddp %st(0),%st(6)
+fsubl crypto_scalarmult_curve25519_athlon_alpha128
+fxch %st(1)
+fstpl 24(%ecx)
+fldl 16(%edx)
+fmull 40(%ebp)
+faddp %st(0),%st(3)
+fldl 32(%edx)
+fmull 16(%ebp)
+faddp %st(0),%st(4)
+fadd %st(0),%st(4)
+fsubrp %st(0),%st(1)
+fstpl 32(%ecx)
+fldl 24(%edx)
+fmull 32(%ebp)
+faddp %st(0),%st(1)
+fldl 40(%edx)
+fmull 8(%ebp)
+faddp %st(0),%st(2)
+fldl crypto_scalarmult_curve25519_athlon_alpha153
+fadd %st(3),%st(0)
+fldl 32(%edx)
+fmull 24(%ebp)
+faddp %st(0),%st(2)
+fldl 48(%edx)
+fmull 0(%ebp)
+faddp %st(0),%st(3)
+fsubl crypto_scalarmult_curve25519_athlon_alpha153
+fldl 40(%edx)
+fmull 16(%ebp)
+faddp %st(0),%st(2)
+fadd %st(0),%st(2)
+fsubrp %st(0),%st(3)
+fxch %st(2)
+fstpl 40(%ecx)
+fldl 48(%edx)
+fmull 8(%ebp)
+faddp %st(0),%st(2)
+fldl crypto_scalarmult_curve25519_athlon_alpha179
+fadd %st(1),%st(0)
+fldl 56(%edx)
+fmull 0(%ebp)
+faddp %st(0),%st(3)
+fsubl crypto_scalarmult_curve25519_athlon_alpha179
+fldl 64(%ecx)
+fldl 72(%ecx)
+fxch %st(2)
+fadd %st(0),%st(4)
+fsubrp %st(0),%st(3)
+fldl crypto_scalarmult_curve25519_athlon_alpha204
+fadd %st(4),%st(0)
+fsubl crypto_scalarmult_curve25519_athlon_alpha204
+fadd %st(0),%st(1)
+fsubrp %st(0),%st(4)
+fldl crypto_scalarmult_curve25519_athlon_alpha230
+fadd %st(1),%st(0)
+fsubl crypto_scalarmult_curve25519_athlon_alpha230
+fsubr %st(0),%st(1)
+faddp %st(0),%st(2)
+fxch %st(2)
+fstpl 48(%ecx)
+fxch %st(2)
+fstpl 56(%ecx)
+fstpl 64(%ecx)
+fstpl 72(%ecx)
+movl 0(%esp),%ebp
+add %eax,%esp
+ret

data/ext/ruby_nacl/NaCl/crypto_scalarmult/curve25519/athlon/smult.c

@@ -0,0 +1,91 @@
+#include "crypto_scalarmult.h"
+
+#define mult crypto_scalarmult_curve25519_athlon_mult
+#define square crypto_scalarmult_curve25519_athlon_square
+
+void crypto_scalarmult_curve25519_athlon_recip(double out[10],const double z[10])
+{
+  double z2[10];
+  double z9[10];
+  double z11[10];
+  double z2_5_0[10];
+  double z2_10_0[10];
+  double z2_20_0[10];
+  double z2_50_0[10];
+  double z2_100_0[10];
+  double t0[10];
+  double t1[10];
+  int i;
+
+  /* 2 */ square(z2,z);
+  /* 4 */ square(t1,z2);
+  /* 8 */ square(t0,t1);
+  /* 9 */ mult(z9,t0,z);
+  /* 11 */ mult(z11,z9,z2);
+  /* 22 */ square(t0,z11);
+  /* 2^5 - 2^0 = 31 */ mult(z2_5_0,t0,z9);
+
+  /* 2^6 - 2^1 */ square(t0,z2_5_0);
+  /* 2^7 - 2^2 */ square(t1,t0);
+  /* 2^8 - 2^3 */ square(t0,t1);
+  /* 2^9 - 2^4 */ square(t1,t0);
+  /* 2^10 - 2^5 */ square(t0,t1);
+  /* 2^10 - 2^0 */ mult(z2_10_0,t0,z2_5_0);
+
+  /* 2^11 - 2^1 */ square(t0,z2_10_0);
+  /* 2^12 - 2^2 */ square(t1,t0);
+  /* 2^20 - 2^10 */ for (i = 2;i < 10;i += 2) { square(t0,t1); square(t1,t0); }
+  /* 2^20 - 2^0 */ mult(z2_20_0,t1,z2_10_0);
+
+  /* 2^21 - 2^1 */ square(t0,z2_20_0);
+  /* 2^22 - 2^2 */ square(t1,t0);
+  /* 2^40 - 2^20 */ for (i = 2;i < 20;i += 2) { square(t0,t1); square(t1,t0); }
+  /* 2^40 - 2^0 */ mult(t0,t1,z2_20_0);
+
+  /* 2^41 - 2^1 */ square(t1,t0);
+  /* 2^42 - 2^2 */ square(t0,t1);
+  /* 2^50 - 2^10 */ for (i = 2;i < 10;i += 2) { square(t1,t0); square(t0,t1); }
+  /* 2^50 - 2^0 */ mult(z2_50_0,t0,z2_10_0);
+
+  /* 2^51 - 2^1 */ square(t0,z2_50_0);
+  /* 2^52 - 2^2 */ square(t1,t0);
+  /* 2^100 - 2^50 */ for (i = 2;i < 50;i += 2) { square(t0,t1); square(t1,t0); }
+  /* 2^100 - 2^0 */ mult(z2_100_0,t1,z2_50_0);
+
+  /* 2^101 - 2^1 */ square(t1,z2_100_0);
+  /* 2^102 - 2^2 */ square(t0,t1);
+  /* 2^200 - 2^100 */ for (i = 2;i < 100;i += 2) { square(t1,t0); square(t0,t1); }
+  /* 2^200 - 2^0 */ mult(t1,t0,z2_100_0);
+
+  /* 2^201 - 2^1 */ square(t0,t1);
+  /* 2^202 - 2^2 */ square(t1,t0);
+  /* 2^250 - 2^50 */ for (i = 2;i < 50;i += 2) { square(t0,t1); square(t1,t0); }
+  /* 2^250 - 2^0 */ mult(t0,t1,z2_50_0);
+
+  /* 2^251 - 2^1 */ square(t1,t0);
+  /* 2^252 - 2^2 */ square(t0,t1);
+  /* 2^253 - 2^3 */ square(t1,t0);
+  /* 2^254 - 2^4 */ square(t0,t1);
+  /* 2^255 - 2^5 */ square(t1,t0);
+  /* 2^255 - 21 */ mult(out,t1,z11);
+}
+
+int crypto_scalarmult(unsigned char *q,
+  const unsigned char *n,
+  const unsigned char *p)
+{
+  double work[30];
+  unsigned char e[32];
+  int i;
+  for (i = 0;i < 32;++i) e[i] = n[i];
+  e[0] &= 248;
+  e[31] &= 127;
+  e[31] |= 64;
+  crypto_scalarmult_curve25519_athlon_init();
+  crypto_scalarmult_curve25519_athlon_todouble(work,p);
+  crypto_scalarmult_curve25519_athlon_mainloop(work,e);
+  crypto_scalarmult_curve25519_athlon_recip(work + 10,work + 10);
+  mult(work + 20,work,work + 10);
+  crypto_scalarmult_curve25519_athlon_fromdouble(q,work + 20);
+  return 0;
+}