| | | 1 | | using System; |
| | | 2 | | |
| | | 3 | | namespace Renci.SshNet.Security.Chaos.NaCl.Internal.Ed25519Ref10 |
| | | 4 | | { |
| | | 5 | | internal static partial class ScalarOperations |
| | | 6 | | { |
| | | 7 | | static Int64 load_3(byte[] input, int offset) |
| | 73 | 8 | | { |
| | | 9 | | Int64 result; |
| | 73 | 10 | | result = (Int64)input[offset + 0]; |
| | 73 | 11 | | result |= ((Int64)input[offset + 1]) << 8; |
| | 73 | 12 | | result |= ((Int64)input[offset + 2]) << 16; |
| | 73 | 13 | | return result; |
| | 73 | 14 | | } |
| | | 15 | | |
| | | 16 | | static Int64 load_4(byte[] input, int offset) |
| | 83 | 17 | | { |
| | | 18 | | Int64 result; |
| | 83 | 19 | | result = (Int64)input[offset + 0]; |
| | 83 | 20 | | result |= ((Int64)input[offset + 1]) << 8; |
| | 83 | 21 | | result |= ((Int64)input[offset + 2]) << 16; |
| | 83 | 22 | | result |= ((Int64)input[offset + 3]) << 24; |
| | 83 | 23 | | return result; |
| | 83 | 24 | | } |
| | | 25 | | |
| | | 26 | | /* |
| | | 27 | | Input: |
| | | 28 | | a[0]+256*a[1]+...+256^31*a[31] = a |
| | | 29 | | b[0]+256*b[1]+...+256^31*b[31] = b |
| | | 30 | | c[0]+256*c[1]+...+256^31*c[31] = c |
| | | 31 | | |
| | | 32 | | Output: |
| | | 33 | | s[0]+256*s[1]+...+256^31*s[31] = (ab+c) mod l |
| | | 34 | | where l = 2^252 + 27742317777372353535851937790883648493. |
| | | 35 | | */ |
| | | 36 | | |
| | | 37 | | internal static void sc_muladd(byte[] s, byte[] a, byte[] b, byte[] c) |
| | 1 | 38 | | { |
| | 1 | 39 | | Int64 a0 = 2097151 & load_3(a, 0); |
| | 1 | 40 | | Int64 a1 = 2097151 & (load_4(a, 2) >> 5); |
| | 1 | 41 | | Int64 a2 = 2097151 & (load_3(a, 5) >> 2); |
| | 1 | 42 | | Int64 a3 = 2097151 & (load_4(a, 7) >> 7); |
| | 1 | 43 | | Int64 a4 = 2097151 & (load_4(a, 10) >> 4); |
| | 1 | 44 | | Int64 a5 = 2097151 & (load_3(a, 13) >> 1); |
| | 1 | 45 | | Int64 a6 = 2097151 & (load_4(a, 15) >> 6); |
| | 1 | 46 | | Int64 a7 = 2097151 & (load_3(a, 18) >> 3); |
| | 1 | 47 | | Int64 a8 = 2097151 & load_3(a, 21); |
| | 1 | 48 | | Int64 a9 = 2097151 & (load_4(a, 23) >> 5); |
| | 1 | 49 | | Int64 a10 = 2097151 & (load_3(a, 26) >> 2); |
| | 1 | 50 | | Int64 a11 = (load_4(a, 28) >> 7); |
| | 1 | 51 | | Int64 b0 = 2097151 & load_3(b, 0); |
| | 1 | 52 | | Int64 b1 = 2097151 & (load_4(b, 2) >> 5); |
| | 1 | 53 | | Int64 b2 = 2097151 & (load_3(b, 5) >> 2); |
| | 1 | 54 | | Int64 b3 = 2097151 & (load_4(b, 7) >> 7); |
| | 1 | 55 | | Int64 b4 = 2097151 & (load_4(b, 10) >> 4); |
| | 1 | 56 | | Int64 b5 = 2097151 & (load_3(b, 13) >> 1); |
| | 1 | 57 | | Int64 b6 = 2097151 & (load_4(b, 15) >> 6); |
| | 1 | 58 | | Int64 b7 = 2097151 & (load_3(b, 18) >> 3); |
| | 1 | 59 | | Int64 b8 = 2097151 & load_3(b, 21); |
| | 1 | 60 | | Int64 b9 = 2097151 & (load_4(b, 23) >> 5); |
| | 1 | 61 | | Int64 b10 = 2097151 & (load_3(b, 26) >> 2); |
| | 1 | 62 | | Int64 b11 = (load_4(b, 28) >> 7); |
| | 1 | 63 | | Int64 c0 = 2097151 & load_3(c, 0); |
| | 1 | 64 | | Int64 c1 = 2097151 & (load_4(c, 2) >> 5); |
| | 1 | 65 | | Int64 c2 = 2097151 & (load_3(c, 5) >> 2); |
| | 1 | 66 | | Int64 c3 = 2097151 & (load_4(c, 7) >> 7); |
| | 1 | 67 | | Int64 c4 = 2097151 & (load_4(c, 10) >> 4); |
| | 1 | 68 | | Int64 c5 = 2097151 & (load_3(c, 13) >> 1); |
| | 1 | 69 | | Int64 c6 = 2097151 & (load_4(c, 15) >> 6); |
| | 1 | 70 | | Int64 c7 = 2097151 & (load_3(c, 18) >> 3); |
| | 1 | 71 | | Int64 c8 = 2097151 & load_3(c, 21); |
| | 1 | 72 | | Int64 c9 = 2097151 & (load_4(c, 23) >> 5); |
| | 1 | 73 | | Int64 c10 = 2097151 & (load_3(c, 26) >> 2); |
| | 1 | 74 | | Int64 c11 = (load_4(c, 28) >> 7); |
| | | 75 | | Int64 s0; |
| | | 76 | | Int64 s1; |
| | | 77 | | Int64 s2; |
| | | 78 | | Int64 s3; |
| | | 79 | | Int64 s4; |
| | | 80 | | Int64 s5; |
| | | 81 | | Int64 s6; |
| | | 82 | | Int64 s7; |
| | | 83 | | Int64 s8; |
| | | 84 | | Int64 s9; |
| | | 85 | | Int64 s10; |
| | | 86 | | Int64 s11; |
| | | 87 | | Int64 s12; |
| | | 88 | | Int64 s13; |
| | | 89 | | Int64 s14; |
| | | 90 | | Int64 s15; |
| | | 91 | | Int64 s16; |
| | | 92 | | Int64 s17; |
| | | 93 | | Int64 s18; |
| | | 94 | | Int64 s19; |
| | | 95 | | Int64 s20; |
| | | 96 | | Int64 s21; |
| | | 97 | | Int64 s22; |
| | | 98 | | Int64 s23; |
| | | 99 | | Int64 carry0; |
| | | 100 | | Int64 carry1; |
| | | 101 | | Int64 carry2; |
| | | 102 | | Int64 carry3; |
| | | 103 | | Int64 carry4; |
| | | 104 | | Int64 carry5; |
| | | 105 | | Int64 carry6; |
| | | 106 | | Int64 carry7; |
| | | 107 | | Int64 carry8; |
| | | 108 | | Int64 carry9; |
| | | 109 | | Int64 carry10; |
| | | 110 | | Int64 carry11; |
| | | 111 | | Int64 carry12; |
| | | 112 | | Int64 carry13; |
| | | 113 | | Int64 carry14; |
| | | 114 | | Int64 carry15; |
| | | 115 | | Int64 carry16; |
| | | 116 | | Int64 carry17; |
| | | 117 | | Int64 carry18; |
| | | 118 | | Int64 carry19; |
| | | 119 | | Int64 carry20; |
| | | 120 | | Int64 carry21; |
| | | 121 | | Int64 carry22; |
| | | 122 | | |
| | 1 | 123 | | s0 = c0 + a0 * b0; |
| | 1 | 124 | | s1 = c1 + a0 * b1 + a1 * b0; |
| | 1 | 125 | | s2 = c2 + a0 * b2 + a1 * b1 + a2 * b0; |
| | 1 | 126 | | s3 = c3 + a0 * b3 + a1 * b2 + a2 * b1 + a3 * b0; |
| | 1 | 127 | | s4 = c4 + a0 * b4 + a1 * b3 + a2 * b2 + a3 * b1 + a4 * b0; |
| | 1 | 128 | | s5 = c5 + a0 * b5 + a1 * b4 + a2 * b3 + a3 * b2 + a4 * b1 + a5 * b0; |
| | 1 | 129 | | s6 = c6 + a0 * b6 + a1 * b5 + a2 * b4 + a3 * b3 + a4 * b2 + a5 * b1 + a6 * b0; |
| | 1 | 130 | | s7 = c7 + a0 * b7 + a1 * b6 + a2 * b5 + a3 * b4 + a4 * b3 + a5 * b2 + a6 * b1 + a7 * b0; |
| | 1 | 131 | | s8 = c8 + a0 * b8 + a1 * b7 + a2 * b6 + a3 * b5 + a4 * b4 + a5 * b3 + a6 * b2 + a7 * b1 + a8 * b0; |
| | 1 | 132 | | s9 = c9 + a0 * b9 + a1 * b8 + a2 * b7 + a3 * b6 + a4 * b5 + a5 * b4 + a6 * b3 + a7 * b2 + a8 * b1 + a9 * b0; |
| | 1 | 133 | | s10 = c10 + a0 * b10 + a1 * b9 + a2 * b8 + a3 * b7 + a4 * b6 + a5 * b5 + a6 * b4 + a7 * b3 + a8 * b2 + a9 * |
| | 1 | 134 | | s11 = c11 + a0 * b11 + a1 * b10 + a2 * b9 + a3 * b8 + a4 * b7 + a5 * b6 + a6 * b5 + a7 * b4 + a8 * b3 + a9 * |
| | 1 | 135 | | s12 = a1 * b11 + a2 * b10 + a3 * b9 + a4 * b8 + a5 * b7 + a6 * b6 + a7 * b5 + a8 * b4 + a9 * b3 + a10 * b2 + |
| | 1 | 136 | | s13 = a2 * b11 + a3 * b10 + a4 * b9 + a5 * b8 + a6 * b7 + a7 * b6 + a8 * b5 + a9 * b4 + a10 * b3 + a11 * b2; |
| | 1 | 137 | | s14 = a3 * b11 + a4 * b10 + a5 * b9 + a6 * b8 + a7 * b7 + a8 * b6 + a9 * b5 + a10 * b4 + a11 * b3; |
| | 1 | 138 | | s15 = a4 * b11 + a5 * b10 + a6 * b9 + a7 * b8 + a8 * b7 + a9 * b6 + a10 * b5 + a11 * b4; |
| | 1 | 139 | | s16 = a5 * b11 + a6 * b10 + a7 * b9 + a8 * b8 + a9 * b7 + a10 * b6 + a11 * b5; |
| | 1 | 140 | | s17 = a6 * b11 + a7 * b10 + a8 * b9 + a9 * b8 + a10 * b7 + a11 * b6; |
| | 1 | 141 | | s18 = a7 * b11 + a8 * b10 + a9 * b9 + a10 * b8 + a11 * b7; |
| | 1 | 142 | | s19 = a8 * b11 + a9 * b10 + a10 * b9 + a11 * b8; |
| | 1 | 143 | | s20 = a9 * b11 + a10 * b10 + a11 * b9; |
| | 1 | 144 | | s21 = a10 * b11 + a11 * b10; |
| | 1 | 145 | | s22 = a11 * b11; |
| | 1 | 146 | | s23 = 0; |
| | | 147 | | |
| | 3 | 148 | | carry0 = (s0 + (1 << 20)) >> 21; s1 += carry0; s0 -= carry0 << 21; |
| | 3 | 149 | | carry2 = (s2 + (1 << 20)) >> 21; s3 += carry2; s2 -= carry2 << 21; |
| | 3 | 150 | | carry4 = (s4 + (1 << 20)) >> 21; s5 += carry4; s4 -= carry4 << 21; |
| | 3 | 151 | | carry6 = (s6 + (1 << 20)) >> 21; s7 += carry6; s6 -= carry6 << 21; |
| | 3 | 152 | | carry8 = (s8 + (1 << 20)) >> 21; s9 += carry8; s8 -= carry8 << 21; |
| | 3 | 153 | | carry10 = (s10 + (1 << 20)) >> 21; s11 += carry10; s10 -= carry10 << 21; |
| | 3 | 154 | | carry12 = (s12 + (1 << 20)) >> 21; s13 += carry12; s12 -= carry12 << 21; |
| | 3 | 155 | | carry14 = (s14 + (1 << 20)) >> 21; s15 += carry14; s14 -= carry14 << 21; |
| | 3 | 156 | | carry16 = (s16 + (1 << 20)) >> 21; s17 += carry16; s16 -= carry16 << 21; |
| | 3 | 157 | | carry18 = (s18 + (1 << 20)) >> 21; s19 += carry18; s18 -= carry18 << 21; |
| | 3 | 158 | | carry20 = (s20 + (1 << 20)) >> 21; s21 += carry20; s20 -= carry20 << 21; |
| | 3 | 159 | | carry22 = (s22 + (1 << 20)) >> 21; s23 += carry22; s22 -= carry22 << 21; |
| | | 160 | | |
| | 3 | 161 | | carry1 = (s1 + (1 << 20)) >> 21; s2 += carry1; s1 -= carry1 << 21; |
| | 3 | 162 | | carry3 = (s3 + (1 << 20)) >> 21; s4 += carry3; s3 -= carry3 << 21; |
| | 3 | 163 | | carry5 = (s5 + (1 << 20)) >> 21; s6 += carry5; s5 -= carry5 << 21; |
| | 3 | 164 | | carry7 = (s7 + (1 << 20)) >> 21; s8 += carry7; s7 -= carry7 << 21; |
| | 3 | 165 | | carry9 = (s9 + (1 << 20)) >> 21; s10 += carry9; s9 -= carry9 << 21; |
| | 3 | 166 | | carry11 = (s11 + (1 << 20)) >> 21; s12 += carry11; s11 -= carry11 << 21; |
| | 3 | 167 | | carry13 = (s13 + (1 << 20)) >> 21; s14 += carry13; s13 -= carry13 << 21; |
| | 3 | 168 | | carry15 = (s15 + (1 << 20)) >> 21; s16 += carry15; s15 -= carry15 << 21; |
| | 3 | 169 | | carry17 = (s17 + (1 << 20)) >> 21; s18 += carry17; s17 -= carry17 << 21; |
| | 3 | 170 | | carry19 = (s19 + (1 << 20)) >> 21; s20 += carry19; s19 -= carry19 << 21; |
| | 3 | 171 | | carry21 = (s21 + (1 << 20)) >> 21; s22 += carry21; s21 -= carry21 << 21; |
| | | 172 | | |
| | 1 | 173 | | s11 += s23 * 666643; |
| | 1 | 174 | | s12 += s23 * 470296; |
| | 1 | 175 | | s13 += s23 * 654183; |
| | 1 | 176 | | s14 -= s23 * 997805; |
| | 1 | 177 | | s15 += s23 * 136657; |
| | 1 | 178 | | s16 -= s23 * 683901; |
| | 1 | 179 | | s23 = 0; |
| | | 180 | | |
| | 1 | 181 | | s10 += s22 * 666643; |
| | 1 | 182 | | s11 += s22 * 470296; |
| | 1 | 183 | | s12 += s22 * 654183; |
| | 1 | 184 | | s13 -= s22 * 997805; |
| | 1 | 185 | | s14 += s22 * 136657; |
| | 1 | 186 | | s15 -= s22 * 683901; |
| | 1 | 187 | | s22 = 0; |
| | | 188 | | |
| | 1 | 189 | | s9 += s21 * 666643; |
| | 1 | 190 | | s10 += s21 * 470296; |
| | 1 | 191 | | s11 += s21 * 654183; |
| | 1 | 192 | | s12 -= s21 * 997805; |
| | 1 | 193 | | s13 += s21 * 136657; |
| | 1 | 194 | | s14 -= s21 * 683901; |
| | 1 | 195 | | s21 = 0; |
| | | 196 | | |
| | 1 | 197 | | s8 += s20 * 666643; |
| | 1 | 198 | | s9 += s20 * 470296; |
| | 1 | 199 | | s10 += s20 * 654183; |
| | 1 | 200 | | s11 -= s20 * 997805; |
| | 1 | 201 | | s12 += s20 * 136657; |
| | 1 | 202 | | s13 -= s20 * 683901; |
| | 1 | 203 | | s20 = 0; |
| | | 204 | | |
| | 1 | 205 | | s7 += s19 * 666643; |
| | 1 | 206 | | s8 += s19 * 470296; |
| | 1 | 207 | | s9 += s19 * 654183; |
| | 1 | 208 | | s10 -= s19 * 997805; |
| | 1 | 209 | | s11 += s19 * 136657; |
| | 1 | 210 | | s12 -= s19 * 683901; |
| | 1 | 211 | | s19 = 0; |
| | | 212 | | |
| | 1 | 213 | | s6 += s18 * 666643; |
| | 1 | 214 | | s7 += s18 * 470296; |
| | 1 | 215 | | s8 += s18 * 654183; |
| | 1 | 216 | | s9 -= s18 * 997805; |
| | 1 | 217 | | s10 += s18 * 136657; |
| | 1 | 218 | | s11 -= s18 * 683901; |
| | 1 | 219 | | s18 = 0; |
| | | 220 | | |
| | 3 | 221 | | carry6 = (s6 + (1 << 20)) >> 21; s7 += carry6; s6 -= carry6 << 21; |
| | 3 | 222 | | carry8 = (s8 + (1 << 20)) >> 21; s9 += carry8; s8 -= carry8 << 21; |
| | 3 | 223 | | carry10 = (s10 + (1 << 20)) >> 21; s11 += carry10; s10 -= carry10 << 21; |
| | 3 | 224 | | carry12 = (s12 + (1 << 20)) >> 21; s13 += carry12; s12 -= carry12 << 21; |
| | 3 | 225 | | carry14 = (s14 + (1 << 20)) >> 21; s15 += carry14; s14 -= carry14 << 21; |
| | 3 | 226 | | carry16 = (s16 + (1 << 20)) >> 21; s17 += carry16; s16 -= carry16 << 21; |
| | | 227 | | |
| | 3 | 228 | | carry7 = (s7 + (1 << 20)) >> 21; s8 += carry7; s7 -= carry7 << 21; |
| | 3 | 229 | | carry9 = (s9 + (1 << 20)) >> 21; s10 += carry9; s9 -= carry9 << 21; |
| | 3 | 230 | | carry11 = (s11 + (1 << 20)) >> 21; s12 += carry11; s11 -= carry11 << 21; |
| | 3 | 231 | | carry13 = (s13 + (1 << 20)) >> 21; s14 += carry13; s13 -= carry13 << 21; |
| | 3 | 232 | | carry15 = (s15 + (1 << 20)) >> 21; s16 += carry15; s15 -= carry15 << 21; |
| | | 233 | | |
| | 1 | 234 | | s5 += s17 * 666643; |
| | 1 | 235 | | s6 += s17 * 470296; |
| | 1 | 236 | | s7 += s17 * 654183; |
| | 1 | 237 | | s8 -= s17 * 997805; |
| | 1 | 238 | | s9 += s17 * 136657; |
| | 1 | 239 | | s10 -= s17 * 683901; |
| | 1 | 240 | | s17 = 0; |
| | | 241 | | |
| | 1 | 242 | | s4 += s16 * 666643; |
| | 1 | 243 | | s5 += s16 * 470296; |
| | 1 | 244 | | s6 += s16 * 654183; |
| | 1 | 245 | | s7 -= s16 * 997805; |
| | 1 | 246 | | s8 += s16 * 136657; |
| | 1 | 247 | | s9 -= s16 * 683901; |
| | 1 | 248 | | s16 = 0; |
| | | 249 | | |
| | 1 | 250 | | s3 += s15 * 666643; |
| | 1 | 251 | | s4 += s15 * 470296; |
| | 1 | 252 | | s5 += s15 * 654183; |
| | 1 | 253 | | s6 -= s15 * 997805; |
| | 1 | 254 | | s7 += s15 * 136657; |
| | 1 | 255 | | s8 -= s15 * 683901; |
| | 1 | 256 | | s15 = 0; |
| | | 257 | | |
| | 1 | 258 | | s2 += s14 * 666643; |
| | 1 | 259 | | s3 += s14 * 470296; |
| | 1 | 260 | | s4 += s14 * 654183; |
| | 1 | 261 | | s5 -= s14 * 997805; |
| | 1 | 262 | | s6 += s14 * 136657; |
| | 1 | 263 | | s7 -= s14 * 683901; |
| | 1 | 264 | | s14 = 0; |
| | | 265 | | |
| | 1 | 266 | | s1 += s13 * 666643; |
| | 1 | 267 | | s2 += s13 * 470296; |
| | 1 | 268 | | s3 += s13 * 654183; |
| | 1 | 269 | | s4 -= s13 * 997805; |
| | 1 | 270 | | s5 += s13 * 136657; |
| | 1 | 271 | | s6 -= s13 * 683901; |
| | 1 | 272 | | s13 = 0; |
| | | 273 | | |
| | 1 | 274 | | s0 += s12 * 666643; |
| | 1 | 275 | | s1 += s12 * 470296; |
| | 1 | 276 | | s2 += s12 * 654183; |
| | 1 | 277 | | s3 -= s12 * 997805; |
| | 1 | 278 | | s4 += s12 * 136657; |
| | 1 | 279 | | s5 -= s12 * 683901; |
| | 1 | 280 | | s12 = 0; |
| | | 281 | | |
| | 3 | 282 | | carry0 = (s0 + (1 << 20)) >> 21; s1 += carry0; s0 -= carry0 << 21; |
| | 3 | 283 | | carry2 = (s2 + (1 << 20)) >> 21; s3 += carry2; s2 -= carry2 << 21; |
| | 3 | 284 | | carry4 = (s4 + (1 << 20)) >> 21; s5 += carry4; s4 -= carry4 << 21; |
| | 3 | 285 | | carry6 = (s6 + (1 << 20)) >> 21; s7 += carry6; s6 -= carry6 << 21; |
| | 3 | 286 | | carry8 = (s8 + (1 << 20)) >> 21; s9 += carry8; s8 -= carry8 << 21; |
| | 3 | 287 | | carry10 = (s10 + (1 << 20)) >> 21; s11 += carry10; s10 -= carry10 << 21; |
| | | 288 | | |
| | 3 | 289 | | carry1 = (s1 + (1 << 20)) >> 21; s2 += carry1; s1 -= carry1 << 21; |
| | 3 | 290 | | carry3 = (s3 + (1 << 20)) >> 21; s4 += carry3; s3 -= carry3 << 21; |
| | 3 | 291 | | carry5 = (s5 + (1 << 20)) >> 21; s6 += carry5; s5 -= carry5 << 21; |
| | 3 | 292 | | carry7 = (s7 + (1 << 20)) >> 21; s8 += carry7; s7 -= carry7 << 21; |
| | 3 | 293 | | carry9 = (s9 + (1 << 20)) >> 21; s10 += carry9; s9 -= carry9 << 21; |
| | 3 | 294 | | carry11 = (s11 + (1 << 20)) >> 21; s12 += carry11; s11 -= carry11 << 21; |
| | | 295 | | |
| | 1 | 296 | | s0 += s12 * 666643; |
| | 1 | 297 | | s1 += s12 * 470296; |
| | 1 | 298 | | s2 += s12 * 654183; |
| | 1 | 299 | | s3 -= s12 * 997805; |
| | 1 | 300 | | s4 += s12 * 136657; |
| | 1 | 301 | | s5 -= s12 * 683901; |
| | 1 | 302 | | s12 = 0; |
| | | 303 | | |
| | 3 | 304 | | carry0 = s0 >> 21; s1 += carry0; s0 -= carry0 << 21; |
| | 3 | 305 | | carry1 = s1 >> 21; s2 += carry1; s1 -= carry1 << 21; |
| | 3 | 306 | | carry2 = s2 >> 21; s3 += carry2; s2 -= carry2 << 21; |
| | 3 | 307 | | carry3 = s3 >> 21; s4 += carry3; s3 -= carry3 << 21; |
| | 3 | 308 | | carry4 = s4 >> 21; s5 += carry4; s4 -= carry4 << 21; |
| | 3 | 309 | | carry5 = s5 >> 21; s6 += carry5; s5 -= carry5 << 21; |
| | 3 | 310 | | carry6 = s6 >> 21; s7 += carry6; s6 -= carry6 << 21; |
| | 3 | 311 | | carry7 = s7 >> 21; s8 += carry7; s7 -= carry7 << 21; |
| | 3 | 312 | | carry8 = s8 >> 21; s9 += carry8; s8 -= carry8 << 21; |
| | 3 | 313 | | carry9 = s9 >> 21; s10 += carry9; s9 -= carry9 << 21; |
| | 3 | 314 | | carry10 = s10 >> 21; s11 += carry10; s10 -= carry10 << 21; |
| | 3 | 315 | | carry11 = s11 >> 21; s12 += carry11; s11 -= carry11 << 21; |
| | | 316 | | |
| | 1 | 317 | | s0 += s12 * 666643; |
| | 1 | 318 | | s1 += s12 * 470296; |
| | 1 | 319 | | s2 += s12 * 654183; |
| | 1 | 320 | | s3 -= s12 * 997805; |
| | 1 | 321 | | s4 += s12 * 136657; |
| | 1 | 322 | | s5 -= s12 * 683901; |
| | 1 | 323 | | s12 = 0; |
| | | 324 | | |
| | 3 | 325 | | carry0 = s0 >> 21; s1 += carry0; s0 -= carry0 << 21; |
| | 3 | 326 | | carry1 = s1 >> 21; s2 += carry1; s1 -= carry1 << 21; |
| | 3 | 327 | | carry2 = s2 >> 21; s3 += carry2; s2 -= carry2 << 21; |
| | 3 | 328 | | carry3 = s3 >> 21; s4 += carry3; s3 -= carry3 << 21; |
| | 3 | 329 | | carry4 = s4 >> 21; s5 += carry4; s4 -= carry4 << 21; |
| | 3 | 330 | | carry5 = s5 >> 21; s6 += carry5; s5 -= carry5 << 21; |
| | 3 | 331 | | carry6 = s6 >> 21; s7 += carry6; s6 -= carry6 << 21; |
| | 3 | 332 | | carry7 = s7 >> 21; s8 += carry7; s7 -= carry7 << 21; |
| | 3 | 333 | | carry8 = s8 >> 21; s9 += carry8; s8 -= carry8 << 21; |
| | 3 | 334 | | carry9 = s9 >> 21; s10 += carry9; s9 -= carry9 << 21; |
| | 3 | 335 | | carry10 = s10 >> 21; s11 += carry10; s10 -= carry10 << 21; |
| | | 336 | | |
| | | 337 | | unchecked |
| | 1 | 338 | | { |
| | 1 | 339 | | s[0] = (byte)(s0 >> 0); |
| | 1 | 340 | | s[1] = (byte)(s0 >> 8); |
| | 1 | 341 | | s[2] = (byte)((s0 >> 16) | (s1 << 5)); |
| | 1 | 342 | | s[3] = (byte)(s1 >> 3); |
| | 1 | 343 | | s[4] = (byte)(s1 >> 11); |
| | 1 | 344 | | s[5] = (byte)((s1 >> 19) | (s2 << 2)); |
| | 1 | 345 | | s[6] = (byte)(s2 >> 6); |
| | 1 | 346 | | s[7] = (byte)((s2 >> 14) | (s3 << 7)); |
| | 1 | 347 | | s[8] = (byte)(s3 >> 1); |
| | 1 | 348 | | s[9] = (byte)(s3 >> 9); |
| | 1 | 349 | | s[10] = (byte)((s3 >> 17) | (s4 << 4)); |
| | 1 | 350 | | s[11] = (byte)(s4 >> 4); |
| | 1 | 351 | | s[12] = (byte)(s4 >> 12); |
| | 1 | 352 | | s[13] = (byte)((s4 >> 20) | (s5 << 1)); |
| | 1 | 353 | | s[14] = (byte)(s5 >> 7); |
| | 1 | 354 | | s[15] = (byte)((s5 >> 15) | (s6 << 6)); |
| | 1 | 355 | | s[16] = (byte)(s6 >> 2); |
| | 1 | 356 | | s[17] = (byte)(s6 >> 10); |
| | 1 | 357 | | s[18] = (byte)((s6 >> 18) | (s7 << 3)); |
| | 1 | 358 | | s[19] = (byte)(s7 >> 5); |
| | 1 | 359 | | s[20] = (byte)(s7 >> 13); |
| | 1 | 360 | | s[21] = (byte)(s8 >> 0); |
| | 1 | 361 | | s[22] = (byte)(s8 >> 8); |
| | 1 | 362 | | s[23] = (byte)((s8 >> 16) | (s9 << 5)); |
| | 1 | 363 | | s[24] = (byte)(s9 >> 3); |
| | 1 | 364 | | s[25] = (byte)(s9 >> 11); |
| | 1 | 365 | | s[26] = (byte)((s9 >> 19) | (s10 << 2)); |
| | 1 | 366 | | s[27] = (byte)(s10 >> 6); |
| | 1 | 367 | | s[28] = (byte)((s10 >> 14) | (s11 << 7)); |
| | 1 | 368 | | s[29] = (byte)(s11 >> 1); |
| | 1 | 369 | | s[30] = (byte)(s11 >> 9); |
| | 1 | 370 | | s[31] = (byte)(s11 >> 17); |
| | 1 | 371 | | } |
| | 1 | 372 | | } |
| | | 373 | | } |
| | | 374 | | } |