4  * Permission is hereby granted, free of charge, to any person obtaining a
16  * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.  IN NO EVENT SHALL
36 #define dfixed_const(A) (u32)(((A) << 12))/*  + ((B + 0.000122)*4096)) */  argument
37 #define dfixed_const_half(A) (u32)(((A) << 12) + 2048)  argument
38 #define dfixed_const_666(A) (u32)(((A) << 12) + 2731)  argument
39 #define dfixed_const_8(A) (u32)(((A) << 12) + 3277)  argument
40 #define dfixed_mul(A, B) ((u64)((u64)(A).full * (B).full + 2048) >> 12)  argument
41 #define dfixed_init(A) { .full = dfixed_const((A)) }  argument
42 #define dfixed_init_half(A) { .full = dfixed_const_half((A)) }  argument
43 #define dfixed_trunc(A) ((A).full >> 12)  argument
44 #define dfixed_frac(A) ((A).full & ((1 << 12) - 1))  argument
46 static inline u32 dfixed_floor(fixed20_12 A)  in dfixed_floor()  argument
48 	u32 non_frac = dfixed_trunc(A);  in dfixed_floor()
53 static inline u32 dfixed_ceil(fixed20_12 A)  in dfixed_ceil()  argument
55 	u32 non_frac = dfixed_trunc(A);  in dfixed_ceil()
57 	if (A.full > dfixed_const(non_frac))  in dfixed_ceil()
58 		return dfixed_const(non_frac + 1);  in dfixed_ceil()
63 static inline u32 dfixed_div(fixed20_12 A, fixed20_12 B)  in dfixed_div()  argument
65 	u64 tmp = ((u64)A.full << 13);  in dfixed_div()
68 	tmp += 1;  in dfixed_div()
74 #define DRM_FIXED_ONE		(1ULL << DRM_FIXED_POINT)
75 #define DRM_FIXED_DECIMAL_MASK	(DRM_FIXED_ONE - 1)
77 #define DRM_FIXED_EPSILON	1LL
80 static inline s64 drm_int2fixp(int a)  in drm_int2fixp()  argument
82 	return ((s64)a) << DRM_FIXED_POINT;  in drm_int2fixp()
85 static inline int drm_fixp2int(s64 a)  in drm_fixp2int()  argument
87 	return ((s64)a) >> DRM_FIXED_POINT;  in drm_fixp2int()
90 static inline int drm_fixp2int_round(s64 a)  in drm_fixp2int_round()  argument
92 	return drm_fixp2int(a + DRM_FIXED_ONE / 2);  in drm_fixp2int_round()
95 static inline int drm_fixp2int_ceil(s64 a)  in drm_fixp2int_ceil()  argument
97 	if (a >= 0)  in drm_fixp2int_ceil()
98 		return drm_fixp2int(a + DRM_FIXED_ALMOST_ONE);  in drm_fixp2int_ceil()
100 		return drm_fixp2int(a - DRM_FIXED_ALMOST_ONE);  in drm_fixp2int_ceil()
103 static inline unsigned drm_fixp_msbset(s64 a)  in drm_fixp_msbset()  argument
105 	unsigned shift, sign = (a >> 63) & 1;  in drm_fixp_msbset()
108 		if (((a >> shift) & 1) != sign)  in drm_fixp_msbset()
114 static inline s64 drm_fixp_mul(s64 a, s64 b)  in drm_fixp_mul()  argument
116 	unsigned shift = drm_fixp_msbset(a) + drm_fixp_msbset(b);  in drm_fixp_mul()
121 		a >>= (shift >> 1) + (shift & 1);  in drm_fixp_mul()
122 		b >>= shift >> 1;  in drm_fixp_mul()
126 	result = a * b;  in drm_fixp_mul()
137 static inline s64 drm_fixp_div(s64 a, s64 b)  in drm_fixp_div()  argument
139 	unsigned shift = 62 - drm_fixp_msbset(a);  in drm_fixp_div()
142 	a <<= shift;  in drm_fixp_div()
147 	result = div64_s64(a, b);  in drm_fixp_div()
155 static inline s64 drm_fixp_from_fraction(s64 a, s64 b)  in drm_fixp_from_fraction()  argument
158 	bool a_neg = a < 0;  in drm_fixp_from_fraction()
160 	u64 a_abs = a_neg ? -a : a;  in drm_fixp_from_fraction()
172 			rem <<= 1;  in drm_fixp_from_fraction()
173 			res_abs <<= 1;  in drm_fixp_from_fraction()
175 				res_abs |= 1;  in drm_fixp_from_fraction()
183 		u64 summand = (rem << 1) >= b_abs;  in drm_fixp_from_fraction()
198 	u64 count = 1;  in drm_fixp_exp()
201 		y = -1 * x;  in drm_fixp_exp()
207 		count = count + 1;  in drm_fixp_exp()