/* * jsimd_powerpc.c * * Copyright 2009 Pierre Ossman for Cendio AB * Copyright (C) 2009-2011, 2014-2016, 2018, D. R. Commander. * Copyright (C) 2015-2016, 2018, Matthieu Darbois. * * Based on the x86 SIMD extension for IJG JPEG library, * Copyright (C) 1999-2006, MIYASAKA Masaru. * For conditions of distribution and use, see copyright notice in jsimdext.inc * * This file contains the interface between the "normal" portions * of the library and the SIMD implementations when running on a * PowerPC architecture. */ #ifdef __amigaos4__ /* This must be defined first as it re-defines GLOBAL otherwise */ #include #endif #define JPEG_INTERNALS #include "../../jinclude.h" #include "../../jpeglib.h" #include "../../jsimd.h" #include "../../jdct.h" #include "../../jsimddct.h" #include "../jsimd.h" #include #if defined(__APPLE__) #include #include #elif defined(__OpenBSD__) #include #include #include #elif defined(__FreeBSD__) #include #include #endif static unsigned int simd_support = ~0; #if !defined(__ALTIVEC__) && (defined(__linux__) || defined(ANDROID) || defined(__ANDROID__)) #define SOMEWHAT_SANE_PROC_CPUINFO_SIZE_LIMIT (1024 * 1024) LOCAL(int) check_feature(char *buffer, char *feature) { char *p; if (*feature == 0) return 0; if (strncmp(buffer, "cpu", 3) != 0) return 0; buffer += 3; while (isspace(*buffer)) buffer++; /* Check if 'feature' is present in the buffer as a separate word */ while ((p = strstr(buffer, feature))) { if (p > buffer && !isspace(*(p - 1))) { buffer++; continue; } p += strlen(feature); if (*p != 0 && !isspace(*p)) { buffer++; continue; } return 1; } return 0; } LOCAL(int) parse_proc_cpuinfo(int bufsize) { char *buffer = (char *)malloc(bufsize); FILE *fd; simd_support = 0; if (!buffer) return 0; fd = fopen("/proc/cpuinfo", "r"); if (fd) { while (fgets(buffer, bufsize, fd)) { if (!strchr(buffer, '\n') && !feof(fd)) { /* "impossible" happened - insufficient size of the buffer! */ fclose(fd); free(buffer); return 0; } if (check_feature(buffer, "altivec")) simd_support |= JSIMD_ALTIVEC; } fclose(fd); } free(buffer); return 1; } #endif /* * Check what SIMD accelerations are supported. * * FIXME: This code is racy under a multi-threaded environment. */ LOCAL(void) init_simd(void) { #ifndef NO_GETENV char *env = NULL; #endif #if !defined(__ALTIVEC__) && (defined(__linux__) || defined(ANDROID) || defined(__ANDROID__)) int bufsize = 1024; /* an initial guess for the line buffer size limit */ #elif defined(__amigaos4__) uint32 altivec = 0; #elif defined(__APPLE__) int mib[2] = { CTL_HW, HW_VECTORUNIT }; int altivec; size_t len = sizeof(altivec); #elif defined(__OpenBSD__) int mib[2] = { CTL_MACHDEP, CPU_ALTIVEC }; int altivec; size_t len = sizeof(altivec); #elif defined(__FreeBSD__) unsigned long cpufeatures = 0; #endif if (simd_support != ~0U) return; simd_support = 0; #if defined(__ALTIVEC__) simd_support |= JSIMD_ALTIVEC; #elif defined(__linux__) || defined(ANDROID) || defined(__ANDROID__) while (!parse_proc_cpuinfo(bufsize)) { bufsize *= 2; if (bufsize > SOMEWHAT_SANE_PROC_CPUINFO_SIZE_LIMIT) break; } #elif defined(__amigaos4__) IExec->GetCPUInfoTags(GCIT_VectorUnit, &altivec, TAG_DONE); if (altivec == VECTORTYPE_ALTIVEC) simd_support |= JSIMD_ALTIVEC; #elif defined(__APPLE__) || defined(__OpenBSD__) if (sysctl(mib, 2, &altivec, &len, NULL, 0) == 0 && altivec != 0) simd_support |= JSIMD_ALTIVEC; #elif defined(__FreeBSD__) elf_aux_info(AT_HWCAP, &cpufeatures, sizeof(cpufeatures)); if (cpufeatures & PPC_FEATURE_HAS_ALTIVEC) simd_support |= JSIMD_ALTIVEC; #endif #ifndef NO_GETENV /* Force different settings through environment variables */ env = getenv("JSIMD_FORCEALTIVEC"); if ((env != NULL) && (strcmp(env, "1") == 0)) simd_support = JSIMD_ALTIVEC; env = getenv("JSIMD_FORCENONE"); if ((env != NULL) && (strcmp(env, "1") == 0)) simd_support = 0; #endif } GLOBAL(int) jsimd_can_rgb_ycc(void) { init_simd(); /* The code is optimised for these values only */ if (BITS_IN_JSAMPLE != 8) return 0; if (sizeof(JDIMENSION) != 4) return 0; if ((RGB_PIXELSIZE != 3) && (RGB_PIXELSIZE != 4)) return 0; if (simd_support & JSIMD_ALTIVEC) return 1; return 0; } GLOBAL(int) jsimd_can_rgb_gray(void) { init_simd(); /* The code is optimised for these values only */ if (BITS_IN_JSAMPLE != 8) return 0; if (sizeof(JDIMENSION) != 4) return 0; if ((RGB_PIXELSIZE != 3) && (RGB_PIXELSIZE != 4)) return 0; if (simd_support & JSIMD_ALTIVEC) return 1; return 0; } GLOBAL(int) jsimd_can_ycc_rgb(void) { init_simd(); /* The code is optimised for these values only */ if (BITS_IN_JSAMPLE != 8) return 0; if (sizeof(JDIMENSION) != 4) return 0; if ((RGB_PIXELSIZE != 3) && (RGB_PIXELSIZE != 4)) return 0; if (simd_support & JSIMD_ALTIVEC) return 1; return 0; } GLOBAL(int) jsimd_can_ycc_rgb565(void) { return 0; } GLOBAL(void) jsimd_rgb_ycc_convert(j_compress_ptr cinfo, JSAMPARRAY input_buf, JSAMPIMAGE output_buf, JDIMENSION output_row, int num_rows) { void (*altivecfct) (JDIMENSION, JSAMPARRAY, JSAMPIMAGE, JDIMENSION, int); switch (cinfo->in_color_space) { case JCS_EXT_RGB: altivecfct = jsimd_extrgb_ycc_convert_altivec; break; case JCS_EXT_RGBX: case JCS_EXT_RGBA: altivecfct = jsimd_extrgbx_ycc_convert_altivec; break; case JCS_EXT_BGR: altivecfct = jsimd_extbgr_ycc_convert_altivec; break; case JCS_EXT_BGRX: case JCS_EXT_BGRA: altivecfct = jsimd_extbgrx_ycc_convert_altivec; break; case JCS_EXT_XBGR: case JCS_EXT_ABGR: altivecfct = jsimd_extxbgr_ycc_convert_altivec; break; case JCS_EXT_XRGB: case JCS_EXT_ARGB: altivecfct = jsimd_extxrgb_ycc_convert_altivec; break; default: altivecfct = jsimd_rgb_ycc_convert_altivec; break; } altivecfct(cinfo->image_width, input_buf, output_buf, output_row, num_rows); } GLOBAL(void) jsimd_rgb_gray_convert(j_compress_ptr cinfo, JSAMPARRAY input_buf, JSAMPIMAGE output_buf, JDIMENSION output_row, int num_rows) { void (*altivecfct) (JDIMENSION, JSAMPARRAY, JSAMPIMAGE, JDIMENSION, int); switch (cinfo->in_color_space) { case JCS_EXT_RGB: altivecfct = jsimd_extrgb_gray_convert_altivec; break; case JCS_EXT_RGBX: case JCS_EXT_RGBA: altivecfct = jsimd_extrgbx_gray_convert_altivec; break; case JCS_EXT_BGR: altivecfct = jsimd_extbgr_gray_convert_altivec; break; case JCS_EXT_BGRX: case JCS_EXT_BGRA: altivecfct = jsimd_extbgrx_gray_convert_altivec; break; case JCS_EXT_XBGR: case JCS_EXT_ABGR: altivecfct = jsimd_extxbgr_gray_convert_altivec; break; case JCS_EXT_XRGB: case JCS_EXT_ARGB: altivecfct = jsimd_extxrgb_gray_convert_altivec; break; default: altivecfct = jsimd_rgb_gray_convert_altivec; break; } altivecfct(cinfo->image_width, input_buf, output_buf, output_row, num_rows); } GLOBAL(void) jsimd_ycc_rgb_convert(j_decompress_ptr cinfo, JSAMPIMAGE input_buf, JDIMENSION input_row, JSAMPARRAY output_buf, int num_rows) { void (*altivecfct) (JDIMENSION, JSAMPIMAGE, JDIMENSION, JSAMPARRAY, int); switch (cinfo->out_color_space) { case JCS_EXT_RGB: altivecfct = jsimd_ycc_extrgb_convert_altivec; break; case JCS_EXT_RGBX: case JCS_EXT_RGBA: altivecfct = jsimd_ycc_extrgbx_convert_altivec; break; case JCS_EXT_BGR: altivecfct = jsimd_ycc_extbgr_convert_altivec; break; case JCS_EXT_BGRX: case JCS_EXT_BGRA: altivecfct = jsimd_ycc_extbgrx_convert_altivec; break; case JCS_EXT_XBGR: case JCS_EXT_ABGR: altivecfct = jsimd_ycc_extxbgr_convert_altivec; break; case JCS_EXT_XRGB: case JCS_EXT_ARGB: altivecfct = jsimd_ycc_extxrgb_convert_altivec; break; default: altivecfct = jsimd_ycc_rgb_convert_altivec; break; } altivecfct(cinfo->output_width, input_buf, input_row, output_buf, num_rows); } GLOBAL(void) jsimd_ycc_rgb565_convert(j_decompress_ptr cinfo, JSAMPIMAGE input_buf, JDIMENSION input_row, JSAMPARRAY output_buf, int num_rows) { } GLOBAL(int) jsimd_can_h2v2_downsample(void) { init_simd(); /* The code is optimised for these values only */ if (BITS_IN_JSAMPLE != 8) return 0; if (sizeof(JDIMENSION) != 4) return 0; if (simd_support & JSIMD_ALTIVEC) return 1; return 0; } GLOBAL(int) jsimd_can_h2v1_downsample(void) { init_simd(); /* The code is optimised for these values only */ if (BITS_IN_JSAMPLE != 8) return 0; if (sizeof(JDIMENSION) != 4) return 0; if (simd_support & JSIMD_ALTIVEC) return 1; return 0; } GLOBAL(void) jsimd_h2v2_downsample(j_compress_ptr cinfo, jpeg_component_info *compptr, JSAMPARRAY input_data, JSAMPARRAY output_data) { jsimd_h2v2_downsample_altivec(cinfo->image_width, cinfo->max_v_samp_factor, compptr->v_samp_factor, compptr->width_in_blocks, input_data, output_data); } GLOBAL(void) jsimd_h2v1_downsample(j_compress_ptr cinfo, jpeg_component_info *compptr, JSAMPARRAY input_data, JSAMPARRAY output_data) { jsimd_h2v1_downsample_altivec(cinfo->image_width, cinfo->max_v_samp_factor, compptr->v_samp_factor, compptr->width_in_blocks, input_data, output_data); } GLOBAL(int) jsimd_can_h2v2_upsample(void) { init_simd(); /* The code is optimised for these values only */ if (BITS_IN_JSAMPLE != 8) return 0; if (sizeof(JDIMENSION) != 4) return 0; if (simd_support & JSIMD_ALTIVEC) return 1; return 0; } GLOBAL(int) jsimd_can_h2v1_upsample(void) { init_simd(); /* The code is optimised for these values only */ if (BITS_IN_JSAMPLE != 8) return 0; if (sizeof(JDIMENSION) != 4) return 0; if (simd_support & JSIMD_ALTIVEC) return 1; return 0; } GLOBAL(void) jsimd_h2v2_upsample(j_decompress_ptr cinfo, jpeg_component_info *compptr, JSAMPARRAY input_data, JSAMPARRAY *output_data_ptr) { jsimd_h2v2_upsample_altivec(cinfo->max_v_samp_factor, cinfo->output_width, input_data, output_data_ptr); } GLOBAL(void) jsimd_h2v1_upsample(j_decompress_ptr cinfo, jpeg_component_info *compptr, JSAMPARRAY input_data, JSAMPARRAY *output_data_ptr) { jsimd_h2v1_upsample_altivec(cinfo->max_v_samp_factor, cinfo->output_width, input_data, output_data_ptr); } GLOBAL(int) jsimd_can_h2v2_fancy_upsample(void) { init_simd(); /* The code is optimised for these values only */ if (BITS_IN_JSAMPLE != 8) return 0; if (sizeof(JDIMENSION) != 4) return 0; if (simd_support & JSIMD_ALTIVEC) return 1; return 0; } GLOBAL(int) jsimd_can_h2v1_fancy_upsample(void) { init_simd(); /* The code is optimised for these values only */ if (BITS_IN_JSAMPLE != 8) return 0; if (sizeof(JDIMENSION) != 4) return 0; if (simd_support & JSIMD_ALTIVEC) return 1; return 0; } GLOBAL(void) jsimd_h2v2_fancy_upsample(j_decompress_ptr cinfo, jpeg_component_info *compptr, JSAMPARRAY input_data, JSAMPARRAY *output_data_ptr) { jsimd_h2v2_fancy_upsample_altivec(cinfo->max_v_samp_factor, compptr->downsampled_width, input_data, output_data_ptr); } GLOBAL(void) jsimd_h2v1_fancy_upsample(j_decompress_ptr cinfo, jpeg_component_info *compptr, JSAMPARRAY input_data, JSAMPARRAY *output_data_ptr) { jsimd_h2v1_fancy_upsample_altivec(cinfo->max_v_samp_factor, compptr->downsampled_width, input_data, output_data_ptr); } GLOBAL(int) jsimd_can_h2v2_merged_upsample(void) { init_simd(); /* The code is optimised for these values only */ if (BITS_IN_JSAMPLE != 8) return 0; if (sizeof(JDIMENSION) != 4) return 0; if (simd_support & JSIMD_ALTIVEC) return 1; return 0; } GLOBAL(int) jsimd_can_h2v1_merged_upsample(void) { init_simd(); /* The code is optimised for these values only */ if (BITS_IN_JSAMPLE != 8) return 0; if (sizeof(JDIMENSION) != 4) return 0; if (simd_support & JSIMD_ALTIVEC) return 1; return 0; } GLOBAL(void) jsimd_h2v2_merged_upsample(j_decompress_ptr cinfo, JSAMPIMAGE input_buf, JDIMENSION in_row_group_ctr, JSAMPARRAY output_buf) { void (*altivecfct) (JDIMENSION, JSAMPIMAGE, JDIMENSION, JSAMPARRAY); switch (cinfo->out_color_space) { case JCS_EXT_RGB: altivecfct = jsimd_h2v2_extrgb_merged_upsample_altivec; break; case JCS_EXT_RGBX: case JCS_EXT_RGBA: altivecfct = jsimd_h2v2_extrgbx_merged_upsample_altivec; break; case JCS_EXT_BGR: altivecfct = jsimd_h2v2_extbgr_merged_upsample_altivec; break; case JCS_EXT_BGRX: case JCS_EXT_BGRA: altivecfct = jsimd_h2v2_extbgrx_merged_upsample_altivec; break; case JCS_EXT_XBGR: case JCS_EXT_ABGR: altivecfct = jsimd_h2v2_extxbgr_merged_upsample_altivec; break; case JCS_EXT_XRGB: case JCS_EXT_ARGB: altivecfct = jsimd_h2v2_extxrgb_merged_upsample_altivec; break; default: altivecfct = jsimd_h2v2_merged_upsample_altivec; break; } altivecfct(cinfo->output_width, input_buf, in_row_group_ctr, output_buf); } GLOBAL(void) jsimd_h2v1_merged_upsample(j_decompress_ptr cinfo, JSAMPIMAGE input_buf, JDIMENSION in_row_group_ctr, JSAMPARRAY output_buf) { void (*altivecfct) (JDIMENSION, JSAMPIMAGE, JDIMENSION, JSAMPARRAY); switch (cinfo->out_color_space) { case JCS_EXT_RGB: altivecfct = jsimd_h2v1_extrgb_merged_upsample_altivec; break; case JCS_EXT_RGBX: case JCS_EXT_RGBA: altivecfct = jsimd_h2v1_extrgbx_merged_upsample_altivec; break; case JCS_EXT_BGR: altivecfct = jsimd_h2v1_extbgr_merged_upsample_altivec; break; case JCS_EXT_BGRX: case JCS_EXT_BGRA: altivecfct = jsimd_h2v1_extbgrx_merged_upsample_altivec; break; case JCS_EXT_XBGR: case JCS_EXT_ABGR: altivecfct = jsimd_h2v1_extxbgr_merged_upsample_altivec; break; case JCS_EXT_XRGB: case JCS_EXT_ARGB: altivecfct = jsimd_h2v1_extxrgb_merged_upsample_altivec; break; default: altivecfct = jsimd_h2v1_merged_upsample_altivec; break; } altivecfct(cinfo->output_width, input_buf, in_row_group_ctr, output_buf); } GLOBAL(int) jsimd_can_convsamp(void) { init_simd(); /* The code is optimised for these values only */ if (DCTSIZE != 8) return 0; if (BITS_IN_JSAMPLE != 8) return 0; if (sizeof(JDIMENSION) != 4) return 0; if (sizeof(DCTELEM) != 2) return 0; if (simd_support & JSIMD_ALTIVEC) return 1; return 0; } GLOBAL(int) jsimd_can_convsamp_float(void) { return 0; } GLOBAL(void) jsimd_convsamp(JSAMPARRAY sample_data, JDIMENSION start_col, DCTELEM *workspace) { jsimd_convsamp_altivec(sample_data, start_col, workspace); } GLOBAL(void) jsimd_convsamp_float(JSAMPARRAY sample_data, JDIMENSION start_col, FAST_FLOAT *workspace) { } GLOBAL(int) jsimd_can_fdct_islow(void) { init_simd(); /* The code is optimised for these values only */ if (DCTSIZE != 8) return 0; if (sizeof(DCTELEM) != 2) return 0; if (simd_support & JSIMD_ALTIVEC) return 1; return 0; } GLOBAL(int) jsimd_can_fdct_ifast(void) { init_simd(); /* The code is optimised for these values only */ if (DCTSIZE != 8) return 0; if (sizeof(DCTELEM) != 2) return 0; if (simd_support & JSIMD_ALTIVEC) return 1; return 0; } GLOBAL(int) jsimd_can_fdct_float(void) { return 0; } GLOBAL(void) jsimd_fdct_islow(DCTELEM *data) { jsimd_fdct_islow_altivec(data); } GLOBAL(void) jsimd_fdct_ifast(DCTELEM *data) { jsimd_fdct_ifast_altivec(data); } GLOBAL(void) jsimd_fdct_float(FAST_FLOAT *data) { } GLOBAL(int) jsimd_can_quantize(void) { init_simd(); /* The code is optimised for these values only */ if (DCTSIZE != 8) return 0; if (sizeof(JCOEF) != 2) return 0; if (sizeof(DCTELEM) != 2) return 0; if (simd_support & JSIMD_ALTIVEC) return 1; return 0; } GLOBAL(int) jsimd_can_quantize_float(void) { return 0; } GLOBAL(void) jsimd_quantize(JCOEFPTR coef_block, DCTELEM *divisors, DCTELEM *workspace) { jsimd_quantize_altivec(coef_block, divisors, workspace); } GLOBAL(void) jsimd_quantize_float(JCOEFPTR coef_block, FAST_FLOAT *divisors, FAST_FLOAT *workspace) { } GLOBAL(int) jsimd_can_idct_2x2(void) { return 0; } GLOBAL(int) jsimd_can_idct_4x4(void) { return 0; } GLOBAL(void) jsimd_idct_2x2(j_decompress_ptr cinfo, jpeg_component_info *compptr, JCOEFPTR coef_block, JSAMPARRAY output_buf, JDIMENSION output_col) { } GLOBAL(void) jsimd_idct_4x4(j_decompress_ptr cinfo, jpeg_component_info *compptr, JCOEFPTR coef_block, JSAMPARRAY output_buf, JDIMENSION output_col) { } GLOBAL(int) jsimd_can_idct_islow(void) { init_simd(); /* The code is optimised for these values only */ if (DCTSIZE != 8) return 0; if (sizeof(JCOEF) != 2) return 0; if (simd_support & JSIMD_ALTIVEC) return 1; return 0; } GLOBAL(int) jsimd_can_idct_ifast(void) { init_simd(); /* The code is optimised for these values only */ if (DCTSIZE != 8) return 0; if (sizeof(JCOEF) != 2) return 0; if (simd_support & JSIMD_ALTIVEC) return 1; return 0; } GLOBAL(int) jsimd_can_idct_float(void) { return 0; } GLOBAL(void) jsimd_idct_islow(j_decompress_ptr cinfo, jpeg_component_info *compptr, JCOEFPTR coef_block, JSAMPARRAY output_buf, JDIMENSION output_col) { jsimd_idct_islow_altivec(compptr->dct_table, coef_block, output_buf, output_col); } GLOBAL(void) jsimd_idct_ifast(j_decompress_ptr cinfo, jpeg_component_info *compptr, JCOEFPTR coef_block, JSAMPARRAY output_buf, JDIMENSION output_col) { jsimd_idct_ifast_altivec(compptr->dct_table, coef_block, output_buf, output_col); } GLOBAL(void) jsimd_idct_float(j_decompress_ptr cinfo, jpeg_component_info *compptr, JCOEFPTR coef_block, JSAMPARRAY output_buf, JDIMENSION output_col) { } GLOBAL(int) jsimd_can_huff_encode_one_block(void) { return 0; } GLOBAL(JOCTET *) jsimd_huff_encode_one_block(void *state, JOCTET *buffer, JCOEFPTR block, int last_dc_val, c_derived_tbl *dctbl, c_derived_tbl *actbl) { return NULL; } GLOBAL(int) jsimd_can_encode_mcu_AC_first_prepare(void) { return 0; } GLOBAL(void) jsimd_encode_mcu_AC_first_prepare(const JCOEF *block, const int *jpeg_natural_order_start, int Sl, int Al, JCOEF *values, size_t *zerobits) { } GLOBAL(int) jsimd_can_encode_mcu_AC_refine_prepare(void) { return 0; } GLOBAL(int) jsimd_encode_mcu_AC_refine_prepare(const JCOEF *block, const int *jpeg_natural_order_start, int Sl, int Al, JCOEF *absvalues, size_t *bits) { return 0; }