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