1
2 #include "simint/ostei/gen/ostei_generated.h"
3 //////////////////////////////////////////////
4 // VRR: ( d s | f s )
5 //////////////////////////////////////////////
VRR_K_d_s_f_s(SIMINT_DBLTYPE * const restrict PRIM_INT__d_s_f_s,const SIMINT_DBLTYPE * const restrict PRIM_INT__d_s_d_s,const SIMINT_DBLTYPE * const restrict PRIM_INT__d_s_p_s,const SIMINT_DBLTYPE * const restrict PRIM_INT__p_s_d_s,const SIMINT_DBLTYPE * Q_PA,const SIMINT_DBLTYPE a_over_q,const SIMINT_DBLTYPE * aoq_PQ,const SIMINT_DBLTYPE one_over_2pq,const SIMINT_DBLTYPE one_over_2q,const int num_n)6 void VRR_K_d_s_f_s(
7 SIMINT_DBLTYPE * const restrict PRIM_INT__d_s_f_s,
8 const SIMINT_DBLTYPE * const restrict PRIM_INT__d_s_d_s,
9 const SIMINT_DBLTYPE * const restrict PRIM_INT__d_s_p_s,
10 const SIMINT_DBLTYPE * const restrict PRIM_INT__p_s_d_s,
11 const SIMINT_DBLTYPE * Q_PA,
12 const SIMINT_DBLTYPE a_over_q,
13 const SIMINT_DBLTYPE * aoq_PQ,
14 const SIMINT_DBLTYPE one_over_2pq,
15 const SIMINT_DBLTYPE one_over_2q,
16 const int num_n)
17 {
18 int n = 0;
19 const SIMINT_DBLTYPE vrr_const_2_over_2q = SIMINT_MUL(SIMINT_DBLSET1(2), one_over_2q);
20 const SIMINT_DBLTYPE vrr_const_1_over_2pq = SIMINT_MUL(SIMINT_DBLSET1(1), one_over_2pq);
21 const SIMINT_DBLTYPE vrr_const_2_over_2pq = SIMINT_MUL(SIMINT_DBLSET1(2), one_over_2pq);
22
23 // Forming PRIM_INT__d_s_f_s[num_n * 60];
24 for(n = 0; n < num_n; ++n) // loop over orders of auxiliary function
25 {
26
27 PRIM_INT__d_s_f_s[n * 60 + 0] = SIMINT_MUL(Q_PA[0], PRIM_INT__d_s_d_s[n * 36 + 0]);
28 PRIM_INT__d_s_f_s[n * 60 + 0] = SIMINT_FMADD( aoq_PQ[0], PRIM_INT__d_s_d_s[(n+1) * 36 + 0], PRIM_INT__d_s_f_s[n * 60 + 0]);
29 PRIM_INT__d_s_f_s[n * 60 + 0] = SIMINT_FMADD( vrr_const_2_over_2q, SIMINT_FMADD(a_over_q, PRIM_INT__d_s_p_s[(n+1) * 18 + 0], PRIM_INT__d_s_p_s[n * 18 + 0]), PRIM_INT__d_s_f_s[n * 60 + 0]);
30 PRIM_INT__d_s_f_s[n * 60 + 0] = SIMINT_FMADD( vrr_const_2_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 0], PRIM_INT__d_s_f_s[n * 60 + 0]);
31
32 PRIM_INT__d_s_f_s[n * 60 + 1] = SIMINT_MUL(Q_PA[1], PRIM_INT__d_s_d_s[n * 36 + 0]);
33 PRIM_INT__d_s_f_s[n * 60 + 1] = SIMINT_FMADD( aoq_PQ[1], PRIM_INT__d_s_d_s[(n+1) * 36 + 0], PRIM_INT__d_s_f_s[n * 60 + 1]);
34
35 PRIM_INT__d_s_f_s[n * 60 + 2] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 0]);
36 PRIM_INT__d_s_f_s[n * 60 + 2] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 0], PRIM_INT__d_s_f_s[n * 60 + 2]);
37
38 PRIM_INT__d_s_f_s[n * 60 + 3] = SIMINT_MUL(Q_PA[0], PRIM_INT__d_s_d_s[n * 36 + 3]);
39 PRIM_INT__d_s_f_s[n * 60 + 3] = SIMINT_FMADD( aoq_PQ[0], PRIM_INT__d_s_d_s[(n+1) * 36 + 3], PRIM_INT__d_s_f_s[n * 60 + 3]);
40 PRIM_INT__d_s_f_s[n * 60 + 3] = SIMINT_FMADD( vrr_const_2_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 3], PRIM_INT__d_s_f_s[n * 60 + 3]);
41
42 PRIM_INT__d_s_f_s[n * 60 + 4] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 1]);
43 PRIM_INT__d_s_f_s[n * 60 + 4] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 1], PRIM_INT__d_s_f_s[n * 60 + 4]);
44
45 PRIM_INT__d_s_f_s[n * 60 + 5] = SIMINT_MUL(Q_PA[0], PRIM_INT__d_s_d_s[n * 36 + 5]);
46 PRIM_INT__d_s_f_s[n * 60 + 5] = SIMINT_FMADD( aoq_PQ[0], PRIM_INT__d_s_d_s[(n+1) * 36 + 5], PRIM_INT__d_s_f_s[n * 60 + 5]);
47 PRIM_INT__d_s_f_s[n * 60 + 5] = SIMINT_FMADD( vrr_const_2_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 5], PRIM_INT__d_s_f_s[n * 60 + 5]);
48
49 PRIM_INT__d_s_f_s[n * 60 + 6] = SIMINT_MUL(Q_PA[1], PRIM_INT__d_s_d_s[n * 36 + 3]);
50 PRIM_INT__d_s_f_s[n * 60 + 6] = SIMINT_FMADD( aoq_PQ[1], PRIM_INT__d_s_d_s[(n+1) * 36 + 3], PRIM_INT__d_s_f_s[n * 60 + 6]);
51 PRIM_INT__d_s_f_s[n * 60 + 6] = SIMINT_FMADD( vrr_const_2_over_2q, SIMINT_FMADD(a_over_q, PRIM_INT__d_s_p_s[(n+1) * 18 + 1], PRIM_INT__d_s_p_s[n * 18 + 1]), PRIM_INT__d_s_f_s[n * 60 + 6]);
52
53 PRIM_INT__d_s_f_s[n * 60 + 7] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 3]);
54 PRIM_INT__d_s_f_s[n * 60 + 7] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 3], PRIM_INT__d_s_f_s[n * 60 + 7]);
55
56 PRIM_INT__d_s_f_s[n * 60 + 8] = SIMINT_MUL(Q_PA[1], PRIM_INT__d_s_d_s[n * 36 + 5]);
57 PRIM_INT__d_s_f_s[n * 60 + 8] = SIMINT_FMADD( aoq_PQ[1], PRIM_INT__d_s_d_s[(n+1) * 36 + 5], PRIM_INT__d_s_f_s[n * 60 + 8]);
58
59 PRIM_INT__d_s_f_s[n * 60 + 9] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 5]);
60 PRIM_INT__d_s_f_s[n * 60 + 9] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 5], PRIM_INT__d_s_f_s[n * 60 + 9]);
61 PRIM_INT__d_s_f_s[n * 60 + 9] = SIMINT_FMADD( vrr_const_2_over_2q, SIMINT_FMADD(a_over_q, PRIM_INT__d_s_p_s[(n+1) * 18 + 2], PRIM_INT__d_s_p_s[n * 18 + 2]), PRIM_INT__d_s_f_s[n * 60 + 9]);
62
63 PRIM_INT__d_s_f_s[n * 60 + 10] = SIMINT_MUL(Q_PA[0], PRIM_INT__d_s_d_s[n * 36 + 6]);
64 PRIM_INT__d_s_f_s[n * 60 + 10] = SIMINT_FMADD( aoq_PQ[0], PRIM_INT__d_s_d_s[(n+1) * 36 + 6], PRIM_INT__d_s_f_s[n * 60 + 10]);
65 PRIM_INT__d_s_f_s[n * 60 + 10] = SIMINT_FMADD( vrr_const_2_over_2q, SIMINT_FMADD(a_over_q, PRIM_INT__d_s_p_s[(n+1) * 18 + 3], PRIM_INT__d_s_p_s[n * 18 + 3]), PRIM_INT__d_s_f_s[n * 60 + 10]);
66 PRIM_INT__d_s_f_s[n * 60 + 10] = SIMINT_FMADD( vrr_const_1_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 6], PRIM_INT__d_s_f_s[n * 60 + 10]);
67
68 PRIM_INT__d_s_f_s[n * 60 + 11] = SIMINT_MUL(Q_PA[1], PRIM_INT__d_s_d_s[n * 36 + 6]);
69 PRIM_INT__d_s_f_s[n * 60 + 11] = SIMINT_FMADD( aoq_PQ[1], PRIM_INT__d_s_d_s[(n+1) * 36 + 6], PRIM_INT__d_s_f_s[n * 60 + 11]);
70 PRIM_INT__d_s_f_s[n * 60 + 11] = SIMINT_FMADD( vrr_const_1_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 0], PRIM_INT__d_s_f_s[n * 60 + 11]);
71
72 PRIM_INT__d_s_f_s[n * 60 + 12] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 6]);
73 PRIM_INT__d_s_f_s[n * 60 + 12] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 6], PRIM_INT__d_s_f_s[n * 60 + 12]);
74
75 PRIM_INT__d_s_f_s[n * 60 + 13] = SIMINT_MUL(Q_PA[0], PRIM_INT__d_s_d_s[n * 36 + 9]);
76 PRIM_INT__d_s_f_s[n * 60 + 13] = SIMINT_FMADD( aoq_PQ[0], PRIM_INT__d_s_d_s[(n+1) * 36 + 9], PRIM_INT__d_s_f_s[n * 60 + 13]);
77 PRIM_INT__d_s_f_s[n * 60 + 13] = SIMINT_FMADD( vrr_const_1_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 9], PRIM_INT__d_s_f_s[n * 60 + 13]);
78
79 PRIM_INT__d_s_f_s[n * 60 + 14] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 7]);
80 PRIM_INT__d_s_f_s[n * 60 + 14] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 7], PRIM_INT__d_s_f_s[n * 60 + 14]);
81
82 PRIM_INT__d_s_f_s[n * 60 + 15] = SIMINT_MUL(Q_PA[0], PRIM_INT__d_s_d_s[n * 36 + 11]);
83 PRIM_INT__d_s_f_s[n * 60 + 15] = SIMINT_FMADD( aoq_PQ[0], PRIM_INT__d_s_d_s[(n+1) * 36 + 11], PRIM_INT__d_s_f_s[n * 60 + 15]);
84 PRIM_INT__d_s_f_s[n * 60 + 15] = SIMINT_FMADD( vrr_const_1_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 11], PRIM_INT__d_s_f_s[n * 60 + 15]);
85
86 PRIM_INT__d_s_f_s[n * 60 + 16] = SIMINT_MUL(Q_PA[1], PRIM_INT__d_s_d_s[n * 36 + 9]);
87 PRIM_INT__d_s_f_s[n * 60 + 16] = SIMINT_FMADD( aoq_PQ[1], PRIM_INT__d_s_d_s[(n+1) * 36 + 9], PRIM_INT__d_s_f_s[n * 60 + 16]);
88 PRIM_INT__d_s_f_s[n * 60 + 16] = SIMINT_FMADD( vrr_const_2_over_2q, SIMINT_FMADD(a_over_q, PRIM_INT__d_s_p_s[(n+1) * 18 + 4], PRIM_INT__d_s_p_s[n * 18 + 4]), PRIM_INT__d_s_f_s[n * 60 + 16]);
89 PRIM_INT__d_s_f_s[n * 60 + 16] = SIMINT_FMADD( vrr_const_1_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 3], PRIM_INT__d_s_f_s[n * 60 + 16]);
90
91 PRIM_INT__d_s_f_s[n * 60 + 17] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 9]);
92 PRIM_INT__d_s_f_s[n * 60 + 17] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 9], PRIM_INT__d_s_f_s[n * 60 + 17]);
93
94 PRIM_INT__d_s_f_s[n * 60 + 18] = SIMINT_MUL(Q_PA[1], PRIM_INT__d_s_d_s[n * 36 + 11]);
95 PRIM_INT__d_s_f_s[n * 60 + 18] = SIMINT_FMADD( aoq_PQ[1], PRIM_INT__d_s_d_s[(n+1) * 36 + 11], PRIM_INT__d_s_f_s[n * 60 + 18]);
96 PRIM_INT__d_s_f_s[n * 60 + 18] = SIMINT_FMADD( vrr_const_1_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 5], PRIM_INT__d_s_f_s[n * 60 + 18]);
97
98 PRIM_INT__d_s_f_s[n * 60 + 19] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 11]);
99 PRIM_INT__d_s_f_s[n * 60 + 19] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 11], PRIM_INT__d_s_f_s[n * 60 + 19]);
100 PRIM_INT__d_s_f_s[n * 60 + 19] = SIMINT_FMADD( vrr_const_2_over_2q, SIMINT_FMADD(a_over_q, PRIM_INT__d_s_p_s[(n+1) * 18 + 5], PRIM_INT__d_s_p_s[n * 18 + 5]), PRIM_INT__d_s_f_s[n * 60 + 19]);
101
102 PRIM_INT__d_s_f_s[n * 60 + 20] = SIMINT_MUL(Q_PA[0], PRIM_INT__d_s_d_s[n * 36 + 12]);
103 PRIM_INT__d_s_f_s[n * 60 + 20] = SIMINT_FMADD( aoq_PQ[0], PRIM_INT__d_s_d_s[(n+1) * 36 + 12], PRIM_INT__d_s_f_s[n * 60 + 20]);
104 PRIM_INT__d_s_f_s[n * 60 + 20] = SIMINT_FMADD( vrr_const_2_over_2q, SIMINT_FMADD(a_over_q, PRIM_INT__d_s_p_s[(n+1) * 18 + 6], PRIM_INT__d_s_p_s[n * 18 + 6]), PRIM_INT__d_s_f_s[n * 60 + 20]);
105 PRIM_INT__d_s_f_s[n * 60 + 20] = SIMINT_FMADD( vrr_const_1_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 12], PRIM_INT__d_s_f_s[n * 60 + 20]);
106
107 PRIM_INT__d_s_f_s[n * 60 + 21] = SIMINT_MUL(Q_PA[1], PRIM_INT__d_s_d_s[n * 36 + 12]);
108 PRIM_INT__d_s_f_s[n * 60 + 21] = SIMINT_FMADD( aoq_PQ[1], PRIM_INT__d_s_d_s[(n+1) * 36 + 12], PRIM_INT__d_s_f_s[n * 60 + 21]);
109
110 PRIM_INT__d_s_f_s[n * 60 + 22] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 12]);
111 PRIM_INT__d_s_f_s[n * 60 + 22] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 12], PRIM_INT__d_s_f_s[n * 60 + 22]);
112 PRIM_INT__d_s_f_s[n * 60 + 22] = SIMINT_FMADD( vrr_const_1_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 0], PRIM_INT__d_s_f_s[n * 60 + 22]);
113
114 PRIM_INT__d_s_f_s[n * 60 + 23] = SIMINT_MUL(Q_PA[0], PRIM_INT__d_s_d_s[n * 36 + 15]);
115 PRIM_INT__d_s_f_s[n * 60 + 23] = SIMINT_FMADD( aoq_PQ[0], PRIM_INT__d_s_d_s[(n+1) * 36 + 15], PRIM_INT__d_s_f_s[n * 60 + 23]);
116 PRIM_INT__d_s_f_s[n * 60 + 23] = SIMINT_FMADD( vrr_const_1_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 15], PRIM_INT__d_s_f_s[n * 60 + 23]);
117
118 PRIM_INT__d_s_f_s[n * 60 + 24] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 13]);
119 PRIM_INT__d_s_f_s[n * 60 + 24] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 13], PRIM_INT__d_s_f_s[n * 60 + 24]);
120 PRIM_INT__d_s_f_s[n * 60 + 24] = SIMINT_FMADD( vrr_const_1_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 1], PRIM_INT__d_s_f_s[n * 60 + 24]);
121
122 PRIM_INT__d_s_f_s[n * 60 + 25] = SIMINT_MUL(Q_PA[0], PRIM_INT__d_s_d_s[n * 36 + 17]);
123 PRIM_INT__d_s_f_s[n * 60 + 25] = SIMINT_FMADD( aoq_PQ[0], PRIM_INT__d_s_d_s[(n+1) * 36 + 17], PRIM_INT__d_s_f_s[n * 60 + 25]);
124 PRIM_INT__d_s_f_s[n * 60 + 25] = SIMINT_FMADD( vrr_const_1_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 17], PRIM_INT__d_s_f_s[n * 60 + 25]);
125
126 PRIM_INT__d_s_f_s[n * 60 + 26] = SIMINT_MUL(Q_PA[1], PRIM_INT__d_s_d_s[n * 36 + 15]);
127 PRIM_INT__d_s_f_s[n * 60 + 26] = SIMINT_FMADD( aoq_PQ[1], PRIM_INT__d_s_d_s[(n+1) * 36 + 15], PRIM_INT__d_s_f_s[n * 60 + 26]);
128 PRIM_INT__d_s_f_s[n * 60 + 26] = SIMINT_FMADD( vrr_const_2_over_2q, SIMINT_FMADD(a_over_q, PRIM_INT__d_s_p_s[(n+1) * 18 + 7], PRIM_INT__d_s_p_s[n * 18 + 7]), PRIM_INT__d_s_f_s[n * 60 + 26]);
129
130 PRIM_INT__d_s_f_s[n * 60 + 27] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 15]);
131 PRIM_INT__d_s_f_s[n * 60 + 27] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 15], PRIM_INT__d_s_f_s[n * 60 + 27]);
132 PRIM_INT__d_s_f_s[n * 60 + 27] = SIMINT_FMADD( vrr_const_1_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 3], PRIM_INT__d_s_f_s[n * 60 + 27]);
133
134 PRIM_INT__d_s_f_s[n * 60 + 28] = SIMINT_MUL(Q_PA[1], PRIM_INT__d_s_d_s[n * 36 + 17]);
135 PRIM_INT__d_s_f_s[n * 60 + 28] = SIMINT_FMADD( aoq_PQ[1], PRIM_INT__d_s_d_s[(n+1) * 36 + 17], PRIM_INT__d_s_f_s[n * 60 + 28]);
136
137 PRIM_INT__d_s_f_s[n * 60 + 29] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 17]);
138 PRIM_INT__d_s_f_s[n * 60 + 29] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 17], PRIM_INT__d_s_f_s[n * 60 + 29]);
139 PRIM_INT__d_s_f_s[n * 60 + 29] = SIMINT_FMADD( vrr_const_2_over_2q, SIMINT_FMADD(a_over_q, PRIM_INT__d_s_p_s[(n+1) * 18 + 8], PRIM_INT__d_s_p_s[n * 18 + 8]), PRIM_INT__d_s_f_s[n * 60 + 29]);
140 PRIM_INT__d_s_f_s[n * 60 + 29] = SIMINT_FMADD( vrr_const_1_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 5], PRIM_INT__d_s_f_s[n * 60 + 29]);
141
142 PRIM_INT__d_s_f_s[n * 60 + 30] = SIMINT_MUL(Q_PA[0], PRIM_INT__d_s_d_s[n * 36 + 18]);
143 PRIM_INT__d_s_f_s[n * 60 + 30] = SIMINT_FMADD( aoq_PQ[0], PRIM_INT__d_s_d_s[(n+1) * 36 + 18], PRIM_INT__d_s_f_s[n * 60 + 30]);
144 PRIM_INT__d_s_f_s[n * 60 + 30] = SIMINT_FMADD( vrr_const_2_over_2q, SIMINT_FMADD(a_over_q, PRIM_INT__d_s_p_s[(n+1) * 18 + 9], PRIM_INT__d_s_p_s[n * 18 + 9]), PRIM_INT__d_s_f_s[n * 60 + 30]);
145
146 PRIM_INT__d_s_f_s[n * 60 + 31] = SIMINT_MUL(Q_PA[1], PRIM_INT__d_s_d_s[n * 36 + 18]);
147 PRIM_INT__d_s_f_s[n * 60 + 31] = SIMINT_FMADD( aoq_PQ[1], PRIM_INT__d_s_d_s[(n+1) * 36 + 18], PRIM_INT__d_s_f_s[n * 60 + 31]);
148 PRIM_INT__d_s_f_s[n * 60 + 31] = SIMINT_FMADD( vrr_const_2_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 6], PRIM_INT__d_s_f_s[n * 60 + 31]);
149
150 PRIM_INT__d_s_f_s[n * 60 + 32] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 18]);
151 PRIM_INT__d_s_f_s[n * 60 + 32] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 18], PRIM_INT__d_s_f_s[n * 60 + 32]);
152
153 PRIM_INT__d_s_f_s[n * 60 + 33] = SIMINT_MUL(Q_PA[0], PRIM_INT__d_s_d_s[n * 36 + 21]);
154 PRIM_INT__d_s_f_s[n * 60 + 33] = SIMINT_FMADD( aoq_PQ[0], PRIM_INT__d_s_d_s[(n+1) * 36 + 21], PRIM_INT__d_s_f_s[n * 60 + 33]);
155
156 PRIM_INT__d_s_f_s[n * 60 + 34] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 19]);
157 PRIM_INT__d_s_f_s[n * 60 + 34] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 19], PRIM_INT__d_s_f_s[n * 60 + 34]);
158
159 PRIM_INT__d_s_f_s[n * 60 + 35] = SIMINT_MUL(Q_PA[0], PRIM_INT__d_s_d_s[n * 36 + 23]);
160 PRIM_INT__d_s_f_s[n * 60 + 35] = SIMINT_FMADD( aoq_PQ[0], PRIM_INT__d_s_d_s[(n+1) * 36 + 23], PRIM_INT__d_s_f_s[n * 60 + 35]);
161
162 PRIM_INT__d_s_f_s[n * 60 + 36] = SIMINT_MUL(Q_PA[1], PRIM_INT__d_s_d_s[n * 36 + 21]);
163 PRIM_INT__d_s_f_s[n * 60 + 36] = SIMINT_FMADD( aoq_PQ[1], PRIM_INT__d_s_d_s[(n+1) * 36 + 21], PRIM_INT__d_s_f_s[n * 60 + 36]);
164 PRIM_INT__d_s_f_s[n * 60 + 36] = SIMINT_FMADD( vrr_const_2_over_2q, SIMINT_FMADD(a_over_q, PRIM_INT__d_s_p_s[(n+1) * 18 + 10], PRIM_INT__d_s_p_s[n * 18 + 10]), PRIM_INT__d_s_f_s[n * 60 + 36]);
165 PRIM_INT__d_s_f_s[n * 60 + 36] = SIMINT_FMADD( vrr_const_2_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 9], PRIM_INT__d_s_f_s[n * 60 + 36]);
166
167 PRIM_INT__d_s_f_s[n * 60 + 37] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 21]);
168 PRIM_INT__d_s_f_s[n * 60 + 37] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 21], PRIM_INT__d_s_f_s[n * 60 + 37]);
169
170 PRIM_INT__d_s_f_s[n * 60 + 38] = SIMINT_MUL(Q_PA[1], PRIM_INT__d_s_d_s[n * 36 + 23]);
171 PRIM_INT__d_s_f_s[n * 60 + 38] = SIMINT_FMADD( aoq_PQ[1], PRIM_INT__d_s_d_s[(n+1) * 36 + 23], PRIM_INT__d_s_f_s[n * 60 + 38]);
172 PRIM_INT__d_s_f_s[n * 60 + 38] = SIMINT_FMADD( vrr_const_2_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 11], PRIM_INT__d_s_f_s[n * 60 + 38]);
173
174 PRIM_INT__d_s_f_s[n * 60 + 39] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 23]);
175 PRIM_INT__d_s_f_s[n * 60 + 39] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 23], PRIM_INT__d_s_f_s[n * 60 + 39]);
176 PRIM_INT__d_s_f_s[n * 60 + 39] = SIMINT_FMADD( vrr_const_2_over_2q, SIMINT_FMADD(a_over_q, PRIM_INT__d_s_p_s[(n+1) * 18 + 11], PRIM_INT__d_s_p_s[n * 18 + 11]), PRIM_INT__d_s_f_s[n * 60 + 39]);
177
178 PRIM_INT__d_s_f_s[n * 60 + 40] = SIMINT_MUL(Q_PA[0], PRIM_INT__d_s_d_s[n * 36 + 24]);
179 PRIM_INT__d_s_f_s[n * 60 + 40] = SIMINT_FMADD( aoq_PQ[0], PRIM_INT__d_s_d_s[(n+1) * 36 + 24], PRIM_INT__d_s_f_s[n * 60 + 40]);
180 PRIM_INT__d_s_f_s[n * 60 + 40] = SIMINT_FMADD( vrr_const_2_over_2q, SIMINT_FMADD(a_over_q, PRIM_INT__d_s_p_s[(n+1) * 18 + 12], PRIM_INT__d_s_p_s[n * 18 + 12]), PRIM_INT__d_s_f_s[n * 60 + 40]);
181
182 PRIM_INT__d_s_f_s[n * 60 + 41] = SIMINT_MUL(Q_PA[1], PRIM_INT__d_s_d_s[n * 36 + 24]);
183 PRIM_INT__d_s_f_s[n * 60 + 41] = SIMINT_FMADD( aoq_PQ[1], PRIM_INT__d_s_d_s[(n+1) * 36 + 24], PRIM_INT__d_s_f_s[n * 60 + 41]);
184 PRIM_INT__d_s_f_s[n * 60 + 41] = SIMINT_FMADD( vrr_const_1_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 12], PRIM_INT__d_s_f_s[n * 60 + 41]);
185
186 PRIM_INT__d_s_f_s[n * 60 + 42] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 24]);
187 PRIM_INT__d_s_f_s[n * 60 + 42] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 24], PRIM_INT__d_s_f_s[n * 60 + 42]);
188 PRIM_INT__d_s_f_s[n * 60 + 42] = SIMINT_FMADD( vrr_const_1_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 6], PRIM_INT__d_s_f_s[n * 60 + 42]);
189
190 PRIM_INT__d_s_f_s[n * 60 + 43] = SIMINT_MUL(Q_PA[0], PRIM_INT__d_s_d_s[n * 36 + 27]);
191 PRIM_INT__d_s_f_s[n * 60 + 43] = SIMINT_FMADD( aoq_PQ[0], PRIM_INT__d_s_d_s[(n+1) * 36 + 27], PRIM_INT__d_s_f_s[n * 60 + 43]);
192
193 PRIM_INT__d_s_f_s[n * 60 + 44] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 25]);
194 PRIM_INT__d_s_f_s[n * 60 + 44] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 25], PRIM_INT__d_s_f_s[n * 60 + 44]);
195 PRIM_INT__d_s_f_s[n * 60 + 44] = SIMINT_FMADD( vrr_const_1_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 7], PRIM_INT__d_s_f_s[n * 60 + 44]);
196
197 PRIM_INT__d_s_f_s[n * 60 + 45] = SIMINT_MUL(Q_PA[0], PRIM_INT__d_s_d_s[n * 36 + 29]);
198 PRIM_INT__d_s_f_s[n * 60 + 45] = SIMINT_FMADD( aoq_PQ[0], PRIM_INT__d_s_d_s[(n+1) * 36 + 29], PRIM_INT__d_s_f_s[n * 60 + 45]);
199
200 PRIM_INT__d_s_f_s[n * 60 + 46] = SIMINT_MUL(Q_PA[1], PRIM_INT__d_s_d_s[n * 36 + 27]);
201 PRIM_INT__d_s_f_s[n * 60 + 46] = SIMINT_FMADD( aoq_PQ[1], PRIM_INT__d_s_d_s[(n+1) * 36 + 27], PRIM_INT__d_s_f_s[n * 60 + 46]);
202 PRIM_INT__d_s_f_s[n * 60 + 46] = SIMINT_FMADD( vrr_const_2_over_2q, SIMINT_FMADD(a_over_q, PRIM_INT__d_s_p_s[(n+1) * 18 + 13], PRIM_INT__d_s_p_s[n * 18 + 13]), PRIM_INT__d_s_f_s[n * 60 + 46]);
203 PRIM_INT__d_s_f_s[n * 60 + 46] = SIMINT_FMADD( vrr_const_1_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 15], PRIM_INT__d_s_f_s[n * 60 + 46]);
204
205 PRIM_INT__d_s_f_s[n * 60 + 47] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 27]);
206 PRIM_INT__d_s_f_s[n * 60 + 47] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 27], PRIM_INT__d_s_f_s[n * 60 + 47]);
207 PRIM_INT__d_s_f_s[n * 60 + 47] = SIMINT_FMADD( vrr_const_1_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 9], PRIM_INT__d_s_f_s[n * 60 + 47]);
208
209 PRIM_INT__d_s_f_s[n * 60 + 48] = SIMINT_MUL(Q_PA[1], PRIM_INT__d_s_d_s[n * 36 + 29]);
210 PRIM_INT__d_s_f_s[n * 60 + 48] = SIMINT_FMADD( aoq_PQ[1], PRIM_INT__d_s_d_s[(n+1) * 36 + 29], PRIM_INT__d_s_f_s[n * 60 + 48]);
211 PRIM_INT__d_s_f_s[n * 60 + 48] = SIMINT_FMADD( vrr_const_1_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 17], PRIM_INT__d_s_f_s[n * 60 + 48]);
212
213 PRIM_INT__d_s_f_s[n * 60 + 49] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 29]);
214 PRIM_INT__d_s_f_s[n * 60 + 49] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 29], PRIM_INT__d_s_f_s[n * 60 + 49]);
215 PRIM_INT__d_s_f_s[n * 60 + 49] = SIMINT_FMADD( vrr_const_2_over_2q, SIMINT_FMADD(a_over_q, PRIM_INT__d_s_p_s[(n+1) * 18 + 14], PRIM_INT__d_s_p_s[n * 18 + 14]), PRIM_INT__d_s_f_s[n * 60 + 49]);
216 PRIM_INT__d_s_f_s[n * 60 + 49] = SIMINT_FMADD( vrr_const_1_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 11], PRIM_INT__d_s_f_s[n * 60 + 49]);
217
218 PRIM_INT__d_s_f_s[n * 60 + 50] = SIMINT_MUL(Q_PA[0], PRIM_INT__d_s_d_s[n * 36 + 30]);
219 PRIM_INT__d_s_f_s[n * 60 + 50] = SIMINT_FMADD( aoq_PQ[0], PRIM_INT__d_s_d_s[(n+1) * 36 + 30], PRIM_INT__d_s_f_s[n * 60 + 50]);
220 PRIM_INT__d_s_f_s[n * 60 + 50] = SIMINT_FMADD( vrr_const_2_over_2q, SIMINT_FMADD(a_over_q, PRIM_INT__d_s_p_s[(n+1) * 18 + 15], PRIM_INT__d_s_p_s[n * 18 + 15]), PRIM_INT__d_s_f_s[n * 60 + 50]);
221
222 PRIM_INT__d_s_f_s[n * 60 + 51] = SIMINT_MUL(Q_PA[1], PRIM_INT__d_s_d_s[n * 36 + 30]);
223 PRIM_INT__d_s_f_s[n * 60 + 51] = SIMINT_FMADD( aoq_PQ[1], PRIM_INT__d_s_d_s[(n+1) * 36 + 30], PRIM_INT__d_s_f_s[n * 60 + 51]);
224
225 PRIM_INT__d_s_f_s[n * 60 + 52] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 30]);
226 PRIM_INT__d_s_f_s[n * 60 + 52] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 30], PRIM_INT__d_s_f_s[n * 60 + 52]);
227 PRIM_INT__d_s_f_s[n * 60 + 52] = SIMINT_FMADD( vrr_const_2_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 12], PRIM_INT__d_s_f_s[n * 60 + 52]);
228
229 PRIM_INT__d_s_f_s[n * 60 + 53] = SIMINT_MUL(Q_PA[0], PRIM_INT__d_s_d_s[n * 36 + 33]);
230 PRIM_INT__d_s_f_s[n * 60 + 53] = SIMINT_FMADD( aoq_PQ[0], PRIM_INT__d_s_d_s[(n+1) * 36 + 33], PRIM_INT__d_s_f_s[n * 60 + 53]);
231
232 PRIM_INT__d_s_f_s[n * 60 + 54] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 31]);
233 PRIM_INT__d_s_f_s[n * 60 + 54] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 31], PRIM_INT__d_s_f_s[n * 60 + 54]);
234 PRIM_INT__d_s_f_s[n * 60 + 54] = SIMINT_FMADD( vrr_const_2_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 13], PRIM_INT__d_s_f_s[n * 60 + 54]);
235
236 PRIM_INT__d_s_f_s[n * 60 + 55] = SIMINT_MUL(Q_PA[0], PRIM_INT__d_s_d_s[n * 36 + 35]);
237 PRIM_INT__d_s_f_s[n * 60 + 55] = SIMINT_FMADD( aoq_PQ[0], PRIM_INT__d_s_d_s[(n+1) * 36 + 35], PRIM_INT__d_s_f_s[n * 60 + 55]);
238
239 PRIM_INT__d_s_f_s[n * 60 + 56] = SIMINT_MUL(Q_PA[1], PRIM_INT__d_s_d_s[n * 36 + 33]);
240 PRIM_INT__d_s_f_s[n * 60 + 56] = SIMINT_FMADD( aoq_PQ[1], PRIM_INT__d_s_d_s[(n+1) * 36 + 33], PRIM_INT__d_s_f_s[n * 60 + 56]);
241 PRIM_INT__d_s_f_s[n * 60 + 56] = SIMINT_FMADD( vrr_const_2_over_2q, SIMINT_FMADD(a_over_q, PRIM_INT__d_s_p_s[(n+1) * 18 + 16], PRIM_INT__d_s_p_s[n * 18 + 16]), PRIM_INT__d_s_f_s[n * 60 + 56]);
242
243 PRIM_INT__d_s_f_s[n * 60 + 57] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 33]);
244 PRIM_INT__d_s_f_s[n * 60 + 57] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 33], PRIM_INT__d_s_f_s[n * 60 + 57]);
245 PRIM_INT__d_s_f_s[n * 60 + 57] = SIMINT_FMADD( vrr_const_2_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 15], PRIM_INT__d_s_f_s[n * 60 + 57]);
246
247 PRIM_INT__d_s_f_s[n * 60 + 58] = SIMINT_MUL(Q_PA[1], PRIM_INT__d_s_d_s[n * 36 + 35]);
248 PRIM_INT__d_s_f_s[n * 60 + 58] = SIMINT_FMADD( aoq_PQ[1], PRIM_INT__d_s_d_s[(n+1) * 36 + 35], PRIM_INT__d_s_f_s[n * 60 + 58]);
249
250 PRIM_INT__d_s_f_s[n * 60 + 59] = SIMINT_MUL(Q_PA[2], PRIM_INT__d_s_d_s[n * 36 + 35]);
251 PRIM_INT__d_s_f_s[n * 60 + 59] = SIMINT_FMADD( aoq_PQ[2], PRIM_INT__d_s_d_s[(n+1) * 36 + 35], PRIM_INT__d_s_f_s[n * 60 + 59]);
252 PRIM_INT__d_s_f_s[n * 60 + 59] = SIMINT_FMADD( vrr_const_2_over_2q, SIMINT_FMADD(a_over_q, PRIM_INT__d_s_p_s[(n+1) * 18 + 17], PRIM_INT__d_s_p_s[n * 18 + 17]), PRIM_INT__d_s_f_s[n * 60 + 59]);
253 PRIM_INT__d_s_f_s[n * 60 + 59] = SIMINT_FMADD( vrr_const_2_over_2pq, PRIM_INT__p_s_d_s[(n+1) * 18 + 17], PRIM_INT__d_s_f_s[n * 60 + 59]);
254
255 }
256
257 }
258
259
260