Initial revision
[blender.git] / source / blender / imbuf / intern / hamx.c
1 /**
2  *
3  * ***** BEGIN GPL/BL DUAL LICENSE BLOCK *****
4  *
5  * This program is free software; you can redistribute it and/or
6  * modify it under the terms of the GNU General Public License
7  * as published by the Free Software Foundation; either version 2
8  * of the License, or (at your option) any later version. The Blender
9  * Foundation also sells licenses for use in proprietary software under
10  * the Blender License.  See http://www.blender.org/BL/ for information
11  * about this.
12  *
13  * This program is distributed in the hope that it will be useful,
14  * but WITHOUT ANY WARRANTY; without even the implied warranty of
15  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
16  * GNU General Public License for more details.
17  *
18  * You should have received a copy of the GNU General Public License
19  * along with this program; if not, write to the Free Software Foundation,
20  * Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.
21  *
22  * The Original Code is Copyright (C) 2001-2002 by NaN Holding BV.
23  * All rights reserved.
24  *
25  * The Original Code is: all of this file.
26  *
27  * Contributor(s): none yet.
28  *
29  * ***** END GPL/BL DUAL LICENSE BLOCK *****
30  * hamx.c
31  *
32  * $Id$
33  */
34
35 #include "BLI_blenlib.h"
36 #ifdef WIN32
37 #include "BLI_winstuff.h"
38 #include <io.h>
39 #endif
40
41
42 #include "imbuf.h"
43 #include "imbuf_patch.h"
44 #include "IMB_imbuf_types.h"
45 #include "IMB_imbuf.h"
46
47 #include "IMB_allocimbuf.h"
48 #include "IMB_filter.h"
49 #include "IMB_ham.h"
50 #include "IMB_hamx.h"
51
52 /* actually hard coded endianness */
53 #define GET_BIG_LONG(x) (((uchar *) (x))[0] << 24 | ((uchar *) (x))[1] << 16 | ((uchar *) (x))[2] << 8 | ((uchar *) (x))[3])
54 #define GET_LITTLE_LONG(x) (((uchar *) (x))[3] << 24 | ((uchar *) (x))[2] << 16 | ((uchar *) (x))[1] << 8 | ((uchar *) (x))[0])
55 #define SWAP_L(x) (((x << 24) & 0xff000000) | ((x << 8) & 0xff0000) | ((x >> 8) & 0xff00) | ((x >> 24) & 0xff))
56 #define SWAP_S(x) (((x << 8) & 0xff00) | ((x >> 8) & 0xff))
57
58 /* more endianness... should move to a separate file... */
59 #if defined(__sgi) || defined (__sparc) || defined (__PPC__) || defined (__ppc__) || defined (__BIG_ENDIAN__)
60 #define GET_ID GET_BIG_LONG
61 #define LITTLE_LONG SWAP_LONG
62 #else
63 #define GET_ID GET_LITTLE_LONG
64 #define LITTLE_LONG ENDIAN_NOP
65 #endif
66
67 #ifndef ABS
68 #define ABS(x)  ((x) < 0 ? -(x) : (x))
69 #endif
70
71 static uchar hamx_array_char[] = {
72         0x00,0x00,0xFF,0xFF, 0x00,0x00,0xFF,0xFF, 0x00,0x00,0xFF,0xFF, 0x00,0x00,0xFF,0xFF, 0x00,0x00,0xFF,0xFF, 0x00,0x00,0xFF,0xFF, 0x00,0x00,0xFF,0xFF, 0x00,0x00,0xFF,0xFF,
73         0x00,0x00,0xFF,0xFF, 0x00,0x00,0xFF,0xFF, 0x00,0x00,0xFF,0xFF, 0x00,0x00,0xFF,0xFF, 0x00,0x00,0xFF,0xFF, 0x00,0x00,0xFF,0xFF, 0x00,0x00,0xFF,0xFF, 0x00,0x00,0xFF,0xFF,
74         0x00,0xFF,0x00,0xFF, 0x00,0xFF,0x00,0xFF, 0x00,0xFF,0x00,0xFF, 0x00,0xFF,0x00,0xFF, 0x00,0xFF,0x00,0xFF, 0x00,0xFF,0x00,0xFF, 0x00,0xFF,0x00,0xFF, 0x00,0xFF,0x00,0xFF,
75         0x00,0xFF,0x00,0xFF, 0x00,0xFF,0x00,0xFF, 0x00,0xFF,0x00,0xFF, 0x00,0xFF,0x00,0xFF, 0x00,0xFF,0x00,0xFF, 0x00,0xFF,0x00,0xFF, 0x00,0xFF,0x00,0xFF, 0x00,0xFF,0x00,0xFF,
76         0x00,0xFF,0xFF,0x00, 0x00,0xFF,0xFF,0x00, 0x00,0xFF,0xFF,0x00, 0x00,0xFF,0xFF,0x00, 0x00,0xFF,0xFF,0x00, 0x00,0xFF,0xFF,0x00, 0x00,0xFF,0xFF,0x00, 0x00,0xFF,0xFF,0x00,
77         0x00,0xFF,0xFF,0x00, 0x00,0xFF,0xFF,0x00, 0x00,0xFF,0xFF,0x00, 0x00,0xFF,0xFF,0x00, 0x00,0xFF,0xFF,0x00, 0x00,0xFF,0xFF,0x00, 0x00,0xFF,0xFF,0x00, 0x00,0xFF,0xFF,0x00,
78         
79         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
80         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
81         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
82         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
83         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
84         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
85         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
86         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
87         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
88         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
89         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
90         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
91         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
92         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
93         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
94         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
95         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
96         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
97         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
98         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
99         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
100         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
101         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
102         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
103         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
104         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x00,
105         
106         0x00,0x00,0x00,0x00, 0x00,0x10,0x00,0x00, 0x00,0x20,0x00,0x00, 0x00,0x30,0x00,0x00, 0x00,0x40,0x00,0x00, 0x00,0x50,0x00,0x00, 0x00,0x60,0x00,0x00, 0x00,0x70,0x00,0x00,
107         0x00,0x80,0x00,0x00, 0x00,0x90,0x00,0x00, 0x00,0xA0,0x00,0x00, 0x00,0xB0,0x00,0x00, 0x00,0xC0,0x00,0x00, 0x00,0xD0,0x00,0x00, 0x00,0xE0,0x00,0x00, 0x00,0xF0,0x00,0x00,
108         0x00,0x00,0x00,0x00, 0x00,0x00,0x10,0x00, 0x00,0x00,0x20,0x00, 0x00,0x00,0x30,0x00, 0x00,0x00,0x40,0x00, 0x00,0x00,0x50,0x00, 0x00,0x00,0x60,0x00, 0x00,0x00,0x70,0x00,
109         0x00,0x00,0x80,0x00, 0x00,0x00,0x90,0x00, 0x00,0x00,0xA0,0x00, 0x00,0x00,0xB0,0x00, 0x00,0x00,0xC0,0x00, 0x00,0x00,0xD0,0x00, 0x00,0x00,0xE0,0x00, 0x00,0x00,0xF0,0x00,
110         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x10, 0x00,0x00,0x00,0x20, 0x00,0x00,0x00,0x30, 0x00,0x00,0x00,0x40, 0x00,0x00,0x00,0x50, 0x00,0x00,0x00,0x60, 0x00,0x00,0x00,0x70,
111         0x00,0x00,0x00,0x80, 0x00,0x00,0x00,0x90, 0x00,0x00,0x00,0xA0, 0x00,0x00,0x00,0xB0, 0x00,0x00,0x00,0xC0, 0x00,0x00,0x00,0xD0, 0x00,0x00,0x00,0xE0, 0x00,0x00,0x00,0xF0,
112         
113         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x30, 0x00,0x00,0x00,0x60, 0x00,0x00,0x00,0x90, 0x00,0x00,0x00,0xC0, 0x00,0x00,0x00,0xF0,
114         0x00,0x00,0x20,0x00, 0x00,0x00,0x20,0x30, 0x00,0x00,0x20,0x60, 0x00,0x00,0x20,0x90, 0x00,0x00,0x20,0xC0, 0x00,0x00,0x20,0xF0,
115         0x00,0x00,0x40,0x00, 0x00,0x00,0x40,0x30, 0x00,0x00,0x40,0x60, 0x00,0x00,0x40,0x90, 0x00,0x00,0x40,0xC0, 0x00,0x00,0x40,0xF0,
116         0x00,0x00,0x60,0x00, 0x00,0x00,0x60,0x30, 0x00,0x00,0x60,0x60, 0x00,0x00,0x60,0x90, 0x00,0x00,0x60,0xC0, 0x00,0x00,0x60,0xF0,
117         0x00,0x00,0x90,0x00, 0x00,0x00,0x90,0x30, 0x00,0x00,0x90,0x60, 0x00,0x00,0x90,0x90, 0x00,0x00,0x90,0xC0, 0x00,0x00,0x90,0xF0,
118         0x00,0x00,0xB0,0x00, 0x00,0x00,0xB0,0x30, 0x00,0x00,0xB0,0x60, 0x00,0x00,0xB0,0x90, 0x00,0x00,0xB0,0xC0, 0x00,0x00,0xB0,0xF0,
119         0x00,0x00,0xD0,0x00, 0x00,0x00,0xD0,0x30, 0x00,0x00,0xD0,0x60, 0x00,0x00,0xD0,0x90, 0x00,0x00,0xD0,0xC0, 0x00,0x00,0xD0,0xF0,
120         0x00,0x00,0xF0,0x00, 0x00,0x00,0xF0,0x30, 0x00,0x00,0xF0,0x60, 0x00,0x00,0xF0,0x90, 0x00,0x00,0xF0,0xC0, 0x00,0x00,0xF0,0xF0,
121         0x00,0x50,0x00,0x00, 0x00,0x50,0x00,0x30, 0x00,0x50,0x00,0x60, 0x00,0x50,0x00,0x90, 0x00,0x50,0x00,0xC0, 0x00,0x50,0x00,0xF0,
122         0x00,0x50,0x20,0x00, 0x00,0x50,0x20,0x30, 0x00,0x50,0x20,0x60, 0x00,0x50,0x20,0x90, 0x00,0x50,0x20,0xC0, 0x00,0x50,0x20,0xF0,
123         0x00,0x50,0x40,0x00, 0x00,0x50,0x40,0x30, 0x00,0x50,0x40,0x60, 0x00,0x50,0x40,0x90, 0x00,0x50,0x40,0xC0, 0x00,0x50,0x40,0xF0,
124         0x00,0x50,0x60,0x00, 0x00,0x50,0x60,0x30, 0x00,0x50,0x60,0x60, 0x00,0x50,0x60,0x90, 0x00,0x50,0x60,0xC0, 0x00,0x50,0x60,0xF0,
125         0x00,0x50,0x90,0x00, 0x00,0x50,0x90,0x30, 0x00,0x50,0x90,0x60, 0x00,0x50,0x90,0x90, 0x00,0x50,0x90,0xC0, 0x00,0x50,0x90,0xF0,
126         0x00,0x50,0xB0,0x00, 0x00,0x50,0xB0,0x30, 0x00,0x50,0xB0,0x60, 0x00,0x50,0xB0,0x90, 0x00,0x50,0xB0,0xC0, 0x00,0x50,0xB0,0xF0,
127         0x00,0x50,0xD0,0x00, 0x00,0x50,0xD0,0x30, 0x00,0x50,0xD0,0x60, 0x00,0x50,0xD0,0x90, 0x00,0x50,0xD0,0xC0, 0x00,0x50,0xD0,0xF0,
128         0x00,0x50,0xF0,0x00, 0x00,0x50,0xF0,0x30, 0x00,0x50,0xF0,0x60, 0x00,0x50,0xF0,0x90, 0x00,0x50,0xF0,0xC0, 0x00,0x50,0xF0,0xF0,
129         0x00,0xA0,0x00,0x00, 0x00,0xA0,0x00,0x30, 0x00,0xA0,0x00,0x60, 0x00,0xA0,0x00,0x90, 0x00,0xA0,0x00,0xC0, 0x00,0xA0,0x00,0xF0,
130         0x00,0xA0,0x20,0x00, 0x00,0xA0,0x20,0x30, 0x00,0xA0,0x20,0x60, 0x00,0xA0,0x20,0x90, 0x00,0xA0,0x20,0xC0, 0x00,0xA0,0x20,0xF0,
131         0x00,0xA0,0x40,0x00, 0x00,0xA0,0x40,0x30, 0x00,0xA0,0x40,0x60, 0x00,0xA0,0x40,0x90, 0x00,0xA0,0x40,0xC0, 0x00,0xA0,0x40,0xF0,
132         0x00,0xA0,0x60,0x00, 0x00,0xA0,0x60,0x30, 0x00,0xA0,0x60,0x60, 0x00,0xA0,0x60,0x90, 0x00,0xA0,0x60,0xC0, 0x00,0xA0,0x60,0xF0,
133         0x00,0xA0,0x90,0x00, 0x00,0xA0,0x90,0x30, 0x00,0xA0,0x90,0x60, 0x00,0xA0,0x90,0x90, 0x00,0xA0,0x90,0xC0, 0x00,0xA0,0x90,0xF0,
134         0x00,0xA0,0xB0,0x00, 0x00,0xA0,0xB0,0x30, 0x00,0xA0,0xB0,0x60, 0x00,0xA0,0xB0,0x90, 0x00,0xA0,0xB0,0xC0, 0x00,0xA0,0xB0,0xF0,
135         0x00,0xA0,0xD0,0x00, 0x00,0xA0,0xD0,0x30, 0x00,0xA0,0xD0,0x60, 0x00,0xA0,0xD0,0x90, 0x00,0xA0,0xD0,0xC0, 0x00,0xA0,0xD0,0xF0,
136         0x00,0xA0,0xF0,0x00, 0x00,0xA0,0xF0,0x30, 0x00,0xA0,0xF0,0x60, 0x00,0xA0,0xF0,0x90, 0x00,0xA0,0xF0,0xC0, 0x00,0xA0,0xF0,0xF0,
137         0x00,0xF0,0x00,0x00, 0x00,0xF0,0x00,0x30, 0x00,0xF0,0x00,0x60, 0x00,0xF0,0x00,0x90, 0x00,0xF0,0x00,0xC0, 0x00,0xF0,0x00,0xF0,
138         0x00,0xF0,0x20,0x00, 0x00,0xF0,0x20,0x30, 0x00,0xF0,0x20,0x60, 0x00,0xF0,0x20,0x90, 0x00,0xF0,0x20,0xC0, 0x00,0xF0,0x20,0xF0,
139         0x00,0xF0,0x40,0x00, 0x00,0xF0,0x40,0x30, 0x00,0xF0,0x40,0x60, 0x00,0xF0,0x40,0x90, 0x00,0xF0,0x40,0xC0, 0x00,0xF0,0x40,0xF0,
140         0x00,0xF0,0x60,0x00, 0x00,0xF0,0x60,0x30, 0x00,0xF0,0x60,0x60, 0x00,0xF0,0x60,0x90, 0x00,0xF0,0x60,0xC0, 0x00,0xF0,0x60,0xF0,
141         0x00,0xF0,0x90,0x00, 0x00,0xF0,0x90,0x30, 0x00,0xF0,0x90,0x60, 0x00,0xF0,0x90,0x90, 0x00,0xF0,0x90,0xC0, 0x00,0xF0,0x90,0xF0,
142         0x00,0xF0,0xB0,0x00, 0x00,0xF0,0xB0,0x30, 0x00,0xF0,0xB0,0x60, 0x00,0xF0,0xB0,0x90, 0x00,0xF0,0xB0,0xC0, 0x00,0xF0,0xB0,0xF0,
143         0x00,0xF0,0xD0,0x00, 0x00,0xF0,0xD0,0x30, 0x00,0xF0,0xD0,0x60, 0x00,0xF0,0xD0,0x90, 0x00,0xF0,0xD0,0xC0, 0x00,0xF0,0xD0,0xF0,
144         0x00,0xF0,0xF0,0x00, 0x00,0xF0,0xF0,0x30, 0x00,0xF0,0xF0,0x60, 0x00,0xF0,0xF0,0x90, 0x00,0xF0,0xF0,0xC0, 0x00,0xF0,0xF0,0xF0,
145         
146         0x00,0x10,0x10,0x10, 0x00,0x20,0x20,0x20, 0x00,0x30,0x30,0x30, 0x00,0x40,0x40,0x40,
147         0x00,0x50,0x50,0x50, 0x00,0x60,0x60,0x60, 0x00,0x70,0x70,0x70, 0x00,0x80,0x80,0x80,
148         0x00,0x90,0x90,0x90, 0x00,0xA0,0xA0,0xA0, 0x00,0xB0,0xB0,0xB0, 0x00,0xC0,0xC0,0xC0,
149         0x00,0xD0,0xD0,0xD0, 0x00,0xE0,0xE0,0xE0
150 };
151
152 static int * hamx_array = (int *) (hamx_array_char);
153
154 static uchar cmap_hamx[] = {
155         0x00,0x00,0x00,0x00, 0x00,0x00,0x00,0x30, 0x00,0x00,0x00,0x60, 0x00,0x00,0x00,0x90, 0x00,0x00,0x00,0xC0, 0x00,0x00,0x00,0xF0,
156         0x00,0x00,0x20,0x00, 0x00,0x00,0x20,0x30, 0x00,0x00,0x20,0x60, 0x00,0x00,0x20,0x90, 0x00,0x00,0x20,0xC0, 0x00,0x00,0x20,0xF0,
157         0x00,0x00,0x40,0x00, 0x00,0x00,0x40,0x30, 0x00,0x00,0x40,0x60, 0x00,0x00,0x40,0x90, 0x00,0x00,0x40,0xC0, 0x00,0x00,0x40,0xF0,
158         0x00,0x00,0x60,0x00, 0x00,0x00,0x60,0x30, 0x00,0x00,0x60,0x60, 0x00,0x00,0x60,0x90, 0x00,0x00,0x60,0xC0, 0x00,0x00,0x60,0xF0,
159         0x00,0x00,0x90,0x00, 0x00,0x00,0x90,0x30, 0x00,0x00,0x90,0x60, 0x00,0x00,0x90,0x90, 0x00,0x00,0x90,0xC0, 0x00,0x00,0x90,0xF0,
160         0x00,0x00,0xB0,0x00, 0x00,0x00,0xB0,0x30, 0x00,0x00,0xB0,0x60, 0x00,0x00,0xB0,0x90, 0x00,0x00,0xB0,0xC0, 0x00,0x00,0xB0,0xF0,
161         0x00,0x00,0xD0,0x00, 0x00,0x00,0xD0,0x30, 0x00,0x00,0xD0,0x60, 0x00,0x00,0xD0,0x90, 0x00,0x00,0xD0,0xC0, 0x00,0x00,0xD0,0xF0,
162         0x00,0x00,0xF0,0x00, 0x00,0x00,0xF0,0x30, 0x00,0x00,0xF0,0x60, 0x00,0x00,0xF0,0x90, 0x00,0x00,0xF0,0xC0, 0x00,0x00,0xF0,0xF0,
163         0x00,0x50,0x00,0x00, 0x00,0x50,0x00,0x30, 0x00,0x50,0x00,0x60, 0x00,0x50,0x00,0x90, 0x00,0x50,0x00,0xC0, 0x00,0x50,0x00,0xF0,
164         0x00,0x50,0x20,0x00, 0x00,0x50,0x20,0x30, 0x00,0x50,0x20,0x60, 0x00,0x50,0x20,0x90, 0x00,0x50,0x20,0xC0, 0x00,0x50,0x20,0xF0,
165         0x00,0x50,0x40,0x00, 0x00,0x50,0x40,0x30, 0x00,0x50,0x40,0x60, 0x00,0x50,0x40,0x90, 0x00,0x50,0x40,0xC0, 0x00,0x50,0x40,0xF0,
166         0x00,0x50,0x60,0x00, 0x00,0x50,0x60,0x30, 0x00,0x50,0x60,0x60, 0x00,0x50,0x60,0x90, 0x00,0x50,0x60,0xC0, 0x00,0x50,0x60,0xF0,
167         0x00,0x50,0x90,0x00, 0x00,0x50,0x90,0x30, 0x00,0x50,0x90,0x60, 0x00,0x50,0x90,0x90, 0x00,0x50,0x90,0xC0, 0x00,0x50,0x90,0xF0,
168         0x00,0x50,0xB0,0x00, 0x00,0x50,0xB0,0x30, 0x00,0x50,0xB0,0x60, 0x00,0x50,0xB0,0x90, 0x00,0x50,0xB0,0xC0, 0x00,0x50,0xB0,0xF0,
169         0x00,0x50,0xD0,0x00, 0x00,0x50,0xD0,0x30, 0x00,0x50,0xD0,0x60, 0x00,0x50,0xD0,0x90, 0x00,0x50,0xD0,0xC0, 0x00,0x50,0xD0,0xF0,
170         0x00,0x50,0xF0,0x00, 0x00,0x50,0xF0,0x30, 0x00,0x50,0xF0,0x60, 0x00,0x50,0xF0,0x90, 0x00,0x50,0xF0,0xC0, 0x00,0x50,0xF0,0xF0,
171         0x00,0xA0,0x00,0x00, 0x00,0xA0,0x00,0x30, 0x00,0xA0,0x00,0x60, 0x00,0xA0,0x00,0x90, 0x00,0xA0,0x00,0xC0, 0x00,0xA0,0x00,0xF0,
172         0x00,0xA0,0x20,0x00, 0x00,0xA0,0x20,0x30, 0x00,0xA0,0x20,0x60, 0x00,0xA0,0x20,0x90, 0x00,0xA0,0x20,0xC0, 0x00,0xA0,0x20,0xF0,
173         0x00,0xA0,0x40,0x00, 0x00,0xA0,0x40,0x30, 0x00,0xA0,0x40,0x60, 0x00,0xA0,0x40,0x90, 0x00,0xA0,0x40,0xC0, 0x00,0xA0,0x40,0xF0,
174         0x00,0xA0,0x60,0x00, 0x00,0xA0,0x60,0x30, 0x00,0xA0,0x60,0x60, 0x00,0xA0,0x60,0x90, 0x00,0xA0,0x60,0xC0, 0x00,0xA0,0x60,0xF0,
175         0x00,0xA0,0x90,0x00, 0x00,0xA0,0x90,0x30, 0x00,0xA0,0x90,0x60, 0x00,0xA0,0x90,0x90, 0x00,0xA0,0x90,0xC0, 0x00,0xA0,0x90,0xF0,
176         0x00,0xA0,0xB0,0x00, 0x00,0xA0,0xB0,0x30, 0x00,0xA0,0xB0,0x60, 0x00,0xA0,0xB0,0x90, 0x00,0xA0,0xB0,0xC0, 0x00,0xA0,0xB0,0xF0,
177         0x00,0xA0,0xD0,0x00, 0x00,0xA0,0xD0,0x30, 0x00,0xA0,0xD0,0x60, 0x00,0xA0,0xD0,0x90, 0x00,0xA0,0xD0,0xC0, 0x00,0xA0,0xD0,0xF0,
178         0x00,0xA0,0xF0,0x00, 0x00,0xA0,0xF0,0x30, 0x00,0xA0,0xF0,0x60, 0x00,0xA0,0xF0,0x90, 0x00,0xA0,0xF0,0xC0, 0x00,0xA0,0xF0,0xF0,
179         0x00,0xF0,0x00,0x00, 0x00,0xF0,0x00,0x30, 0x00,0xF0,0x00,0x60, 0x00,0xF0,0x00,0x90, 0x00,0xF0,0x00,0xC0, 0x00,0xF0,0x00,0xF0,
180         0x00,0xF0,0x20,0x00, 0x00,0xF0,0x20,0x30, 0x00,0xF0,0x20,0x60, 0x00,0xF0,0x20,0x90, 0x00,0xF0,0x20,0xC0, 0x00,0xF0,0x20,0xF0,
181         0x00,0xF0,0x40,0x00, 0x00,0xF0,0x40,0x30, 0x00,0xF0,0x40,0x60, 0x00,0xF0,0x40,0x90, 0x00,0xF0,0x40,0xC0, 0x00,0xF0,0x40,0xF0,
182         0x00,0xF0,0x60,0x00, 0x00,0xF0,0x60,0x30, 0x00,0xF0,0x60,0x60, 0x00,0xF0,0x60,0x90, 0x00,0xF0,0x60,0xC0, 0x00,0xF0,0x60,0xF0,
183         0x00,0xF0,0x90,0x00, 0x00,0xF0,0x90,0x30, 0x00,0xF0,0x90,0x60, 0x00,0xF0,0x90,0x90, 0x00,0xF0,0x90,0xC0, 0x00,0xF0,0x90,0xF0,
184         0x00,0xF0,0xB0,0x00, 0x00,0xF0,0xB0,0x30, 0x00,0xF0,0xB0,0x60, 0x00,0xF0,0xB0,0x90, 0x00,0xF0,0xB0,0xC0, 0x00,0xF0,0xB0,0xF0,
185         0x00,0xF0,0xD0,0x00, 0x00,0xF0,0xD0,0x30, 0x00,0xF0,0xD0,0x60, 0x00,0xF0,0xD0,0x90, 0x00,0xF0,0xD0,0xC0, 0x00,0xF0,0xD0,0xF0,
186         0x00,0xF0,0xF0,0x00, 0x00,0xF0,0xF0,0x30, 0x00,0xF0,0xF0,0x60, 0x00,0xF0,0xF0,0x90, 0x00,0xF0,0xF0,0xC0, 0x00,0xF0,0xF0,0xF0,
187         0x00,0x10,0x10,0x10, 0x00,0x20,0x20,0x20, 0x00,0x30,0x30,0x30, 0x00,0x40,0x40,0x40, 0x00,0x50,0x50,0x50, 0x00,0x60,0x60,0x60,
188         0x00,0x70,0x70,0x70, 0x00,0x80,0x80,0x80, 0x00,0x90,0x90,0x90, 0x00,0xA0,0xA0,0xA0, 0x00,0xB0,0xB0,0xB0, 0x00,0xC0,0xC0,0xC0,
189         0x00,0xD0,0xD0,0xD0, 0x00,0xE0,0xE0,0xE0
190 };
191
192
193 float adat_gamma = 1.0;
194 float adat_distort = 1.0;
195
196 /*
197  * 
198  * Nieuwe versie:
199  * 
200  * 32 helderheden       Y                        15 direct te bereiken (zwart & wit zijn specials)
201  * 16 kleuren           H ue            
202  * 7 intensiteiten      S aturation
203  * 
204  * Totaal 3584 'verschillende' kleuren. Eerste 512 kleuren vrij.
205  * 
206  * 
207  */
208
209 void imb_convhamx(struct ImBuf *ibuf, unsigned char coltab[][4], short *deltab)
210 /*  struct ImBuf *ibuf; */
211 /*  uchar coltab[][4]; */
212 /*  short *deltab; */
213 {
214         short r,g,b,lr,lg,lb,dr,dg,db,col,fout,type,step;
215         int i;
216         uchar *rect;
217
218         /*
219                 b       = 0000 xxxx
220                 g       = 0001 xxxx
221                 r       = 0010 xxxx
222                 cmap  >=  48
223         */
224
225         for (step = 0 ; step < 2 ; step ++){
226                 rect =  (uchar *) ibuf->rect;
227                 rect += 4*step;
228                 i = ((ibuf->x * ibuf->y) + 2 - step - 1) / 2;
229
230                 lb = coltab[0][1];
231                 lg = coltab[0][2];
232                 lr = coltab[0][3];
233                 type = col = 0;
234
235                 for ( ;i>0;i--){
236                         b = rect[2] >> 4;
237                         g = rect[1] >> 4;
238                         r = rect[0] >> 4;
239
240                         if ((b-lb) | (g-lg) | (r-lr)){
241                                 col = ((b<<8) + (g<<4) + r) << 1;
242                                 fout = deltab[col + 1];
243                                 col = deltab[col];
244                                 type = 0;
245                                 dr = quadr[lr-r] ;
246                                 dg = quadr[lg-g] ;
247                                 db = quadr[lb-b];
248
249                                 if ((dr+dg)<=fout) {
250                                         fout = dr+dg ;
251                                         type = 1;
252                                 }
253                                 if ((dg+db)<=fout) {
254                                         fout = dg+db;
255                                         type = 2;
256                                 }
257                                 if ((dr+db)<=fout) {
258                                         fout = dr+db;
259                                         type = 4;
260                                 }
261
262                                 switch(type){
263                                 case 1:
264                                         lb = b ;
265                                         col = b;
266                                         break;
267                                 case 4:
268                                         lg = g ;
269                                         col = g+16;
270                                         break;
271                                 case 2:
272                                         lr = r ;
273                                         col = r + 32;
274                                         break;
275                                 default:
276                                         /*printf("%04x %5d %5d  ", (b<<8) + (g<<4) + r, col, fout);*/
277                                         
278                                         lb = coltab[col][1];
279                                         lg = coltab[col][2];
280                                         lr = coltab[col][3];
281                                         /*printf("%01x%01x%01x %01x%01x%01x\n", b, g, r, lb, lg, lr);*/
282                                         col += 48;
283                                 }
284                         }
285                         rect[3] = col;
286                         rect += 8;
287                 }
288         }
289 }
290
291 static short dec_hamx(struct ImBuf * ibuf, unsigned char *body, int cmap[])
292 /*  struct ImBuf * ibuf; */
293 /*  uchar *body; */
294 /*  int cmap[]; */
295 {
296         int todo,i;
297         int j,step,col;
298         unsigned int *rect;
299
300         for (step = 0 ; step < 2 ; step ++){
301                 rect = ibuf->rect;
302                 rect += step;
303                 todo = (ibuf->x * ibuf->y + 2 - step - 1) / 2;
304                 col = cmap[0];
305                 while (todo>0){
306                         i = *body++;
307
308                         if (i & 128){                   /* fill */
309
310                                 i = 257-i;
311                                 todo -= i;
312                                 j = *(body++);
313
314                                 col = ((col & hamx_array[j]) | hamx_array[j + 256]);
315
316                                 do{
317                                         *rect = col;
318                                         rect += 2;
319                                 }while (--i);
320                         } else{                         /* copy */
321                                 i++;
322                                 todo-=i;
323
324                                 do{
325                                         j = *(body++);
326                                         *rect = col = ((col & hamx_array[j]) | hamx_array[j + 256]);
327                                         rect += 2;
328                                 }while (--i);
329                         }
330                 }
331                 if (todo) return (0);
332         }
333         return(1);
334 }
335
336
337 struct ImBuf *imb_loadanim(int *iffmem, int flags)
338 {
339         int chunk, totlen, len, *mem, cmaplen = 0;
340         unsigned int *cmap;
341         uchar *body = 0;
342         struct Adat adat;
343         struct ImBuf *ibuf=0;
344         static int is_flipped = FALSE;
345         
346         mem=iffmem;
347         if (GET_ID(mem) != FORM) return (0);
348         if (GET_ID(mem + 2) != ANIM) return (0);
349         totlen= (GET_BIG_LONG(mem + 1) + 1) & ~1;
350         mem += 3;
351         totlen -= 4;
352         adat.w = 0;
353         adat.xorig = 0;
354         adat.yorig = 0;
355         adat.gamma = adat_gamma;
356         adat.distort = adat_distort;
357         
358         while(totlen > 0){
359                 chunk = GET_ID(mem);
360                 len = (GET_BIG_LONG(mem + 1) + 1) & ~1;
361                 mem += 2;
362
363                 totlen -= len+8;
364                 switch (chunk){
365                 case ADAT:
366                         if (len > sizeof(struct Adat)){
367                                 memcpy(&adat,mem,sizeof(struct Adat));
368                         } else{
369                                 memcpy(&adat,mem,len);
370                         }
371                         adat.w = BIG_SHORT(adat.w);
372                         adat.h = BIG_SHORT(adat.h);
373                         adat.type = BIG_SHORT(adat.type);
374                         adat.xorig = BIG_SHORT(adat.xorig);
375                         adat.yorig = BIG_SHORT(adat.yorig);
376                         break;
377                 case CMAP:
378                         cmap = (unsigned int *) mem;
379                         cmaplen = len;
380                         break;
381                 case BODY:
382                         body = (uchar *) mem;
383                         break;
384                 }
385                 mem = (int *)((uchar *)mem +len);
386         }
387
388         if (body == 0) return (0);
389         if (adat.w == 0) return (0);
390
391         adat_gamma = adat.gamma;
392         adat_distort = adat.distort;
393
394         if (flags & IB_test) ibuf=IMB_allocImBuf(adat.w, adat.h, 24, 0, 0);
395         else ibuf=IMB_allocImBuf(adat.w, adat.h, 24, IB_rect, 0);
396         if (ibuf==0) return (0);
397
398         ibuf->ftype = (Anim | adat.type);
399         ibuf->xorig = adat.xorig;
400         ibuf->yorig = adat.yorig;
401         ibuf->flags = flags;
402         
403         if (cmaplen){
404                 ibuf->cmap = malloc(cmaplen);
405                 memcpy(ibuf->cmap, cmap, cmaplen);
406                 ibuf->maxcol = cmaplen >> 2;
407         }
408
409         if (flags & IB_test){
410                 if (flags & IB_freem) free(iffmem);
411                 return(ibuf);
412         }
413
414         switch (adat.type){
415         case HAMX:
416                 if (flags & IB_rect){
417                         if (!is_flipped) {
418                                 int i;
419                                 unsigned int * t;
420                                 t = (unsigned int *) hamx_array_char;
421                                 for (i = 0; i < sizeof(hamx_array_char) / sizeof(int) ; i++) {
422                                         t[i] = SWAP_LONG(t[i]);
423                                 }
424                                 
425                                 t = (unsigned int *) cmap_hamx;
426                         
427                                 for (i = 0; i < sizeof(cmap_hamx) / sizeof(int) ; i++) {
428                                         t[i] = SWAP_LONG(t[i]);
429                                 }
430                                 
431                                 is_flipped= TRUE;
432                         }
433
434                         if (dec_hamx(ibuf,body,(int*) cmap_hamx) == 0){
435                                 IMB_freeImBuf(ibuf);
436                                 ibuf = 0;
437                         }
438                         if (flags & IB_ttob) IMB_flipy(ibuf);
439                 }
440                 break;
441         default:
442                 IMB_freeImBuf(ibuf);
443                 ibuf = 0;
444         }
445         
446         if (flags & IB_freem) free(iffmem);
447
448         return (ibuf);
449 }
450
451
452 static unsigned char *makebody_anim(int bytes,
453                                                          unsigned char *buf,
454                                                          unsigned char *rect)
455 /*  register uchar *buf; */
456 /*  register uchar *rect; */
457 /*  int bytes; */
458 {
459         register uchar last,this;
460         register int copy;
461         register uchar *rectstart,*temp;
462
463         bytes--;
464         rectstart = rect;
465         last = *rect++;
466         this = *rect++;
467         copy = last^this;
468         while (bytes>0){
469                 if (copy){
470                         do{
471                                 last = this;
472                                 this = *rect++;
473                                 if (last == this){
474                                         if (this == rect[-3]){          /* drie dezelfde? */
475                                                 bytes --;                                       /* bytes goed zetten */
476                                                 break;
477                                         }
478                                 }
479                         }while (--bytes != 0);
480
481                         copy = rect-rectstart;
482                         copy --;
483                         if (bytes) copy -= 2;
484
485                         temp = rect;
486                         rect = rectstart;
487
488                         while (copy){
489                                 last = copy;
490                                 if (copy>MAXDAT) last = MAXDAT;
491                                 copy -= last;
492                                 *buf++ = last-1;
493                                 do{
494                                         *buf++ = *rect++;
495                                 }while(--last != 0);
496                         }
497                         rectstart = rect;
498                         rect = temp;
499                         last = this;
500
501                         copy = FALSE;
502                 } else {
503                         while (*rect++ == this){                /* zoek naar eerste afwijkende byte */
504                                 if (--bytes == 0) break;        /* of einde regel */
505                         }
506                         rect --;
507                         copy = rect-rectstart;
508                         rectstart = rect;
509                         bytes --;
510                         this = *rect++;
511
512                         while (copy){
513                                 if (copy>MAXRUN){
514                                         *buf++ = -(MAXRUN-1);
515                                         *buf++ = last;
516                                         copy -= MAXRUN;
517                                 } else {
518                                         *buf++ = -(copy-1);
519                                         *buf++ = last;
520                                         break;
521                                 }
522                         }
523                         copy=TRUE;
524                 }
525         }
526         return (buf);
527 }
528
529
530 short imb_enc_anim(struct ImBuf *ibuf, int file)
531 /*  struct ImBuf *ibuf; */
532 /*  int file; */
533 {
534         int step, steps, size, i, skip;
535         uchar *buf1, *crect, *_buf1, *_buf2, *bufend;
536         short ok = TRUE;
537
538         if (ibuf == 0) return (0);
539         if (file < 0 ) return (0);
540         if (ibuf->rect == 0) return(0);
541
542         /* dither toevoegen */
543
544         switch(ibuf->ftype){
545         case AN_hamx:
546                 ibuf->cmap = (unsigned int *) cmap_hamx;
547                 ibuf->mincol = 0;
548                 ibuf->maxcol = sizeof(cmap_hamx) / 4;
549                 imb_converttoham(ibuf);
550                 steps = 2;
551                 break;
552         }
553
554         size = ((ibuf->x + 1)* (ibuf->y + 1)) / steps + 1024;
555         if ((_buf1  = malloc(size)) == 0) return(0);
556         if ((_buf2  = malloc(size)) == 0){
557                 free(_buf1);
558                 return(0);
559         }
560
561         skip = 4 * steps;
562         for (step = 0 ; step < steps ; step ++){
563                 crect = (uchar *) ibuf->rect;
564                 crect += 4 * step;
565                 size = (ibuf->x * ibuf->y + steps - step - 1) / steps;
566                 buf1 = _buf1;
567                 if ((ibuf->ftype == AN_hamx) || (ibuf->ftype == AN_yuvx)){
568                         crect += 3;
569                         for (i = size ; i>0 ; i--){
570                                 *(buf1 ++) = *crect;
571                                 crect += skip;
572                         }
573                 } else{
574                         for (i = size ; i>0 ; i--){
575                                 *(buf1 ++) = crect[1] + (crect[2] >> 2) + (crect[3] >> 5);
576                                 crect += skip;
577                         }
578                 }
579                 bufend = makebody_anim(size,_buf2,_buf1);
580                 if (bufend == 0){
581                         ok = FALSE;
582                         break;
583                 }
584                 size = bufend - _buf2;
585                 if (write(file, _buf2, size) != size){
586                         ok = FALSE;
587                         break;
588                 }
589         }
590         free(_buf1);
591         free(_buf2);
592         return (ok);
593 }