2012-10-31 22:40:12 +03:00
|
|
|
;******************************************************************************
|
|
|
|
|
;* x86 optimized dithering format conversion
|
|
|
|
|
;* Copyright (c) 2012 Justin Ruggles <justin.ruggles@gmail.com>
|
|
|
|
|
;*
|
|
|
|
|
;* This file is part of Libav.
|
|
|
|
|
;*
|
|
|
|
|
;* Libav 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.
|
|
|
|
|
;*
|
|
|
|
|
;* Libav 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 Libav; if not, write to the Free Software
|
|
|
|
|
;* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
|
|
|
|
|
;******************************************************************************
|
|
|
|
|
|
|
|
|
|
%include "libavutil/x86/x86util.asm"
|
|
|
|
|
|
|
|
|
|
SECTION_RODATA 32
|
|
|
|
|
|
2012-10-31 22:40:12 +03:00
|
|
|
; 1.0f / (2.0f * INT32_MAX)
|
|
|
|
|
pf_dither_scale: times 8 dd 2.32830643762e-10
|
|
|
|
|
|
2012-10-31 22:40:12 +03:00
|
|
|
pf_s16_scale: times 4 dd 32753.0
|
|
|
|
|
|
|
|
|
|
SECTION_TEXT
|
|
|
|
|
|
|
|
|
|
;------------------------------------------------------------------------------
|
|
|
|
|
; void ff_quantize(int16_t *dst, float *src, float *dither, int len);
|
|
|
|
|
;------------------------------------------------------------------------------
|
|
|
|
|
|
|
|
|
|
INIT_XMM sse2
|
|
|
|
|
cglobal quantize, 4,4,3, dst, src, dither, len
|
|
|
|
|
lea lenq, [2*lend]
|
|
|
|
|
add dstq, lenq
|
|
|
|
|
lea srcq, [srcq+2*lenq]
|
|
|
|
|
lea ditherq, [ditherq+2*lenq]
|
|
|
|
|
neg lenq
|
|
|
|
|
mova m2, [pf_s16_scale]
|
|
|
|
|
.loop:
|
|
|
|
|
mulps m0, m2, [srcq+2*lenq]
|
|
|
|
|
mulps m1, m2, [srcq+2*lenq+mmsize]
|
|
|
|
|
addps m0, [ditherq+2*lenq]
|
|
|
|
|
addps m1, [ditherq+2*lenq+mmsize]
|
|
|
|
|
cvtps2dq m0, m0
|
|
|
|
|
cvtps2dq m1, m1
|
|
|
|
|
packssdw m0, m1
|
|
|
|
|
mova [dstq+lenq], m0
|
|
|
|
|
add lenq, mmsize
|
|
|
|
|
jl .loop
|
|
|
|
|
REP_RET
|
2012-10-31 22:40:12 +03:00
|
|
|
|
|
|
|
|
;------------------------------------------------------------------------------
|
|
|
|
|
; void ff_dither_int_to_float_rectangular(float *dst, int *src, int len)
|
|
|
|
|
;------------------------------------------------------------------------------
|
|
|
|
|
|
|
|
|
|
%macro DITHER_INT_TO_FLOAT_RECTANGULAR 0
|
|
|
|
|
cglobal dither_int_to_float_rectangular, 3,3,3, dst, src, len
|
|
|
|
|
lea lenq, [4*lend]
|
|
|
|
|
add srcq, lenq
|
|
|
|
|
add dstq, lenq
|
|
|
|
|
neg lenq
|
|
|
|
|
mova m0, [pf_dither_scale]
|
|
|
|
|
.loop:
|
|
|
|
|
cvtdq2ps m1, [srcq+lenq]
|
|
|
|
|
cvtdq2ps m2, [srcq+lenq+mmsize]
|
|
|
|
|
mulps m1, m1, m0
|
|
|
|
|
mulps m2, m2, m0
|
|
|
|
|
mova [dstq+lenq], m1
|
|
|
|
|
mova [dstq+lenq+mmsize], m2
|
|
|
|
|
add lenq, 2*mmsize
|
|
|
|
|
jl .loop
|
|
|
|
|
REP_RET
|
|
|
|
|
%endmacro
|
|
|
|
|
|
|
|
|
|
INIT_XMM sse2
|
|
|
|
|
DITHER_INT_TO_FLOAT_RECTANGULAR
|
|
|
|
|
INIT_YMM avx
|
|
|
|
|
DITHER_INT_TO_FLOAT_RECTANGULAR
|
|
|
|
|
|
|
|
|
|
;------------------------------------------------------------------------------
|
|
|
|
|
; void ff_dither_int_to_float_triangular(float *dst, int *src0, int len)
|
|
|
|
|
;------------------------------------------------------------------------------
|
|
|
|
|
|
|
|
|
|
%macro DITHER_INT_TO_FLOAT_TRIANGULAR 0
|
|
|
|
|
cglobal dither_int_to_float_triangular, 3,4,5, dst, src0, len, src1
|
|
|
|
|
lea lenq, [4*lend]
|
|
|
|
|
lea src1q, [src0q+2*lenq]
|
|
|
|
|
add src0q, lenq
|
|
|
|
|
add dstq, lenq
|
|
|
|
|
neg lenq
|
|
|
|
|
mova m0, [pf_dither_scale]
|
|
|
|
|
.loop:
|
|
|
|
|
cvtdq2ps m1, [src0q+lenq]
|
|
|
|
|
cvtdq2ps m2, [src0q+lenq+mmsize]
|
|
|
|
|
cvtdq2ps m3, [src1q+lenq]
|
|
|
|
|
cvtdq2ps m4, [src1q+lenq+mmsize]
|
|
|
|
|
addps m1, m1, m3
|
|
|
|
|
addps m2, m2, m4
|
|
|
|
|
mulps m1, m1, m0
|
|
|
|
|
mulps m2, m2, m0
|
|
|
|
|
mova [dstq+lenq], m1
|
|
|
|
|
mova [dstq+lenq+mmsize], m2
|
|
|
|
|
add lenq, 2*mmsize
|
|
|
|
|
jl .loop
|
|
|
|
|
REP_RET
|
|
|
|
|
%endmacro
|
|
|
|
|
|
|
|
|
|
INIT_XMM sse2
|
|
|
|
|
DITHER_INT_TO_FLOAT_TRIANGULAR
|
|
|
|
|
INIT_YMM avx
|
|
|
|
|
DITHER_INT_TO_FLOAT_TRIANGULAR
|
