pd.if.org Git - zpackage/blob - crypto/ref10/fe_sub.c

   1 #include <stdint.h>
   2
   3 /*
   4 h = f - g
   5 Can overlap h with f or g.
   6
   7 Preconditions:
   8    |f| bounded by 1.1*2^25,1.1*2^24,1.1*2^25,1.1*2^24,etc.
   9    |g| bounded by 1.1*2^25,1.1*2^24,1.1*2^25,1.1*2^24,etc.
  10
  11 Postconditions:
  12    |h| bounded by 1.1*2^26,1.1*2^25,1.1*2^26,1.1*2^25,etc.
  13 */
  14
  15 void fe_sub(int32_t h[10], int32_t f[10], int32_t g[10]) {
  16         int32_t f0 = f[0];
  17         int32_t f1 = f[1];
  18         int32_t f2 = f[2];
  19         int32_t f3 = f[3];
  20         int32_t f4 = f[4];
  21         int32_t f5 = f[5];
  22         int32_t f6 = f[6];
  23         int32_t f7 = f[7];
  24         int32_t f8 = f[8];
  25         int32_t f9 = f[9];
  26         int32_t g0 = g[0];
  27         int32_t g1 = g[1];
  28         int32_t g2 = g[2];
  29         int32_t g3 = g[3];
  30         int32_t g4 = g[4];
  31         int32_t g5 = g[5];
  32         int32_t g6 = g[6];
  33         int32_t g7 = g[7];
  34         int32_t g8 = g[8];
  35         int32_t g9 = g[9];
  36         int32_t h0 = f0 - g0;
  37         int32_t h1 = f1 - g1;
  38         int32_t h2 = f2 - g2;
  39         int32_t h3 = f3 - g3;
  40         int32_t h4 = f4 - g4;
  41         int32_t h5 = f5 - g5;
  42         int32_t h6 = f6 - g6;
  43         int32_t h7 = f7 - g7;
  44         int32_t h8 = f8 - g8;
  45         int32_t h9 = f9 - g9;
  46         h[0] = h0;
  47         h[1] = h1;
  48         h[2] = h2;
  49         h[3] = h3;
  50         h[4] = h4;
  51         h[5] = h5;
  52         h[6] = h6;
  53         h[7] = h7;
  54         h[8] = h8;
  55         h[9] = h9;
  56 }