summaryrefslogtreecommitdiff
path: root/simd/powerpc/jidctfst-altivec.c
diff options
context:
space:
mode:
Diffstat (limited to 'simd/powerpc/jidctfst-altivec.c')
-rw-r--r--simd/powerpc/jidctfst-altivec.c255
1 files changed, 0 insertions, 255 deletions
diff --git a/simd/powerpc/jidctfst-altivec.c b/simd/powerpc/jidctfst-altivec.c
deleted file mode 100644
index 456c6c6..0000000
--- a/simd/powerpc/jidctfst-altivec.c
+++ /dev/null
@@ -1,255 +0,0 @@
-/*
- * AltiVec optimizations for libjpeg-turbo
- *
- * Copyright (C) 2014-2015, D. R. Commander. All Rights Reserved.
- *
- * This software is provided 'as-is', without any express or implied
- * warranty. In no event will the authors be held liable for any damages
- * arising from the use of this software.
- *
- * Permission is granted to anyone to use this software for any purpose,
- * including commercial applications, and to alter it and redistribute it
- * freely, subject to the following restrictions:
- *
- * 1. The origin of this software must not be misrepresented; you must not
- * claim that you wrote the original software. If you use this software
- * in a product, an acknowledgment in the product documentation would be
- * appreciated but is not required.
- * 2. Altered source versions must be plainly marked as such, and must not be
- * misrepresented as being the original software.
- * 3. This notice may not be removed or altered from any source distribution.
- */
-
-/* FAST INTEGER INVERSE DCT
- *
- * This is similar to the SSE2 implementation, except that we left-shift the
- * constants by 1 less bit (the -1 in CONST_SHIFT.) This is because
- * vec_madds(arg1, arg2, arg3) generates the 16-bit saturated sum of:
- * the elements in arg3 + the most significant 17 bits of
- * (the elements in arg1 * the elements in arg2).
- */
-
-#include "jsimd_altivec.h"
-
-
-#define F_1_082 277 /* FIX(1.082392200) */
-#define F_1_414 362 /* FIX(1.414213562) */
-#define F_1_847 473 /* FIX(1.847759065) */
-#define F_2_613 669 /* FIX(2.613125930) */
-#define F_1_613 (F_2_613 - 256) /* FIX(2.613125930) - FIX(1) */
-
-#define CONST_BITS 8
-#define PASS1_BITS 2
-#define PRE_MULTIPLY_SCALE_BITS 2
-#define CONST_SHIFT (16 - PRE_MULTIPLY_SCALE_BITS - CONST_BITS - 1)
-
-
-#define DO_IDCT(in) { \
- /* Even part */ \
- \
- tmp10 = vec_add(in##0, in##4); \
- tmp11 = vec_sub(in##0, in##4); \
- tmp13 = vec_add(in##2, in##6); \
- \
- tmp12 = vec_sub(in##2, in##6); \
- tmp12 = vec_sl(tmp12, pre_multiply_scale_bits); \
- tmp12 = vec_madds(tmp12, pw_F1414, pw_zero); \
- tmp12 = vec_sub(tmp12, tmp13); \
- \
- tmp0 = vec_add(tmp10, tmp13); \
- tmp3 = vec_sub(tmp10, tmp13); \
- tmp1 = vec_add(tmp11, tmp12); \
- tmp2 = vec_sub(tmp11, tmp12); \
- \
- /* Odd part */ \
- \
- z13 = vec_add(in##5, in##3); \
- z10 = vec_sub(in##5, in##3); \
- z10s = vec_sl(z10, pre_multiply_scale_bits); \
- z11 = vec_add(in##1, in##7); \
- z12s = vec_sub(in##1, in##7); \
- z12s = vec_sl(z12s, pre_multiply_scale_bits); \
- \
- tmp11 = vec_sub(z11, z13); \
- tmp11 = vec_sl(tmp11, pre_multiply_scale_bits); \
- tmp11 = vec_madds(tmp11, pw_F1414, pw_zero); \
- \
- tmp7 = vec_add(z11, z13); \
- \
- /* To avoid overflow... \
- * \
- * (Original) \
- * tmp12 = -2.613125930 * z10 + z5; \
- * \
- * (This implementation) \
- * tmp12 = (-1.613125930 - 1) * z10 + z5; \
- * = -1.613125930 * z10 - z10 + z5; \
- */ \
- \
- z5 = vec_add(z10s, z12s); \
- z5 = vec_madds(z5, pw_F1847, pw_zero); \
- \
- tmp10 = vec_madds(z12s, pw_F1082, pw_zero); \
- tmp10 = vec_sub(tmp10, z5); \
- tmp12 = vec_madds(z10s, pw_MF1613, z5); \
- tmp12 = vec_sub(tmp12, z10); \
- \
- tmp6 = vec_sub(tmp12, tmp7); \
- tmp5 = vec_sub(tmp11, tmp6); \
- tmp4 = vec_add(tmp10, tmp5); \
- \
- out0 = vec_add(tmp0, tmp7); \
- out1 = vec_add(tmp1, tmp6); \
- out2 = vec_add(tmp2, tmp5); \
- out3 = vec_sub(tmp3, tmp4); \
- out4 = vec_add(tmp3, tmp4); \
- out5 = vec_sub(tmp2, tmp5); \
- out6 = vec_sub(tmp1, tmp6); \
- out7 = vec_sub(tmp0, tmp7); \
-}
-
-
-void jsimd_idct_ifast_altivec(void *dct_table_, JCOEFPTR coef_block,
- JSAMPARRAY output_buf, JDIMENSION output_col)
-{
- short *dct_table = (short *)dct_table_;
- int *outptr;
-
- __vector short row0, row1, row2, row3, row4, row5, row6, row7,
- col0, col1, col2, col3, col4, col5, col6, col7,
- quant0, quant1, quant2, quant3, quant4, quant5, quant6, quant7,
- tmp0, tmp1, tmp2, tmp3, tmp4, tmp5, tmp6, tmp7, tmp10, tmp11, tmp12, tmp13,
- z5, z10, z10s, z11, z12s, z13,
- out0, out1, out2, out3, out4, out5, out6, out7;
- __vector signed char outb;
-
- /* Constants */
- __vector short pw_zero = { __8X(0) },
- pw_F1414 = { __8X(F_1_414 << CONST_SHIFT) },
- pw_F1847 = { __8X(F_1_847 << CONST_SHIFT) },
- pw_MF1613 = { __8X(-F_1_613 << CONST_SHIFT) },
- pw_F1082 = { __8X(F_1_082 << CONST_SHIFT) };
- __vector unsigned short
- pre_multiply_scale_bits = { __8X(PRE_MULTIPLY_SCALE_BITS) },
- pass1_bits3 = { __8X(PASS1_BITS + 3) };
- __vector signed char pb_centerjsamp = { __16X(CENTERJSAMPLE) };
-
- /* Pass 1: process columns */
-
- col0 = vec_ld(0, coef_block);
- col1 = vec_ld(16, coef_block);
- col2 = vec_ld(32, coef_block);
- col3 = vec_ld(48, coef_block);
- col4 = vec_ld(64, coef_block);
- col5 = vec_ld(80, coef_block);
- col6 = vec_ld(96, coef_block);
- col7 = vec_ld(112, coef_block);
-
- tmp1 = vec_or(col1, col2);
- tmp2 = vec_or(col3, col4);
- tmp1 = vec_or(tmp1, tmp2);
- tmp3 = vec_or(col5, col6);
- tmp3 = vec_or(tmp3, col7);
- tmp1 = vec_or(tmp1, tmp3);
-
- quant0 = vec_ld(0, dct_table);
- col0 = vec_mladd(col0, quant0, pw_zero);
-
- if (vec_all_eq(tmp1, pw_zero)) {
- /* AC terms all zero */
-
- row0 = vec_splat(col0, 0);
- row1 = vec_splat(col0, 1);
- row2 = vec_splat(col0, 2);
- row3 = vec_splat(col0, 3);
- row4 = vec_splat(col0, 4);
- row5 = vec_splat(col0, 5);
- row6 = vec_splat(col0, 6);
- row7 = vec_splat(col0, 7);
-
- } else {
-
- quant1 = vec_ld(16, dct_table);
- quant2 = vec_ld(32, dct_table);
- quant3 = vec_ld(48, dct_table);
- quant4 = vec_ld(64, dct_table);
- quant5 = vec_ld(80, dct_table);
- quant6 = vec_ld(96, dct_table);
- quant7 = vec_ld(112, dct_table);
-
- col1 = vec_mladd(col1, quant1, pw_zero);
- col2 = vec_mladd(col2, quant2, pw_zero);
- col3 = vec_mladd(col3, quant3, pw_zero);
- col4 = vec_mladd(col4, quant4, pw_zero);
- col5 = vec_mladd(col5, quant5, pw_zero);
- col6 = vec_mladd(col6, quant6, pw_zero);
- col7 = vec_mladd(col7, quant7, pw_zero);
-
- DO_IDCT(col);
-
- TRANSPOSE(out, row);
- }
-
- /* Pass 2: process rows */
-
- DO_IDCT(row);
-
- out0 = vec_sra(out0, pass1_bits3);
- out1 = vec_sra(out1, pass1_bits3);
- out2 = vec_sra(out2, pass1_bits3);
- out3 = vec_sra(out3, pass1_bits3);
- out4 = vec_sra(out4, pass1_bits3);
- out5 = vec_sra(out5, pass1_bits3);
- out6 = vec_sra(out6, pass1_bits3);
- out7 = vec_sra(out7, pass1_bits3);
-
- TRANSPOSE(out, col);
-
- outb = vec_packs(col0, col0);
- outb = vec_add(outb, pb_centerjsamp);
- outptr = (int *)(output_buf[0] + output_col);
- vec_ste((__vector int)outb, 0, outptr);
- vec_ste((__vector int)outb, 4, outptr);
-
- outb = vec_packs(col1, col1);
- outb = vec_add(outb, pb_centerjsamp);
- outptr = (int *)(output_buf[1] + output_col);
- vec_ste((__vector int)outb, 0, outptr);
- vec_ste((__vector int)outb, 4, outptr);
-
- outb = vec_packs(col2, col2);
- outb = vec_add(outb, pb_centerjsamp);
- outptr = (int *)(output_buf[2] + output_col);
- vec_ste((__vector int)outb, 0, outptr);
- vec_ste((__vector int)outb, 4, outptr);
-
- outb = vec_packs(col3, col3);
- outb = vec_add(outb, pb_centerjsamp);
- outptr = (int *)(output_buf[3] + output_col);
- vec_ste((__vector int)outb, 0, outptr);
- vec_ste((__vector int)outb, 4, outptr);
-
- outb = vec_packs(col4, col4);
- outb = vec_add(outb, pb_centerjsamp);
- outptr = (int *)(output_buf[4] + output_col);
- vec_ste((__vector int)outb, 0, outptr);
- vec_ste((__vector int)outb, 4, outptr);
-
- outb = vec_packs(col5, col5);
- outb = vec_add(outb, pb_centerjsamp);
- outptr = (int *)(output_buf[5] + output_col);
- vec_ste((__vector int)outb, 0, outptr);
- vec_ste((__vector int)outb, 4, outptr);
-
- outb = vec_packs(col6, col6);
- outb = vec_add(outb, pb_centerjsamp);
- outptr = (int *)(output_buf[6] + output_col);
- vec_ste((__vector int)outb, 0, outptr);
- vec_ste((__vector int)outb, 4, outptr);
-
- outb = vec_packs(col7, col7);
- outb = vec_add(outb, pb_centerjsamp);
- outptr = (int *)(output_buf[7] + output_col);
- vec_ste((__vector int)outb, 0, outptr);
- vec_ste((__vector int)outb, 4, outptr);
-}
generated by cgit v1.2.3 (git 2.25.1) at 2025-10-03 08:01:15 +0000