35 static const double TwoTo16 = 65536.0;
36 static const double TwoToMinus16 = 1.0/65536.0;
37 static const double Zero = 0.0;
38 static const double TwoTo32 = 65536.0 * 65536.0;
39 static const double TwoToMinus32 = 1.0 / (65536.0 * 65536.0);
43 double upper32(double);
44 double lower32(double, double);
45 double mod(double, double, double);
49 static double
50 upper32(double x)  in upper32()
56 static double
57 lower32(double x, double y)  in lower32()
62 static double
63 mod(double x, double oneoverm, double m)  in mod()
72 cleanup(double *dt, int from, int tlen)  in cleanup()
75 	double tmp, tmp1, x, x1;  in cleanup()
91 conv_d16_to_i32(uint32_t *i32, double *d16, int64_t *tmp, int ilen)  in conv_d16_to_i32()
120 conv_i32_to_d32(double *d32, uint32_t *i32, int len)  in conv_i32_to_d32()
126 		d32[i] = (double)(i32[i]);  in conv_i32_to_d32()
131 conv_i32_to_d16(double *d16, uint32_t *i32, int len)  in conv_i32_to_d16()
139 		d16[2 * i] = (double)(a & 0xffff);  in conv_i32_to_d16()
140 		d16[2 * i + 1] = (double)(a >> 16);  in conv_i32_to_d16()
147 i16_to_d16_and_d32x4(const double *,	/* 1/(2^16) */
148 			const double *,	/* 2^16 */
149 			const double *,	/* 0 */
150 			double *,	/* result16 */
151 			double *,	/* result32 */
159 i16_to_d16_and_d32x4(const double *dummy1,	/* 1/(2^16) */  in i16_to_d16_and_d32x4()
160 			const double *dummy2,	/* 2^16 */  in i16_to_d16_and_d32x4()
161 			const double *dummy3,	/* 0 */  in i16_to_d16_and_d32x4()
162 			double *result16,  in i16_to_d16_and_d32x4()
163 			double *result32,  in i16_to_d16_and_d32x4()
175 	result16[0] = (double)(a & 0xffff);  in i16_to_d16_and_d32x4()
176 	result16[1] = (double)(a >> 16);  in i16_to_d16_and_d32x4()
177 	result32[0] = (double)a;  in i16_to_d16_and_d32x4()
178 	result16[2] = (double)(b & 0xffff);  in i16_to_d16_and_d32x4()
179 	result16[3] = (double)(b >> 16);  in i16_to_d16_and_d32x4()
180 	result32[1] = (double)b;  in i16_to_d16_and_d32x4()
181 	result16[4] = (double)(c & 0xffff);  in i16_to_d16_and_d32x4()
182 	result16[5] = (double)(c >> 16);  in i16_to_d16_and_d32x4()
183 	result32[2] = (double)c;  in i16_to_d16_and_d32x4()
184 	result16[6] = (double)(d & 0xffff);  in i16_to_d16_and_d32x4()
185 	result16[7] = (double)(d >> 16);  in i16_to_d16_and_d32x4()
186 	result32[3] = (double)d;  in i16_to_d16_and_d32x4()
193 conv_i32_to_d32_and_d16(double *d32, double *d16, uint32_t *i32, int len)  in conv_i32_to_d32_and_d16()
205 		d32[i] = (double)(i32[i]);  in conv_i32_to_d32_and_d16()
206 		d16[2 * i] = (double)(a & 0xffff);  in conv_i32_to_d32_and_d16()
207 		d16[2 * i + 1] = (double)(a >> 16);  in conv_i32_to_d32_and_d16()
242 			double *dm1, double *dm2, double *dt,  in mont_mulf_noconv()
243 			double *dn, uint32_t *nint,  in mont_mulf_noconv()
244 			int nlen, double dn0)  in mont_mulf_noconv()
247 	double digit, m2j, a, b;  in mont_mulf_noconv()
248 	double *pdm1, *pdm2, *pdn, *pdtj, pdn_0, pdm1_0;  in mont_mulf_noconv()