|
|
|
|
/*
|
|
|
|
|
* Copyright (c) 2016 Clément Bœsch <clement stupeflix.com>
|
|
|
|
|
*
|
|
|
|
|
* This file is part of FFmpeg.
|
|
|
|
|
*
|
|
|
|
|
* FFmpeg is free software; you can redistribute it and/or
|
|
|
|
|
* modify it under the terms of the GNU Lesser General Public
|
|
|
|
|
* License as published by the Free Software Foundation; either
|
|
|
|
|
* version 2.1 of the License, or (at your option) any later version.
|
|
|
|
|
*
|
|
|
|
|
* FFmpeg is distributed in the hope that it will be useful,
|
|
|
|
|
* but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
|
|
|
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
|
|
|
|
|
* Lesser General Public License for more details.
|
|
|
|
|
*
|
|
|
|
|
* You should have received a copy of the GNU Lesser General Public
|
|
|
|
|
* License along with FFmpeg; if not, write to the Free Software
|
|
|
|
|
* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
|
|
|
|
|
*/
|
|
|
|
|
|
|
|
|
|
#include "libavutil/aarch64/asm.S"
|
|
|
|
|
|
|
|
|
|
function ff_hscale_8_to_15_neon, export=1
|
|
|
|
|
add x10, x4, w6, UXTW #1 // filter2 = filter + filterSize*2 (x2 because int16)
|
|
|
|
|
1: ldr w8, [x5], #4 // filterPos[0]
|
|
|
|
|
ldr w9, [x5], #4 // filterPos[1]
|
|
|
|
|
movi v4.4S, #0 // val sum part 1 (for dst[0])
|
|
|
|
|
movi v5.4S, #0 // val sum part 2 (for dst[1])
|
|
|
|
|
mov w7, w6 // filterSize counter
|
|
|
|
|
mov x13, x3 // srcp = src
|
|
|
|
|
2: add x11, x13, w8, UXTW // srcp + filterPos[0]
|
|
|
|
|
add x12, x13, w9, UXTW // srcp + filterPos[1]
|
|
|
|
|
ld1 {v0.8B}, [x11] // srcp[filterPos[0] + {0..7}]
|
|
|
|
|
ld1 {v1.8B}, [x12] // srcp[filterPos[1] + {0..7}]
|
|
|
|
|
ld1 {v2.8H}, [x4], #16 // load 8x16-bit filter values, part 1
|
|
|
|
|
ld1 {v3.8H}, [x10], #16 // ditto at filter+filterSize for part 2
|
|
|
|
|
uxtl v0.8H, v0.8B // unpack part 1 to 16-bit
|
|
|
|
|
uxtl v1.8H, v1.8B // unpack part 2 to 16-bit
|
|
|
|
|
smull v16.4S, v0.4H, v2.4H // v16.i32{0..3} = part 1 of: srcp[filterPos[0] + {0..7}] * filter[{0..7}]
|
|
|
|
|
smull v18.4S, v1.4H, v3.4H // v18.i32{0..3} = part 1 of: srcp[filterPos[1] + {0..7}] * filter[{0..7}]
|
|
|
|
|
smull2 v17.4S, v0.8H, v2.8H // v17.i32{0..3} = part 2 of: srcp[filterPos[0] + {0..7}] * filter[{0..7}]
|
|
|
|
|
smull2 v19.4S, v1.8H, v3.8H // v19.i32{0..3} = part 2 of: srcp[filterPos[1] + {0..7}] * filter[{0..7}]
|
|
|
|
|
addp v16.4S, v16.4S, v17.4S // horizontal pair adding of the 8x32-bit multiplied values for part 1 into 4x32-bit
|
|
|
|
|
addp v18.4S, v18.4S, v19.4S // horizontal pair adding of the 8x32-bit multiplied values for part 2 into 4x32-bit
|
|
|
|
|
add v4.4S, v4.4S, v16.4S // update val accumulator for part 1
|
|
|
|
|
add v5.4S, v5.4S, v18.4S // update val accumulator for part 2
|
|
|
|
|
add x13, x13, #8 // srcp += 8
|
|
|
|
|
subs w7, w7, #8 // processed 8/filterSize
|
|
|
|
|
b.gt 2b // inner loop if filterSize not consumed completely
|
|
|
|
|
mov x4, x10 // filter = filter2
|
|
|
|
|
add x10, x10, w6, UXTW #1 // filter2 += filterSize*2
|
|
|
|
|
addp v4.4S, v4.4S, v5.4S // horizontal pair adding of the 8x32-bit sums into 4x32-bit
|
|
|
|
|
addp v4.4S, v4.4S, v4.4S // horizontal pair adding of the 4x32-bit sums into 2x32-bit
|
|
|
|
|
sqshrn v4.4H, v4.4S, #7 // shift and clip the 2x16-bit final values
|
|
|
|
|
st1 {v4.S}[0], [x1], #4 // write to destination
|
|
|
|
|
subs w2, w2, #2 // dstW -= 2
|
|
|
|
|
b.gt 1b // loop until end of line
|
|
|
|
|
ret
|
|
|
|
|
endfunc
|
