rlm@1: /*
rlm@1:  * This file is part of the Advance project.
rlm@1:  *
rlm@1:  * Copyright (C) 2003 Andrea Mazzoleni
rlm@1:  *
rlm@1:  * This program is free software; you can redistribute it and/or modify
rlm@1:  * it under the terms of the GNU General Public License as published by
rlm@1:  * the Free Software Foundation; either version 2 of the License, or
rlm@1:  * (at your option) any later version.
rlm@1:  *
rlm@1:  * This program is distributed in the hope that it will be useful,
rlm@1:  * but WITHOUT ANY WARRANTY; without even the implied warranty of
rlm@1:  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
rlm@1:  * GNU General Public License for more details.
rlm@1:  *
rlm@1:  * You should have received a copy of the GNU General Public License
rlm@1:  * along with this program; if not, write to the Free Software
rlm@1:  * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
rlm@1:  *
rlm@1:  * In addition, as a special exception, Andrea Mazzoleni
rlm@1:  * gives permission to link the code of this program with
rlm@1:  * the MAME library (or with modified versions of MAME that use the
rlm@1:  * same license as MAME), and distribute linked combinations including
rlm@1:  * the two.  You must obey the GNU General Public License in all
rlm@1:  * respects for all of the code used other than MAME.  If you modify
rlm@1:  * this file, you may extend this exception to your version of the
rlm@1:  * file, but you are not obligated to do so.  If you do not wish to
rlm@1:  * do so, delete this exception statement from your version.
rlm@1:  */
rlm@1: 
rlm@1: #ifndef __INTERP_H
rlm@1: #define __INTERP_H
rlm@1: 
rlm@1: /***************************************************************************/
rlm@1: /* Basic types */
rlm@1: 
rlm@1: /***************************************************************************/
rlm@1: /* interpolation */
rlm@1: 
rlm@1: extern unsigned interp_mask[2];
rlm@1: extern unsigned interp_bits_per_pixel;
rlm@1: 
rlm@1: #define INTERP_16_MASK_1(v) (v & interp_mask[0])
rlm@1: #define INTERP_16_MASK_2(v) (v & interp_mask[1])
rlm@1: 
rlm@1: static inline u16 interp_16_521(u16 p1, u16 p2, u16 p3)
rlm@1: {
rlm@1:   return INTERP_16_MASK_1((INTERP_16_MASK_1(p1)*5 + INTERP_16_MASK_1(p2)*2 + INTERP_16_MASK_1(p3)*1) / 8)
rlm@1:     | INTERP_16_MASK_2((INTERP_16_MASK_2(p1)*5 + INTERP_16_MASK_2(p2)*2 + INTERP_16_MASK_2(p3)*1) / 8);
rlm@1: }
rlm@1: 
rlm@1: static inline u16 interp_16_332(u16 p1, u16 p2, u16 p3)
rlm@1: {
rlm@1:   return INTERP_16_MASK_1((INTERP_16_MASK_1(p1)*3 + INTERP_16_MASK_1(p2)*3 + INTERP_16_MASK_1(p3)*2) / 8)
rlm@1:     | INTERP_16_MASK_2((INTERP_16_MASK_2(p1)*3 + INTERP_16_MASK_2(p2)*3 + INTERP_16_MASK_2(p3)*2) / 8);
rlm@1: }
rlm@1: 
rlm@1: static inline u16 interp_16_611(u16 p1, u16 p2, u16 p3)
rlm@1: {
rlm@1:   return INTERP_16_MASK_1((INTERP_16_MASK_1(p1)*6 + INTERP_16_MASK_1(p2) + INTERP_16_MASK_1(p3)) / 8)
rlm@1:     | INTERP_16_MASK_2((INTERP_16_MASK_2(p1)*6 + INTERP_16_MASK_2(p2) + INTERP_16_MASK_2(p3)) / 8);
rlm@1: }
rlm@1: 
rlm@1: static inline u16 interp_16_71(u16 p1, u16 p2)
rlm@1: {
rlm@1:   return INTERP_16_MASK_1((INTERP_16_MASK_1(p1)*7 + INTERP_16_MASK_1(p2)) / 8)
rlm@1:     | INTERP_16_MASK_2((INTERP_16_MASK_2(p1)*7 + INTERP_16_MASK_2(p2)) / 8);
rlm@1: }
rlm@1: 
rlm@1: static inline u16 interp_16_211(u16 p1, u16 p2, u16 p3)
rlm@1: {
rlm@1:   return INTERP_16_MASK_1((INTERP_16_MASK_1(p1)*2 + INTERP_16_MASK_1(p2) + INTERP_16_MASK_1(p3)) / 4)
rlm@1:     | INTERP_16_MASK_2((INTERP_16_MASK_2(p1)*2 + INTERP_16_MASK_2(p2) + INTERP_16_MASK_2(p3)) / 4);
rlm@1: }
rlm@1: 
rlm@1: static inline u16 interp_16_772(u16 p1, u16 p2, u16 p3)
rlm@1: {
rlm@1:   return INTERP_16_MASK_1(((INTERP_16_MASK_1(p1) + INTERP_16_MASK_1(p2))*7 + INTERP_16_MASK_1(p3)*2) / 16)
rlm@1:     | INTERP_16_MASK_2(((INTERP_16_MASK_2(p1) + INTERP_16_MASK_2(p2))*7 + INTERP_16_MASK_2(p3)*2) / 16);
rlm@1: }
rlm@1: 
rlm@1: static inline u16 interp_16_11(u16 p1, u16 p2)
rlm@1: {
rlm@1:   return INTERP_16_MASK_1((INTERP_16_MASK_1(p1) + INTERP_16_MASK_1(p2)) / 2)
rlm@1:     | INTERP_16_MASK_2((INTERP_16_MASK_2(p1) + INTERP_16_MASK_2(p2)) / 2);
rlm@1: }
rlm@1: 
rlm@1: static inline u16 interp_16_31(u16 p1, u16 p2)
rlm@1: {
rlm@1:   return INTERP_16_MASK_1((INTERP_16_MASK_1(p1)*3 + INTERP_16_MASK_1(p2)) / 4)
rlm@1:     | INTERP_16_MASK_2((INTERP_16_MASK_2(p1)*3 + INTERP_16_MASK_2(p2)) / 4);
rlm@1: }
rlm@1: 
rlm@1: static inline u16 interp_16_1411(u16 p1, u16 p2, u16 p3)
rlm@1: {
rlm@1:   return INTERP_16_MASK_1((INTERP_16_MASK_1(p1)*14 + INTERP_16_MASK_1(p2) + INTERP_16_MASK_1(p3)) / 16)
rlm@1:     | INTERP_16_MASK_2((INTERP_16_MASK_2(p1)*14 + INTERP_16_MASK_2(p2) + INTERP_16_MASK_2(p3)) / 16);
rlm@1: }
rlm@1: 
rlm@1: static inline u16 interp_16_431(u16 p1, u16 p2, u16 p3)
rlm@1: {
rlm@1:   return INTERP_16_MASK_1((INTERP_16_MASK_1(p1)*4 + INTERP_16_MASK_1(p2)*3 + INTERP_16_MASK_1(p3)) / 8)
rlm@1:     | INTERP_16_MASK_2((INTERP_16_MASK_2(p1)*4 + INTERP_16_MASK_2(p2)*3 + INTERP_16_MASK_2(p3)) / 8);
rlm@1: }
rlm@1: 
rlm@1: static inline u16 interp_16_53(u16 p1, u16 p2)
rlm@1: {
rlm@1:   return INTERP_16_MASK_1((INTERP_16_MASK_1(p1)*5 + INTERP_16_MASK_1(p2)*3) / 8)
rlm@1:     | INTERP_16_MASK_2((INTERP_16_MASK_2(p1)*5 + INTERP_16_MASK_2(p2)*3) / 8);
rlm@1: }
rlm@1: 
rlm@1: static inline u16 interp_16_151(u16 p1, u16 p2)
rlm@1: {
rlm@1:   return INTERP_16_MASK_1((INTERP_16_MASK_1(p1)*15 + INTERP_16_MASK_1(p2)) / 16)
rlm@1:     | INTERP_16_MASK_2((INTERP_16_MASK_2(p1)*15 + INTERP_16_MASK_2(p2)) / 16);
rlm@1: }
rlm@1: 
rlm@1: static inline u16 interp_16_97(u16 p1, u16 p2)
rlm@1: {
rlm@1:   return INTERP_16_MASK_1((INTERP_16_MASK_1(p1)*9 + INTERP_16_MASK_1(p2)*7) / 16)
rlm@1:     | INTERP_16_MASK_2((INTERP_16_MASK_2(p1)*9 + INTERP_16_MASK_2(p2)*7) / 16);
rlm@1: }
rlm@1: 
rlm@1: #define INTERP_32_MASK_1(v) (v & 0xFF00FF)
rlm@1: #define INTERP_32_MASK_2(v) (v & 0x00FF00)
rlm@1: 
rlm@1: static inline u32 interp_32_521(u32 p1, u32 p2, u32 p3)
rlm@1: {
rlm@1:   return INTERP_32_MASK_1((INTERP_32_MASK_1(p1)*5 + INTERP_32_MASK_1(p2)*2 + INTERP_32_MASK_1(p3)*1) / 8)
rlm@1:     | INTERP_32_MASK_2((INTERP_32_MASK_2(p1)*5 + INTERP_32_MASK_2(p2)*2 + INTERP_32_MASK_2(p3)*1) / 8);
rlm@1: }
rlm@1: 
rlm@1: static inline u32 interp_32_332(u32 p1, u32 p2, u32 p3)
rlm@1: {
rlm@1:   return INTERP_32_MASK_1((INTERP_32_MASK_1(p1)*3 + INTERP_32_MASK_1(p2)*3 + INTERP_32_MASK_1(p3)*2) / 8)
rlm@1:     | INTERP_32_MASK_2((INTERP_32_MASK_2(p1)*3 + INTERP_32_MASK_2(p2)*3 + INTERP_32_MASK_2(p3)*2) / 8);
rlm@1: }
rlm@1: 
rlm@1: static inline u32 interp_32_211(u32 p1, u32 p2, u32 p3)
rlm@1: {
rlm@1:   return INTERP_32_MASK_1((INTERP_32_MASK_1(p1)*2 + INTERP_32_MASK_1(p2) + INTERP_32_MASK_1(p3)) / 4)
rlm@1:     | INTERP_32_MASK_2((INTERP_32_MASK_2(p1)*2 + INTERP_32_MASK_2(p2) + INTERP_32_MASK_2(p3)) / 4);
rlm@1: }
rlm@1: 
rlm@1: static inline u32 interp_32_611(u32 p1, u32 p2, u32 p3)
rlm@1: {
rlm@1:   return INTERP_32_MASK_1((INTERP_32_MASK_1(p1)*6 + INTERP_32_MASK_1(p2) + INTERP_32_MASK_1(p3)) / 8)
rlm@1:     | INTERP_32_MASK_2((INTERP_32_MASK_2(p1)*6 + INTERP_32_MASK_2(p2) + INTERP_32_MASK_2(p3)) / 8);
rlm@1: }
rlm@1: 
rlm@1: static inline u32 interp_32_71(u32 p1, u32 p2)
rlm@1: {
rlm@1:   return INTERP_32_MASK_1((INTERP_32_MASK_1(p1)*7 + INTERP_32_MASK_1(p2)) / 8)
rlm@1:     | INTERP_32_MASK_2((INTERP_32_MASK_2(p1)*7 + INTERP_32_MASK_2(p2)) / 8);
rlm@1: }
rlm@1: 
rlm@1: static inline u32 interp_32_772(u32 p1, u32 p2, u32 p3)
rlm@1: {
rlm@1:   return INTERP_32_MASK_1(((INTERP_32_MASK_1(p1) + INTERP_32_MASK_1(p2))*7 + INTERP_32_MASK_1(p3)*2) / 16)
rlm@1:     | INTERP_32_MASK_2(((INTERP_32_MASK_2(p1) + INTERP_32_MASK_2(p2))*7 + INTERP_32_MASK_2(p3)*2) / 16);
rlm@1: }
rlm@1: 
rlm@1: static inline u32 interp_32_11(u32 p1, u32 p2)
rlm@1: {
rlm@1:   return INTERP_32_MASK_1((INTERP_32_MASK_1(p1) + INTERP_32_MASK_1(p2)) / 2)
rlm@1:     | INTERP_32_MASK_2((INTERP_32_MASK_2(p1) + INTERP_32_MASK_2(p2)) / 2);
rlm@1: }
rlm@1: 
rlm@1: static inline u32 interp_32_31(u32 p1, u32 p2)
rlm@1: {
rlm@1:   return INTERP_32_MASK_1((INTERP_32_MASK_1(p1)*3 + INTERP_32_MASK_1(p2)) / 4)
rlm@1:     | INTERP_32_MASK_2((INTERP_32_MASK_2(p1)*3 + INTERP_32_MASK_2(p2)) / 4);
rlm@1: }
rlm@1: 
rlm@1: static inline u32 interp_32_1411(u32 p1, u32 p2, u32 p3)
rlm@1: {
rlm@1:   return INTERP_32_MASK_1((INTERP_32_MASK_1(p1)*14 + INTERP_32_MASK_1(p2) + INTERP_32_MASK_1(p3)) / 16)
rlm@1:     | INTERP_32_MASK_2((INTERP_32_MASK_2(p1)*14 + INTERP_32_MASK_2(p2) + INTERP_32_MASK_2(p3)) / 16);
rlm@1: }
rlm@1: 
rlm@1: static inline u32 interp_32_431(u32 p1, u32 p2, u32 p3)
rlm@1: {
rlm@1:   return INTERP_32_MASK_1((INTERP_32_MASK_1(p1)*4 + INTERP_32_MASK_1(p2)*3 + INTERP_32_MASK_1(p3)) / 8)
rlm@1:     | INTERP_32_MASK_2((INTERP_32_MASK_2(p1)*4 + INTERP_32_MASK_2(p2)*3 + INTERP_32_MASK_2(p3)) / 8);
rlm@1: }
rlm@1: 
rlm@1: static inline u32 interp_32_53(u32 p1, u32 p2)
rlm@1: {
rlm@1:   return INTERP_32_MASK_1((INTERP_32_MASK_1(p1)*5 + INTERP_32_MASK_1(p2)*3) / 8)
rlm@1:     | INTERP_32_MASK_2((INTERP_32_MASK_2(p1)*5 + INTERP_32_MASK_2(p2)*3) / 8);
rlm@1: }
rlm@1: 
rlm@1: static inline u32 interp_32_151(u32 p1, u32 p2)
rlm@1: {
rlm@1:   return INTERP_32_MASK_1((INTERP_32_MASK_1(p1)*15 + INTERP_32_MASK_1(p2)) / 16)
rlm@1:     | INTERP_32_MASK_2((INTERP_32_MASK_2(p1)*15 + INTERP_32_MASK_2(p2)) / 16);
rlm@1: }
rlm@1: 
rlm@1: static inline u32 interp_32_97(u32 p1, u32 p2)
rlm@1: {
rlm@1:   return INTERP_32_MASK_1((INTERP_32_MASK_1(p1)*9 + INTERP_32_MASK_1(p2)*7) / 16)
rlm@1:     | INTERP_32_MASK_2((INTERP_32_MASK_2(p1)*9 + INTERP_32_MASK_2(p2)*7) / 16);
rlm@1: }
rlm@1: 
rlm@1: /***************************************************************************/
rlm@1: /* diff */
rlm@1: 
rlm@1: #define INTERP_Y_LIMIT (0x30*4)
rlm@1: #define INTERP_U_LIMIT (0x07*4)
rlm@1: #define INTERP_V_LIMIT (0x06*8)
rlm@1: 
rlm@1: static int interp_16_diff(u16 p1, u16 p2)
rlm@1: {
rlm@1:   int r, g, b;
rlm@1:   int y, u, v;
rlm@1: 
rlm@1:   if (p1 == p2)
rlm@1:     return 0;
rlm@1: 
rlm@1:   if (interp_bits_per_pixel == 16) {
rlm@1:     b = (int)((p1 & 0x1F) - (p2 & 0x1F)) << 3;
rlm@1:     g = (int)((p1 & 0x7E0) - (p2 & 0x7E0)) >> 3;
rlm@1:     r = (int)((p1 & 0xF800) - (p2 & 0xF800)) >> 8;
rlm@1:   } else {
rlm@1:     b = (int)((p1 & 0x1F) - (p2 & 0x1F)) << 3;
rlm@1:     g = (int)((p1 & 0x3E0) - (p2 & 0x3E0)) >> 2;
rlm@1:     r = (int)((p1 & 0x7C00) - (p2 & 0x7C00)) >> 7;
rlm@1:   }
rlm@1: 
rlm@1:   y = r + g + b;
rlm@1:   u = r - b;
rlm@1:   v = -r + 2*g - b;
rlm@1: 
rlm@1:   if (y < -INTERP_Y_LIMIT || y > INTERP_Y_LIMIT)
rlm@1:     return 1;
rlm@1: 
rlm@1:   if (u < -INTERP_U_LIMIT || u > INTERP_U_LIMIT)
rlm@1:     return 1;
rlm@1: 
rlm@1:   if (v < -INTERP_V_LIMIT || v > INTERP_V_LIMIT)
rlm@1:     return 1;
rlm@1: 
rlm@1: return 0;
rlm@1: }
rlm@1: 
rlm@1: static int interp_32_diff(u32 p1, u32 p2)
rlm@1: {
rlm@1:   int r, g, b;
rlm@1:   int y, u, v;
rlm@1: 
rlm@1:   if ((p1 & 0xF8F8F8) == (p2 & 0xF8F8F8))
rlm@1:     return 0;
rlm@1: 
rlm@1:   b = (int)((p1 & 0xFF) - (p2 & 0xFF));
rlm@1:   g = (int)((p1 & 0xFF00) - (p2 & 0xFF00)) >> 8;
rlm@1:   r = (int)((p1 & 0xFF0000) - (p2 & 0xFF0000)) >> 16;
rlm@1: 
rlm@1:   y = r + g + b;
rlm@1:   u = r - b;
rlm@1:   v = -r + 2*g - b;
rlm@1: 
rlm@1:   if (y < -INTERP_Y_LIMIT || y > INTERP_Y_LIMIT)
rlm@1:     return 1;
rlm@1: 
rlm@1:   if (u < -INTERP_U_LIMIT || u > INTERP_U_LIMIT)
rlm@1:     return 1;
rlm@1: 
rlm@1:   if (v < -INTERP_V_LIMIT || v > INTERP_V_LIMIT)
rlm@1:     return 1;
rlm@1: 
rlm@1:   return 0;
rlm@1: }
rlm@1: 
rlm@1: 
rlm@1: #define INTERP_LIMIT2 (96000)
rlm@1: #define ABS(x) ((x) < 0 ? -(x) : (x))
rlm@1: #define MAX(x,y) ((x) > (y) ? (x) : (y))
rlm@1: #define MIN(x,y) ((x) < (y) ? (x) : (y))
rlm@1: 
rlm@1: static int interp_16_diff2(u16 p1, u16 p2)
rlm@1: {
rlm@1:   int r, g, b;
rlm@1:   int y, u, v;
rlm@1: 
rlm@1:   if ((p1 & 0xF79E) == (p2 & 0xF79E))
rlm@1:     return 0;
rlm@1: 
rlm@1:   if (interp_bits_per_pixel == 16) {
rlm@1:     b = (int)((p1 & 0x1F) - (p2 & 0x1F)) << 3;
rlm@1:     g = (int)((p1 & 0x7E0) - (p2 & 0x7E0)) >> 3;
rlm@1:     r = (int)((p1 & 0xF800) - (p2 & 0xF800)) >> 8;
rlm@1:   } else {
rlm@1:     b = (int)((p1 & 0x1F) - (p2 & 0x1F)) << 3;
rlm@1:     g = (int)((p1 & 0x3E0) - (p2 & 0x3E0)) >> 2;
rlm@1:     r = (int)((p1 & 0x7C00) - (p2 & 0x7C00)) >> 7;
rlm@1:   }
rlm@1: 
rlm@1: //  yb =  30*r + 58*g + 12*b;
rlm@1:   y =  33*r + 36*g + 31*b;
rlm@1:   u = -14*r - 29*g + 44*b; 
rlm@1:   v =  62*r - 51*g - 10*b;
rlm@1: 
rlm@1:   if (11*ABS(y) + 8*ABS(u) + 6*ABS(v) > INTERP_LIMIT2)
rlm@1:     return 1;
rlm@1:   return 0;
rlm@1: }
rlm@1: 
rlm@1: static int interp_32_diff2(u32 p1, u32 p2)
rlm@1: {
rlm@1:   int r, g, b;
rlm@1:   int y, u, v;
rlm@1: 
rlm@1:   if ((p1 & 0xF0F0F0) == (p2 & 0xF0F0F0))
rlm@1:     return 0;
rlm@1: 
rlm@1:   b = (int)((p1 & 0xF8) - (p2 & 0xF8));
rlm@1:   g = (int)((p1 & 0xF800) - (p2 & 0xF800)) >> 8;
rlm@1:   r = (int)((p1 & 0xF80000) - (p2 & 0xF80000)) >> 16;
rlm@1: 
rlm@1: //  y =  30*r + 58*g + 12*b;
rlm@1:   y =  33*r + 36*g + 31*b;
rlm@1:   u = -14*r - 29*g + 44*b; 
rlm@1:   v =  62*r - 51*g - 10*b;
rlm@1: 
rlm@1:   if (11*ABS(y) + 8*ABS(u) + 6*ABS(v) > INTERP_LIMIT2)
rlm@1:     return 1;
rlm@1: 
rlm@1:   return 0;
rlm@1: }
rlm@1: 
rlm@1: static void interp_set(unsigned bits_per_pixel)
rlm@1: {
rlm@1:   interp_bits_per_pixel = bits_per_pixel;
rlm@1: 
rlm@1:   switch (bits_per_pixel) {
rlm@1:   case 15 :
rlm@1:     interp_mask[0] = 0x7C1F;
rlm@1:     interp_mask[1] = 0x03E0;
rlm@1:     break;
rlm@1:   case 16 :
rlm@1:     interp_mask[0] = 0xF81F;
rlm@1:     interp_mask[1] = 0x07E0;
rlm@1:     break;
rlm@1:   case 32 :
rlm@1:     interp_mask[0] = 0xFF00FF;
rlm@1:     interp_mask[1] = 0x00FF00;
rlm@1:     break;
rlm@1:   }
rlm@1: }
rlm@1: 
rlm@1: #endif