-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathfp_x4.ptxh
65 lines (63 loc) · 2.41 KB
/
fp_x4.ptxh
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
// bls12_381: Arithmetic for BLS12-381
// Copyright 2022-2023 Dag Arne Osvik
// Copyright 2022-2023 Luan Cardoso dos Santos
#ifndef FP_X4
/**
* @brief PTX macro for multiplication by 4. Z ← 4*X
* Z and X may be the same. Modifies T6.
*
* For inclusion directly into ptx code by the C preprocessor.
*/
#define FP_X4(Z, X, T) \
/* z = x + x */ newline\
\
add.u64.cc Z##0, X##0, X##0; newline\
addc.u64.cc Z##1, X##1, X##1; newline\
addc.u64.cc Z##2, X##2, X##2; newline\
addc.u64.cc Z##3, X##3, X##3; newline\
addc.u64.cc Z##4, X##4, X##4; newline\
addc.u64.cc Z##5, X##5, X##5; newline\
addc.u32 T##6, 0, 0; newline\
\
/* z = z + z */ newline\
\
add.u64.cc Z##0, Z##0, Z##0; newline\
addc.u64.cc Z##1, Z##1, Z##1; newline\
addc.u64.cc Z##2, Z##2, Z##2; newline\
addc.u64.cc Z##3, Z##3, Z##3; newline\
addc.u64.cc Z##4, Z##4, Z##4; newline\
addc.u64.cc Z##5, Z##5, Z##5; newline\
addc.u32 T##6, T##6, T##6; newline\
\
/* if z >= 2^385 then z -= 19m */ newline\
\
setp.gt.u32 gt, T##6, 1; newline\
@gt sub.u64.cc Z##0, Z##0, 0xCDECFFFFFFF9AAB1U; newline\
@gt subc.u64.cc Z##1, Z##1, 0x46C3FFE7293BFFFAU; newline\
@gt subc.u64.cc Z##2, Z##2, 0xA89FA1F24F2244AEU; newline\
@gt subc.u64.cc Z##3, Z##3, 0x74DA9ADE12E06434U; newline\
@gt subc.u64.cc Z##4, Z##4, 0x930D7286FE9DD3FCU; newline\
@gt subc.u64.cc Z##5, Z##5, 0xEE145462447E1D73U; newline\
@gt subc.u32 T##6, T##6, 1; newline\
\
/* if z >= 2^384 then z -= 9m */ newline\
\
setp.gt.u32 gt, T##6, 0; newline\
@gt sub.u64.cc Z##0, Z##0, 0x89F6FFFFFFFD0003U; newline\
@gt subc.u64.cc Z##1, Z##1, 0x140BFFF43BF3FFFDU; newline\
@gt subc.u64.cc Z##2, Z##2, 0xA0B767A8AC38A745U; newline\
@gt subc.u64.cc Z##3, Z##3, 0x8831A7AC8FADA8BAU; newline\
@gt subc.u64.cc Z##4, Z##4, 0xA3F8E5685DA91392U; newline\
@gt subc.u64.cc Z##5, Z##5, 0xEA09A13C057F1B6CU; newline\
@gt subc.u32 T##6, T##6, 0; newline\
\
/* if z >= 2^384 then z -= 9m */ newline\
\
setp.gt.u32 gt, T##6, 0; newline\
@gt sub.u64.cc Z##0, Z##0, 0x89F6FFFFFFFD0003U; newline\
@gt subc.u64.cc Z##1, Z##1, 0x140BFFF43BF3FFFDU; newline\
@gt subc.u64.cc Z##2, Z##2, 0xA0B767A8AC38A745U; newline\
@gt subc.u64.cc Z##3, Z##3, 0x8831A7AC8FADA8BAU; newline\
@gt subc.u64.cc Z##4, Z##4, 0xA3F8E5685DA91392U; newline\
@gt subc.u64.cc Z##5, Z##5, 0xEA09A13C057F1B6CU
#endif