Context Navigation

fflas_bounds.inl

Timestamp:

03/01/08 14:02:46 (10 months ago)

Author:

pernet

Message:

Merge MAbshoff and Stefan Kaspar's patches (sun arch and inlines)

Files:

r44	r56
13	13	#else
14	14	#define FFLAS_INT_TYPE long unsigned int
	15	#endif
	16
	17	/* From M Abshoff */
	18	#if defined(__sun)
	19	#define lround(x) my_lround(x)
	20	static long my_lround(double x)
	21	{
	22	return (long) ((x) >= 0 ? (x) + 0.5 : (x) - 0.5);
	23	}
15	24	#endif
16	25
…	…
97	106	// Computes nmax s.t. (p-1)/2*(p^{nmax-1} + (p-2)^{nmax-1}) < 2^53
98	107	//---------------------------------------------------------------------
99		size_t bound_compute_double(const long long pi) {
	108	inline size_t bound_compute_double(const long long pi) {
100	109
101	110	long long p=pi,p1=1,p2=1;
…	…
110	119	return nmax;
111	120	}
112		size_t bound_compute_double_balanced(const long long pi) {
	121	inline size_t bound_compute_double_balanced(const long long pi) {
113	122
114	123	long long p=(pi+1)/2,p1=1;
…	…
121	130	return nmax;
122	131	}
123		size_t bound_compute_float(const long long pi) {
	132	inline size_t bound_compute_float(const long long pi) {
124	133
125	134	long long p=pi,p1=1,p2=1;
…	…
134	143	return nmax;
135	144	}
136		size_t bound_compute_float_balanced(const long long pi) {
	145	inline size_t bound_compute_float_balanced(const long long pi) {
137	146
138	147	long long p=(pi+1)/2,p1=1;