0N/A/*
2362N/A * Copyright (c) 1998, 2003, Oracle and/or its affiliates. All rights reserved.
0N/A * DO NOT ALTER OR REMOVE COPYRIGHT NOTICES OR THIS FILE HEADER.
0N/A *
0N/A * This code is free software; you can redistribute it and/or modify it
0N/A * under the terms of the GNU General Public License version 2 only, as
2362N/A * published by the Free Software Foundation. Oracle designates this
0N/A * particular file as subject to the "Classpath" exception as provided
2362N/A * by Oracle in the LICENSE file that accompanied this code.
0N/A *
0N/A * This code is distributed in the hope that it will be useful, but WITHOUT
0N/A * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
0N/A * FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
0N/A * version 2 for more details (a copy is included in the LICENSE file that
0N/A * accompanied this code).
0N/A *
0N/A * You should have received a copy of the GNU General Public License version
0N/A * 2 along with this work; if not, write to the Free Software Foundation,
0N/A * Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA.
0N/A *
2362N/A * Please contact Oracle, 500 Oracle Parkway, Redwood Shores, CA 94065 USA
2362N/A * or visit www.oracle.com if you need additional information or have any
2362N/A * questions.
0N/A */
0N/A
0N/A
0N/A#include "mlib_image.h"
0N/A#include "mlib_ImageAffine.h"
0N/A#include "mlib_ImageColormap.h"
0N/A
0N/A/***************************************************************/
0N/A#define MLIB_LIMIT 512
0N/A#define MLIB_SHIFT 16
0N/A#define MLIB_PREC (1 << MLIB_SHIFT)
0N/A#define MLIB_MASK (MLIB_PREC - 1)
0N/A
0N/A/***************************************************************/
0N/A#define DTYPE MLIB_TYPE
0N/A
0N/A/***************************************************************/
0N/A#define DECLAREVAR_IND() \
0N/A DECLAREVAR0(); \
0N/A mlib_s32 *warp_tbl = param -> warp_tbl; \
0N/A mlib_s32 xSrc, ySrc; \
0N/A mlib_s32 srcYStride = param -> srcYStride; \
0N/A mlib_s32 max_xsize = param -> max_xsize; \
0N/A mlib_filter filter = param -> filter; \
0N/A MLIB_TYPE *sp, *dl; \
0N/A mlib_d64 xf0, xf1, xf2, xf3; \
0N/A mlib_d64 yf0, yf1, yf2, yf3; \
0N/A mlib_d64 c0, c1, c2, c3, val0; \
0N/A mlib_s32 filterpos; \
0N/A mlib_f32 *fptr; \
0N/A mlib_d64 s0, s1, s2, s3; \
0N/A mlib_s32 i, size
0N/A
0N/A/***************************************************************/
0N/A#define GET_FILTERS_KOEF() \
0N/A filterpos = (X >> FILTER_SHIFT) & FILTER_MASK; \
0N/A fptr = (mlib_f32 *) ((mlib_u8 *)mlib_filters_table + filterpos); \
0N/A \
0N/A xf0 = fptr[0]; \
0N/A xf1 = fptr[1]; \
0N/A xf2 = fptr[2]; \
0N/A xf3 = fptr[3]; \
0N/A \
0N/A filterpos = (Y >> FILTER_SHIFT) & FILTER_MASK; \
0N/A fptr = (mlib_f32 *) ((mlib_u8 *)mlib_filters_table + filterpos); \
0N/A \
0N/A yf0 = fptr[0]; \
0N/A yf1 = fptr[1]; \
0N/A yf2 = fptr[2]; \
0N/A yf3 = fptr[3]
0N/A
0N/A/***************************************************************/
0N/A#define GET_POINTER() \
0N/A xSrc = (X >> MLIB_SHIFT)-1; \
0N/A ySrc = (Y >> MLIB_SHIFT)-1; \
0N/A sp = ((MLIB_TYPE **)lineAddr)[ySrc] + xSrc
0N/A
0N/A/***************************************************************/
0N/A#define LOAD_FIRST_ROW(nchan, chan) \
0N/A s0 = *(lut + sp[0]*nchan + chan); \
0N/A s1 = *(lut + sp[1]*nchan + chan); \
0N/A s2 = *(lut + sp[2]*nchan + chan); \
0N/A s3 = *(lut + sp[3]*nchan + chan)
0N/A
0N/A/***************************************************************/
0N/A#define COUNT_NEXT_ROW(dst, nchan, chan) \
0N/A sp = (MLIB_TYPE*)((mlib_addr)sp + srcYStride); \
0N/A dst = ((*(lut + sp[0]*nchan + chan))*xf0 + \
0N/A (*(lut + sp[1]*nchan + chan))*xf1 + \
0N/A (*(lut + sp[2]*nchan + chan))*xf2 + \
0N/A (*(lut + sp[3]*nchan + chan))*xf3)
0N/A
0N/A/***************************************************************/
0N/A#ifdef MLIB_USE_FTOI_CLAMPING
0N/A
0N/A/***********/
0N/A#define STORE_SAT_VALUE_U8(ind) \
0N/A dp[ind] = ((mlib_s32)(val0 - (mlib_d64)0x7F800000) >> 24) ^ 0x80
0N/A
0N/A/***********/
0N/A#define STORE_SAT_VALUE_S16(ind) \
0N/A dp[ind] = ((mlib_s32)(val0)) >> 16
0N/A
0N/A#else
0N/A
0N/A/***********/
0N/A#define STORE_SAT_VALUE_U8(ind) \
0N/A val0 -= (mlib_d64)0x7F800000; \
0N/A if (val0 >= MLIB_S32_MAX) \
0N/A dp[ind] = MLIB_U8_MAX; \
0N/A else if (val0 <= MLIB_S32_MIN) \
0N/A dp[ind] = MLIB_U8_MIN; \
0N/A else \
0N/A dp[ind] = ((mlib_s32)val0 >> 24) ^ 0x80
0N/A
0N/A/***********/
0N/A#define STORE_SAT_VALUE_S16(ind) \
0N/A if (val0 >= MLIB_S32_MAX) \
0N/A dp[ind] = MLIB_S16_MAX; \
0N/A else if (val0 <= MLIB_S32_MIN) \
0N/A dp[ind] = MLIB_S16_MIN; \
0N/A else \
0N/A dp[ind] = (mlib_s32)val0 >> 16
0N/A
0N/A#endif /* MLIB_USE_FTOI_CLAMPING */
0N/A
0N/A/***************************************************************/
0N/A#define MAKE_BC_3CH(lut_format) \
0N/A X += dX; \
0N/A Y += dY; \
0N/A c0 = (s0*xf0 + s1*xf1 + s2*xf2 + s3*xf3); \
0N/A COUNT_NEXT_ROW(c1, 3, 0); \
0N/A COUNT_NEXT_ROW(c2, 3, 0); \
0N/A COUNT_NEXT_ROW(c3, 3, 0); \
0N/A val0 = (c0*yf0 + c1*yf1 + c2*yf2 + c3*yf3); \
0N/A sp = (MLIB_TYPE*)((mlib_addr)sp - 3*srcYStride); \
0N/A LOAD_FIRST_ROW(3, 1); \
0N/A STORE_SAT_VALUE_##lut_format(0); \
0N/A c0 = (s0*xf0 + s1*xf1 + s2*xf2 + s3*xf3); \
0N/A COUNT_NEXT_ROW(c1, 3, 1); \
0N/A COUNT_NEXT_ROW(c2, 3, 1); \
0N/A COUNT_NEXT_ROW(c3, 3, 1); \
0N/A val0 = (c0*yf0 + c1*yf1 + c2*yf2 + c3*yf3); \
0N/A sp = (MLIB_TYPE*)((mlib_addr)sp - 3*srcYStride); \
0N/A LOAD_FIRST_ROW(3, 2); \
0N/A STORE_SAT_VALUE_##lut_format(1); \
0N/A c0 = (s0*xf0 + s1*xf1 + s2*xf2 + s3*xf3); \
0N/A COUNT_NEXT_ROW(c1, 3, 2); \
0N/A COUNT_NEXT_ROW(c2, 3, 2); \
0N/A COUNT_NEXT_ROW(c3, 3, 2); \
0N/A val0 = (c0*yf0 + c1*yf1 + c2*yf2 + c3*yf3); \
0N/A GET_FILTERS_KOEF(); \
0N/A GET_POINTER(); \
0N/A LOAD_FIRST_ROW(3, 0); \
0N/A STORE_SAT_VALUE_##lut_format(2);
0N/A
0N/A/***************************************************************/
0N/A#define MAKE_LAST_PIXEL_BC_3CH(lut_format) \
0N/A c0 = (s0*xf0 + s1*xf1 + s2*xf2 + s3*xf3); \
0N/A COUNT_NEXT_ROW(c1, 3, 0); \
0N/A COUNT_NEXT_ROW(c2, 3, 0); \
0N/A COUNT_NEXT_ROW(c3, 3, 0); \
0N/A val0 = (c0*yf0 + c1*yf1 + c2*yf2 + c3*yf3); \
0N/A sp = (MLIB_TYPE*)((mlib_addr)sp - 3*srcYStride); \
0N/A LOAD_FIRST_ROW(3, 1); \
0N/A STORE_SAT_VALUE_##lut_format(0); \
0N/A c0 = (s0*xf0 + s1*xf1 + s2*xf2 + s3*xf3); \
0N/A COUNT_NEXT_ROW(c1, 3, 1); \
0N/A COUNT_NEXT_ROW(c2, 3, 1); \
0N/A COUNT_NEXT_ROW(c3, 3, 1); \
0N/A val0 = (c0*yf0 + c1*yf1 + c2*yf2 + c3*yf3); \
0N/A sp = (MLIB_TYPE*)((mlib_addr)sp - 3*srcYStride); \
0N/A LOAD_FIRST_ROW(3, 2); \
0N/A STORE_SAT_VALUE_##lut_format(1); \
0N/A c0 = (s0*xf0 + s1*xf1 + s2*xf2 + s3*xf3); \
0N/A COUNT_NEXT_ROW(c1, 3, 2); \
0N/A COUNT_NEXT_ROW(c2, 3, 2); \
0N/A COUNT_NEXT_ROW(c3, 3, 2); \
0N/A val0 = (c0*yf0 + c1*yf1 + c2*yf2 + c3*yf3); \
0N/A STORE_SAT_VALUE_##lut_format(2);
0N/A
0N/A/***************************************************************/
0N/A#define MAKE_BC_4CH(lut_format) \
0N/A X += dX; \
0N/A Y += dY; \
0N/A c0 = (s0*xf0 + s1*xf1 + s2*xf2 + s3*xf3); \
0N/A COUNT_NEXT_ROW(c1, 4, 0); \
0N/A COUNT_NEXT_ROW(c2, 4, 0); \
0N/A COUNT_NEXT_ROW(c3, 4, 0); \
0N/A val0 = (c0*yf0 + c1*yf1 + c2*yf2 + c3*yf3); \
0N/A sp = (MLIB_TYPE*)((mlib_addr)sp - 3*srcYStride); \
0N/A LOAD_FIRST_ROW(4, 1); \
0N/A STORE_SAT_VALUE_##lut_format(0); \
0N/A c0 = (s0*xf0 + s1*xf1 + s2*xf2 + s3*xf3); \
0N/A COUNT_NEXT_ROW(c1, 4, 1); \
0N/A COUNT_NEXT_ROW(c2, 4, 1); \
0N/A COUNT_NEXT_ROW(c3, 4, 1); \
0N/A val0 = (c0*yf0 + c1*yf1 + c2*yf2 + c3*yf3); \
0N/A sp = (MLIB_TYPE*)((mlib_addr)sp - 3*srcYStride); \
0N/A LOAD_FIRST_ROW(4, 2); \
0N/A STORE_SAT_VALUE_##lut_format(1); \
0N/A c0 = (s0*xf0 + s1*xf1 + s2*xf2 + s3*xf3); \
0N/A COUNT_NEXT_ROW(c1, 4, 2); \
0N/A COUNT_NEXT_ROW(c2, 4, 2); \
0N/A COUNT_NEXT_ROW(c3, 4, 2); \
0N/A val0 = (c0*yf0 + c1*yf1 + c2*yf2 + c3*yf3); \
0N/A sp = (MLIB_TYPE*)((mlib_addr)sp - 3*srcYStride); \
0N/A LOAD_FIRST_ROW(4, 3); \
0N/A STORE_SAT_VALUE_##lut_format(2); \
0N/A c0 = (s0*xf0 + s1*xf1 + s2*xf2 + s3*xf3); \
0N/A COUNT_NEXT_ROW(c1, 4, 3); \
0N/A COUNT_NEXT_ROW(c2, 4, 3); \
0N/A COUNT_NEXT_ROW(c3, 4, 3); \
0N/A val0 = (c0*yf0 + c1*yf1 + c2*yf2 + c3*yf3); \
0N/A GET_FILTERS_KOEF(); \
0N/A GET_POINTER(); \
0N/A LOAD_FIRST_ROW(4, 0); \
0N/A STORE_SAT_VALUE_##lut_format(3);
0N/A
0N/A/***************************************************************/
0N/A#define MAKE_LAST_PIXEL_BC_4CH(lut_format) \
0N/A c0 = (s0*xf0 + s1*xf1 + s2*xf2 + s3*xf3); \
0N/A COUNT_NEXT_ROW(c1, 4, 0); \
0N/A COUNT_NEXT_ROW(c2, 4, 0); \
0N/A COUNT_NEXT_ROW(c3, 4, 0); \
0N/A val0 = (c0*yf0 + c1*yf1 + c2*yf2 + c3*yf3); \
0N/A sp = (MLIB_TYPE*)((mlib_addr)sp - 3*srcYStride); \
0N/A LOAD_FIRST_ROW(4, 1); \
0N/A STORE_SAT_VALUE_##lut_format(0); \
0N/A c0 = (s0*xf0 + s1*xf1 + s2*xf2 + s3*xf3); \
0N/A COUNT_NEXT_ROW(c1, 4, 1); \
0N/A COUNT_NEXT_ROW(c2, 4, 1); \
0N/A COUNT_NEXT_ROW(c3, 4, 1); \
0N/A val0 = (c0*yf0 + c1*yf1 + c2*yf2 + c3*yf3); \
0N/A sp = (MLIB_TYPE*)((mlib_addr)sp - 3*srcYStride); \
0N/A LOAD_FIRST_ROW(4, 2); \
0N/A STORE_SAT_VALUE_##lut_format(1); \
0N/A c0 = (s0*xf0 + s1*xf1 + s2*xf2 + s3*xf3); \
0N/A COUNT_NEXT_ROW(c1, 4, 2); \
0N/A COUNT_NEXT_ROW(c2, 4, 2); \
0N/A COUNT_NEXT_ROW(c3, 4, 2); \
0N/A val0 = (c0*yf0 + c1*yf1 + c2*yf2 + c3*yf3); \
0N/A sp = (MLIB_TYPE*)((mlib_addr)sp - 3*srcYStride); \
0N/A LOAD_FIRST_ROW(4, 3); \
0N/A STORE_SAT_VALUE_##lut_format(2); \
0N/A c0 = (s0*xf0 + s1*xf1 + s2*xf2 + s3*xf3); \
0N/A COUNT_NEXT_ROW(c1, 4, 3); \
0N/A COUNT_NEXT_ROW(c2, 4, 3); \
0N/A COUNT_NEXT_ROW(c3, 4, 3); \
0N/A val0 = (c0*yf0 + c1*yf1 + c2*yf2 + c3*yf3); \
0N/A STORE_SAT_VALUE_##lut_format(3);
0N/A
0N/A/***************************************************************/
0N/A#define FILTER_U8 ((filter == MLIB_BICUBIC) ? mlib_filters_u8f_bc : mlib_filters_u8f_bc2)
0N/A#define FILTER_S16 ((filter == MLIB_BICUBIC) ? mlib_filters_s16f_bc : mlib_filters_s16f_bc2)
0N/A
0N/A/***************************************************************/
0N/A#define mlib_U8 mlib_u8
0N/A#define mlib_S16 mlib_s16
0N/A
0N/A/***************************************************************/
0N/A#define FUNC_AFFINEINDEX_BC_0(ITYPE, LTYPE, NCHAN) \
0N/A mlib_status mlib_ImageAffineIndex_##ITYPE##_##LTYPE##_##NCHAN##CH_BC(mlib_affine_param *param, \
0N/A const void *colormap) \
0N/A { \
0N/A DECLAREVAR_IND(); \
0N/A mlib_##LTYPE buff_lcl[NCHAN * MLIB_LIMIT], *pbuff = buff_lcl, *dp; \
0N/A mlib_d64 *lut = ((mlib_d64*)mlib_ImageGetLutDoubleData(colormap) - \
0N/A NCHAN * mlib_ImageGetLutOffset(colormap)); \
0N/A const mlib_f32 *mlib_filters_table = FILTER_##LTYPE; \
0N/A \
0N/A if (max_xsize > MLIB_LIMIT) { \
0N/A pbuff = mlib_malloc(NCHAN * sizeof(mlib_##LTYPE) * max_xsize); \
0N/A if (pbuff == NULL) return MLIB_FAILURE; \
0N/A } \
0N/A \
0N/A for (j = yStart; j <= yFinish; j++) { \
0N/A \
0N/A NEW_LINE(1); \
0N/A dp = pbuff; \
0N/A \
0N/A GET_FILTERS_KOEF(); \
0N/A GET_POINTER(); \
0N/A LOAD_FIRST_ROW(NCHAN, 0);
0N/A
0N/A /* pragma pipeloop(0) must be here */
0N/A
0N/A/***************************************************************/
0N/A#define FUNC_AFFINEINDEX_BC_1(ITYPE, LTYPE, NCHAN) \
0N/A \
0N/A for (i = 0; i < (xRight - xLeft); i++, dp += NCHAN) { \
0N/A MAKE_BC_##NCHAN##CH(LTYPE); \
0N/A } \
0N/A \
0N/A MAKE_LAST_PIXEL_BC_##NCHAN##CH(LTYPE); \
0N/A \
0N/A mlib_ImageColorTrue2IndexLine_##LTYPE##_##ITYPE##_##NCHAN \
0N/A (pbuff, dl, xRight - xLeft + 1, colormap); \
0N/A } \
0N/A \
0N/A if (pbuff != buff_lcl) mlib_free(pbuff); \
0N/A \
0N/A return MLIB_SUCCESS; \
0N/A }
0N/A
0N/A/***************************************************************/
0N/A#undef MLIB_TYPE
0N/A#define MLIB_TYPE mlib_u8
0N/A
0N/A/***************************************************************/
0N/A#undef FILTER_SHIFT
0N/A#define FILTER_SHIFT 4
0N/A#undef FILTER_MASK
0N/A#define FILTER_MASK (((1 << 8) - 1) << 4)
0N/A
0N/AFUNC_AFFINEINDEX_BC_0(U8, U8, 3)
0N/A#ifdef __SUNPRO_C
0N/A#pragma pipeloop(0)
0N/A#endif /* __SUNPRO_C */
0N/AFUNC_AFFINEINDEX_BC_1(U8, U8, 3)
0N/A
0N/A/***************************************************************/
0N/A#undef FILTER_SHIFT
0N/A#define FILTER_SHIFT 3
0N/A#undef FILTER_MASK
0N/A#define FILTER_MASK (((1 << 9) - 1) << 4)
0N/A
0N/AFUNC_AFFINEINDEX_BC_0(U8, S16, 3)
0N/A#ifdef __SUNPRO_C
0N/A#pragma pipeloop(0)
0N/A#endif /* __SUNPRO_C */
0N/AFUNC_AFFINEINDEX_BC_1(U8, S16, 3)
0N/A
0N/A/***************************************************************/
0N/A#undef FILTER_SHIFT
0N/A#define FILTER_SHIFT 4
0N/A#undef FILTER_MASK
0N/A#define FILTER_MASK (((1 << 8) - 1) << 4)
0N/A
0N/AFUNC_AFFINEINDEX_BC_0(U8, U8, 4)
0N/A#ifdef __SUNPRO_C
0N/A#pragma pipeloop(0)
0N/A#endif /* __SUNPRO_C */
0N/AFUNC_AFFINEINDEX_BC_1(U8, U8, 4)
0N/A
0N/A/***************************************************************/
0N/A#undef FILTER_SHIFT
0N/A#define FILTER_SHIFT 3
0N/A#undef FILTER_MASK
0N/A#define FILTER_MASK (((1 << 9) - 1) << 4)
0N/A
0N/AFUNC_AFFINEINDEX_BC_0(U8, S16, 4)
0N/A#ifdef __SUNPRO_C
0N/A#pragma pipeloop(0)
0N/A#endif /* __SUNPRO_C */
0N/AFUNC_AFFINEINDEX_BC_1(U8, S16, 4)
0N/A
0N/A/***************************************************************/
0N/A#undef MLIB_TYPE
0N/A#define MLIB_TYPE mlib_s16
0N/A
0N/A/***************************************************************/
0N/A#undef FILTER_SHIFT
0N/A#define FILTER_SHIFT 4
0N/A#undef FILTER_MASK
0N/A#define FILTER_MASK (((1 << 8) - 1) << 4)
0N/A
0N/AFUNC_AFFINEINDEX_BC_0(S16, U8, 3)
0N/A#ifdef __SUNPRO_C
0N/A#pragma pipeloop(0)
0N/A#endif /* __SUNPRO_C */
0N/AFUNC_AFFINEINDEX_BC_1(S16, U8, 3)
0N/A
0N/A/***************************************************************/
0N/A#undef FILTER_SHIFT
0N/A#define FILTER_SHIFT 3
0N/A#undef FILTER_MASK
0N/A#define FILTER_MASK (((1 << 9) - 1) << 4)
0N/A
0N/AFUNC_AFFINEINDEX_BC_0(S16, S16, 3)
0N/A#ifdef __SUNPRO_C
0N/A#pragma pipeloop(0)
0N/A#endif /* __SUNPRO_C */
0N/AFUNC_AFFINEINDEX_BC_1(S16, S16, 3)
0N/A
0N/A/***************************************************************/
0N/A#undef FILTER_SHIFT
0N/A#define FILTER_SHIFT 4
0N/A#undef FILTER_MASK
0N/A#define FILTER_MASK (((1 << 8) - 1) << 4)
0N/A
0N/AFUNC_AFFINEINDEX_BC_0(S16, U8, 4)
0N/A#ifdef __SUNPRO_C
0N/A#pragma pipeloop(0)
0N/A#endif /* __SUNPRO_C */
0N/AFUNC_AFFINEINDEX_BC_1(S16, U8, 4)
0N/A
0N/A/***************************************************************/
0N/A#undef FILTER_SHIFT
0N/A#define FILTER_SHIFT 3
0N/A#undef FILTER_MASK
0N/A#define FILTER_MASK (((1 << 9) - 1) << 4)
0N/A
0N/AFUNC_AFFINEINDEX_BC_0(S16, S16, 4)
0N/A#ifdef __SUNPRO_C
0N/A#pragma pipeloop(0)
0N/A#endif /* __SUNPRO_C */
0N/AFUNC_AFFINEINDEX_BC_1(S16, S16, 4)
0N/A
0N/A/***************************************************************/