2010-07-12 15:41:10 +00:00
|
|
|
#ifndef _WMAPRO_MATH_H_
|
|
|
|
|
#define _WMAPRO_MATH_H_
|
|
|
|
|
|
2010-07-05 22:33:37 +00:00
|
|
|
#include <inttypes.h>
|
|
|
|
|
|
2010-07-26 21:43:07 +00:00
|
|
|
/* rockbox: not used
|
2010-07-05 22:33:37 +00:00
|
|
|
#define fixtof16(x) (float)((float)(x) / (float)(1 << 16))
|
2010-07-12 15:41:10 +00:00
|
|
|
#define fixtof31(x) (float)((float)(x) / (float)(1 << 31))
|
2010-07-05 22:33:37 +00:00
|
|
|
#define ftofix16(x) ((int32_t)((x) * (float)(1 << 16) + ((x) < 0 ? -0.5:0.5)))
|
2010-07-12 15:41:10 +00:00
|
|
|
#define ftofix31(x) ((int32_t)((x) * (float)(1 << 31) + ((x) < 0 ? -0.5:0.5)))
|
2010-07-26 21:43:07 +00:00
|
|
|
*/
|
2010-07-05 22:33:37 +00:00
|
|
|
|
2010-07-26 21:43:07 +00:00
|
|
|
#if defined(CPU_ARM)
|
|
|
|
|
/* Calculates: result = (X*Y)>>Z */
|
|
|
|
|
#define fixmulshift(X,Y,Z) \
|
|
|
|
|
({ \
|
|
|
|
|
int32_t lo; \
|
|
|
|
|
int32_t hi; \
|
|
|
|
|
asm volatile ( \
|
|
|
|
|
"smull %[lo], %[hi], %[x], %[y] \n\t" /* multiply */ \
|
|
|
|
|
"mov %[lo], %[lo], lsr %[shr] \n\t" /* lo >>= Z */ \
|
|
|
|
|
"orr %[lo], %[lo], %[hi], lsl %[shl]" /* lo |= (hi << (32-Z)) */ \
|
|
|
|
|
: [lo]"=&r"(lo), [hi]"=&r"(hi) \
|
|
|
|
|
: [x]"r"(X), [y]"r"(Y), [shr]"r"(Z), [shl]"r"(32-Z)); \
|
|
|
|
|
lo; \
|
|
|
|
|
})
|
|
|
|
|
|
|
|
|
|
/* Calculates: result = (X*Y)>>16 */
|
|
|
|
|
#define fixmul16(X,Y) \
|
|
|
|
|
({ \
|
|
|
|
|
int32_t lo; \
|
|
|
|
|
int32_t hi; \
|
|
|
|
|
asm volatile ( \
|
|
|
|
|
"smull %[lo], %[hi], %[x], %[y] \n\t" /* multiply */ \
|
|
|
|
|
"mov %[lo], %[lo], lsr #16 \n\t" /* lo >>= 16 */ \
|
|
|
|
|
"orr %[lo], %[lo], %[hi], lsl #16" /* lo |= (hi << 16) */ \
|
|
|
|
|
: [lo]"=&r"(lo), [hi]"=&r"(hi) \
|
|
|
|
|
: [x]"r"(X), [y]"r"(Y)); \
|
|
|
|
|
lo; \
|
|
|
|
|
})
|
|
|
|
|
|
|
|
|
|
/* Calculates: result = (X*Y)>>24 */
|
|
|
|
|
#define fixmul24(X,Y) \
|
|
|
|
|
({ \
|
|
|
|
|
int32_t lo; \
|
|
|
|
|
int32_t hi; \
|
|
|
|
|
asm volatile ( \
|
|
|
|
|
"smull %[lo], %[hi], %[x], %[y] \n\t" /* multiply */ \
|
|
|
|
|
"mov %[lo], %[lo], lsr #24 \n\t" /* lo >>= 24 */ \
|
|
|
|
|
"orr %[lo], %[lo], %[hi], lsl #8" /* lo |= (hi << 8) */ \
|
|
|
|
|
: [lo]"=&r"(lo), [hi]"=&r"(hi) \
|
|
|
|
|
: [x]"r"(X), [y]"r"(Y)); \
|
|
|
|
|
lo; \
|
|
|
|
|
})
|
|
|
|
|
|
|
|
|
|
/* Calculates: result = (X*Y)>>31 */
|
|
|
|
|
#define fixmul31(X,Y) \
|
|
|
|
|
({ \
|
|
|
|
|
int32_t lo; \
|
|
|
|
|
int32_t hi; \
|
|
|
|
|
asm volatile ( \
|
|
|
|
|
"smull %[lo], %[hi], %[x], %[y] \n\t" /* multiply */ \
|
|
|
|
|
"mov %[lo], %[lo], lsr #31 \n\t" /* lo >>= 31 */ \
|
|
|
|
|
"orr %[lo], %[lo], %[hi], lsl #1" /* lo |= (hi << 1) */ \
|
|
|
|
|
: [lo]"=&r"(lo), [hi]"=&r"(hi) \
|
|
|
|
|
: [x]"r"(X), [y]"r"(Y)); \
|
|
|
|
|
lo; \
|
|
|
|
|
})
|
|
|
|
|
#elif defined(CPU_COLDFIRE)
|
|
|
|
|
/* Calculates: result = (X*Y)>>Z */
|
|
|
|
|
#define fixmulshift(X,Y,Z) \
|
|
|
|
|
({ \
|
|
|
|
|
int32_t t1; \
|
|
|
|
|
int32_t t2; \
|
|
|
|
|
asm volatile ( \
|
|
|
|
|
"mac.l %[x],%[y],%%acc0\n\t" /* multiply */ \
|
|
|
|
|
"mulu.l %[y],%[x] \n\t" /* get lower half, avoid emac stall */ \
|
|
|
|
|
"movclr.l %%acc0,%[t1] \n\t" /* get higher half */ \
|
|
|
|
|
"moveq.l #31,%[t2] \n\t" \
|
|
|
|
|
"sub.l %[sh],%[t2] \n\t" /* t2 = 31 - shift */ \
|
|
|
|
|
"ble.s 1f \n\t" \
|
|
|
|
|
"asl.l %[t2],%[t1] \n\t" /* hi <<= 31 - shift */ \
|
|
|
|
|
"lsr.l %[sh],%[x] \n\t" /* (unsigned)lo >>= shift */ \
|
|
|
|
|
"or.l %[x],%[t1] \n\t" /* combine result */ \
|
|
|
|
|
"bra.s 2f \n\t" \
|
|
|
|
|
"1: \n\t" \
|
|
|
|
|
"neg.l %[t2] \n\t" /* t2 = shift - 31 */ \
|
|
|
|
|
"asr.l %[t2],%[t1] \n\t" /* hi >>= t2 */ \
|
|
|
|
|
"2: \n" \
|
|
|
|
|
: [t1]"=&d"(t1), [t2]"=&d"(t2) \
|
|
|
|
|
: [x] "d"((X)), [y] "d"((Y)), [sh]"d"((Z))); \
|
|
|
|
|
t1; \
|
|
|
|
|
})
|
2010-07-05 22:33:37 +00:00
|
|
|
|
2010-07-26 21:43:07 +00:00
|
|
|
/* Calculates: result = (X*Y)>>16 */
|
|
|
|
|
#define fixmul16(X,Y) \
|
|
|
|
|
({ \
|
2010-07-26 23:00:22 +00:00
|
|
|
int32_t t, x = (X); \
|
2010-07-26 21:43:07 +00:00
|
|
|
asm volatile ( \
|
2010-07-26 23:00:22 +00:00
|
|
|
"mac.l %[x],%[y],%%acc0\n\t" /* multiply */ \
|
|
|
|
|
"mulu.l %[y],%[x] \n\t" /* get lower half, avoid emac stall */ \
|
|
|
|
|
"movclr.l %%acc0,%[t] \n\t" /* get higher half */ \
|
|
|
|
|
"lsr.l #1,%[t] \n\t" /* hi >>= 1 to compensate emac shift */ \
|
|
|
|
|
"move.w %[t],%[x] \n\t" /* combine halfwords */\
|
|
|
|
|
"swap %[x] \n\t" \
|
|
|
|
|
: [t]"=&d"(t), [x] "+d" (x) \
|
|
|
|
|
: [y] "d" ((Y))); \
|
|
|
|
|
x; \
|
2010-07-26 21:43:07 +00:00
|
|
|
})
|
2010-07-05 22:33:37 +00:00
|
|
|
|
2010-08-05 10:35:06 +00:00
|
|
|
/* Calculates: result = (X*Y)>>31 (may lose msb to overflow) */
|
2010-07-26 21:43:07 +00:00
|
|
|
#define fixmul31(X,Y) \
|
|
|
|
|
({ \
|
|
|
|
|
int32_t t; \
|
|
|
|
|
asm volatile ( \
|
|
|
|
|
"mac.l %[x], %[y], %%acc0\n\t" /* multiply */ \
|
|
|
|
|
"movclr.l %%acc0, %[t]\n\t" /* get higher half as result */ \
|
|
|
|
|
: [t] "=d" (t) \
|
|
|
|
|
: [x] "r" ((X)), [y] "r" ((Y))); \
|
|
|
|
|
t; \
|
|
|
|
|
})
|
|
|
|
|
#else
|
|
|
|
|
static inline int32_t fixmulshift(int32_t x, int32_t y, int shamt)
|
|
|
|
|
{
|
|
|
|
|
int64_t temp;
|
|
|
|
|
temp = x;
|
|
|
|
|
temp *= y;
|
|
|
|
|
|
|
|
|
|
temp >>= shamt;
|
|
|
|
|
|
|
|
|
|
return (int32_t)temp;
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
static inline int32_t fixmul31(int32_t x, int32_t y)
|
|
|
|
|
{
|
|
|
|
|
int64_t temp;
|
|
|
|
|
temp = x;
|
|
|
|
|
temp *= y;
|
|
|
|
|
|
|
|
|
|
temp >>= 31;
|
|
|
|
|
|
|
|
|
|
return (int32_t)temp;
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
static inline int32_t fixmul24(int32_t x, int32_t y)
|
|
|
|
|
{
|
|
|
|
|
int64_t temp;
|
|
|
|
|
temp = x;
|
|
|
|
|
temp *= y;
|
|
|
|
|
|
|
|
|
|
temp >>= 24;
|
|
|
|
|
|
|
|
|
|
return (int32_t)temp;
|
|
|
|
|
}
|
|
|
|
|
|
|
|
|
|
static inline int32_t fixmul16(int32_t x, int32_t y)
|
|
|
|
|
{
|
|
|
|
|
int64_t temp;
|
|
|
|
|
temp = x;
|
|
|
|
|
temp *= y;
|
|
|
|
|
|
|
|
|
|
temp >>= 16;
|
|
|
|
|
|
|
|
|
|
return (int32_t)temp;
|
|
|
|
|
}
|
|
|
|
|
#endif /* CPU_COLDFIRE, CPU_ARM */
|
2010-07-05 22:33:37 +00:00
|
|
|
|
2010-07-28 19:36:15 +00:00
|
|
|
#if defined(CPU_COLDFIRE)
|
2010-07-28 20:46:51 +00:00
|
|
|
#define VECT_MUL_WIN_KERNEL(i, j, s0, s1, wi, wj) \
|
|
|
|
|
asm volatile ( \
|
|
|
|
|
"mac.l %[s0], %[wj], %%acc0 \n\t" \
|
|
|
|
|
"msac.l %[s1], %[wi], %%acc0 \n\t" \
|
|
|
|
|
"mac.l %[s0], %[wi], %%acc1 \n\t" \
|
|
|
|
|
"mac.l %[s1], %[wj], %%acc1 \n\t" \
|
|
|
|
|
"movclr.l %%acc0, %[s0] \n\t" \
|
|
|
|
|
"move.l %[s0], (%[dst_i]) \n\t" \
|
|
|
|
|
"movclr.l %%acc1, %[s0] \n\t" \
|
|
|
|
|
"move.l %[s0], (%[dst_j]) \n\t" \
|
|
|
|
|
: [s0] "+r" (s0) /* register is clobbered so specify it as an input */ \
|
|
|
|
|
: [dst_i] "a" (&dst[i]), [dst_j] "a" (&dst[j]), \
|
|
|
|
|
[s1] "r" (s1), [wi] "r" (wi), [wj] "r" (wj) \
|
|
|
|
|
: "cc", "memory");
|
2010-07-26 11:15:25 +00:00
|
|
|
#else
|
2010-07-28 20:46:51 +00:00
|
|
|
#define VECT_MUL_WIN_KERNEL(i, j, s0, s1, wi, wj) \
|
|
|
|
|
dst[i] = fixmul31(s0, wj) - fixmul31(s1, wi); \
|
|
|
|
|
dst[j] = fixmul31(s0, wi) + fixmul31(s1, wj);
|
|
|
|
|
#endif /* CPU_COLDFIRE */
|
|
|
|
|
|
2010-07-17 08:00:13 +00:00
|
|
|
static inline void vector_fixmul_window(int32_t *dst, const int32_t *src0,
|
|
|
|
|
const int32_t *src1, const int32_t *win,
|
2010-07-26 11:15:25 +00:00
|
|
|
int len)
|
2010-07-05 22:33:37 +00:00
|
|
|
{
|
|
|
|
|
int i, j;
|
|
|
|
|
dst += len;
|
|
|
|
|
win += len;
|
|
|
|
|
src0+= len;
|
2010-07-28 18:15:53 +00:00
|
|
|
for(i=-len, j=len-1; i<0; i++, j--) {
|
|
|
|
|
int32_t s0 = src0[i]; /* s0 = src0[ 0 ... len-1] */
|
|
|
|
|
int32_t s1 = src1[j]; /* s1 = src1[2*len-1 ... len] */
|
|
|
|
|
int32_t wi = -win[i]; /* wi = -win[ 0 ... len-1] */
|
|
|
|
|
int32_t wj = -win[j]; /* wj = -win[2*len-1 ... len] */
|
2010-07-28 20:46:51 +00:00
|
|
|
VECT_MUL_WIN_KERNEL(i, j, s0, s1, wi, wj);
|
2010-07-26 11:15:25 +00:00
|
|
|
}
|
2010-07-05 22:33:37 +00:00
|
|
|
}
|
|
|
|
|
|
2010-07-28 19:36:15 +00:00
|
|
|
#if defined(CPU_ARM)
|
2010-07-28 20:46:51 +00:00
|
|
|
#define VECT_MUL_SCALAR_KERNEL(dst, src, mul) \
|
|
|
|
|
asm volatile ( \
|
|
|
|
|
"ldmia %[src]!, {r1-r4} \n\t" \
|
|
|
|
|
"smull r0, r5, r1, %[mul] \n\t" \
|
2010-08-04 23:21:35 +00:00
|
|
|
"mov r0, r0, lsr #16 \n\t" \
|
|
|
|
|
"orr r0, r0, r5, lsl #16\n\t" \
|
2010-07-28 20:46:51 +00:00
|
|
|
"smull r1, r5, r2, %[mul] \n\t" \
|
2010-08-04 23:21:35 +00:00
|
|
|
"mov r1, r1, lsr #16 \n\t" \
|
|
|
|
|
"orr r1, r1, r5, lsl #16\n\t" \
|
2010-07-28 20:46:51 +00:00
|
|
|
"smull r2, r5, r3, %[mul] \n\t" \
|
2010-08-04 23:21:35 +00:00
|
|
|
"mov r2, r2, lsr #16 \n\t" \
|
|
|
|
|
"orr r2, r2, r5, lsl #16\n\t" \
|
2010-07-28 20:46:51 +00:00
|
|
|
"smull r3, r5, r4, %[mul] \n\t" \
|
2010-08-04 23:21:35 +00:00
|
|
|
"mov r3, r3, lsr #16 \n\t" \
|
|
|
|
|
"orr r3, r3, r5, lsl #16\n\t" \
|
2010-07-28 20:46:51 +00:00
|
|
|
"stmia %[dst]!, {r0-r3} \n" \
|
|
|
|
|
: [dst]"+r"(dst), [src]"+r"(src) \
|
|
|
|
|
: [mul]"r"(mul) \
|
2010-07-28 19:36:15 +00:00
|
|
|
: "r0", "r1", "r2", "r3", "r4", "r5", "memory");
|
2010-08-04 23:21:35 +00:00
|
|
|
/* Disable ColdFire version until a correct version is written
|
2010-07-28 23:25:45 +00:00
|
|
|
#elif defined (CPU_COLDFIRE)
|
|
|
|
|
#define VECT_MUL_SCALAR_KERNEL(dst, src, mul) \
|
|
|
|
|
int32_t tmp; \
|
|
|
|
|
asm volatile ( \
|
|
|
|
|
"movem.l (%[src]), %%d0-%%d3 \n\t" \
|
|
|
|
|
"mac.l %[mul], %%d0, %%acc0 \n\t" \
|
|
|
|
|
"mac.l %[mul], %%d1, %%acc1 \n\t" \
|
|
|
|
|
"mac.l %[mul], %%d2, %%acc2 \n\t" \
|
|
|
|
|
"mac.l %[mul], %%d3, %%acc3 \n\t" \
|
|
|
|
|
"move.l %%accext01, %[tmp] \n\t" \
|
|
|
|
|
"movclr.l %%acc0, %%d0 \n\t" \
|
|
|
|
|
"movclr.l %%acc1, %%d1 \n\t" \
|
|
|
|
|
"lsl.l #7, %%d0 \n\t" \
|
|
|
|
|
"move.b %[tmp], %%d0 \n\t" \
|
|
|
|
|
"swap %[tmp] \n\t" \
|
|
|
|
|
"lsl.l #7, %%d1 \n\t" \
|
|
|
|
|
"move.b %[tmp], %%d1 \n\t" \
|
|
|
|
|
"move.l %%accext23, %[tmp] \n\t" \
|
|
|
|
|
"movclr.l %%acc2, %%d2 \n\t" \
|
|
|
|
|
"movclr.l %%acc3, %%d3 \n\t" \
|
|
|
|
|
"lsl.l #7, %%d2 \n\t" \
|
|
|
|
|
"move.b %[tmp], %%d2 \n\t" \
|
|
|
|
|
"swap %[tmp] \n\t" \
|
|
|
|
|
"lsl.l #7, %%d3 \n\t" \
|
|
|
|
|
"move.b %[tmp], %%d3 \n\t" \
|
|
|
|
|
"movem.l %%d0-%%d3, (%[dst]) \n\t" \
|
|
|
|
|
"lea.l (4*4, %[src]), %[src]\n\t" \
|
|
|
|
|
"lea.l (4*4, %[dst]), %[dst]\n\t" \
|
|
|
|
|
: [dst] "+a" (dst), [src] "+a" (src),\
|
|
|
|
|
[tmp] "=d" (tmp) \
|
|
|
|
|
: [mul] "r" (mul) \
|
|
|
|
|
: "d0", "d1", "d2", "d3", "memory", "cc");
|
2010-08-04 23:21:35 +00:00
|
|
|
*/
|
2010-07-28 19:36:15 +00:00
|
|
|
#else
|
2010-07-28 20:46:51 +00:00
|
|
|
#define VECT_MUL_SCALAR_KERNEL(dst, src, mul) \
|
2010-08-04 23:21:35 +00:00
|
|
|
dst[i ] = fixmul16(src[i ], mul); \
|
|
|
|
|
dst[i+1] = fixmul16(src[i+1], mul); \
|
|
|
|
|
dst[i+2] = fixmul16(src[i+2], mul); \
|
|
|
|
|
dst[i+3] = fixmul16(src[i+3], mul);
|
2010-07-29 22:18:04 +00:00
|
|
|
#endif /* CPU_ARM, CPU_COLDFIRE */
|
2010-07-28 20:46:51 +00:00
|
|
|
|
2010-07-26 21:43:07 +00:00
|
|
|
static inline void vector_fixmul_scalar(int32_t *dst, const int32_t *src,
|
|
|
|
|
int32_t mul, int len)
|
2010-07-05 22:33:37 +00:00
|
|
|
{
|
2010-07-28 18:15:53 +00:00
|
|
|
/* len is _always_ a multiple of 4, because len is the difference of sfb's
|
|
|
|
|
* which themselves are always a multiple of 4. */
|
2010-07-05 22:33:37 +00:00
|
|
|
int i;
|
2010-07-28 18:15:53 +00:00
|
|
|
for (i=0; i<len; i+=4) {
|
2010-07-28 20:46:51 +00:00
|
|
|
VECT_MUL_SCALAR_KERNEL(dst, src, mul);
|
2010-07-28 18:15:53 +00:00
|
|
|
}
|
2010-07-05 22:33:37 +00:00
|
|
|
}
|
2010-07-12 15:41:10 +00:00
|
|
|
|
2010-07-15 05:38:09 +00:00
|
|
|
static inline int av_clip(int a, int amin, int amax)
|
|
|
|
|
{
|
|
|
|
|
if (a < amin) return amin;
|
|
|
|
|
else if (a > amax) return amax;
|
|
|
|
|
else return a;
|
|
|
|
|
}
|
2010-07-12 15:41:10 +00:00
|
|
|
#endif /* _WMAPRO_MATH_H_ */
|
