← Back to branch summary

~ubuntu-branches/ubuntu/jaunty/xvidcap/jaunty-proposed

« back to all changes in this revision

Viewing changes to ffmpeg/libavcodec/armv4l/jrevdct_arm.S

  • Committer: Bazaar Package Importer
  • Author(s): John Dong
  • Date: 2008-02-25 15:47:12 UTC
  • mfrom: (1.1.1 upstream)
  • Revision ID: james.westby@ubuntu.com-20080225154712-qvr11ekcea4c9ry8
Tags: 1.1.6-0.1ubuntu1
https://launchpad.net/bugs/120003
* Merge from debian-multimedia (LP: #120003), Ubuntu Changes:
 - For ffmpeg-related build-deps, remove cvs from package names.
 - Standards-Version 3.7.3
 - Maintainer Spec

Show diffs side-by-side

added added

removed removed

Lines of Context:
ffmpeg/libavcodec/armv4l/jrevdct_arm.S
1
 
/* 
 
1
/*
2
2
   C-like prototype :
3
 
        void j_rev_dct_ARM(DCTBLOCK data)
 
3
        void j_rev_dct_ARM(DCTBLOCK data)
4
4
 
5
5
   With DCTBLOCK being a pointer to an array of 64 'signed shorts'
6
6
 
22
22
   COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER
23
23
   IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN
24
24
   CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
25
 
        
 
25
 
26
26
*/
27
27
#define FIX_0_298631336 2446
28
28
#define FIX_0_541196100 4433
36
36
#define FIX_M_1_847759065 -15137
37
37
#define FIX_M_1_961570560 -16069
38
38
#define FIX_M_2_562915447 -20995
39
 
#define FIX_0xFFFF 0xFFFF       
40
 
                
 
39
#define FIX_0xFFFF 0xFFFF
 
40
 
41
41
#define FIX_0_298631336_ID      0
42
42
#define FIX_0_541196100_ID      4
43
43
#define FIX_0_765366865_ID      8
51
51
#define FIX_M_1_961570560_ID   40
52
52
#define FIX_M_2_562915447_ID   44
53
53
#define FIX_0xFFFF_ID          48
54
 
        .text
55
 
        .align
56
 
        
57
 
        .global j_rev_dct_ARM
 
54
        .text
 
55
        .align
 
56
 
 
57
        .global j_rev_dct_ARM
58
58
j_rev_dct_ARM:
59
 
        stmdb   sp!, { r4 - r12, lr }   @ all callee saved regs
60
 
 
61
 
        sub sp, sp, #4                  @ reserve some space on the stack
62
 
        str r0, [ sp ]                  @ save the DCT pointer to the stack
63
 
 
64
 
        mov lr, r0                      @ lr = pointer to the current row
65
 
        mov r12, #8                     @ r12 = row-counter
66
 
        add r11, pc, #(const_array-.-8) @ r11 = base pointer to the constants array     
 
59
        stmdb   sp!, { r4 - r12, lr }   @ all callee saved regs
 
60
 
 
61
        sub sp, sp, #4                  @ reserve some space on the stack
 
62
        str r0, [ sp ]                  @ save the DCT pointer to the stack
 
63
 
 
64
        mov lr, r0                      @ lr = pointer to the current row
 
65
        mov r12, #8                     @ r12 = row-counter
 
66
        add r11, pc, #(const_array-.-8) @ r11 = base pointer to the constants array
67
67
row_loop:
68
 
        ldrsh r0, [lr, # 0]             @ r0 = 'd0'
69
 
        ldrsh r1, [lr, # 8]             @ r1 = 'd1'
70
 
 
71
 
        @ Optimization for row that have all items except the first set to 0
72
 
        @ (this works as the DCTELEMS are always 4-byte aligned)
73
 
        ldr r5, [lr, # 0]
74
 
        ldr r2, [lr, # 4]
75
 
        ldr r3, [lr, # 8]
76
 
        ldr r4, [lr, #12]
77
 
        orr r3, r3, r4
78
 
        orr r3, r3, r2
79
 
        orrs r5, r3, r5
80
 
        beq end_of_row_loop             @ nothing to be done as ALL of them are '0'
81
 
        orrs r2, r3, r1
82
 
        beq empty_row
83
 
        
84
 
        ldrsh r2, [lr, # 2]             @ r2 = 'd2'
85
 
        ldrsh r4, [lr, # 4]             @ r4 = 'd4'
86
 
        ldrsh r6, [lr, # 6]             @ r6 = 'd6'
87
 
        
88
 
        ldr r3, [r11, #FIX_0_541196100_ID]
89
 
        add r7, r2, r6
90
 
        ldr r5, [r11, #FIX_M_1_847759065_ID]
91
 
        mul r7, r3, r7                      @ r7 = z1
92
 
        ldr r3, [r11, #FIX_0_765366865_ID]
93
 
        mla r6, r5, r6, r7                  @ r6 = tmp2
94
 
        add r5, r0, r4                      @ r5 = tmp0
95
 
        mla r2, r3, r2, r7                  @ r2 = tmp3
96
 
        sub r3, r0, r4                      @ r3 = tmp1
97
 
 
98
 
        add r0, r2, r5, lsl #13             @ r0 = tmp10
99
 
        rsb r2, r2, r5, lsl #13             @ r2 = tmp13
100
 
        add r4, r6, r3, lsl #13             @ r4 = tmp11
101
 
        rsb r3, r6, r3, lsl #13             @ r3 = tmp12
102
 
 
103
 
        stmdb   sp!, { r0, r2, r3, r4 } @ save on the stack tmp10, tmp13, tmp12, tmp11
104
 
        
105
 
        ldrsh r3, [lr, #10]             @ r3 = 'd3'
106
 
        ldrsh r5, [lr, #12]             @ r5 = 'd5'
107
 
        ldrsh r7, [lr, #14]             @ r7 = 'd7'
108
 
 
109
 
        add r0, r3, r5                  @ r0 = 'z2'
110
 
        add r2, r1, r7                  @ r2 = 'z1'
111
 
        add r4, r3, r7                  @ r4 = 'z3'
112
 
        add r6, r1, r5                  @ r6 = 'z4'
113
 
        ldr r9, [r11, #FIX_1_175875602_ID]
114
 
        add r8, r4, r6                  @ r8 = z3 + z4
115
 
        ldr r10, [r11, #FIX_M_0_899976223_ID]
116
 
        mul r8, r9, r8                  @ r8 = 'z5'
117
 
        ldr r9, [r11, #FIX_M_2_562915447_ID]
118
 
        mul r2, r10, r2                 @ r2 = 'z1'
119
 
        ldr r10, [r11, #FIX_M_1_961570560_ID]
120
 
        mul r0, r9, r0                  @ r0 = 'z2'
121
 
        ldr r9, [r11, #FIX_M_0_390180644_ID]
122
 
        mla r4, r10, r4, r8             @ r4 = 'z3'
123
 
        ldr r10, [r11, #FIX_0_298631336_ID]
124
 
        mla r6, r9, r6, r8              @ r6 = 'z4'
125
 
        ldr r9, [r11, #FIX_2_053119869_ID]
126
 
        mla r7, r10, r7, r2             @ r7 = tmp0 + z1
127
 
        ldr r10, [r11, #FIX_3_072711026_ID]
128
 
        mla r5, r9, r5, r0              @ r5 = tmp1 + z2
129
 
        ldr r9, [r11, #FIX_1_501321110_ID]
130
 
        mla r3, r10, r3, r0             @ r3 = tmp2 + z2
131
 
        add r7, r7, r4                  @ r7 = tmp0
132
 
        mla r1, r9, r1, r2              @ r1 = tmp3 + z1
133
 
        add r5, r5, r6                  @ r5 = tmp1
134
 
        add r3, r3, r4                  @ r3 = tmp2
135
 
        add r1, r1, r6                  @ r1 = tmp3
136
 
 
137
 
        ldmia sp!, { r0, r2, r4, r6 } @ r0 = tmp10 / r2 = tmp13 / r4 = tmp12 / r6 = tmp11
138
 
                                      @ r1 = tmp3  / r3 = tmp2  / r5 = tmp1  / r7 = tmp0
139
 
        
140
 
        @ Compute DESCALE(tmp10 + tmp3, CONST_BITS-PASS1_BITS)
141
 
        add r8, r0, r1
142
 
        add r8, r8, #(1<<10)
143
 
        mov r8, r8, asr #11
144
 
        strh r8, [lr, # 0]
145
 
        
146
 
        @ Compute DESCALE(tmp10 - tmp3, CONST_BITS-PASS1_BITS)
147
 
        sub r8, r0, r1
148
 
        add r8, r8, #(1<<10)
149
 
        mov r8, r8, asr #11
150
 
        strh r8, [lr, #14]
151
 
        
152
 
        @ Compute DESCALE(tmp11 + tmp2, CONST_BITS-PASS1_BITS)
153
 
        add r8, r6, r3
154
 
        add r8, r8, #(1<<10)
155
 
        mov r8, r8, asr #11
156
 
        strh r8, [lr, # 2]
157
 
        
158
 
        @ Compute DESCALE(tmp11 - tmp2, CONST_BITS-PASS1_BITS)
159
 
        sub r8, r6, r3
160
 
        add r8, r8, #(1<<10)
161
 
        mov r8, r8, asr #11
162
 
        strh r8, [lr, #12]
163
 
        
164
 
        @ Compute DESCALE(tmp12 + tmp1, CONST_BITS-PASS1_BITS)
165
 
        add r8, r4, r5
166
 
        add r8, r8, #(1<<10)
167
 
        mov r8, r8, asr #11
168
 
        strh r8, [lr, # 4]
169
 
        
170
 
        @ Compute DESCALE(tmp12 - tmp1, CONST_BITS-PASS1_BITS)
171
 
        sub r8, r4, r5
172
 
        add r8, r8, #(1<<10)
173
 
        mov r8, r8, asr #11
174
 
        strh r8, [lr, #10]
175
 
        
176
 
        @ Compute DESCALE(tmp13 + tmp0, CONST_BITS-PASS1_BITS)
177
 
        add r8, r2, r7
178
 
        add r8, r8, #(1<<10)
179
 
        mov r8, r8, asr #11
180
 
        strh r8, [lr, # 6]
181
 
        
182
 
        @ Compute DESCALE(tmp13 - tmp0, CONST_BITS-PASS1_BITS)
183
 
        sub r8, r2, r7
184
 
        add r8, r8, #(1<<10)
185
 
        mov r8, r8, asr #11
186
 
        strh r8, [lr, # 8]
187
 
 
188
 
        @ End of row loop
189
 
        add lr, lr, #16
190
 
        subs r12, r12, #1
191
 
        bne row_loop
192
 
        beq start_column_loop
193
 
        
 
68
        ldrsh r0, [lr, # 0]             @ r0 = 'd0'
 
69
        ldrsh r1, [lr, # 8]             @ r1 = 'd1'
 
70
 
 
71
        @ Optimization for row that have all items except the first set to 0
 
72
        @ (this works as the DCTELEMS are always 4-byte aligned)
 
73
        ldr r5, [lr, # 0]
 
74
        ldr r2, [lr, # 4]
 
75
        ldr r3, [lr, # 8]
 
76
        ldr r4, [lr, #12]
 
77
        orr r3, r3, r4
 
78
        orr r3, r3, r2
 
79
        orrs r5, r3, r5
 
80
        beq end_of_row_loop             @ nothing to be done as ALL of them are '0'
 
81
        orrs r2, r3, r1
 
82
        beq empty_row
 
83
 
 
84
        ldrsh r2, [lr, # 2]             @ r2 = 'd2'
 
85
        ldrsh r4, [lr, # 4]             @ r4 = 'd4'
 
86
        ldrsh r6, [lr, # 6]             @ r6 = 'd6'
 
87
 
 
88
        ldr r3, [r11, #FIX_0_541196100_ID]
 
89
        add r7, r2, r6
 
90
        ldr r5, [r11, #FIX_M_1_847759065_ID]
 
91
        mul r7, r3, r7                      @ r7 = z1
 
92
        ldr r3, [r11, #FIX_0_765366865_ID]
 
93
        mla r6, r5, r6, r7                  @ r6 = tmp2
 
94
        add r5, r0, r4                      @ r5 = tmp0
 
95
        mla r2, r3, r2, r7                  @ r2 = tmp3
 
96
        sub r3, r0, r4                      @ r3 = tmp1
 
97
 
 
98
        add r0, r2, r5, lsl #13             @ r0 = tmp10
 
99
        rsb r2, r2, r5, lsl #13             @ r2 = tmp13
 
100
        add r4, r6, r3, lsl #13             @ r4 = tmp11
 
101
        rsb r3, r6, r3, lsl #13             @ r3 = tmp12
 
102
 
 
103
        stmdb   sp!, { r0, r2, r3, r4 } @ save on the stack tmp10, tmp13, tmp12, tmp11
 
104
 
 
105
        ldrsh r3, [lr, #10]             @ r3 = 'd3'
 
106
        ldrsh r5, [lr, #12]             @ r5 = 'd5'
 
107
        ldrsh r7, [lr, #14]             @ r7 = 'd7'
 
108
 
 
109
        add r0, r3, r5                        @ r0 = 'z2'
 
110
        add r2, r1, r7                  @ r2 = 'z1'
 
111
        add r4, r3, r7                  @ r4 = 'z3'
 
112
        add r6, r1, r5                  @ r6 = 'z4'
 
113
        ldr r9, [r11, #FIX_1_175875602_ID]
 
114
        add r8, r4, r6                  @ r8 = z3 + z4
 
115
        ldr r10, [r11, #FIX_M_0_899976223_ID]
 
116
        mul r8, r9, r8                  @ r8 = 'z5'
 
117
        ldr r9, [r11, #FIX_M_2_562915447_ID]
 
118
        mul r2, r10, r2                 @ r2 = 'z1'
 
119
        ldr r10, [r11, #FIX_M_1_961570560_ID]
 
120
        mul r0, r9, r0                  @ r0 = 'z2'
 
121
        ldr r9, [r11, #FIX_M_0_390180644_ID]
 
122
        mla r4, r10, r4, r8             @ r4 = 'z3'
 
123
        ldr r10, [r11, #FIX_0_298631336_ID]
 
124
        mla r6, r9, r6, r8              @ r6 = 'z4'
 
125
        ldr r9, [r11, #FIX_2_053119869_ID]
 
126
        mla r7, r10, r7, r2             @ r7 = tmp0 + z1
 
127
        ldr r10, [r11, #FIX_3_072711026_ID]
 
128
        mla r5, r9, r5, r0              @ r5 = tmp1 + z2
 
129
        ldr r9, [r11, #FIX_1_501321110_ID]
 
130
        mla r3, r10, r3, r0             @ r3 = tmp2 + z2
 
131
        add r7, r7, r4                  @ r7 = tmp0
 
132
        mla r1, r9, r1, r2              @ r1 = tmp3 + z1
 
133
        add r5,        r5, r6                  @ r5 = tmp1
 
134
        add r3, r3, r4                  @ r3 = tmp2
 
135
        add r1, r1, r6                  @ r1 = tmp3
 
136
 
 
137
        ldmia sp!, { r0, r2, r4, r6 } @ r0 = tmp10 / r2 = tmp13 / r4 = tmp12 / r6 = tmp11
 
138
                                      @ r1 = tmp3  / r3 = tmp2  / r5 = tmp1  / r7 = tmp0
 
139
 
 
140
        @ Compute DESCALE(tmp10 + tmp3, CONST_BITS-PASS1_BITS)
 
141
        add r8, r0, r1
 
142
        add r8, r8, #(1<<10)
 
143
        mov r8, r8, asr #11
 
144
        strh r8, [lr, # 0]
 
145
 
 
146
        @ Compute DESCALE(tmp10 - tmp3, CONST_BITS-PASS1_BITS)
 
147
        sub r8, r0, r1
 
148
        add r8, r8, #(1<<10)
 
149
        mov r8, r8, asr #11
 
150
        strh r8, [lr, #14]
 
151
 
 
152
        @ Compute DESCALE(tmp11 + tmp2, CONST_BITS-PASS1_BITS)
 
153
        add r8, r6, r3
 
154
        add r8, r8, #(1<<10)
 
155
        mov r8, r8, asr #11
 
156
        strh r8, [lr, # 2]
 
157
 
 
158
        @ Compute DESCALE(tmp11 - tmp2, CONST_BITS-PASS1_BITS)
 
159
        sub r8, r6, r3
 
160
        add r8, r8, #(1<<10)
 
161
        mov r8, r8, asr #11
 
162
        strh r8, [lr, #12]
 
163
 
 
164
        @ Compute DESCALE(tmp12 + tmp1, CONST_BITS-PASS1_BITS)
 
165
        add r8, r4, r5
 
166
        add r8, r8, #(1<<10)
 
167
        mov r8, r8, asr #11
 
168
        strh r8, [lr, # 4]
 
169
 
 
170
        @ Compute DESCALE(tmp12 - tmp1, CONST_BITS-PASS1_BITS)
 
171
        sub r8, r4, r5
 
172
        add r8, r8, #(1<<10)
 
173
        mov r8, r8, asr #11
 
174
        strh r8, [lr, #10]
 
175
 
 
176
        @ Compute DESCALE(tmp13 + tmp0, CONST_BITS-PASS1_BITS)
 
177
        add r8, r2, r7
 
178
        add r8, r8, #(1<<10)
 
179
        mov r8, r8, asr #11
 
180
        strh r8, [lr, # 6]
 
181
 
 
182
        @ Compute DESCALE(tmp13 - tmp0, CONST_BITS-PASS1_BITS)
 
183
        sub r8, r2, r7
 
184
        add r8, r8, #(1<<10)
 
185
        mov r8, r8, asr #11
 
186
        strh r8, [lr, # 8]
 
187
 
 
188
        @ End of row loop
 
189
        add lr, lr, #16
 
190
        subs r12, r12, #1
 
191
        bne row_loop
 
192
        beq start_column_loop
 
193
 
194
194
empty_row:
195
 
        ldr r1, [r11, #FIX_0xFFFF_ID]
196
 
        mov r0, r0, lsl #2
197
 
        and r0, r0, r1
198
 
        add r0, r0, r0, lsl #16
199
 
        str r0, [lr, # 0]
200
 
        str r0, [lr, # 4]
201
 
        str r0, [lr, # 8]
202
 
        str r0, [lr, #12]
 
195
        ldr r1, [r11, #FIX_0xFFFF_ID]
 
196
        mov r0, r0, lsl #2
 
197
        and r0, r0, r1
 
198
        add r0, r0, r0, lsl #16
 
199
        str r0, [lr, # 0]
 
200
        str r0, [lr, # 4]
 
201
        str r0, [lr, # 8]
 
202
        str r0, [lr, #12]
203
203
 
204
204
end_of_row_loop:
205
 
        @ End of loop
206
 
        add lr, lr, #16
207
 
        subs r12, r12, #1
208
 
        bne row_loop
 
205
        @ End of loop
 
206
        add lr, lr, #16
 
207
        subs r12, r12, #1
 
208
        bne row_loop
209
209
 
210
210
start_column_loop:
211
 
        @ Start of column loop
212
 
        ldr lr, [ sp ]
213
 
        mov r12, #8
 
211
        @ Start of column loop
 
212
        ldr lr, [ sp ]
 
213
        mov r12, #8
214
214
column_loop:
215
 
        ldrsh r0, [lr, #( 0*8)]             @ r0 = 'd0'
216
 
        ldrsh r2, [lr, #( 4*8)]             @ r2 = 'd2'
217
 
        ldrsh r4, [lr, #( 8*8)]             @ r4 = 'd4'
218
 
        ldrsh r6, [lr, #(12*8)]             @ r6 = 'd6'
219
 
 
220
 
        ldr r3, [r11, #FIX_0_541196100_ID]
221
 
        add r1, r2, r6
222
 
        ldr r5, [r11, #FIX_M_1_847759065_ID]
223
 
        mul r1, r3, r1                      @ r1 = z1
224
 
        ldr r3, [r11, #FIX_0_765366865_ID]
225
 
        mla r6, r5, r6, r1                  @ r6 = tmp2
226
 
        add r5, r0, r4                      @ r5 = tmp0
227
 
        mla r2, r3, r2, r1                  @ r2 = tmp3
228
 
        sub r3, r0, r4                      @ r3 = tmp1
229
 
 
230
 
        add r0, r2, r5, lsl #13             @ r0 = tmp10
231
 
        rsb r2, r2, r5, lsl #13             @ r2 = tmp13
232
 
        add r4, r6, r3, lsl #13             @ r4 = tmp11
233
 
        rsb r6, r6, r3, lsl #13             @ r6 = tmp12
234
 
 
235
 
        ldrsh r1, [lr, #( 2*8)]             @ r1 = 'd1'
236
 
        ldrsh r3, [lr, #( 6*8)]             @ r3 = 'd3'
237
 
        ldrsh r5, [lr, #(10*8)]             @ r5 = 'd5'
238
 
        ldrsh r7, [lr, #(14*8)]             @ r7 = 'd7'
239
 
 
240
 
        @ Check for empty odd column (happens about 20 to 25 % of the time according to my stats)
241
 
        orr r9, r1, r3
242
 
        orr r10, r5, r7
243
 
        orrs r10, r9, r10
244
 
        beq empty_odd_column
245
 
 
246
 
        stmdb   sp!, { r0, r2, r4, r6 } @ save on the stack tmp10, tmp13, tmp12, tmp11
247
 
                
248
 
        add r0, r3, r5                  @ r0 = 'z2'
249
 
        add r2, r1, r7                  @ r2 = 'z1'
250
 
        add r4, r3, r7                  @ r4 = 'z3'
251
 
        add r6, r1, r5                  @ r6 = 'z4'
252
 
        ldr r9, [r11, #FIX_1_175875602_ID]
253
 
        add r8, r4, r6
254
 
        ldr r10, [r11, #FIX_M_0_899976223_ID]
255
 
        mul r8, r9, r8                  @ r8 = 'z5'
256
 
        ldr r9, [r11, #FIX_M_2_562915447_ID]
257
 
        mul r2, r10, r2                 @ r2 = 'z1'
258
 
        ldr r10, [r11, #FIX_M_1_961570560_ID]
259
 
        mul r0, r9, r0                  @ r0 = 'z2'
260
 
        ldr r9, [r11, #FIX_M_0_390180644_ID]
261
 
        mla r4, r10, r4, r8             @ r4 = 'z3'
262
 
        ldr r10, [r11, #FIX_0_298631336_ID]
263
 
        mla r6, r9, r6, r8              @ r6 = 'z4'
264
 
        ldr r9, [r11, #FIX_2_053119869_ID]
265
 
        mla r7, r10, r7, r2             @ r7 = tmp0 + z1
266
 
        ldr r10, [r11, #FIX_3_072711026_ID]
267
 
        mla r5, r9, r5, r0              @ r5 = tmp1 + z2
268
 
        ldr r9, [r11, #FIX_1_501321110_ID]
269
 
        mla r3, r10, r3, r0             @ r3 = tmp2 + z2
270
 
        add r7, r7, r4                  @ r7 = tmp0
271
 
        mla r1, r9, r1, r2              @ r1 = tmp3 + z1
272
 
        add r5, r5, r6                  @ r5 = tmp1
273
 
        add r3, r3, r4                  @ r3 = tmp2
274
 
        add r1, r1, r6                  @ r1 = tmp3     
275
 
        
276
 
        ldmia sp!, { r0, r2, r4, r6 } @ r0 = tmp10 / r2 = tmp13 / r4 = tmp11 / r6 = tmp12
277
 
                                      @ r1 = tmp3  / r3 = tmp2  / r5 = tmp1  / r7 = tmp0        
278
 
 
279
 
        @ Compute DESCALE(tmp10 + tmp3, CONST_BITS+PASS1_BITS+3)
280
 
        add r8, r0, r1
281
 
        add r8, r8, #(1<<17)
282
 
        mov r8, r8, asr #18
283
 
        strh r8, [lr, #( 0*8)]
284
 
        
285
 
        @ Compute DESCALE(tmp10 - tmp3, CONST_BITS+PASS1_BITS+3)
286
 
        sub r8, r0, r1
287
 
        add r8, r8, #(1<<17)
288
 
        mov r8, r8, asr #18
289
 
        strh r8, [lr, #(14*8)]
290
 
        
291
 
        @ Compute DESCALE(tmp11 + tmp2, CONST_BITS+PASS1_BITS+3)
292
 
        add r8, r4, r3
293
 
        add r8, r8, #(1<<17)
294
 
        mov r8, r8, asr #18
295
 
        strh r8, [lr, #( 2*8)]
296
 
        
297
 
        @ Compute DESCALE(tmp11 - tmp2, CONST_BITS+PASS1_BITS+3)
298
 
        sub r8, r4, r3
299
 
        add r8, r8, #(1<<17)
300
 
        mov r8, r8, asr #18
301
 
        strh r8, [lr, #(12*8)]
302
 
        
303
 
        @ Compute DESCALE(tmp12 + tmp1, CONST_BITS+PASS1_BITS+3)
304
 
        add r8, r6, r5
305
 
        add r8, r8, #(1<<17)
306
 
        mov r8, r8, asr #18
307
 
        strh r8, [lr, #( 4*8)]
308
 
        
309
 
        @ Compute DESCALE(tmp12 - tmp1, CONST_BITS+PASS1_BITS+3)
310
 
        sub r8, r6, r5
311
 
        add r8, r8, #(1<<17)
312
 
        mov r8, r8, asr #18
313
 
        strh r8, [lr, #(10*8)]
314
 
        
315
 
        @ Compute DESCALE(tmp13 + tmp0, CONST_BITS+PASS1_BITS+3)
316
 
        add r8, r2, r7
317
 
        add r8, r8, #(1<<17)
318
 
        mov r8, r8, asr #18
319
 
        strh r8, [lr, #( 6*8)]
320
 
        
321
 
        @ Compute DESCALE(tmp13 - tmp0, CONST_BITS+PASS1_BITS+3)
322
 
        sub r8, r2, r7
323
 
        add r8, r8, #(1<<17)
324
 
        mov r8, r8, asr #18
325
 
        strh r8, [lr, #( 8*8)]
326
 
 
327
 
        @ End of row loop
328
 
        add lr, lr, #2
329
 
        subs r12, r12, #1
330
 
        bne column_loop
331
 
        beq the_end
332
 
        
 
215
        ldrsh r0, [lr, #( 0*8)]             @ r0 = 'd0'
 
216
        ldrsh r2, [lr, #( 4*8)]             @ r2 = 'd2'
 
217
        ldrsh r4, [lr, #( 8*8)]             @ r4 = 'd4'
 
218
        ldrsh r6, [lr, #(12*8)]             @ r6 = 'd6'
 
219
 
 
220
        ldr r3, [r11, #FIX_0_541196100_ID]
 
221
        add r1, r2, r6
 
222
        ldr r5, [r11, #FIX_M_1_847759065_ID]
 
223
        mul r1, r3, r1                      @ r1 = z1
 
224
        ldr r3, [r11, #FIX_0_765366865_ID]
 
225
        mla r6, r5, r6, r1                  @ r6 = tmp2
 
226
        add r5, r0, r4                      @ r5 = tmp0
 
227
        mla r2, r3, r2, r1                  @ r2 = tmp3
 
228
        sub r3, r0, r4                      @ r3 = tmp1
 
229
 
 
230
        add r0, r2, r5, lsl #13             @ r0 = tmp10
 
231
        rsb r2, r2, r5, lsl #13             @ r2 = tmp13
 
232
        add r4, r6, r3, lsl #13             @ r4 = tmp11
 
233
        rsb r6, r6, r3, lsl #13             @ r6 = tmp12
 
234
 
 
235
        ldrsh r1, [lr, #( 2*8)]             @ r1 = 'd1'
 
236
        ldrsh r3, [lr, #( 6*8)]             @ r3 = 'd3'
 
237
        ldrsh r5, [lr, #(10*8)]             @ r5 = 'd5'
 
238
        ldrsh r7, [lr, #(14*8)]             @ r7 = 'd7'
 
239
 
 
240
        @ Check for empty odd column (happens about 20 to 25 % of the time according to my stats)
 
241
        orr r9, r1, r3
 
242
        orr r10, r5, r7
 
243
        orrs r10, r9, r10
 
244
        beq empty_odd_column
 
245
 
 
246
        stmdb   sp!, { r0, r2, r4, r6 } @ save on the stack tmp10, tmp13, tmp12, tmp11
 
247
 
 
248
        add r0, r3, r5                  @ r0 = 'z2'
 
249
        add r2, r1, r7                  @ r2 = 'z1'
 
250
        add r4, r3, r7                  @ r4 = 'z3'
 
251
        add r6, r1, r5                  @ r6 = 'z4'
 
252
        ldr r9, [r11, #FIX_1_175875602_ID]
 
253
        add r8, r4, r6
 
254
        ldr r10, [r11, #FIX_M_0_899976223_ID]
 
255
        mul r8, r9, r8                  @ r8 = 'z5'
 
256
        ldr r9, [r11, #FIX_M_2_562915447_ID]
 
257
        mul r2, r10, r2                 @ r2 = 'z1'
 
258
        ldr r10, [r11, #FIX_M_1_961570560_ID]
 
259
        mul r0, r9, r0                  @ r0 = 'z2'
 
260
        ldr r9, [r11, #FIX_M_0_390180644_ID]
 
261
        mla r4, r10, r4, r8             @ r4 = 'z3'
 
262
        ldr r10, [r11, #FIX_0_298631336_ID]
 
263
        mla r6, r9, r6, r8              @ r6 = 'z4'
 
264
        ldr r9, [r11, #FIX_2_053119869_ID]
 
265
        mla r7, r10, r7, r2             @ r7 = tmp0 + z1
 
266
        ldr r10, [r11, #FIX_3_072711026_ID]
 
267
        mla r5, r9, r5, r0              @ r5 = tmp1 + z2
 
268
        ldr r9, [r11, #FIX_1_501321110_ID]
 
269
        mla r3, r10, r3, r0             @ r3 = tmp2 + z2
 
270
        add r7, r7, r4                  @ r7 = tmp0
 
271
        mla r1, r9, r1, r2              @ r1 = tmp3 + z1
 
272
        add r5,        r5, r6                  @ r5 = tmp1
 
273
        add r3, r3, r4                  @ r3 = tmp2
 
274
        add r1, r1, r6                  @ r1 = tmp3
 
275
 
 
276
        ldmia sp!, { r0, r2, r4, r6 } @ r0 = tmp10 / r2 = tmp13 / r4 = tmp11 / r6 = tmp12
 
277
                                      @ r1 = tmp3  / r3 = tmp2  / r5 = tmp1  / r7 = tmp0
 
278
 
 
279
        @ Compute DESCALE(tmp10 + tmp3, CONST_BITS+PASS1_BITS+3)
 
280
        add r8, r0, r1
 
281
        add r8, r8, #(1<<17)
 
282
        mov r8, r8, asr #18
 
283
        strh r8, [lr, #( 0*8)]
 
284
 
 
285
        @ Compute DESCALE(tmp10 - tmp3, CONST_BITS+PASS1_BITS+3)
 
286
        sub r8, r0, r1
 
287
        add r8, r8, #(1<<17)
 
288
        mov r8, r8, asr #18
 
289
        strh r8, [lr, #(14*8)]
 
290
 
 
291
        @ Compute DESCALE(tmp11 + tmp2, CONST_BITS+PASS1_BITS+3)
 
292
        add r8, r4, r3
 
293
        add r8, r8, #(1<<17)
 
294
        mov r8, r8, asr #18
 
295
        strh r8, [lr, #( 2*8)]
 
296
 
 
297
        @ Compute DESCALE(tmp11 - tmp2, CONST_BITS+PASS1_BITS+3)
 
298
        sub r8, r4, r3
 
299
        add r8, r8, #(1<<17)
 
300
        mov r8, r8, asr #18
 
301
        strh r8, [lr, #(12*8)]
 
302
 
 
303
        @ Compute DESCALE(tmp12 + tmp1, CONST_BITS+PASS1_BITS+3)
 
304
        add r8, r6, r5
 
305
        add r8, r8, #(1<<17)
 
306
        mov r8, r8, asr #18
 
307
        strh r8, [lr, #( 4*8)]
 
308
 
 
309
        @ Compute DESCALE(tmp12 - tmp1, CONST_BITS+PASS1_BITS+3)
 
310
        sub r8, r6, r5
 
311
        add r8, r8, #(1<<17)
 
312
        mov r8, r8, asr #18
 
313
        strh r8, [lr, #(10*8)]
 
314
 
 
315
        @ Compute DESCALE(tmp13 + tmp0, CONST_BITS+PASS1_BITS+3)
 
316
        add r8, r2, r7
 
317
        add r8, r8, #(1<<17)
 
318
        mov r8, r8, asr #18
 
319
        strh r8, [lr, #( 6*8)]
 
320
 
 
321
        @ Compute DESCALE(tmp13 - tmp0, CONST_BITS+PASS1_BITS+3)
 
322
        sub r8, r2, r7
 
323
        add r8, r8, #(1<<17)
 
324
        mov r8, r8, asr #18
 
325
        strh r8, [lr, #( 8*8)]
 
326
 
 
327
        @ End of row loop
 
328
        add lr, lr, #2
 
329
        subs r12, r12, #1
 
330
        bne column_loop
 
331
        beq the_end
 
332
 
333
333
empty_odd_column:
334
 
        @ Compute DESCALE(tmp10 + tmp3, CONST_BITS+PASS1_BITS+3)
335
 
        @ Compute DESCALE(tmp10 - tmp3, CONST_BITS+PASS1_BITS+3)
336
 
        add r0, r0, #(1<<17)
337
 
        mov r0, r0, asr #18
338
 
        strh r0, [lr, #( 0*8)]
339
 
        strh r0, [lr, #(14*8)]
340
 
        
341
 
        @ Compute DESCALE(tmp11 + tmp2, CONST_BITS+PASS1_BITS+3)
342
 
        @ Compute DESCALE(tmp11 - tmp2, CONST_BITS+PASS1_BITS+3)
343
 
        add r4, r4, #(1<<17)
344
 
        mov r4, r4, asr #18
345
 
        strh r4, [lr, #( 2*8)]
346
 
        strh r4, [lr, #(12*8)]
347
 
        
348
 
        @ Compute DESCALE(tmp12 + tmp1, CONST_BITS+PASS1_BITS+3)
349
 
        @ Compute DESCALE(tmp12 - tmp1, CONST_BITS+PASS1_BITS+3)
350
 
        add r6, r6, #(1<<17)
351
 
        mov r6, r6, asr #18
352
 
        strh r6, [lr, #( 4*8)]
353
 
        strh r6, [lr, #(10*8)]
354
 
        
355
 
        @ Compute DESCALE(tmp13 + tmp0, CONST_BITS+PASS1_BITS+3)
356
 
        @ Compute DESCALE(tmp13 - tmp0, CONST_BITS+PASS1_BITS+3)
357
 
        add r2, r2, #(1<<17)
358
 
        mov r2, r2, asr #18
359
 
        strh r2, [lr, #( 6*8)]
360
 
        strh r2, [lr, #( 8*8)]
361
 
 
362
 
        @ End of row loop
363
 
        add lr, lr, #2
364
 
        subs r12, r12, #1
365
 
        bne column_loop
366
 
                
367
 
the_end:        
368
 
        @ The end....
369
 
        add sp, sp, #4
370
 
        ldmia   sp!, { r4 - r12, pc }   @ restore callee saved regs and return
 
334
        @ Compute DESCALE(tmp10 + tmp3, CONST_BITS+PASS1_BITS+3)
 
335
        @ Compute DESCALE(tmp10 - tmp3, CONST_BITS+PASS1_BITS+3)
 
336
        add r0, r0, #(1<<17)
 
337
        mov r0, r0, asr #18
 
338
        strh r0, [lr, #( 0*8)]
 
339
        strh r0, [lr, #(14*8)]
 
340
 
 
341
        @ Compute DESCALE(tmp11 + tmp2, CONST_BITS+PASS1_BITS+3)
 
342
        @ Compute DESCALE(tmp11 - tmp2, CONST_BITS+PASS1_BITS+3)
 
343
        add r4, r4, #(1<<17)
 
344
        mov r4, r4, asr #18
 
345
        strh r4, [lr, #( 2*8)]
 
346
        strh r4, [lr, #(12*8)]
 
347
 
 
348
        @ Compute DESCALE(tmp12 + tmp1, CONST_BITS+PASS1_BITS+3)
 
349
        @ Compute DESCALE(tmp12 - tmp1, CONST_BITS+PASS1_BITS+3)
 
350
        add r6, r6, #(1<<17)
 
351
        mov r6, r6, asr #18
 
352
        strh r6, [lr, #( 4*8)]
 
353
        strh r6, [lr, #(10*8)]
 
354
 
 
355
        @ Compute DESCALE(tmp13 + tmp0, CONST_BITS+PASS1_BITS+3)
 
356
        @ Compute DESCALE(tmp13 - tmp0, CONST_BITS+PASS1_BITS+3)
 
357
        add r2, r2, #(1<<17)
 
358
        mov r2, r2, asr #18
 
359
        strh r2, [lr, #( 6*8)]
 
360
        strh r2, [lr, #( 8*8)]
 
361
 
 
362
        @ End of row loop
 
363
        add lr, lr, #2
 
364
        subs r12, r12, #1
 
365
        bne column_loop
 
366
 
 
367
the_end:
 
368
        @ The end....
 
369
        add sp, sp, #4
 
370
        ldmia   sp!, { r4 - r12, pc }   @ restore callee saved regs and return
371
371
 
372
372
const_array:
373
 
        .align
374
 
        .word FIX_0_298631336
375
 
        .word FIX_0_541196100
376
 
        .word FIX_0_765366865
377
 
        .word FIX_1_175875602
378
 
        .word FIX_1_501321110
379
 
        .word FIX_2_053119869
380
 
        .word FIX_3_072711026
381
 
        .word FIX_M_0_390180644
382
 
        .word FIX_M_0_899976223
383
 
        .word FIX_M_1_847759065
384
 
        .word FIX_M_1_961570560
385
 
        .word FIX_M_2_562915447
386
 
        .word FIX_0xFFFF
 
373
        .align
 
374
        .word FIX_0_298631336
 
375
        .word FIX_0_541196100
 
376
        .word FIX_0_765366865
 
377
        .word FIX_1_175875602
 
378
        .word FIX_1_501321110
 
379
        .word FIX_2_053119869
 
380
        .word FIX_3_072711026
 
381
        .word FIX_M_0_390180644
 
382
        .word FIX_M_0_899976223
 
383
        .word FIX_M_1_847759065
 
384
        .word FIX_M_1_961570560
 
385
        .word FIX_M_2_562915447
 
386
        .word FIX_0xFFFF