; ; Copyright (c) 2010 The WebM project authors. All Rights Reserved. ; ; Use of this source code is governed by a BSD-style license ; that can be found in the LICENSE file in the root of the source ; tree. An additional intellectual property rights grant can be found ; in the file PATENTS. All contributing project authors may ; be found in the AUTHORS file in the root of the source tree. ; EXPORT |vp8_short_idct4x4llm_v6_dual| AREA |.text|, CODE, READONLY ; void vp8_short_idct4x4llm_c(short *input, unsigned char *pred, int pitch, ; unsigned char *dst, int stride) ; r0 short* input ; r1 unsigned char* pred ; r2 int pitch ; r3 unsigned char* dst ; sp int stride |vp8_short_idct4x4llm_v6_dual| PROC stmdb sp!, {r4-r11, lr} sub sp, sp, #4 mov r4, #0x00008A00 ; sin orr r4, r4, #0x0000008C ; sinpi8sqrt2 mov r5, #0x00004E00 ; cos orr r5, r5, #0x0000007B ; cospi8sqrt2minus1 orr r5, r5, #1<<31 ; loop counter on top bit loop1_dual ldr r6, [r0, #(4*2)] ; i5 | i4 ldr r12, [r0, #(12*2)] ; i13|i12 ldr r14, [r0, #(8*2)] ; i9 | i8 smulbt r9, r5, r6 ; (ip[5] * cospi8sqrt2minus1) >> 16 smulbb r7, r5, r6 ; (ip[4] * cospi8sqrt2minus1) >> 16 smulwt r10, r4, r6 ; (ip[5] * sinpi8sqrt2) >> 16 smulwb r8, r4, r6 ; (ip[4] * sinpi8sqrt2) >> 16 smulbt r11, r5, r12 ; (ip[13] * cospi8sqrt2minus1) >> 16 pkhtb r7, r9, r7, asr #16 ; 5c | 4c pkhbt r8, r8, r10, lsl #16 ; 5s | 4s uadd16 r6, r6, r7 ; 5c+5 | 4c+4 smulwt r7, r4, r12 ; (ip[13] * sinpi8sqrt2) >> 16 smulbb r9, r5, r12 ; (ip[12] * cospi8sqrt2minus1) >> 16 smulwb r10, r4, r12 ; (ip[12] * sinpi8sqrt2) >> 16 subs r5, r5, #1<<31 ; i-- pkhtb r9, r11, r9, asr #16 ; 13c | 12c ldr r11, [r0] ; i1 | i0 pkhbt r10, r10, r7, lsl #16 ; 13s | 12s uadd16 r7, r12, r9 ; 13c+13 | 12c+12 usub16 r7, r8, r7 ; c uadd16 r6, r6, r10 ; d uadd16 r10, r11, r14 ; a usub16 r8, r11, r14 ; b uadd16 r9, r10, r6 ; a+d usub16 r10, r10, r6 ; a-d uadd16 r6, r8, r7 ; b+c usub16 r7, r8, r7 ; b-c ; use input buffer to store intermediate results str r6, [r0, #(4*2)] ; o5 | o4 str r7, [r0, #(8*2)] ; o9 | o8 str r10,[r0, #(12*2)] ; o13|o12 str r9, [r0], #4 ; o1 | o0 bcs loop1_dual sub r0, r0, #8 ; reset input/output str r0, [sp] loop2_dual ldr r6, [r0, #(4*2)] ; i5 | i4 ldr r12,[r0, #(2*2)] ; i3 | i2 ldr r14,[r0, #(6*2)] ; i7 | i6 ldr r0, [r0, #(0*2)] ; i1 | i0 smulbt r9, r5, r6 ; (ip[5] * cospi8sqrt2minus1) >> 16 smulbt r7, r5, r0 ; (ip[1] * cospi8sqrt2minus1) >> 16 smulwt r10, r4, r6 ; (ip[5] * sinpi8sqrt2) >> 16 smulwt r8, r4, r0 ; (ip[1] * sinpi8sqrt2) >> 16 pkhbt r11, r6, r0, lsl #16 ; i0 | i4 pkhtb r7, r7, r9, asr #16 ; 1c | 5c pkhtb r0, r0, r6, asr #16 ; i1 | i5 pkhbt r8, r10, r8, lsl #16 ; 1s | 5s = temp1 uadd16 r0, r7, r0 ; 1c+1 | 5c+5 = temp2 pkhbt r9, r14, r12, lsl #16 ; i2 | i6 uadd16 r10, r11, r9 ; a usub16 r9, r11, r9 ; b pkhtb r6, r12, r14, asr #16 ; i3 | i7 subs r5, r5, #1<<31 ; i-- smulbt r7, r5, r6 ; (ip[3] * cospi8sqrt2minus1) >> 16 smulwt r11, r4, r6 ; (ip[3] * sinpi8sqrt2) >> 16 smulbb r12, r5, r6 ; (ip[7] * cospi8sqrt2minus1) >> 16 smulwb r14, r4, r6 ; (ip[7] * sinpi8sqrt2) >> 16 pkhtb r7, r7, r12, asr #16 ; 3c | 7c pkhbt r11, r14, r11, lsl #16 ; 3s | 7s = temp1 uadd16 r6, r7, r6 ; 3c+3 | 7c+7 = temp2 usub16 r12, r8, r6 ; c (o1 | o5) uadd16 r6, r11, r0 ; d (o3 | o7) uadd16 r7, r10, r6 ; a+d mov r8, #4 ; set up 4's orr r8, r8, #0x40000 ; 4|4 usub16 r6, r10, r6 ; a-d uadd16 r6, r6, r8 ; a-d+4, 3|7 uadd16 r7, r7, r8 ; a+d+4, 0|4 uadd16 r10, r9, r12 ; b+c usub16 r0, r9, r12 ; b-c uadd16 r10, r10, r8 ; b+c+4, 1|5 uadd16 r8, r0, r8 ; b-c+4, 2|6 ldr lr, [sp, #40] ; dst stride ldrb r0, [r1] ; pred p0 ldrb r11, [r1, #1] ; pred p1 ldrb r12, [r1, #2] ; pred p2 add r0, r0, r7, asr #19 ; p0 + o0 add r11, r11, r10, asr #19 ; p1 + o1 add r12, r12, r8, asr #19 ; p2 + o2 usat r0, #8, r0 ; d0 = clip8(p0 + o0) usat r11, #8, r11 ; d1 = clip8(p1 + o1) usat r12, #8, r12 ; d2 = clip8(p2 + o2) add r0, r0, r11, lsl #8 ; |--|--|d1|d0| ldrb r11, [r1, #3] ; pred p3 add r0, r0, r12, lsl #16 ; |--|d2|d1|d0| add r11, r11, r6, asr #19 ; p3 + o3 sxth r7, r7 ; sxth r10, r10 ; usat r11, #8, r11 ; d3 = clip8(p3 + o3) sxth r8, r8 ; sxth r6, r6 ; add r0, r0, r11, lsl #24 ; |d3|d2|d1|d0| ldrb r12, [r1, r2]! ; pred p4 str r0, [r3], lr ldrb r11, [r1, #1] ; pred p5 add r12, r12, r7, asr #3 ; p4 + o4 add r11, r11, r10, asr #3 ; p5 + o5 usat r12, #8, r12 ; d4 = clip8(p4 + o4) usat r11, #8, r11 ; d5 = clip8(p5 + o5) ldrb r7, [r1, #2] ; pred p6 ldrb r10, [r1, #3] ; pred p6 add r12, r12, r11, lsl #8 ; |--|--|d5|d4| add r7, r7, r8, asr #3 ; p6 + o6 add r10, r10, r6, asr #3 ; p7 + o7 ldr r0, [sp] ; load input pointer usat r7, #8, r7 ; d6 = clip8(p6 + o6) usat r10, #8, r10 ; d7 = clip8(p7 + o7) add r12, r12, r7, lsl #16 ; |--|d6|d5|d4| add r12, r12, r10, lsl #24 ; |d7|d6|d5|d4| str r12, [r3], lr add r0, r0, #16 add r1, r1, r2 ; pred + pitch bcs loop2_dual add sp, sp, #4 ; idct_output buffer ldmia sp!, {r4 - r11, pc} ENDP END