~ppsspp/ppsspp/ppsspp_1.3.0

« back to all changes in this revision

Viewing changes to ext/native/ext/rg_etc1/rg_etc1.cpp

  • Committer: Sérgio Benjamim
  • Date: 2017-01-02 00:12:05 UTC
  • Revision ID: sergio_br2@yahoo.com.br-20170102001205-cxbta9za203nmjwm
1.3.0 source (from ppsspp_1.3.0-r160.p5.l1762.a165.t83~56~ubuntu16.04.1.tar.xz).

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
// File: rg_etc1.cpp - Fast, high quality ETC1 block packer/unpacker - Rich Geldreich <richgel99@gmail.com>
 
2
// Please see ZLIB license at the end of rg_etc1.h.
 
3
//
 
4
// For more information Ericsson Texture Compression (ETC/ETC1), see:
 
5
// http://www.khronos.org/registry/gles/extensions/OES/OES_compressed_ETC1_RGB8_texture.txt
 
6
//
 
7
// v1.04 - 5/15/14 - Fix signed vs. unsigned subtraction problem (noticed when compiled with gcc) in pack_etc1_block_init().
 
8
//         This issue would cause an assert when this func. was called in debug. (Note this module was developed/testing with MSVC,
 
9
//         I still need to test it throughly when compiled with gcc.)
 
10
//
 
11
// v1.03 - 5/12/13 - Initial public release
 
12
#include "rg_etc1.h"
 
13
 
 
14
#include <stdlib.h>
 
15
#include <memory.h>
 
16
#include <assert.h>
 
17
//#include <stdio.h>
 
18
#include <math.h>
 
19
 
 
20
#ifdef _MSC_VER
 
21
#pragma warning (disable: 4201) //  nonstandard extension used : nameless struct/union
 
22
#endif
 
23
 
 
24
#if defined(_DEBUG) || defined(DEBUG)
 
25
#define RG_ETC1_BUILD_DEBUG
 
26
#endif
 
27
 
 
28
#define RG_ETC1_ASSERT assert
 
29
 
 
30
namespace rg_etc1
 
31
{
 
32
   typedef unsigned char uint8;
 
33
   typedef unsigned short uint16;
 
34
   typedef unsigned int uint;
 
35
   typedef unsigned int uint32;
 
36
   typedef long long int64;
 
37
   typedef unsigned long long uint64;
 
38
 
 
39
   const uint32 cUINT32_MAX = 0xFFFFFFFFU;
 
40
   const uint64 cUINT64_MAX = 0xFFFFFFFFFFFFFFFFULL; //0xFFFFFFFFFFFFFFFFui64;
 
41
   
 
42
   template<typename T> inline T minimum(T a, T b) { return (a < b) ? a : b; }
 
43
   template<typename T> inline T minimum(T a, T b, T c) { return minimum(minimum(a, b), c); }
 
44
   template<typename T> inline T maximum(T a, T b) { return (a > b) ? a : b; }
 
45
   template<typename T> inline T maximum(T a, T b, T c) { return maximum(maximum(a, b), c); }
 
46
   template<typename T> inline T clamp(T value, T low, T high) { return (value < low) ? low : ((value > high) ? high : value); }
 
47
   template<typename T> inline T square(T value) { return value * value; }
 
48
   template<typename T> inline void zero_object(T& obj) { memset((void*)&obj, 0, sizeof(obj)); }
 
49
   template<typename T> inline void zero_this(T* pObj) { memset((void*)pObj, 0, sizeof(*pObj)); }
 
50
 
 
51
   template<class T, size_t N> T decay_array_to_subtype(T (&a)[N]);   
 
52
 
 
53
#define RG_ETC1_ARRAY_SIZE(X) (sizeof(X) / sizeof(decay_array_to_subtype(X)))
 
54
 
 
55
   enum eNoClamp { cNoClamp };
 
56
 
 
57
   struct color_quad_u8
 
58
   {
 
59
      static inline int clamp(int v) { if (v & 0xFFFFFF00U) v = (~(static_cast<int>(v) >> 31)) & 0xFF; return v; }
 
60
 
 
61
      struct component_traits { enum { cSigned = false, cFloat = false, cMin = 0U, cMax = 255U }; };
 
62
 
 
63
   public:
 
64
      typedef unsigned char component_t;
 
65
      typedef int parameter_t;
 
66
 
 
67
      enum { cNumComps = 4 };
 
68
 
 
69
      union
 
70
      {
 
71
         struct
 
72
         {
 
73
            component_t r;
 
74
            component_t g;
 
75
            component_t b;
 
76
            component_t a;
 
77
         };
 
78
 
 
79
         component_t c[cNumComps];
 
80
 
 
81
         uint32 m_u32;
 
82
      };
 
83
 
 
84
      inline color_quad_u8()
 
85
      {
 
86
      }
 
87
 
 
88
      inline color_quad_u8(const color_quad_u8& other) : m_u32(other.m_u32)
 
89
      {
 
90
      }
 
91
 
 
92
      explicit inline color_quad_u8(parameter_t y, parameter_t alpha = component_traits::cMax)
 
93
      {
 
94
         set(y, alpha);
 
95
      }
 
96
 
 
97
      inline color_quad_u8(parameter_t red, parameter_t green, parameter_t blue, parameter_t alpha = component_traits::cMax)
 
98
      {
 
99
         set(red, green, blue, alpha);
 
100
      }
 
101
 
 
102
      explicit inline color_quad_u8(eNoClamp, parameter_t y, parameter_t alpha = component_traits::cMax)
 
103
      {
 
104
         set_noclamp_y_alpha(y, alpha);
 
105
      }
 
106
 
 
107
      inline color_quad_u8(eNoClamp, parameter_t red, parameter_t green, parameter_t blue, parameter_t alpha = component_traits::cMax)
 
108
      {
 
109
         set_noclamp_rgba(red, green, blue, alpha);
 
110
      }
 
111
 
 
112
      inline void clear()
 
113
      {
 
114
         m_u32 = 0;
 
115
      }
 
116
 
 
117
      inline color_quad_u8& operator= (const color_quad_u8& other)
 
118
      {
 
119
         m_u32 = other.m_u32;
 
120
         return *this;
 
121
      }
 
122
 
 
123
      inline color_quad_u8& set_rgb(const color_quad_u8& other)
 
124
      {
 
125
         r = other.r;
 
126
         g = other.g;
 
127
         b = other.b;
 
128
         return *this;
 
129
      }
 
130
 
 
131
      inline color_quad_u8& operator= (parameter_t y)
 
132
      {
 
133
         set(y, component_traits::cMax);
 
134
         return *this;
 
135
      }
 
136
 
 
137
      inline color_quad_u8& set(parameter_t y, parameter_t alpha = component_traits::cMax)
 
138
      {
 
139
         y = clamp(y);
 
140
         alpha = clamp(alpha);
 
141
         r = static_cast<component_t>(y);
 
142
         g = static_cast<component_t>(y);
 
143
         b = static_cast<component_t>(y);
 
144
         a = static_cast<component_t>(alpha);
 
145
         return *this;
 
146
      }
 
147
 
 
148
      inline color_quad_u8& set_noclamp_y_alpha(parameter_t y, parameter_t alpha = component_traits::cMax)
 
149
      {
 
150
         RG_ETC1_ASSERT( (y >= component_traits::cMin) && (y <= component_traits::cMax) );
 
151
         RG_ETC1_ASSERT( (alpha >= component_traits::cMin) && (alpha <= component_traits::cMax) );
 
152
 
 
153
         r = static_cast<component_t>(y);
 
154
         g = static_cast<component_t>(y);
 
155
         b = static_cast<component_t>(y);
 
156
         a = static_cast<component_t>(alpha);
 
157
         return *this;
 
158
      }
 
159
 
 
160
      inline color_quad_u8& set(parameter_t red, parameter_t green, parameter_t blue, parameter_t alpha = component_traits::cMax)
 
161
      {
 
162
         r = static_cast<component_t>(clamp(red));
 
163
         g = static_cast<component_t>(clamp(green));
 
164
         b = static_cast<component_t>(clamp(blue));
 
165
         a = static_cast<component_t>(clamp(alpha));
 
166
         return *this;
 
167
      }
 
168
 
 
169
      inline color_quad_u8& set_noclamp_rgba(parameter_t red, parameter_t green, parameter_t blue, parameter_t alpha)
 
170
      {
 
171
         RG_ETC1_ASSERT( (red >= component_traits::cMin) && (red <= component_traits::cMax) );
 
172
         RG_ETC1_ASSERT( (green >= component_traits::cMin) && (green <= component_traits::cMax) );
 
173
         RG_ETC1_ASSERT( (blue >= component_traits::cMin) && (blue <= component_traits::cMax) );
 
174
         RG_ETC1_ASSERT( (alpha >= component_traits::cMin) && (alpha <= component_traits::cMax) );
 
175
 
 
176
         r = static_cast<component_t>(red);
 
177
         g = static_cast<component_t>(green);
 
178
         b = static_cast<component_t>(blue);
 
179
         a = static_cast<component_t>(alpha);
 
180
         return *this;
 
181
      }
 
182
 
 
183
      inline color_quad_u8& set_noclamp_rgb(parameter_t red, parameter_t green, parameter_t blue)
 
184
      {
 
185
         RG_ETC1_ASSERT( (red >= component_traits::cMin) && (red <= component_traits::cMax) );
 
186
         RG_ETC1_ASSERT( (green >= component_traits::cMin) && (green <= component_traits::cMax) );
 
187
         RG_ETC1_ASSERT( (blue >= component_traits::cMin) && (blue <= component_traits::cMax) );
 
188
 
 
189
         r = static_cast<component_t>(red);
 
190
         g = static_cast<component_t>(green);
 
191
         b = static_cast<component_t>(blue);
 
192
         return *this;
 
193
      }
 
194
 
 
195
      static inline parameter_t get_min_comp() { return component_traits::cMin; }
 
196
      static inline parameter_t get_max_comp() { return component_traits::cMax; }
 
197
      static inline bool get_comps_are_signed() { return component_traits::cSigned; }
 
198
 
 
199
      inline component_t operator[] (uint i) const { RG_ETC1_ASSERT(i < cNumComps); return c[i]; }
 
200
      inline component_t& operator[] (uint i) { RG_ETC1_ASSERT(i < cNumComps); return c[i]; }
 
201
 
 
202
      inline color_quad_u8& set_component(uint i, parameter_t f)
 
203
      {
 
204
         RG_ETC1_ASSERT(i < cNumComps);
 
205
 
 
206
         c[i] = static_cast<component_t>(clamp(f));
 
207
 
 
208
         return *this;
 
209
      }
 
210
 
 
211
      inline color_quad_u8& set_grayscale(parameter_t l)
 
212
      {
 
213
         component_t x = static_cast<component_t>(clamp(l));
 
214
         c[0] = x;
 
215
         c[1] = x;
 
216
         c[2] = x;
 
217
         return *this;
 
218
      }
 
219
 
 
220
      inline color_quad_u8& clamp(const color_quad_u8& l, const color_quad_u8& h)
 
221
      {
 
222
         for (uint i = 0; i < cNumComps; i++)
 
223
            c[i] = static_cast<component_t>(rg_etc1::clamp<parameter_t>(c[i], l[i], h[i]));
 
224
         return *this;
 
225
      }
 
226
 
 
227
      inline color_quad_u8& clamp(parameter_t l, parameter_t h)
 
228
      {
 
229
         for (uint i = 0; i < cNumComps; i++)
 
230
            c[i] = static_cast<component_t>(rg_etc1::clamp<parameter_t>(c[i], l, h));
 
231
         return *this;
 
232
      }
 
233
 
 
234
      // Returns CCIR 601 luma (consistent with color_utils::RGB_To_Y).
 
235
      inline parameter_t get_luma() const
 
236
      {
 
237
         return static_cast<parameter_t>((19595U * r + 38470U * g + 7471U * b + 32768U) >> 16U);
 
238
      }
 
239
 
 
240
      // Returns REC 709 luma.
 
241
      inline parameter_t get_luma_rec709() const
 
242
      {
 
243
         return static_cast<parameter_t>((13938U * r + 46869U * g + 4729U * b + 32768U) >> 16U);
 
244
      }
 
245
 
 
246
      inline uint squared_distance_rgb(const color_quad_u8& c) const
 
247
      {
 
248
         return rg_etc1::square(r - c.r) + rg_etc1::square(g - c.g) + rg_etc1::square(b - c.b);
 
249
      }
 
250
 
 
251
      inline uint squared_distance_rgba(const color_quad_u8& c) const
 
252
      {
 
253
         return rg_etc1::square(r - c.r) + rg_etc1::square(g - c.g) + rg_etc1::square(b - c.b) + rg_etc1::square(a - c.a);
 
254
      }
 
255
 
 
256
      inline bool rgb_equals(const color_quad_u8& rhs) const
 
257
      {
 
258
         return (r == rhs.r) && (g == rhs.g) && (b == rhs.b);
 
259
      }
 
260
 
 
261
      inline bool operator== (const color_quad_u8& rhs) const
 
262
      {
 
263
         return m_u32 == rhs.m_u32;
 
264
      }
 
265
 
 
266
      color_quad_u8& operator+= (const color_quad_u8& other)
 
267
      {
 
268
         for (uint i = 0; i < 4; i++)
 
269
            c[i] = static_cast<component_t>(clamp(c[i] + other.c[i]));
 
270
         return *this;
 
271
      }
 
272
 
 
273
      color_quad_u8& operator-= (const color_quad_u8& other)
 
274
      {
 
275
         for (uint i = 0; i < 4; i++)
 
276
            c[i] = static_cast<component_t>(clamp(c[i] - other.c[i]));
 
277
         return *this;
 
278
      }
 
279
 
 
280
      friend color_quad_u8 operator+ (const color_quad_u8& lhs, const color_quad_u8& rhs)
 
281
      {
 
282
         color_quad_u8 result(lhs);
 
283
         result += rhs;
 
284
         return result;
 
285
      }
 
286
 
 
287
      friend color_quad_u8 operator- (const color_quad_u8& lhs, const color_quad_u8& rhs)
 
288
      {
 
289
         color_quad_u8 result(lhs);
 
290
         result -= rhs;
 
291
         return result;
 
292
      }
 
293
   }; // class color_quad_u8
 
294
 
 
295
   struct vec3F
 
296
   {
 
297
      float m_s[3];
 
298
      
 
299
      inline vec3F() { }
 
300
      inline vec3F(float s) { m_s[0] = s; m_s[1] = s; m_s[2] = s; }
 
301
      inline vec3F(float x, float y, float z) { m_s[0] = x; m_s[1] = y; m_s[2] = z; }
 
302
      
 
303
      inline float operator[] (uint i) const { RG_ETC1_ASSERT(i < 3); return m_s[i]; }
 
304
 
 
305
      inline vec3F& operator += (const vec3F& other) { for (uint i = 0; i < 3; i++) m_s[i] += other.m_s[i]; return *this; }
 
306
 
 
307
      inline vec3F& operator *= (float s) { for (uint i = 0; i < 3; i++) m_s[i] *= s; return *this; }
 
308
   };
 
309
     
 
310
   enum etc_constants
 
311
   {
 
312
      cETC1BytesPerBlock = 8U,
 
313
 
 
314
      cETC1SelectorBits = 2U,
 
315
      cETC1SelectorValues = 1U << cETC1SelectorBits,
 
316
      cETC1SelectorMask = cETC1SelectorValues - 1U,
 
317
 
 
318
      cETC1BlockShift = 2U,
 
319
      cETC1BlockSize = 1U << cETC1BlockShift,
 
320
 
 
321
      cETC1LSBSelectorIndicesBitOffset = 0,
 
322
      cETC1MSBSelectorIndicesBitOffset = 16,
 
323
 
 
324
      cETC1FlipBitOffset = 32,
 
325
      cETC1DiffBitOffset = 33,
 
326
 
 
327
      cETC1IntenModifierNumBits = 3,
 
328
      cETC1IntenModifierValues = 1 << cETC1IntenModifierNumBits,
 
329
      cETC1RightIntenModifierTableBitOffset = 34,
 
330
      cETC1LeftIntenModifierTableBitOffset = 37,
 
331
 
 
332
      // Base+Delta encoding (5 bit bases, 3 bit delta)
 
333
      cETC1BaseColorCompNumBits = 5,
 
334
      cETC1BaseColorCompMax = 1 << cETC1BaseColorCompNumBits,
 
335
 
 
336
      cETC1DeltaColorCompNumBits = 3,
 
337
      cETC1DeltaColorComp = 1 << cETC1DeltaColorCompNumBits,
 
338
      cETC1DeltaColorCompMax = 1 << cETC1DeltaColorCompNumBits,
 
339
 
 
340
      cETC1BaseColor5RBitOffset = 59,
 
341
      cETC1BaseColor5GBitOffset = 51,
 
342
      cETC1BaseColor5BBitOffset = 43,
 
343
 
 
344
      cETC1DeltaColor3RBitOffset = 56,
 
345
      cETC1DeltaColor3GBitOffset = 48,
 
346
      cETC1DeltaColor3BBitOffset = 40,
 
347
 
 
348
      // Absolute (non-delta) encoding (two 4-bit per component bases)
 
349
      cETC1AbsColorCompNumBits = 4,
 
350
      cETC1AbsColorCompMax = 1 << cETC1AbsColorCompNumBits,
 
351
 
 
352
      cETC1AbsColor4R1BitOffset = 60,
 
353
      cETC1AbsColor4G1BitOffset = 52,
 
354
      cETC1AbsColor4B1BitOffset = 44,
 
355
 
 
356
      cETC1AbsColor4R2BitOffset = 56,
 
357
      cETC1AbsColor4G2BitOffset = 48,
 
358
      cETC1AbsColor4B2BitOffset = 40,
 
359
 
 
360
      cETC1ColorDeltaMin = -4,
 
361
      cETC1ColorDeltaMax = 3,
 
362
 
 
363
      // Delta3:
 
364
      // 0   1   2   3   4   5   6   7
 
365
      // 000 001 010 011 100 101 110 111
 
366
      // 0   1   2   3   -4  -3  -2  -1
 
367
   };
 
368
   
 
369
   static uint8 g_quant5_tab[256+16];
 
370
 
 
371
   static const int g_etc1_inten_tables[cETC1IntenModifierValues][cETC1SelectorValues] = 
 
372
   { 
 
373
      { -8,  -2,   2,   8 }, { -17,  -5,  5,  17 }, { -29,  -9,   9,  29 }, {  -42, -13, 13,  42 }, 
 
374
      { -60, -18, 18,  60 }, { -80, -24, 24,  80 }, { -106, -33, 33, 106 }, { -183, -47, 47, 183 } 
 
375
   };
 
376
 
 
377
   static const uint8 g_etc1_to_selector_index[cETC1SelectorValues] = { 2, 3, 1, 0 };
 
378
   static const uint8 g_selector_index_to_etc1[cETC1SelectorValues] = { 3, 2, 0, 1 };
 
379
      
 
380
   // Given an ETC1 diff/inten_table/selector, and an 8-bit desired color, this table encodes the best packed_color in the low byte, and the abs error in the high byte.
 
381
   static uint16 g_etc1_inverse_lookup[2*8*4][256];      // [diff/inten_table/selector][desired_color]
 
382
 
 
383
   // g_color8_to_etc_block_config[color][table_index] = Supplies for each 8-bit color value a list of packed ETC1 diff/intensity table/selectors/packed_colors that map to that color.
 
384
   // To pack: diff | (inten << 1) | (selector << 4) | (packed_c << 8)
 
385
   static const uint16 g_color8_to_etc_block_config_0_255[2][33] =
 
386
   {
 
387
      { 0x0000,  0x0010,  0x0002,  0x0012,  0x0004,  0x0014,  0x0006,  0x0016,  0x0008,  0x0018,  0x000A,  0x001A,  0x000C,  0x001C,  0x000E,  0x001E,
 
388
        0x0001,  0x0011,  0x0003,  0x0013,  0x0005,  0x0015,  0x0007,  0x0017,  0x0009,  0x0019,  0x000B,  0x001B,  0x000D,  0x001D,  0x000F,  0x001F, 0xFFFF },
 
389
      { 0x0F20,  0x0F30,  0x0E32,  0x0F22,  0x0E34,  0x0F24,  0x0D36,  0x0F26,  0x0C38,  0x0E28,  0x0B3A,  0x0E2A,  0x093C,  0x0E2C,  0x053E,  0x0D2E,
 
390
        0x1E31,  0x1F21,  0x1D33,  0x1F23,  0x1C35,  0x1E25,  0x1A37,  0x1E27,  0x1839,  0x1D29,  0x163B,  0x1C2B,  0x133D,  0x1B2D,  0x093F,  0x1A2F, 0xFFFF },
 
391
   };
 
392
 
 
393
   // Really only [254][11].
 
394
   static const uint16 g_color8_to_etc_block_config_1_to_254[254][12] = 
 
395
   {
 
396
      { 0x021C, 0x0D0D, 0xFFFF }, { 0x0020, 0x0021, 0x0A0B, 0x061F, 0xFFFF }, { 0x0113, 0x0217, 0xFFFF }, { 0x0116, 0x031E,
 
397
      0x0B0E, 0x0405, 0xFFFF }, { 0x0022, 0x0204, 0x050A, 0x0023, 0xFFFF }, { 0x0111, 0x0319, 0x0809, 0x170F, 0xFFFF }, {
 
398
      0x0303, 0x0215, 0x0607, 0xFFFF }, { 0x0030, 0x0114, 0x0408, 0x0031, 0x0201, 0x051D, 0xFFFF }, { 0x0100, 0x0024, 0x0306,
 
399
      0x0025, 0x041B, 0x0E0D, 0xFFFF }, { 0x021A, 0x0121, 0x0B0B, 0x071F, 0xFFFF }, { 0x0213, 0x0317, 0xFFFF }, { 0x0112,
 
400
      0x0505, 0xFFFF }, { 0x0026, 0x070C, 0x0123, 0x0027, 0xFFFF }, { 0x0211, 0x0909, 0xFFFF }, { 0x0110, 0x0315, 0x0707,
 
401
      0x0419, 0x180F, 0xFFFF }, { 0x0218, 0x0131, 0x0301, 0x0403, 0x061D, 0xFFFF }, { 0x0032, 0x0202, 0x0033, 0x0125, 0x051B,
 
402
      0x0F0D, 0xFFFF }, { 0x0028, 0x031C, 0x0221, 0x0029, 0xFFFF }, { 0x0120, 0x0313, 0x0C0B, 0x081F, 0xFFFF }, { 0x0605,
 
403
      0x0417, 0xFFFF }, { 0x0216, 0x041E, 0x0C0E, 0x0223, 0x0127, 0xFFFF }, { 0x0122, 0x0304, 0x060A, 0x0311, 0x0A09, 0xFFFF
 
404
      }, { 0x0519, 0x190F, 0xFFFF }, { 0x002A, 0x0231, 0x0503, 0x0415, 0x0807, 0x002B, 0x071D, 0xFFFF }, { 0x0130, 0x0214,
 
405
      0x0508, 0x0401, 0x0133, 0x0225, 0x061B, 0xFFFF }, { 0x0200, 0x0124, 0x0406, 0x0321, 0x0129, 0x100D, 0xFFFF }, { 0x031A,
 
406
      0x0D0B, 0x091F, 0xFFFF }, { 0x0413, 0x0705, 0x0517, 0xFFFF }, { 0x0212, 0x0034, 0x0323, 0x0035, 0x0227, 0xFFFF }, {
 
407
      0x0126, 0x080C, 0x0B09, 0xFFFF }, { 0x0411, 0x0619, 0x1A0F, 0xFFFF }, { 0x0210, 0x0331, 0x0603, 0x0515, 0x0907, 0x012B,
 
408
      0xFFFF }, { 0x0318, 0x002C, 0x0501, 0x0233, 0x0325, 0x071B, 0x002D, 0x081D, 0xFFFF }, { 0x0132, 0x0302, 0x0229, 0x110D,
 
409
      0xFFFF }, { 0x0128, 0x041C, 0x0421, 0x0E0B, 0x0A1F, 0xFFFF }, { 0x0220, 0x0513, 0x0617, 0xFFFF }, { 0x0135, 0x0805,
 
410
      0x0327, 0xFFFF }, { 0x0316, 0x051E, 0x0D0E, 0x0423, 0xFFFF }, { 0x0222, 0x0404, 0x070A, 0x0511, 0x0719, 0x0C09, 0x1B0F,
 
411
      0xFFFF }, { 0x0703, 0x0615, 0x0A07, 0x022B, 0xFFFF }, { 0x012A, 0x0431, 0x0601, 0x0333, 0x012D, 0x091D, 0xFFFF }, {
 
412
      0x0230, 0x0314, 0x0036, 0x0608, 0x0425, 0x0037, 0x0329, 0x081B, 0x120D, 0xFFFF }, { 0x0300, 0x0224, 0x0506, 0x0521,
 
413
      0x0F0B, 0x0B1F, 0xFFFF }, { 0x041A, 0x0613, 0x0717, 0xFFFF }, { 0x0235, 0x0905, 0xFFFF }, { 0x0312, 0x0134, 0x0523,
 
414
      0x0427, 0xFFFF }, { 0x0226, 0x090C, 0x002E, 0x0611, 0x0D09, 0x002F, 0xFFFF }, { 0x0715, 0x0B07, 0x0819, 0x032B, 0x1C0F,
 
415
      0xFFFF }, { 0x0310, 0x0531, 0x0701, 0x0803, 0x022D, 0x0A1D, 0xFFFF }, { 0x0418, 0x012C, 0x0433, 0x0525, 0x0137, 0x091B,
 
416
      0x130D, 0xFFFF }, { 0x0232, 0x0402, 0x0621, 0x0429, 0xFFFF }, { 0x0228, 0x051C, 0x0713, 0x100B, 0x0C1F, 0xFFFF }, {
 
417
      0x0320, 0x0335, 0x0A05, 0x0817, 0xFFFF }, { 0x0623, 0x0527, 0xFFFF }, { 0x0416, 0x061E, 0x0E0E, 0x0711, 0x0E09, 0x012F,
 
418
      0xFFFF }, { 0x0322, 0x0504, 0x080A, 0x0919, 0x1D0F, 0xFFFF }, { 0x0631, 0x0903, 0x0815, 0x0C07, 0x042B, 0x032D, 0x0B1D,
 
419
      0xFFFF }, { 0x022A, 0x0801, 0x0533, 0x0625, 0x0237, 0x0A1B, 0xFFFF }, { 0x0330, 0x0414, 0x0136, 0x0708, 0x0721, 0x0529,
 
420
      0x140D, 0xFFFF }, { 0x0400, 0x0324, 0x0606, 0x0038, 0x0039, 0x110B, 0x0D1F, 0xFFFF }, { 0x051A, 0x0813, 0x0B05, 0x0917,
 
421
      0xFFFF }, { 0x0723, 0x0435, 0x0627, 0xFFFF }, { 0x0412, 0x0234, 0x0F09, 0x022F, 0xFFFF }, { 0x0326, 0x0A0C, 0x012E,
 
422
      0x0811, 0x0A19, 0x1E0F, 0xFFFF }, { 0x0731, 0x0A03, 0x0915, 0x0D07, 0x052B, 0xFFFF }, { 0x0410, 0x0901, 0x0633, 0x0725,
 
423
      0x0337, 0x0B1B, 0x042D, 0x0C1D, 0xFFFF }, { 0x0518, 0x022C, 0x0629, 0x150D, 0xFFFF }, { 0x0332, 0x0502, 0x0821, 0x0139,
 
424
      0x120B, 0x0E1F, 0xFFFF }, { 0x0328, 0x061C, 0x0913, 0x0A17, 0xFFFF }, { 0x0420, 0x0535, 0x0C05, 0x0727, 0xFFFF }, {
 
425
      0x0823, 0x032F, 0xFFFF }, { 0x0516, 0x071E, 0x0F0E, 0x0911, 0x0B19, 0x1009, 0x1F0F, 0xFFFF }, { 0x0422, 0x0604, 0x090A,
 
426
      0x0B03, 0x0A15, 0x0E07, 0x062B, 0xFFFF }, { 0x0831, 0x0A01, 0x0733, 0x052D, 0x0D1D, 0xFFFF }, { 0x032A, 0x0825, 0x0437,
 
427
      0x0729, 0x0C1B, 0x160D, 0xFFFF }, { 0x0430, 0x0514, 0x0236, 0x0808, 0x0921, 0x0239, 0x130B, 0x0F1F, 0xFFFF }, { 0x0500,
 
428
      0x0424, 0x0706, 0x0138, 0x0A13, 0x0B17, 0xFFFF }, { 0x061A, 0x0635, 0x0D05, 0xFFFF }, { 0x0923, 0x0827, 0xFFFF }, {
 
429
      0x0512, 0x0334, 0x003A, 0x0A11, 0x1109, 0x003B, 0x042F, 0xFFFF }, { 0x0426, 0x0B0C, 0x022E, 0x0B15, 0x0F07, 0x0C19,
 
430
      0x072B, 0xFFFF }, { 0x0931, 0x0B01, 0x0C03, 0x062D, 0x0E1D, 0xFFFF }, { 0x0510, 0x0833, 0x0925, 0x0537, 0x0D1B, 0x170D,
 
431
      0xFFFF }, { 0x0618, 0x032C, 0x0A21, 0x0339, 0x0829, 0xFFFF }, { 0x0432, 0x0602, 0x0B13, 0x140B, 0x101F, 0xFFFF }, {
 
432
      0x0428, 0x071C, 0x0735, 0x0E05, 0x0C17, 0xFFFF }, { 0x0520, 0x0A23, 0x0927, 0xFFFF }, { 0x0B11, 0x1209, 0x013B, 0x052F,
 
433
      0xFFFF }, { 0x0616, 0x081E, 0x0D19, 0xFFFF }, { 0x0522, 0x0704, 0x0A0A, 0x0A31, 0x0D03, 0x0C15, 0x1007, 0x082B, 0x072D,
 
434
      0x0F1D, 0xFFFF }, { 0x0C01, 0x0933, 0x0A25, 0x0637, 0x0E1B, 0xFFFF }, { 0x042A, 0x0B21, 0x0929, 0x180D, 0xFFFF }, {
 
435
      0x0530, 0x0614, 0x0336, 0x0908, 0x0439, 0x150B, 0x111F, 0xFFFF }, { 0x0600, 0x0524, 0x0806, 0x0238, 0x0C13, 0x0F05,
 
436
      0x0D17, 0xFFFF }, { 0x071A, 0x0B23, 0x0835, 0x0A27, 0xFFFF }, { 0x1309, 0x023B, 0x062F, 0xFFFF }, { 0x0612, 0x0434,
 
437
      0x013A, 0x0C11, 0x0E19, 0xFFFF }, { 0x0526, 0x0C0C, 0x032E, 0x0B31, 0x0E03, 0x0D15, 0x1107, 0x092B, 0xFFFF }, { 0x0D01,
 
438
      0x0A33, 0x0B25, 0x0737, 0x0F1B, 0x082D, 0x101D, 0xFFFF }, { 0x0610, 0x0A29, 0x190D, 0xFFFF }, { 0x0718, 0x042C, 0x0C21,
 
439
      0x0539, 0x160B, 0x121F, 0xFFFF }, { 0x0532, 0x0702, 0x0D13, 0x0E17, 0xFFFF }, { 0x0528, 0x081C, 0x0935, 0x1005, 0x0B27,
 
440
      0xFFFF }, { 0x0620, 0x0C23, 0x033B, 0x072F, 0xFFFF }, { 0x0D11, 0x0F19, 0x1409, 0xFFFF }, { 0x0716, 0x003C, 0x091E,
 
441
      0x0F03, 0x0E15, 0x1207, 0x0A2B, 0x003D, 0xFFFF }, { 0x0622, 0x0804, 0x0B0A, 0x0C31, 0x0E01, 0x0B33, 0x092D, 0x111D,
 
442
      0xFFFF }, { 0x0C25, 0x0837, 0x0B29, 0x101B, 0x1A0D, 0xFFFF }, { 0x052A, 0x0D21, 0x0639, 0x170B, 0x131F, 0xFFFF }, {
 
443
      0x0630, 0x0714, 0x0436, 0x0A08, 0x0E13, 0x0F17, 0xFFFF }, { 0x0700, 0x0624, 0x0906, 0x0338, 0x0A35, 0x1105, 0xFFFF }, {
 
444
      0x081A, 0x0D23, 0x0C27, 0xFFFF }, { 0x0E11, 0x1509, 0x043B, 0x082F, 0xFFFF }, { 0x0712, 0x0534, 0x023A, 0x0F15, 0x1307,
 
445
      0x1019, 0x0B2B, 0x013D, 0xFFFF }, { 0x0626, 0x0D0C, 0x042E, 0x0D31, 0x0F01, 0x1003, 0x0A2D, 0x121D, 0xFFFF }, { 0x0C33,
 
446
      0x0D25, 0x0937, 0x111B, 0x1B0D, 0xFFFF }, { 0x0710, 0x0E21, 0x0739, 0x0C29, 0xFFFF }, { 0x0818, 0x052C, 0x0F13, 0x180B,
 
447
      0x141F, 0xFFFF }, { 0x0632, 0x0802, 0x0B35, 0x1205, 0x1017, 0xFFFF }, { 0x0628, 0x091C, 0x0E23, 0x0D27, 0xFFFF }, {
 
448
      0x0720, 0x0F11, 0x1609, 0x053B, 0x092F, 0xFFFF }, { 0x1119, 0x023D, 0xFFFF }, { 0x0816, 0x013C, 0x0A1E, 0x0E31, 0x1103,
 
449
      0x1015, 0x1407, 0x0C2B, 0x0B2D, 0x131D, 0xFFFF }, { 0x0722, 0x0904, 0x0C0A, 0x1001, 0x0D33, 0x0E25, 0x0A37, 0x121B,
 
450
      0xFFFF }, { 0x0F21, 0x0D29, 0x1C0D, 0xFFFF }, { 0x062A, 0x0839, 0x190B, 0x151F, 0xFFFF }, { 0x0730, 0x0814, 0x0536,
 
451
      0x0B08, 0x1013, 0x1305, 0x1117, 0xFFFF }, { 0x0800, 0x0724, 0x0A06, 0x0438, 0x0F23, 0x0C35, 0x0E27, 0xFFFF }, { 0x091A,
 
452
      0x1709, 0x063B, 0x0A2F, 0xFFFF }, { 0x1011, 0x1219, 0x033D, 0xFFFF }, { 0x0812, 0x0634, 0x033A, 0x0F31, 0x1203, 0x1115,
 
453
      0x1507, 0x0D2B, 0xFFFF }, { 0x0726, 0x0E0C, 0x052E, 0x1101, 0x0E33, 0x0F25, 0x0B37, 0x131B, 0x0C2D, 0x141D, 0xFFFF }, {
 
454
      0x0E29, 0x1D0D, 0xFFFF }, { 0x0810, 0x1021, 0x0939, 0x1A0B, 0x161F, 0xFFFF }, { 0x0918, 0x062C, 0x1113, 0x1217, 0xFFFF
 
455
      }, { 0x0732, 0x0902, 0x0D35, 0x1405, 0x0F27, 0xFFFF }, { 0x0728, 0x0A1C, 0x1023, 0x073B, 0x0B2F, 0xFFFF }, { 0x0820,
 
456
      0x1111, 0x1319, 0x1809, 0xFFFF }, { 0x1303, 0x1215, 0x1607, 0x0E2B, 0x043D, 0xFFFF }, { 0x0916, 0x023C, 0x0B1E, 0x1031,
 
457
      0x1201, 0x0F33, 0x0D2D, 0x151D, 0xFFFF }, { 0x0822, 0x0A04, 0x0D0A, 0x1025, 0x0C37, 0x0F29, 0x141B, 0x1E0D, 0xFFFF }, {
 
458
      0x1121, 0x0A39, 0x1B0B, 0x171F, 0xFFFF }, { 0x072A, 0x1213, 0x1317, 0xFFFF }, { 0x0830, 0x0914, 0x0636, 0x0C08, 0x0E35,
 
459
      0x1505, 0xFFFF }, { 0x0900, 0x0824, 0x0B06, 0x0538, 0x1123, 0x1027, 0xFFFF }, { 0x0A1A, 0x1211, 0x1909, 0x083B, 0x0C2F,
 
460
      0xFFFF }, { 0x1315, 0x1707, 0x1419, 0x0F2B, 0x053D, 0xFFFF }, { 0x0912, 0x0734, 0x043A, 0x1131, 0x1301, 0x1403, 0x0E2D,
 
461
      0x161D, 0xFFFF }, { 0x0826, 0x0F0C, 0x062E, 0x1033, 0x1125, 0x0D37, 0x151B, 0x1F0D, 0xFFFF }, { 0x1221, 0x0B39, 0x1029,
 
462
      0xFFFF }, { 0x0910, 0x1313, 0x1C0B, 0x181F, 0xFFFF }, { 0x0A18, 0x072C, 0x0F35, 0x1605, 0x1417, 0xFFFF }, { 0x0832,
 
463
      0x0A02, 0x1223, 0x1127, 0xFFFF }, { 0x0828, 0x0B1C, 0x1311, 0x1A09, 0x093B, 0x0D2F, 0xFFFF }, { 0x0920, 0x1519, 0x063D,
 
464
      0xFFFF }, { 0x1231, 0x1503, 0x1415, 0x1807, 0x102B, 0x0F2D, 0x171D, 0xFFFF }, { 0x0A16, 0x033C, 0x0C1E, 0x1401, 0x1133,
 
465
      0x1225, 0x0E37, 0x161B, 0xFFFF }, { 0x0922, 0x0B04, 0x0E0A, 0x1321, 0x1129, 0xFFFF }, { 0x0C39, 0x1D0B, 0x191F, 0xFFFF
 
466
      }, { 0x082A, 0x1413, 0x1705, 0x1517, 0xFFFF }, { 0x0930, 0x0A14, 0x0736, 0x0D08, 0x1323, 0x1035, 0x1227, 0xFFFF }, {
 
467
      0x0A00, 0x0924, 0x0C06, 0x0638, 0x1B09, 0x0A3B, 0x0E2F, 0xFFFF }, { 0x0B1A, 0x1411, 0x1619, 0x073D, 0xFFFF }, { 0x1331,
 
468
      0x1603, 0x1515, 0x1907, 0x112B, 0xFFFF }, { 0x0A12, 0x0834, 0x053A, 0x1501, 0x1233, 0x1325, 0x0F37, 0x171B, 0x102D,
 
469
      0x181D, 0xFFFF }, { 0x0926, 0x072E, 0x1229, 0xFFFF }, { 0x1421, 0x0D39, 0x1E0B, 0x1A1F, 0xFFFF }, { 0x0A10, 0x1513,
 
470
      0x1617, 0xFFFF }, { 0x0B18, 0x082C, 0x1135, 0x1805, 0x1327, 0xFFFF }, { 0x0932, 0x0B02, 0x1423, 0x0B3B, 0x0F2F, 0xFFFF
 
471
      }, { 0x0928, 0x0C1C, 0x1511, 0x1719, 0x1C09, 0xFFFF }, { 0x0A20, 0x1703, 0x1615, 0x1A07, 0x122B, 0x083D, 0xFFFF }, {
 
472
      0x1431, 0x1601, 0x1333, 0x112D, 0x191D, 0xFFFF }, { 0x0B16, 0x043C, 0x0D1E, 0x1425, 0x1037, 0x1329, 0x181B, 0xFFFF }, {
 
473
      0x0A22, 0x0C04, 0x0F0A, 0x1521, 0x0E39, 0x1F0B, 0x1B1F, 0xFFFF }, { 0x1613, 0x1717, 0xFFFF }, { 0x092A, 0x1235, 0x1905,
 
474
      0xFFFF }, { 0x0A30, 0x0B14, 0x0836, 0x0E08, 0x1523, 0x1427, 0xFFFF }, { 0x0B00, 0x0A24, 0x0D06, 0x0738, 0x1611, 0x1D09,
 
475
      0x0C3B, 0x102F, 0xFFFF }, { 0x0C1A, 0x1715, 0x1B07, 0x1819, 0x132B, 0x093D, 0xFFFF }, { 0x1531, 0x1701, 0x1803, 0x122D,
 
476
      0x1A1D, 0xFFFF }, { 0x0B12, 0x0934, 0x063A, 0x1433, 0x1525, 0x1137, 0x191B, 0xFFFF }, { 0x0A26, 0x003E, 0x082E, 0x1621,
 
477
      0x0F39, 0x1429, 0x003F, 0xFFFF }, { 0x1713, 0x1C1F, 0xFFFF }, { 0x0B10, 0x1335, 0x1A05, 0x1817, 0xFFFF }, { 0x0C18,
 
478
      0x092C, 0x1623, 0x1527, 0xFFFF }, { 0x0A32, 0x0C02, 0x1711, 0x1E09, 0x0D3B, 0x112F, 0xFFFF }, { 0x0A28, 0x0D1C, 0x1919,
 
479
      0x0A3D, 0xFFFF }, { 0x0B20, 0x1631, 0x1903, 0x1815, 0x1C07, 0x142B, 0x132D, 0x1B1D, 0xFFFF }, { 0x1801, 0x1533, 0x1625,
 
480
      0x1237, 0x1A1B, 0xFFFF }, { 0x0C16, 0x053C, 0x0E1E, 0x1721, 0x1529, 0x013F, 0xFFFF }, { 0x0B22, 0x0D04, 0x1039, 0x1D1F,
 
481
      0xFFFF }, { 0x1813, 0x1B05, 0x1917, 0xFFFF }, { 0x0A2A, 0x1723, 0x1435, 0x1627, 0xFFFF }, { 0x0B30, 0x0C14, 0x0936,
 
482
      0x0F08, 0x1F09, 0x0E3B, 0x122F, 0xFFFF }, { 0x0C00, 0x0B24, 0x0E06, 0x0838, 0x1811, 0x1A19, 0x0B3D, 0xFFFF }, { 0x0D1A,
 
483
      0x1731, 0x1A03, 0x1915, 0x1D07, 0x152B, 0xFFFF }, { 0x1901, 0x1633, 0x1725, 0x1337, 0x1B1B, 0x142D, 0x1C1D, 0xFFFF }, {
 
484
      0x0C12, 0x0A34, 0x073A, 0x1629, 0x023F, 0xFFFF }, { 0x0B26, 0x013E, 0x092E, 0x1821, 0x1139, 0x1E1F, 0xFFFF }, { 0x1913,
 
485
      0x1A17, 0xFFFF }, { 0x0C10, 0x1535, 0x1C05, 0x1727, 0xFFFF }, { 0x0D18, 0x0A2C, 0x1823, 0x0F3B, 0x132F, 0xFFFF }, {
 
486
      0x0B32, 0x0D02, 0x1911, 0x1B19, 0xFFFF }, { 0x0B28, 0x0E1C, 0x1B03, 0x1A15, 0x1E07, 0x162B, 0x0C3D, 0xFFFF }, { 0x0C20,
 
487
      0x1831, 0x1A01, 0x1733, 0x152D, 0x1D1D, 0xFFFF }, { 0x1825, 0x1437, 0x1729, 0x1C1B, 0x033F, 0xFFFF }, { 0x0D16, 0x063C,
 
488
      0x0F1E, 0x1921, 0x1239, 0x1F1F, 0xFFFF }, { 0x0C22, 0x0E04, 0x1A13, 0x1B17, 0xFFFF }, { 0x1635, 0x1D05, 0xFFFF }, {
 
489
      0x0B2A, 0x1923, 0x1827, 0xFFFF }, { 0x0C30, 0x0D14, 0x0A36, 0x1A11, 0x103B, 0x142F, 0xFFFF }, { 0x0D00, 0x0C24, 0x0F06,
 
490
      0x0938, 0x1B15, 0x1F07, 0x1C19, 0x172B, 0x0D3D, 0xFFFF }, { 0x0E1A, 0x1931, 0x1B01, 0x1C03, 0x162D, 0x1E1D, 0xFFFF }, {
 
491
      0x1833, 0x1925, 0x1537, 0x1D1B, 0xFFFF }, { 0x0D12, 0x0B34, 0x083A, 0x1A21, 0x1339, 0x1829, 0x043F, 0xFFFF }, { 0x0C26,
 
492
      0x023E, 0x0A2E, 0x1B13, 0xFFFF }, { 0x1735, 0x1E05, 0x1C17, 0xFFFF }, { 0x0D10, 0x1A23, 0x1927, 0xFFFF }, { 0x0E18,
 
493
      0x0B2C, 0x1B11, 0x113B, 0x152F, 0xFFFF }, { 0x0C32, 0x0E02, 0x1D19, 0x0E3D, 0xFFFF }, { 0x0C28, 0x0F1C, 0x1A31, 0x1D03,
 
494
      0x1C15, 0x182B, 0x172D, 0x1F1D, 0xFFFF }, { 0x0D20, 0x1C01, 0x1933, 0x1A25, 0x1637, 0x1E1B, 0xFFFF }, { 0x1B21, 0x1929,
 
495
      0x053F, 0xFFFF }, { 0x0E16, 0x073C, 0x1439, 0xFFFF }, { 0x0D22, 0x0F04, 0x1C13, 0x1F05, 0x1D17, 0xFFFF }, { 0x1B23,
 
496
      0x1835, 0x1A27, 0xFFFF }, { 0x0C2A, 0x123B, 0x162F, 0xFFFF }, { 0x0D30, 0x0E14, 0x0B36, 0x1C11, 0x1E19, 0x0F3D, 0xFFFF
 
497
      }, { 0x0E00, 0x0D24, 0x0A38, 0x1B31, 0x1E03, 0x1D15, 0x192B, 0xFFFF }, { 0x0F1A, 0x1D01, 0x1A33, 0x1B25, 0x1737, 0x1F1B,
 
498
      0x182D, 0xFFFF }, { 0x1A29, 0x063F, 0xFFFF }, { 0x0E12, 0x0C34, 0x093A, 0x1C21, 0x1539, 0xFFFF }, { 0x0D26, 0x033E,
 
499
      0x0B2E, 0x1D13, 0x1E17, 0xFFFF }, { 0x1935, 0x1B27, 0xFFFF }, { 0x0E10, 0x1C23, 0x133B, 0x172F, 0xFFFF }, { 0x0F18,
 
500
      0x0C2C, 0x1D11, 0x1F19, 0xFFFF }, { 0x0D32, 0x0F02, 0x1F03, 0x1E15, 0x1A2B, 0x103D, 0xFFFF }, { 0x0D28, 0x1C31, 0x1E01,
 
501
      0x1B33, 0x192D, 0xFFFF }, { 0x0E20, 0x1C25, 0x1837, 0x1B29, 0x073F, 0xFFFF }, { 0x1D21, 0x1639, 0xFFFF }, { 0x0F16,
 
502
      0x083C, 0x1E13, 0x1F17, 0xFFFF }, { 0x0E22, 0x1A35, 0xFFFF }, { 0x1D23, 0x1C27, 0xFFFF }, { 0x0D2A, 0x1E11, 0x143B,
 
503
      0x182F, 0xFFFF }, { 0x0E30, 0x0F14, 0x0C36, 0x1F15, 0x1B2B, 0x113D, 0xFFFF }, { 0x0F00, 0x0E24, 0x0B38, 0x1D31, 0x1F01,
 
504
      0x1A2D, 0xFFFF }, { 0x1C33, 0x1D25, 0x1937, 0xFFFF }, { 0x1E21, 0x1739, 0x1C29, 0x083F, 0xFFFF }, { 0x0F12, 0x0D34,
 
505
      0x0A3A, 0x1F13, 0xFFFF }, { 0x0E26, 0x043E, 0x0C2E, 0x1B35, 0xFFFF }, { 0x1E23, 0x1D27, 0xFFFF }, { 0x0F10, 0x1F11,
 
506
      0x153B, 0x192F, 0xFFFF }, { 0x0D2C, 0x123D, 0xFFFF },
 
507
   };
 
508
 
 
509
   struct etc1_block
 
510
   {
 
511
      // big endian uint64:
 
512
      // bit ofs:  56  48  40  32  24  16   8   0
 
513
      // byte ofs: b0, b1, b2, b3, b4, b5, b6, b7 
 
514
      union 
 
515
      {
 
516
         uint64 m_uint64;
 
517
         uint8 m_bytes[8];
 
518
      };
 
519
 
 
520
      uint8 m_low_color[2];
 
521
      uint8 m_high_color[2];
 
522
 
 
523
      enum { cNumSelectorBytes = 4 };
 
524
      uint8 m_selectors[cNumSelectorBytes];
 
525
 
 
526
      inline void clear()
 
527
      {
 
528
         zero_this(this);
 
529
      }
 
530
 
 
531
      inline uint get_byte_bits(uint ofs, uint num) const
 
532
      {
 
533
         RG_ETC1_ASSERT((ofs + num) <= 64U);
 
534
         RG_ETC1_ASSERT(num && (num <= 8U));
 
535
         RG_ETC1_ASSERT((ofs >> 3) == ((ofs + num - 1) >> 3));
 
536
         const uint byte_ofs = 7 - (ofs >> 3);
 
537
         const uint byte_bit_ofs = ofs & 7;
 
538
         return (m_bytes[byte_ofs] >> byte_bit_ofs) & ((1 << num) - 1);
 
539
      }
 
540
 
 
541
      inline void set_byte_bits(uint ofs, uint num, uint bits)
 
542
      {
 
543
         RG_ETC1_ASSERT((ofs + num) <= 64U);
 
544
         RG_ETC1_ASSERT(num && (num < 32U));
 
545
         RG_ETC1_ASSERT((ofs >> 3) == ((ofs + num - 1) >> 3));
 
546
         RG_ETC1_ASSERT(bits < (1U << num));
 
547
         const uint byte_ofs = 7 - (ofs >> 3);
 
548
         const uint byte_bit_ofs = ofs & 7;
 
549
         const uint mask = (1 << num) - 1;
 
550
         m_bytes[byte_ofs] &= ~(mask << byte_bit_ofs);
 
551
         m_bytes[byte_ofs] |= (bits << byte_bit_ofs);
 
552
      }
 
553
 
 
554
      // false = left/right subblocks
 
555
      // true = upper/lower subblocks
 
556
      inline bool get_flip_bit() const 
 
557
      {
 
558
         return (m_bytes[3] & 1) != 0;
 
559
      }   
 
560
 
 
561
      inline void set_flip_bit(bool flip)
 
562
      {
 
563
         m_bytes[3] &= ~1;
 
564
         m_bytes[3] |= static_cast<uint8>(flip);
 
565
      }
 
566
 
 
567
      inline bool get_diff_bit() const
 
568
      {
 
569
         return (m_bytes[3] & 2) != 0;
 
570
      }
 
571
 
 
572
      inline void set_diff_bit(bool diff)
 
573
      {
 
574
         m_bytes[3] &= ~2;
 
575
         m_bytes[3] |= (static_cast<uint>(diff) << 1);
 
576
      }
 
577
 
 
578
      // Returns intensity modifier table (0-7) used by subblock subblock_id.
 
579
      // subblock_id=0 left/top (CW 1), 1=right/bottom (CW 2)
 
580
      inline uint get_inten_table(uint subblock_id) const
 
581
      {
 
582
         RG_ETC1_ASSERT(subblock_id < 2);
 
583
         const uint ofs = subblock_id ? 2 : 5;
 
584
         return (m_bytes[3] >> ofs) & 7;
 
585
      }
 
586
 
 
587
      // Sets intensity modifier table (0-7) used by subblock subblock_id (0 or 1)
 
588
      inline void set_inten_table(uint subblock_id, uint t)
 
589
      {
 
590
         RG_ETC1_ASSERT(subblock_id < 2);
 
591
         RG_ETC1_ASSERT(t < 8);
 
592
         const uint ofs = subblock_id ? 2 : 5;
 
593
         m_bytes[3] &= ~(7 << ofs);
 
594
         m_bytes[3] |= (t << ofs);
 
595
      }
 
596
 
 
597
      // Returned selector value ranges from 0-3 and is a direct index into g_etc1_inten_tables.
 
598
      inline uint get_selector(uint x, uint y) const
 
599
      {
 
600
         RG_ETC1_ASSERT((x | y) < 4);
 
601
 
 
602
         const uint bit_index = x * 4 + y;
 
603
         const uint byte_bit_ofs = bit_index & 7;
 
604
         const uint8 *p = &m_bytes[7 - (bit_index >> 3)];
 
605
         const uint lsb = (p[0] >> byte_bit_ofs) & 1;
 
606
         const uint msb = (p[-2] >> byte_bit_ofs) & 1;
 
607
         const uint val = lsb | (msb << 1);
 
608
 
 
609
         return g_etc1_to_selector_index[val];
 
610
      }
 
611
 
 
612
      // Selector "val" ranges from 0-3 and is a direct index into g_etc1_inten_tables.
 
613
      inline void set_selector(uint x, uint y, uint val)
 
614
      {
 
615
         RG_ETC1_ASSERT((x | y | val) < 4);
 
616
         const uint bit_index = x * 4 + y;
 
617
 
 
618
         uint8 *p = &m_bytes[7 - (bit_index >> 3)];
 
619
 
 
620
         const uint byte_bit_ofs = bit_index & 7;
 
621
         const uint mask = 1 << byte_bit_ofs;
 
622
 
 
623
         const uint etc1_val = g_selector_index_to_etc1[val];
 
624
 
 
625
         const uint lsb = etc1_val & 1;
 
626
         const uint msb = etc1_val >> 1;
 
627
 
 
628
         p[0] &= ~mask;
 
629
         p[0] |= (lsb << byte_bit_ofs);
 
630
 
 
631
         p[-2] &= ~mask;
 
632
         p[-2] |= (msb << byte_bit_ofs);
 
633
      }
 
634
 
 
635
      inline void set_base4_color(uint idx, uint16 c)
 
636
      {
 
637
         if (idx)
 
638
         {
 
639
            set_byte_bits(cETC1AbsColor4R2BitOffset, 4, (c >> 8) & 15);
 
640
            set_byte_bits(cETC1AbsColor4G2BitOffset, 4, (c >> 4) & 15);
 
641
            set_byte_bits(cETC1AbsColor4B2BitOffset, 4, c & 15);
 
642
         }
 
643
         else
 
644
         {
 
645
            set_byte_bits(cETC1AbsColor4R1BitOffset, 4, (c >> 8) & 15);
 
646
            set_byte_bits(cETC1AbsColor4G1BitOffset, 4, (c >> 4) & 15);
 
647
            set_byte_bits(cETC1AbsColor4B1BitOffset, 4, c & 15);
 
648
         }
 
649
      }
 
650
 
 
651
      inline uint16 get_base4_color(uint idx) const
 
652
      {
 
653
         uint r, g, b;
 
654
         if (idx)
 
655
         {
 
656
            r = get_byte_bits(cETC1AbsColor4R2BitOffset, 4);
 
657
            g = get_byte_bits(cETC1AbsColor4G2BitOffset, 4);
 
658
            b = get_byte_bits(cETC1AbsColor4B2BitOffset, 4);
 
659
         }
 
660
         else
 
661
         {
 
662
            r = get_byte_bits(cETC1AbsColor4R1BitOffset, 4);
 
663
            g = get_byte_bits(cETC1AbsColor4G1BitOffset, 4);
 
664
            b = get_byte_bits(cETC1AbsColor4B1BitOffset, 4);
 
665
         }
 
666
         return static_cast<uint16>(b | (g << 4U) | (r << 8U));
 
667
      }
 
668
 
 
669
      inline void set_base5_color(uint16 c)
 
670
      {
 
671
         set_byte_bits(cETC1BaseColor5RBitOffset, 5, (c >> 10) & 31);
 
672
         set_byte_bits(cETC1BaseColor5GBitOffset, 5, (c >> 5) & 31);
 
673
         set_byte_bits(cETC1BaseColor5BBitOffset, 5, c & 31);
 
674
      }
 
675
 
 
676
      inline uint16 get_base5_color() const
 
677
      {
 
678
         const uint r = get_byte_bits(cETC1BaseColor5RBitOffset, 5);
 
679
         const uint g = get_byte_bits(cETC1BaseColor5GBitOffset, 5);
 
680
         const uint b = get_byte_bits(cETC1BaseColor5BBitOffset, 5);
 
681
         return static_cast<uint16>(b | (g << 5U) | (r << 10U));
 
682
      }
 
683
 
 
684
      void set_delta3_color(uint16 c)
 
685
      {
 
686
         set_byte_bits(cETC1DeltaColor3RBitOffset, 3, (c >> 6) & 7);
 
687
         set_byte_bits(cETC1DeltaColor3GBitOffset, 3, (c >> 3) & 7);
 
688
         set_byte_bits(cETC1DeltaColor3BBitOffset, 3, c & 7);
 
689
      }
 
690
 
 
691
      inline uint16 get_delta3_color() const
 
692
      {
 
693
         const uint r = get_byte_bits(cETC1DeltaColor3RBitOffset, 3);
 
694
         const uint g = get_byte_bits(cETC1DeltaColor3GBitOffset, 3);
 
695
         const uint b = get_byte_bits(cETC1DeltaColor3BBitOffset, 3);
 
696
         return static_cast<uint16>(b | (g << 3U) | (r << 6U));
 
697
      }
 
698
 
 
699
      // Base color 5
 
700
      static uint16 pack_color5(const color_quad_u8& color, bool scaled, uint bias = 127U);
 
701
      static uint16 pack_color5(uint r, uint g, uint b, bool scaled, uint bias = 127U);
 
702
 
 
703
      static color_quad_u8 unpack_color5(uint16 packed_color5, bool scaled, uint alpha = 255U);
 
704
      static void unpack_color5(uint& r, uint& g, uint& b, uint16 packed_color, bool scaled);
 
705
 
 
706
      static bool unpack_color5(color_quad_u8& result, uint16 packed_color5, uint16 packed_delta3, bool scaled, uint alpha = 255U);
 
707
      static bool unpack_color5(uint& r, uint& g, uint& b, uint16 packed_color5, uint16 packed_delta3, bool scaled, uint alpha = 255U);
 
708
 
 
709
      // Delta color 3
 
710
      // Inputs range from -4 to 3 (cETC1ColorDeltaMin to cETC1ColorDeltaMax)
 
711
      static uint16 pack_delta3(int r, int g, int b);
 
712
 
 
713
      // Results range from -4 to 3 (cETC1ColorDeltaMin to cETC1ColorDeltaMax)
 
714
      static void unpack_delta3(int& r, int& g, int& b, uint16 packed_delta3);
 
715
 
 
716
      // Abs color 4
 
717
      static uint16 pack_color4(const color_quad_u8& color, bool scaled, uint bias = 127U);
 
718
      static uint16 pack_color4(uint r, uint g, uint b, bool scaled, uint bias = 127U);
 
719
 
 
720
      static color_quad_u8 unpack_color4(uint16 packed_color4, bool scaled, uint alpha = 255U);
 
721
      static void unpack_color4(uint& r, uint& g, uint& b, uint16 packed_color4, bool scaled);
 
722
 
 
723
      // subblock colors
 
724
      static void get_diff_subblock_colors(color_quad_u8* pDst, uint16 packed_color5, uint table_idx);
 
725
      static bool get_diff_subblock_colors(color_quad_u8* pDst, uint16 packed_color5, uint16 packed_delta3, uint table_idx);
 
726
      static void get_abs_subblock_colors(color_quad_u8* pDst, uint16 packed_color4, uint table_idx);
 
727
 
 
728
      static inline void unscaled_to_scaled_color(color_quad_u8& dst, const color_quad_u8& src, bool color4)
 
729
      {
 
730
         if (color4)
 
731
         {
 
732
            dst.r = src.r | (src.r << 4);
 
733
            dst.g = src.g | (src.g << 4);
 
734
            dst.b = src.b | (src.b << 4);
 
735
         }
 
736
         else
 
737
         {
 
738
            dst.r = (src.r >> 2) | (src.r << 3);
 
739
            dst.g = (src.g >> 2) | (src.g << 3);
 
740
            dst.b = (src.b >> 2) | (src.b << 3);
 
741
         }
 
742
         dst.a = src.a;
 
743
      }
 
744
   };
 
745
 
 
746
   // Returns pointer to sorted array.
 
747
   template<typename T, typename Q>
 
748
   T* indirect_radix_sort(uint num_indices, T* pIndices0, T* pIndices1, const Q* pKeys, uint key_ofs, uint key_size, bool init_indices)
 
749
   {  
 
750
      RG_ETC1_ASSERT((key_ofs >= 0) && (key_ofs < sizeof(T)));
 
751
      RG_ETC1_ASSERT((key_size >= 1) && (key_size <= 4));
 
752
 
 
753
      if (init_indices)
 
754
      {
 
755
         T* p = pIndices0;
 
756
         T* q = pIndices0 + (num_indices >> 1) * 2;
 
757
         uint i;
 
758
         for (i = 0; p != q; p += 2, i += 2)
 
759
         {
 
760
            p[0] = static_cast<T>(i);
 
761
            p[1] = static_cast<T>(i + 1); 
 
762
         }
 
763
 
 
764
         if (num_indices & 1)
 
765
            *p = static_cast<T>(i);
 
766
      }
 
767
 
 
768
      uint hist[256 * 4];
 
769
 
 
770
      memset(hist, 0, sizeof(hist[0]) * 256 * key_size);
 
771
 
 
772
#define RG_ETC1_GET_KEY(p) (*(const uint*)((const uint8*)(pKeys + *(p)) + key_ofs))
 
773
#define RG_ETC1_GET_KEY_FROM_INDEX(i) (*(const uint*)((const uint8*)(pKeys + (i)) + key_ofs))
 
774
 
 
775
      if (key_size == 4)
 
776
      {
 
777
         T* p = pIndices0;
 
778
         T* q = pIndices0 + num_indices;
 
779
         for ( ; p != q; p++)
 
780
         {
 
781
            const uint key = RG_ETC1_GET_KEY(p);
 
782
 
 
783
            hist[        key        & 0xFF]++;
 
784
            hist[256 + ((key >>  8) & 0xFF)]++;
 
785
            hist[512 + ((key >> 16) & 0xFF)]++;
 
786
            hist[768 + ((key >> 24) & 0xFF)]++;
 
787
         }
 
788
      }
 
789
      else if (key_size == 3)
 
790
      {
 
791
         T* p = pIndices0;
 
792
         T* q = pIndices0 + num_indices;
 
793
         for ( ; p != q; p++)
 
794
         {
 
795
            const uint key = RG_ETC1_GET_KEY(p);
 
796
 
 
797
            hist[        key        & 0xFF]++;
 
798
            hist[256 + ((key >>  8) & 0xFF)]++;
 
799
            hist[512 + ((key >> 16) & 0xFF)]++;
 
800
         }
 
801
      }   
 
802
      else if (key_size == 2)
 
803
      {
 
804
         T* p = pIndices0;
 
805
         T* q = pIndices0 + (num_indices >> 1) * 2;
 
806
 
 
807
         for ( ; p != q; p += 2)
 
808
         {
 
809
            const uint key0 = RG_ETC1_GET_KEY(p);
 
810
            const uint key1 = RG_ETC1_GET_KEY(p+1);
 
811
 
 
812
            hist[        key0         & 0xFF]++;
 
813
            hist[256 + ((key0 >>  8) & 0xFF)]++;
 
814
 
 
815
            hist[        key1        & 0xFF]++;
 
816
            hist[256 + ((key1 >>  8) & 0xFF)]++;
 
817
         }
 
818
 
 
819
         if (num_indices & 1)
 
820
         {
 
821
            const uint key = RG_ETC1_GET_KEY(p);
 
822
 
 
823
            hist[        key        & 0xFF]++;
 
824
            hist[256 + ((key >>  8) & 0xFF)]++;
 
825
         }
 
826
      }      
 
827
      else
 
828
      {
 
829
         RG_ETC1_ASSERT(key_size == 1);
 
830
         if (key_size != 1)
 
831
            return NULL;
 
832
 
 
833
         T* p = pIndices0;
 
834
         T* q = pIndices0 + (num_indices >> 1) * 2;
 
835
 
 
836
         for ( ; p != q; p += 2)
 
837
         {
 
838
            const uint key0 = RG_ETC1_GET_KEY(p);
 
839
            const uint key1 = RG_ETC1_GET_KEY(p+1);
 
840
 
 
841
            hist[key0 & 0xFF]++;
 
842
            hist[key1 & 0xFF]++;
 
843
         }
 
844
 
 
845
         if (num_indices & 1)
 
846
         {
 
847
            const uint key = RG_ETC1_GET_KEY(p);
 
848
 
 
849
            hist[key & 0xFF]++;
 
850
         }
 
851
      }      
 
852
 
 
853
      T* pCur = pIndices0;
 
854
      T* pNew = pIndices1;
 
855
 
 
856
      for (uint pass = 0; pass < key_size; pass++)
 
857
      {
 
858
         const uint* pHist = &hist[pass << 8];
 
859
 
 
860
         uint offsets[256];
 
861
 
 
862
         uint cur_ofs = 0;
 
863
         for (uint i = 0; i < 256; i += 2)
 
864
         {
 
865
            offsets[i] = cur_ofs;
 
866
            cur_ofs += pHist[i];
 
867
 
 
868
            offsets[i+1] = cur_ofs;
 
869
            cur_ofs += pHist[i+1];
 
870
         }
 
871
 
 
872
         const uint pass_shift = pass << 3;
 
873
 
 
874
         T* p = pCur;
 
875
         T* q = pCur + (num_indices >> 1) * 2;
 
876
 
 
877
         for ( ; p != q; p += 2)
 
878
         {
 
879
            uint index0 = p[0];
 
880
            uint index1 = p[1];
 
881
 
 
882
            uint c0 = (RG_ETC1_GET_KEY_FROM_INDEX(index0) >> pass_shift) & 0xFF;
 
883
            uint c1 = (RG_ETC1_GET_KEY_FROM_INDEX(index1) >> pass_shift) & 0xFF;
 
884
 
 
885
            if (c0 == c1)
 
886
            {
 
887
               uint dst_offset0 = offsets[c0];
 
888
 
 
889
               offsets[c0] = dst_offset0 + 2;
 
890
 
 
891
               pNew[dst_offset0] = static_cast<T>(index0);
 
892
               pNew[dst_offset0 + 1] = static_cast<T>(index1);
 
893
            }
 
894
            else
 
895
            {
 
896
               uint dst_offset0 = offsets[c0]++;
 
897
               uint dst_offset1 = offsets[c1]++;
 
898
 
 
899
               pNew[dst_offset0] = static_cast<T>(index0);
 
900
               pNew[dst_offset1] = static_cast<T>(index1);
 
901
            }
 
902
         }
 
903
 
 
904
         if (num_indices & 1)
 
905
         {
 
906
            uint index = *p;
 
907
            uint c = (RG_ETC1_GET_KEY_FROM_INDEX(index) >> pass_shift) & 0xFF;
 
908
 
 
909
            uint dst_offset = offsets[c];
 
910
            offsets[c] = dst_offset + 1;
 
911
 
 
912
            pNew[dst_offset] = static_cast<T>(index);
 
913
         }
 
914
 
 
915
         T* t = pCur;
 
916
         pCur = pNew;
 
917
         pNew = t;
 
918
      }            
 
919
 
 
920
      return pCur;
 
921
   }
 
922
 
 
923
#undef RG_ETC1_GET_KEY
 
924
#undef RG_ETC1_GET_KEY_FROM_INDEX
 
925
 
 
926
   uint16 etc1_block::pack_color5(const color_quad_u8& color, bool scaled, uint bias)
 
927
   {
 
928
      return pack_color5(color.r, color.g, color.b, scaled, bias);
 
929
   }
 
930
   
 
931
   uint16 etc1_block::pack_color5(uint r, uint g, uint b, bool scaled, uint bias)
 
932
   {
 
933
      if (scaled)
 
934
      {
 
935
         r = (r * 31U + bias) / 255U;
 
936
         g = (g * 31U + bias) / 255U;
 
937
         b = (b * 31U + bias) / 255U;
 
938
      }
 
939
 
 
940
      r = rg_etc1::minimum(r, 31U);
 
941
      g = rg_etc1::minimum(g, 31U);
 
942
      b = rg_etc1::minimum(b, 31U);
 
943
 
 
944
      return static_cast<uint16>(b | (g << 5U) | (r << 10U));
 
945
   }
 
946
 
 
947
   color_quad_u8 etc1_block::unpack_color5(uint16 packed_color5, bool scaled, uint alpha)
 
948
   {
 
949
      uint b = packed_color5 & 31U;
 
950
      uint g = (packed_color5 >> 5U) & 31U;
 
951
      uint r = (packed_color5 >> 10U) & 31U;
 
952
 
 
953
      if (scaled)
 
954
      {
 
955
         b = (b << 3U) | (b >> 2U);
 
956
         g = (g << 3U) | (g >> 2U);
 
957
         r = (r << 3U) | (r >> 2U);
 
958
      }
 
959
 
 
960
      return color_quad_u8(cNoClamp, r, g, b, rg_etc1::minimum(alpha, 255U));
 
961
   }
 
962
 
 
963
   void etc1_block::unpack_color5(uint& r, uint& g, uint& b, uint16 packed_color5, bool scaled)
 
964
   {
 
965
      color_quad_u8 c(unpack_color5(packed_color5, scaled, 0));
 
966
      r = c.r;
 
967
      g = c.g;
 
968
      b = c.b;
 
969
   }
 
970
 
 
971
   bool etc1_block::unpack_color5(color_quad_u8& result, uint16 packed_color5, uint16 packed_delta3, bool scaled, uint alpha)
 
972
   {
 
973
      int dc_r, dc_g, dc_b;
 
974
      unpack_delta3(dc_r, dc_g, dc_b, packed_delta3);
 
975
      
 
976
      int b = (packed_color5 & 31U) + dc_b;
 
977
      int g = ((packed_color5 >> 5U) & 31U) + dc_g;
 
978
      int r = ((packed_color5 >> 10U) & 31U) + dc_r;
 
979
 
 
980
      bool success = true;
 
981
      if (static_cast<uint>(r | g | b) > 31U)
 
982
      {
 
983
         success = false;
 
984
         r = rg_etc1::clamp<int>(r, 0, 31);
 
985
         g = rg_etc1::clamp<int>(g, 0, 31);
 
986
         b = rg_etc1::clamp<int>(b, 0, 31);
 
987
      }
 
988
 
 
989
      if (scaled)
 
990
      {
 
991
         b = (b << 3U) | (b >> 2U);
 
992
         g = (g << 3U) | (g >> 2U);
 
993
         r = (r << 3U) | (r >> 2U);
 
994
      }
 
995
 
 
996
      result.set_noclamp_rgba(r, g, b, rg_etc1::minimum(alpha, 255U));
 
997
      return success;
 
998
   }
 
999
 
 
1000
   bool etc1_block::unpack_color5(uint& r, uint& g, uint& b, uint16 packed_color5, uint16 packed_delta3, bool scaled, uint alpha)
 
1001
   {
 
1002
      color_quad_u8 result;
 
1003
      const bool success = unpack_color5(result, packed_color5, packed_delta3, scaled, alpha);
 
1004
      r = result.r;
 
1005
      g = result.g;
 
1006
      b = result.b;
 
1007
      return success;
 
1008
   }
 
1009
     
 
1010
   uint16 etc1_block::pack_delta3(int r, int g, int b)
 
1011
   {
 
1012
      RG_ETC1_ASSERT((r >= cETC1ColorDeltaMin) && (r <= cETC1ColorDeltaMax));
 
1013
      RG_ETC1_ASSERT((g >= cETC1ColorDeltaMin) && (g <= cETC1ColorDeltaMax));
 
1014
      RG_ETC1_ASSERT((b >= cETC1ColorDeltaMin) && (b <= cETC1ColorDeltaMax));
 
1015
      if (r < 0) r += 8;
 
1016
      if (g < 0) g += 8;
 
1017
      if (b < 0) b += 8;
 
1018
      return static_cast<uint16>(b | (g << 3) | (r << 6));
 
1019
   }
 
1020
   
 
1021
   void etc1_block::unpack_delta3(int& r, int& g, int& b, uint16 packed_delta3)
 
1022
   {
 
1023
      r = (packed_delta3 >> 6) & 7;
 
1024
      g = (packed_delta3 >> 3) & 7;
 
1025
      b = packed_delta3 & 7;
 
1026
      if (r >= 4) r -= 8;
 
1027
      if (g >= 4) g -= 8;
 
1028
      if (b >= 4) b -= 8;
 
1029
   }
 
1030
 
 
1031
   uint16 etc1_block::pack_color4(const color_quad_u8& color, bool scaled, uint bias)
 
1032
   {
 
1033
      return pack_color4(color.r, color.g, color.b, scaled, bias);
 
1034
   }
 
1035
   
 
1036
   uint16 etc1_block::pack_color4(uint r, uint g, uint b, bool scaled, uint bias)
 
1037
   {
 
1038
      if (scaled)
 
1039
      {
 
1040
         r = (r * 15U + bias) / 255U;
 
1041
         g = (g * 15U + bias) / 255U;
 
1042
         b = (b * 15U + bias) / 255U;
 
1043
      }
 
1044
 
 
1045
      r = rg_etc1::minimum(r, 15U);
 
1046
      g = rg_etc1::minimum(g, 15U);
 
1047
      b = rg_etc1::minimum(b, 15U);
 
1048
 
 
1049
      return static_cast<uint16>(b | (g << 4U) | (r << 8U));
 
1050
   }
 
1051
 
 
1052
   color_quad_u8 etc1_block::unpack_color4(uint16 packed_color4, bool scaled, uint alpha)
 
1053
   {
 
1054
      uint b = packed_color4 & 15U;
 
1055
      uint g = (packed_color4 >> 4U) & 15U;
 
1056
      uint r = (packed_color4 >> 8U) & 15U;
 
1057
 
 
1058
      if (scaled)
 
1059
      {
 
1060
         b = (b << 4U) | b;
 
1061
         g = (g << 4U) | g;
 
1062
         r = (r << 4U) | r;
 
1063
      }
 
1064
 
 
1065
      return color_quad_u8(cNoClamp, r, g, b, rg_etc1::minimum(alpha, 255U));
 
1066
   }
 
1067
   
 
1068
   void etc1_block::unpack_color4(uint& r, uint& g, uint& b, uint16 packed_color4, bool scaled)
 
1069
   {
 
1070
      color_quad_u8 c(unpack_color4(packed_color4, scaled, 0));
 
1071
      r = c.r;
 
1072
      g = c.g;
 
1073
      b = c.b;
 
1074
   }
 
1075
 
 
1076
   void etc1_block::get_diff_subblock_colors(color_quad_u8* pDst, uint16 packed_color5, uint table_idx)
 
1077
   {
 
1078
      RG_ETC1_ASSERT(table_idx < cETC1IntenModifierValues);
 
1079
      const int *pInten_modifer_table = &g_etc1_inten_tables[table_idx][0];
 
1080
 
 
1081
      uint r, g, b;
 
1082
      unpack_color5(r, g, b, packed_color5, true);
 
1083
 
 
1084
      const int ir = static_cast<int>(r), ig = static_cast<int>(g), ib = static_cast<int>(b);
 
1085
 
 
1086
      const int y0 = pInten_modifer_table[0];
 
1087
      pDst[0].set(ir + y0, ig + y0, ib + y0);
 
1088
 
 
1089
      const int y1 = pInten_modifer_table[1];
 
1090
      pDst[1].set(ir + y1, ig + y1, ib + y1);
 
1091
 
 
1092
      const int y2 = pInten_modifer_table[2];
 
1093
      pDst[2].set(ir + y2, ig + y2, ib + y2);
 
1094
 
 
1095
      const int y3 = pInten_modifer_table[3];
 
1096
      pDst[3].set(ir + y3, ig + y3, ib + y3);
 
1097
   }
 
1098
   
 
1099
   bool etc1_block::get_diff_subblock_colors(color_quad_u8* pDst, uint16 packed_color5, uint16 packed_delta3, uint table_idx)
 
1100
   {
 
1101
      RG_ETC1_ASSERT(table_idx < cETC1IntenModifierValues);
 
1102
      const int *pInten_modifer_table = &g_etc1_inten_tables[table_idx][0];
 
1103
 
 
1104
      uint r, g, b;
 
1105
      bool success = unpack_color5(r, g, b, packed_color5, packed_delta3, true);
 
1106
 
 
1107
      const int ir = static_cast<int>(r), ig = static_cast<int>(g), ib = static_cast<int>(b);
 
1108
 
 
1109
      const int y0 = pInten_modifer_table[0];
 
1110
      pDst[0].set(ir + y0, ig + y0, ib + y0);
 
1111
 
 
1112
      const int y1 = pInten_modifer_table[1];
 
1113
      pDst[1].set(ir + y1, ig + y1, ib + y1);
 
1114
 
 
1115
      const int y2 = pInten_modifer_table[2];
 
1116
      pDst[2].set(ir + y2, ig + y2, ib + y2);
 
1117
 
 
1118
      const int y3 = pInten_modifer_table[3];
 
1119
      pDst[3].set(ir + y3, ig + y3, ib + y3);
 
1120
 
 
1121
      return success;
 
1122
   }
 
1123
   
 
1124
   void etc1_block::get_abs_subblock_colors(color_quad_u8* pDst, uint16 packed_color4, uint table_idx)
 
1125
   {
 
1126
      RG_ETC1_ASSERT(table_idx < cETC1IntenModifierValues);
 
1127
      const int *pInten_modifer_table = &g_etc1_inten_tables[table_idx][0];
 
1128
 
 
1129
      uint r, g, b;
 
1130
      unpack_color4(r, g, b, packed_color4, true);
 
1131
      
 
1132
      const int ir = static_cast<int>(r), ig = static_cast<int>(g), ib = static_cast<int>(b);
 
1133
 
 
1134
      const int y0 = pInten_modifer_table[0];
 
1135
      pDst[0].set(ir + y0, ig + y0, ib + y0);
 
1136
      
 
1137
      const int y1 = pInten_modifer_table[1];
 
1138
      pDst[1].set(ir + y1, ig + y1, ib + y1);
 
1139
 
 
1140
      const int y2 = pInten_modifer_table[2];
 
1141
      pDst[2].set(ir + y2, ig + y2, ib + y2);
 
1142
 
 
1143
      const int y3 = pInten_modifer_table[3];
 
1144
      pDst[3].set(ir + y3, ig + y3, ib + y3);
 
1145
   }
 
1146
      
 
1147
   bool unpack_etc1_block(const void* pETC1_block, unsigned int* pDst_pixels_rgba, int stride, bool preserve_alpha)
 
1148
   {
 
1149
      color_quad_u8* pDst = reinterpret_cast<color_quad_u8*>(pDst_pixels_rgba);
 
1150
      const etc1_block& block = *static_cast<const etc1_block*>(pETC1_block);
 
1151
 
 
1152
      const bool diff_flag = block.get_diff_bit();
 
1153
      const bool flip_flag = block.get_flip_bit();
 
1154
      const uint table_index0 = block.get_inten_table(0);
 
1155
      const uint table_index1 = block.get_inten_table(1);
 
1156
 
 
1157
      color_quad_u8 subblock_colors0[4];
 
1158
      color_quad_u8 subblock_colors1[4];
 
1159
      bool success = true;
 
1160
 
 
1161
      if (diff_flag)
 
1162
      {
 
1163
         const uint16 base_color5 = block.get_base5_color();
 
1164
         const uint16 delta_color3 = block.get_delta3_color();
 
1165
         etc1_block::get_diff_subblock_colors(subblock_colors0, base_color5, table_index0);
 
1166
            
 
1167
         if (!etc1_block::get_diff_subblock_colors(subblock_colors1, base_color5, delta_color3, table_index1))
 
1168
            success = false;
 
1169
      }
 
1170
      else
 
1171
      {
 
1172
         const uint16 base_color4_0 = block.get_base4_color(0);
 
1173
         etc1_block::get_abs_subblock_colors(subblock_colors0, base_color4_0, table_index0);
 
1174
 
 
1175
         const uint16 base_color4_1 = block.get_base4_color(1);
 
1176
         etc1_block::get_abs_subblock_colors(subblock_colors1, base_color4_1, table_index1);
 
1177
      }
 
1178
 
 
1179
      if (preserve_alpha)
 
1180
      {
 
1181
         if (flip_flag)
 
1182
         {
 
1183
            for (uint y = 0; y < 2; y++)
 
1184
            {
 
1185
               pDst[0].set_rgb(subblock_colors0[block.get_selector(0, y)]);
 
1186
               pDst[1].set_rgb(subblock_colors0[block.get_selector(1, y)]);
 
1187
               pDst[2].set_rgb(subblock_colors0[block.get_selector(2, y)]);
 
1188
               pDst[3].set_rgb(subblock_colors0[block.get_selector(3, y)]);
 
1189
               pDst += stride;
 
1190
            }
 
1191
 
 
1192
            for (uint y = 2; y < 4; y++)
 
1193
            {
 
1194
               pDst[0].set_rgb(subblock_colors1[block.get_selector(0, y)]);
 
1195
               pDst[1].set_rgb(subblock_colors1[block.get_selector(1, y)]);
 
1196
               pDst[2].set_rgb(subblock_colors1[block.get_selector(2, y)]);
 
1197
               pDst[3].set_rgb(subblock_colors1[block.get_selector(3, y)]);
 
1198
               pDst += stride;
 
1199
            }
 
1200
         }
 
1201
         else
 
1202
         {
 
1203
            for (uint y = 0; y < 4; y++)
 
1204
            {
 
1205
               pDst[0].set_rgb(subblock_colors0[block.get_selector(0, y)]);
 
1206
               pDst[1].set_rgb(subblock_colors0[block.get_selector(1, y)]);
 
1207
               pDst[2].set_rgb(subblock_colors1[block.get_selector(2, y)]);
 
1208
               pDst[3].set_rgb(subblock_colors1[block.get_selector(3, y)]);
 
1209
               pDst += stride;
 
1210
            }
 
1211
         }
 
1212
      }
 
1213
      else 
 
1214
      {
 
1215
         if (flip_flag)
 
1216
         {
 
1217
            // 0000
 
1218
            // 0000
 
1219
            // 1111
 
1220
            // 1111
 
1221
            for (uint y = 0; y < 2; y++)
 
1222
            {
 
1223
               pDst[0] = subblock_colors0[block.get_selector(0, y)];
 
1224
               pDst[1] = subblock_colors0[block.get_selector(1, y)];
 
1225
               pDst[2] = subblock_colors0[block.get_selector(2, y)];
 
1226
               pDst[3] = subblock_colors0[block.get_selector(3, y)];
 
1227
               pDst += stride;
 
1228
            }
 
1229
 
 
1230
            for (uint y = 2; y < 4; y++)
 
1231
            {
 
1232
               pDst[0] = subblock_colors1[block.get_selector(0, y)];
 
1233
               pDst[1] = subblock_colors1[block.get_selector(1, y)];
 
1234
               pDst[2] = subblock_colors1[block.get_selector(2, y)];
 
1235
               pDst[3] = subblock_colors1[block.get_selector(3, y)];
 
1236
               pDst += stride;
 
1237
            }
 
1238
         }
 
1239
         else
 
1240
         {
 
1241
            // 0011
 
1242
            // 0011
 
1243
            // 0011
 
1244
            // 0011
 
1245
            for (uint y = 0; y < 4; y++)
 
1246
            {
 
1247
               pDst[0] = subblock_colors0[block.get_selector(0, y)];
 
1248
               pDst[1] = subblock_colors0[block.get_selector(1, y)];
 
1249
               pDst[2] = subblock_colors1[block.get_selector(2, y)];
 
1250
               pDst[3] = subblock_colors1[block.get_selector(3, y)];
 
1251
               pDst += stride;
 
1252
            }
 
1253
         }
 
1254
      }
 
1255
      
 
1256
      return success;
 
1257
   }
 
1258
 
 
1259
   struct etc1_solution_coordinates
 
1260
   {
 
1261
      inline etc1_solution_coordinates() :
 
1262
      m_unscaled_color(0, 0, 0, 0),
 
1263
         m_inten_table(0),
 
1264
         m_color4(false)
 
1265
      {
 
1266
      }
 
1267
 
 
1268
      inline etc1_solution_coordinates(uint r, uint g, uint b, uint inten_table, bool color4) : 
 
1269
      m_unscaled_color(r, g, b, 255),
 
1270
         m_inten_table(inten_table),
 
1271
         m_color4(color4)
 
1272
      {
 
1273
      }
 
1274
 
 
1275
      inline etc1_solution_coordinates(const color_quad_u8& c, uint inten_table, bool color4) : 
 
1276
      m_unscaled_color(c),
 
1277
         m_inten_table(inten_table),
 
1278
         m_color4(color4)
 
1279
      {
 
1280
      }
 
1281
 
 
1282
      inline etc1_solution_coordinates(const etc1_solution_coordinates& other)
 
1283
      {
 
1284
         *this = other;
 
1285
      }
 
1286
 
 
1287
      inline etc1_solution_coordinates& operator= (const etc1_solution_coordinates& rhs)
 
1288
      {
 
1289
         m_unscaled_color = rhs.m_unscaled_color;
 
1290
         m_inten_table = rhs.m_inten_table;
 
1291
         m_color4 = rhs.m_color4;
 
1292
         return *this;
 
1293
      }
 
1294
 
 
1295
      inline void clear()
 
1296
      {
 
1297
         m_unscaled_color.clear();
 
1298
         m_inten_table = 0;
 
1299
         m_color4 = false;
 
1300
      }
 
1301
 
 
1302
      inline color_quad_u8 get_scaled_color() const
 
1303
      {
 
1304
         int br, bg, bb;
 
1305
         if (m_color4)
 
1306
         {
 
1307
            br = m_unscaled_color.r | (m_unscaled_color.r << 4);
 
1308
            bg = m_unscaled_color.g | (m_unscaled_color.g << 4);
 
1309
            bb = m_unscaled_color.b | (m_unscaled_color.b << 4);
 
1310
         }
 
1311
         else
 
1312
         {
 
1313
            br = (m_unscaled_color.r >> 2) | (m_unscaled_color.r << 3);
 
1314
            bg = (m_unscaled_color.g >> 2) | (m_unscaled_color.g << 3);
 
1315
            bb = (m_unscaled_color.b >> 2) | (m_unscaled_color.b << 3);
 
1316
         }
 
1317
         return color_quad_u8(br, bg, bb);
 
1318
      }
 
1319
 
 
1320
      inline void get_block_colors(color_quad_u8* pBlock_colors)
 
1321
      {
 
1322
         int br, bg, bb;
 
1323
         if (m_color4)
 
1324
         {
 
1325
            br = m_unscaled_color.r | (m_unscaled_color.r << 4);
 
1326
            bg = m_unscaled_color.g | (m_unscaled_color.g << 4);
 
1327
            bb = m_unscaled_color.b | (m_unscaled_color.b << 4);
 
1328
         }
 
1329
         else
 
1330
         {
 
1331
            br = (m_unscaled_color.r >> 2) | (m_unscaled_color.r << 3);
 
1332
            bg = (m_unscaled_color.g >> 2) | (m_unscaled_color.g << 3);
 
1333
            bb = (m_unscaled_color.b >> 2) | (m_unscaled_color.b << 3);
 
1334
         }
 
1335
         const int* pInten_table = g_etc1_inten_tables[m_inten_table];
 
1336
         pBlock_colors[0].set(br + pInten_table[0], bg + pInten_table[0], bb + pInten_table[0]);
 
1337
         pBlock_colors[1].set(br + pInten_table[1], bg + pInten_table[1], bb + pInten_table[1]);
 
1338
         pBlock_colors[2].set(br + pInten_table[2], bg + pInten_table[2], bb + pInten_table[2]);
 
1339
         pBlock_colors[3].set(br + pInten_table[3], bg + pInten_table[3], bb + pInten_table[3]);
 
1340
      }
 
1341
 
 
1342
      color_quad_u8 m_unscaled_color;
 
1343
      uint m_inten_table;
 
1344
      bool m_color4;
 
1345
   };
 
1346
 
 
1347
   class etc1_optimizer
 
1348
   {
 
1349
      etc1_optimizer(const etc1_optimizer&);
 
1350
      etc1_optimizer& operator= (const etc1_optimizer&);
 
1351
 
 
1352
   public:
 
1353
      etc1_optimizer()
 
1354
      {
 
1355
         clear();
 
1356
      }
 
1357
 
 
1358
      void clear()
 
1359
      {
 
1360
         m_pParams = NULL;
 
1361
         m_pResult = NULL;
 
1362
         m_pSorted_luma = NULL;
 
1363
         m_pSorted_luma_indices = NULL;
 
1364
      }
 
1365
 
 
1366
      struct params : etc1_pack_params
 
1367
      {
 
1368
         params()
 
1369
         {
 
1370
            clear();
 
1371
         }
 
1372
 
 
1373
         params(const etc1_pack_params& base_params) : 
 
1374
         etc1_pack_params(base_params)
 
1375
         {
 
1376
            clear_optimizer_params();
 
1377
         }
 
1378
 
 
1379
         void clear()
 
1380
         {
 
1381
            etc1_pack_params::clear();
 
1382
            clear_optimizer_params();
 
1383
         }
 
1384
 
 
1385
         void clear_optimizer_params()
 
1386
         {
 
1387
            m_num_src_pixels = 0;
 
1388
            m_pSrc_pixels = 0;
 
1389
 
 
1390
            m_use_color4 = false;
 
1391
            static const int s_default_scan_delta[] = { 0 };
 
1392
            m_pScan_deltas = s_default_scan_delta;
 
1393
            m_scan_delta_size = 1;
 
1394
 
 
1395
            m_base_color5.clear();
 
1396
            m_constrain_against_base_color5 = false;
 
1397
         }
 
1398
 
 
1399
         uint m_num_src_pixels;
 
1400
         const color_quad_u8* m_pSrc_pixels;
 
1401
 
 
1402
         bool m_use_color4;
 
1403
         const int* m_pScan_deltas;
 
1404
         uint m_scan_delta_size;
 
1405
 
 
1406
         color_quad_u8 m_base_color5;
 
1407
         bool m_constrain_against_base_color5;
 
1408
      };
 
1409
 
 
1410
      struct results
 
1411
      {
 
1412
         uint64 m_error;
 
1413
         color_quad_u8 m_block_color_unscaled;
 
1414
         uint m_block_inten_table;
 
1415
         uint m_n;
 
1416
         uint8* m_pSelectors;
 
1417
         bool m_block_color4;
 
1418
 
 
1419
         inline results& operator= (const results& rhs)
 
1420
         {
 
1421
            m_block_color_unscaled = rhs.m_block_color_unscaled;
 
1422
            m_block_color4 = rhs.m_block_color4;
 
1423
            m_block_inten_table = rhs.m_block_inten_table;
 
1424
            m_error = rhs.m_error;
 
1425
            RG_ETC1_ASSERT(m_n == rhs.m_n);
 
1426
            memcpy(m_pSelectors, rhs.m_pSelectors, rhs.m_n);
 
1427
            return *this;
 
1428
         }
 
1429
      };
 
1430
 
 
1431
      void init(const params& params, results& result);
 
1432
      bool compute();
 
1433
 
 
1434
   private:      
 
1435
      struct potential_solution
 
1436
      {
 
1437
         potential_solution() : m_coords(), m_error(cUINT64_MAX), m_valid(false)
 
1438
         {
 
1439
         }
 
1440
 
 
1441
         etc1_solution_coordinates  m_coords;
 
1442
         uint8                      m_selectors[8];
 
1443
         uint64                     m_error;
 
1444
         bool                       m_valid;
 
1445
 
 
1446
         void clear()
 
1447
         {
 
1448
            m_coords.clear();
 
1449
            m_error = cUINT64_MAX;
 
1450
            m_valid = false;
 
1451
         }
 
1452
      };
 
1453
 
 
1454
      const params* m_pParams;
 
1455
      results* m_pResult;
 
1456
 
 
1457
      int m_limit;
 
1458
 
 
1459
      vec3F m_avg_color;
 
1460
      int m_br, m_bg, m_bb;
 
1461
      uint16 m_luma[8];
 
1462
      uint32 m_sorted_luma[2][8];
 
1463
      const uint32* m_pSorted_luma_indices;
 
1464
      uint32* m_pSorted_luma;
 
1465
 
 
1466
      uint8 m_selectors[8];
 
1467
      uint8 m_best_selectors[8];
 
1468
 
 
1469
      potential_solution m_best_solution;
 
1470
      potential_solution m_trial_solution;
 
1471
      uint8 m_temp_selectors[8];
 
1472
 
 
1473
      bool evaluate_solution(const etc1_solution_coordinates& coords, potential_solution& trial_solution, potential_solution* pBest_solution);
 
1474
      bool evaluate_solution_fast(const etc1_solution_coordinates& coords, potential_solution& trial_solution, potential_solution* pBest_solution);
 
1475
   };
 
1476
      
 
1477
   bool etc1_optimizer::compute()
 
1478
   {
 
1479
      const uint n = m_pParams->m_num_src_pixels;
 
1480
      const int scan_delta_size = m_pParams->m_scan_delta_size;
 
1481
      
 
1482
      // Scan through a subset of the 3D lattice centered around the avg block color trying each 3D (555 or 444) lattice point as a potential block color.
 
1483
      // Each time a better solution is found try to refine the current solution's block color based of the current selectors and intensity table index.
 
1484
      for (int zdi = 0; zdi < scan_delta_size; zdi++)
 
1485
      {
 
1486
         const int zd = m_pParams->m_pScan_deltas[zdi];
 
1487
         const int mbb = m_bb + zd;
 
1488
         if (mbb < 0) continue; else if (mbb > m_limit) break;
 
1489
         
 
1490
         for (int ydi = 0; ydi < scan_delta_size; ydi++)
 
1491
         {
 
1492
            const int yd = m_pParams->m_pScan_deltas[ydi];
 
1493
            const int mbg = m_bg + yd;
 
1494
            if (mbg < 0) continue; else if (mbg > m_limit) break;
 
1495
 
 
1496
            for (int xdi = 0; xdi < scan_delta_size; xdi++)
 
1497
            {
 
1498
               const int xd = m_pParams->m_pScan_deltas[xdi];
 
1499
               const int mbr = m_br + xd;
 
1500
               if (mbr < 0) continue; else if (mbr > m_limit) break;
 
1501
      
 
1502
               etc1_solution_coordinates coords(mbr, mbg, mbb, 0, m_pParams->m_use_color4);
 
1503
               if (m_pParams->m_quality == cHighQuality)
 
1504
               {
 
1505
                  if (!evaluate_solution(coords, m_trial_solution, &m_best_solution))
 
1506
                     continue;
 
1507
               }
 
1508
               else
 
1509
               {
 
1510
                  if (!evaluate_solution_fast(coords, m_trial_solution, &m_best_solution))
 
1511
                     continue;
 
1512
               }
 
1513
               
 
1514
               // Now we have the input block, the avg. color of the input pixels, a set of trial selector indices, and the block color+intensity index.
 
1515
               // Now, for each component, attempt to refine the current solution by solving a simple linear equation. For example, for 4 colors:
 
1516
               // The goal is:
 
1517
               // pixel0 - (block_color+inten_table[selector0]) + pixel1 - (block_color+inten_table[selector1]) + pixel2 - (block_color+inten_table[selector2]) + pixel3 - (block_color+inten_table[selector3]) = 0
 
1518
               // Rearranging this:
 
1519
               // (pixel0 + pixel1 + pixel2 + pixel3) - (block_color+inten_table[selector0]) - (block_color+inten_table[selector1]) - (block_color+inten_table[selector2]) - (block_color+inten_table[selector3]) = 0
 
1520
               // (pixel0 + pixel1 + pixel2 + pixel3) - block_color - inten_table[selector0] - block_color-inten_table[selector1] - block_color-inten_table[selector2] - block_color-inten_table[selector3] = 0
 
1521
               // (pixel0 + pixel1 + pixel2 + pixel3) - 4*block_color - inten_table[selector0] - inten_table[selector1] - inten_table[selector2] - inten_table[selector3] = 0
 
1522
               // (pixel0 + pixel1 + pixel2 + pixel3) - 4*block_color - (inten_table[selector0] + inten_table[selector1] + inten_table[selector2] + inten_table[selector3]) = 0
 
1523
               // (pixel0 + pixel1 + pixel2 + pixel3)/4 - block_color - (inten_table[selector0] + inten_table[selector1] + inten_table[selector2] + inten_table[selector3])/4 = 0
 
1524
               // block_color = (pixel0 + pixel1 + pixel2 + pixel3)/4 - (inten_table[selector0] + inten_table[selector1] + inten_table[selector2] + inten_table[selector3])/4
 
1525
               // So what this means:
 
1526
               // optimal_block_color = avg_input - avg_inten_delta
 
1527
               // So the optimal block color can be computed by taking the average block color and subtracting the current average of the intensity delta.
 
1528
               // Unfortunately, optimal_block_color must then be quantized to 555 or 444 so it's not always possible to improve matters using this formula.
 
1529
               // Also, the above formula is for unclamped intensity deltas. The actual implementation takes into account clamping.
 
1530
 
 
1531
               const uint max_refinement_trials = (m_pParams->m_quality == cLowQuality) ? 2 : (((xd | yd | zd) == 0) ? 4 : 2);
 
1532
               for (uint refinement_trial = 0; refinement_trial < max_refinement_trials; refinement_trial++)
 
1533
               {
 
1534
                  const uint8* pSelectors = m_best_solution.m_selectors;
 
1535
                  const int* pInten_table = g_etc1_inten_tables[m_best_solution.m_coords.m_inten_table];
 
1536
 
 
1537
                  int delta_sum_r = 0, delta_sum_g = 0, delta_sum_b = 0;
 
1538
                  const color_quad_u8 base_color(m_best_solution.m_coords.get_scaled_color());
 
1539
                  for (uint r = 0; r < n; r++)
 
1540
                  {
 
1541
                     const uint s = *pSelectors++;
 
1542
                     const int yd = pInten_table[s];
 
1543
                     // Compute actual delta being applied to each pixel, taking into account clamping.
 
1544
                     delta_sum_r += rg_etc1::clamp<int>(base_color.r + yd, 0, 255) - base_color.r;
 
1545
                     delta_sum_g += rg_etc1::clamp<int>(base_color.g + yd, 0, 255) - base_color.g;
 
1546
                     delta_sum_b += rg_etc1::clamp<int>(base_color.b + yd, 0, 255) - base_color.b;
 
1547
                  }
 
1548
                  if ((!delta_sum_r) && (!delta_sum_g) && (!delta_sum_b))
 
1549
                     break;
 
1550
                  const float avg_delta_r_f = static_cast<float>(delta_sum_r) / n;
 
1551
                  const float avg_delta_g_f = static_cast<float>(delta_sum_g) / n;
 
1552
                  const float avg_delta_b_f = static_cast<float>(delta_sum_b) / n;
 
1553
                  const int br1 = rg_etc1::clamp<int>(static_cast<uint>((m_avg_color[0] - avg_delta_r_f) * m_limit / 255.0f + .5f), 0, m_limit);
 
1554
                  const int bg1 = rg_etc1::clamp<int>(static_cast<uint>((m_avg_color[1] - avg_delta_g_f) * m_limit / 255.0f + .5f), 0, m_limit);
 
1555
                  const int bb1 = rg_etc1::clamp<int>(static_cast<uint>((m_avg_color[2] - avg_delta_b_f) * m_limit / 255.0f + .5f), 0, m_limit);
 
1556
                  
 
1557
                  bool skip = false;
 
1558
                  
 
1559
                  if ((mbr == br1) && (mbg == bg1) && (mbb == bb1))
 
1560
                     skip = true;
 
1561
                  else if ((br1 == m_best_solution.m_coords.m_unscaled_color.r) && (bg1 == m_best_solution.m_coords.m_unscaled_color.g) && (bb1 == m_best_solution.m_coords.m_unscaled_color.b))
 
1562
                     skip = true;
 
1563
                  else if ((m_br == br1) && (m_bg == bg1) && (m_bb == bb1))
 
1564
                     skip = true;
 
1565
 
 
1566
                  if (skip)
 
1567
                     break;
 
1568
 
 
1569
                  etc1_solution_coordinates coords1(br1, bg1, bb1, 0, m_pParams->m_use_color4);
 
1570
                  if (m_pParams->m_quality == cHighQuality)
 
1571
                  {
 
1572
                     if (!evaluate_solution(coords1, m_trial_solution, &m_best_solution)) 
 
1573
                        break;
 
1574
                  }
 
1575
                  else
 
1576
                  {
 
1577
                     if (!evaluate_solution_fast(coords1, m_trial_solution, &m_best_solution))
 
1578
                        break;
 
1579
                  }
 
1580
 
 
1581
               }  // refinement_trial
 
1582
 
 
1583
            } // xdi
 
1584
         } // ydi
 
1585
      } // zdi
 
1586
 
 
1587
      if (!m_best_solution.m_valid)
 
1588
      {
 
1589
         m_pResult->m_error = cUINT32_MAX;
 
1590
         return false;
 
1591
      }
 
1592
      
 
1593
      const uint8* pSelectors = m_best_solution.m_selectors;
 
1594
 
 
1595
#ifdef RG_ETC1_BUILD_DEBUG
 
1596
      {
 
1597
         color_quad_u8 block_colors[4];
 
1598
         m_best_solution.m_coords.get_block_colors(block_colors);
 
1599
 
 
1600
         const color_quad_u8* pSrc_pixels = m_pParams->m_pSrc_pixels;
 
1601
         uint64 actual_error = 0;
 
1602
         for (uint i = 0; i < n; i++)
 
1603
            actual_error += pSrc_pixels[i].squared_distance_rgb(block_colors[pSelectors[i]]);
 
1604
         
 
1605
         RG_ETC1_ASSERT(actual_error == m_best_solution.m_error);
 
1606
      }
 
1607
#endif      
 
1608
      
 
1609
      m_pResult->m_error = m_best_solution.m_error;
 
1610
 
 
1611
      m_pResult->m_block_color_unscaled = m_best_solution.m_coords.m_unscaled_color;
 
1612
      m_pResult->m_block_color4 = m_best_solution.m_coords.m_color4;
 
1613
      
 
1614
      m_pResult->m_block_inten_table = m_best_solution.m_coords.m_inten_table;
 
1615
      memcpy(m_pResult->m_pSelectors, pSelectors, n);
 
1616
      m_pResult->m_n = n;
 
1617
 
 
1618
      return true;
 
1619
   }
 
1620
 
 
1621
   void etc1_optimizer::init(const params& p, results& r)
 
1622
   {
 
1623
      // This version is hardcoded for 8 pixel subblocks.
 
1624
      RG_ETC1_ASSERT(p.m_num_src_pixels == 8);
 
1625
      
 
1626
      m_pParams = &p;
 
1627
      m_pResult = &r;
 
1628
                  
 
1629
      const uint n = 8;
 
1630
      
 
1631
      m_limit = m_pParams->m_use_color4 ? 15 : 31;
 
1632
 
 
1633
      vec3F avg_color(0.0f);
 
1634
 
 
1635
      for (uint i = 0; i < n; i++)
 
1636
      {
 
1637
         const color_quad_u8& c = m_pParams->m_pSrc_pixels[i];
 
1638
         const vec3F fc(c.r, c.g, c.b);
 
1639
 
 
1640
         avg_color += fc;
 
1641
 
 
1642
         m_luma[i] = static_cast<uint16>(c.r + c.g + c.b);
 
1643
         m_sorted_luma[0][i] = i;
 
1644
      }
 
1645
      avg_color *= (1.0f / static_cast<float>(n));
 
1646
      m_avg_color = avg_color;
 
1647
 
 
1648
      m_br = rg_etc1::clamp<int>(static_cast<uint>(m_avg_color[0] * m_limit / 255.0f + .5f), 0, m_limit);
 
1649
      m_bg = rg_etc1::clamp<int>(static_cast<uint>(m_avg_color[1] * m_limit / 255.0f + .5f), 0, m_limit);
 
1650
      m_bb = rg_etc1::clamp<int>(static_cast<uint>(m_avg_color[2] * m_limit / 255.0f + .5f), 0, m_limit);
 
1651
 
 
1652
      if (m_pParams->m_quality <= cMediumQuality)
 
1653
      {
 
1654
         m_pSorted_luma_indices = indirect_radix_sort(n, m_sorted_luma[0], m_sorted_luma[1], m_luma, 0, sizeof(m_luma[0]), false);
 
1655
         m_pSorted_luma = m_sorted_luma[0];
 
1656
         if (m_pSorted_luma_indices == m_sorted_luma[0])
 
1657
            m_pSorted_luma = m_sorted_luma[1];
 
1658
      
 
1659
         for (uint i = 0; i < n; i++)
 
1660
            m_pSorted_luma[i] = m_luma[m_pSorted_luma_indices[i]];
 
1661
      }
 
1662
      
 
1663
      m_best_solution.m_coords.clear();
 
1664
      m_best_solution.m_valid = false;
 
1665
      m_best_solution.m_error = cUINT64_MAX;
 
1666
   }
 
1667
 
 
1668
   bool etc1_optimizer::evaluate_solution(const etc1_solution_coordinates& coords, potential_solution& trial_solution, potential_solution* pBest_solution)
 
1669
   {
 
1670
      trial_solution.m_valid = false;
 
1671
 
 
1672
      if (m_pParams->m_constrain_against_base_color5)
 
1673
      {
 
1674
         const int dr = coords.m_unscaled_color.r - m_pParams->m_base_color5.r;
 
1675
         const int dg = coords.m_unscaled_color.g - m_pParams->m_base_color5.g;
 
1676
         const int db = coords.m_unscaled_color.b - m_pParams->m_base_color5.b;
 
1677
 
 
1678
         if ((rg_etc1::minimum(dr, dg, db) < cETC1ColorDeltaMin) || (rg_etc1::maximum(dr, dg, db) > cETC1ColorDeltaMax))
 
1679
            return false;
 
1680
      }
 
1681
 
 
1682
      const color_quad_u8 base_color(coords.get_scaled_color());
 
1683
      
 
1684
      const uint n = 8;
 
1685
            
 
1686
      trial_solution.m_error = cUINT64_MAX;
 
1687
            
 
1688
      for (uint inten_table = 0; inten_table < cETC1IntenModifierValues; inten_table++)
 
1689
      {
 
1690
         const int* pInten_table = g_etc1_inten_tables[inten_table];
 
1691
 
 
1692
         color_quad_u8 block_colors[4];
 
1693
         for (uint s = 0; s < 4; s++)
 
1694
         {
 
1695
            const int yd = pInten_table[s];
 
1696
            block_colors[s].set(base_color.r + yd, base_color.g + yd, base_color.b + yd, 0);
 
1697
         }
 
1698
         
 
1699
         uint64 total_error = 0;
 
1700
         
 
1701
         const color_quad_u8* pSrc_pixels = m_pParams->m_pSrc_pixels;
 
1702
         for (uint c = 0; c < n; c++)
 
1703
         {
 
1704
            const color_quad_u8& src_pixel = *pSrc_pixels++;
 
1705
            
 
1706
            uint best_selector_index = 0;
 
1707
            uint best_error = rg_etc1::square(src_pixel.r - block_colors[0].r) + rg_etc1::square(src_pixel.g - block_colors[0].g) + rg_etc1::square(src_pixel.b - block_colors[0].b);
 
1708
 
 
1709
            uint trial_error = rg_etc1::square(src_pixel.r - block_colors[1].r) + rg_etc1::square(src_pixel.g - block_colors[1].g) + rg_etc1::square(src_pixel.b - block_colors[1].b);
 
1710
            if (trial_error < best_error)
 
1711
            {
 
1712
               best_error = trial_error;
 
1713
               best_selector_index = 1;
 
1714
            }
 
1715
 
 
1716
            trial_error = rg_etc1::square(src_pixel.r - block_colors[2].r) + rg_etc1::square(src_pixel.g - block_colors[2].g) + rg_etc1::square(src_pixel.b - block_colors[2].b);
 
1717
            if (trial_error < best_error)
 
1718
            {
 
1719
               best_error = trial_error;
 
1720
               best_selector_index = 2;
 
1721
            }
 
1722
 
 
1723
            trial_error = rg_etc1::square(src_pixel.r - block_colors[3].r) + rg_etc1::square(src_pixel.g - block_colors[3].g) + rg_etc1::square(src_pixel.b - block_colors[3].b);
 
1724
            if (trial_error < best_error)
 
1725
            {
 
1726
               best_error = trial_error;
 
1727
               best_selector_index = 3;
 
1728
            }
 
1729
 
 
1730
            m_temp_selectors[c] = static_cast<uint8>(best_selector_index);
 
1731
 
 
1732
            total_error += best_error;
 
1733
            if (total_error >= trial_solution.m_error)
 
1734
               break;
 
1735
         }
 
1736
         
 
1737
         if (total_error < trial_solution.m_error)
 
1738
         {
 
1739
            trial_solution.m_error = total_error;
 
1740
            trial_solution.m_coords.m_inten_table = inten_table;
 
1741
            memcpy(trial_solution.m_selectors, m_temp_selectors, 8);
 
1742
            trial_solution.m_valid = true;
 
1743
         }
 
1744
      }
 
1745
      trial_solution.m_coords.m_unscaled_color = coords.m_unscaled_color;
 
1746
      trial_solution.m_coords.m_color4 = m_pParams->m_use_color4;
 
1747
 
 
1748
      bool success = false;
 
1749
      if (pBest_solution)
 
1750
      {
 
1751
         if (trial_solution.m_error < pBest_solution->m_error)
 
1752
         {
 
1753
            *pBest_solution = trial_solution;
 
1754
            success = true;
 
1755
         }
 
1756
      }
 
1757
 
 
1758
      return success;
 
1759
   }
 
1760
 
 
1761
   bool etc1_optimizer::evaluate_solution_fast(const etc1_solution_coordinates& coords, potential_solution& trial_solution, potential_solution* pBest_solution)
 
1762
   {
 
1763
      if (m_pParams->m_constrain_against_base_color5)
 
1764
      {
 
1765
         const int dr = coords.m_unscaled_color.r - m_pParams->m_base_color5.r;
 
1766
         const int dg = coords.m_unscaled_color.g - m_pParams->m_base_color5.g;
 
1767
         const int db = coords.m_unscaled_color.b - m_pParams->m_base_color5.b;
 
1768
 
 
1769
         if ((rg_etc1::minimum(dr, dg, db) < cETC1ColorDeltaMin) || (rg_etc1::maximum(dr, dg, db) > cETC1ColorDeltaMax))
 
1770
         {
 
1771
            trial_solution.m_valid = false;
 
1772
            return false;
 
1773
         }
 
1774
      }
 
1775
 
 
1776
      const color_quad_u8 base_color(coords.get_scaled_color());
 
1777
 
 
1778
      const uint n = 8;
 
1779
      
 
1780
      trial_solution.m_error = cUINT64_MAX;
 
1781
 
 
1782
      for (int inten_table = cETC1IntenModifierValues - 1; inten_table >= 0; --inten_table)
 
1783
      {
 
1784
         const int* pInten_table = g_etc1_inten_tables[inten_table];
 
1785
 
 
1786
         uint block_inten[4];
 
1787
         color_quad_u8 block_colors[4];
 
1788
         for (uint s = 0; s < 4; s++)
 
1789
         {
 
1790
            const int yd = pInten_table[s];
 
1791
            color_quad_u8 block_color(base_color.r + yd, base_color.g + yd, base_color.b + yd, 0);
 
1792
            block_colors[s] = block_color;
 
1793
            block_inten[s] = block_color.r + block_color.g + block_color.b;
 
1794
         }
 
1795
 
 
1796
         // evaluate_solution_fast() enforces/assumesd a total ordering of the input colors along the intensity (1,1,1) axis to more quickly classify the inputs to selectors.
 
1797
         // The inputs colors have been presorted along the projection onto this axis, and ETC1 block colors are always ordered along the intensity axis, so this classification is fast.
 
1798
         // 0   1   2   3
 
1799
         //   01  12  23
 
1800
         const uint block_inten_midpoints[3] = { block_inten[0] + block_inten[1], block_inten[1] + block_inten[2], block_inten[2] + block_inten[3] };
 
1801
 
 
1802
         uint64 total_error = 0;
 
1803
         const color_quad_u8* pSrc_pixels = m_pParams->m_pSrc_pixels;
 
1804
         if ((m_pSorted_luma[n - 1] * 2) < block_inten_midpoints[0])
 
1805
         {
 
1806
            if (block_inten[0] > m_pSorted_luma[n - 1])
 
1807
            {
 
1808
               const uint min_error = (uint)(labs(block_inten[0] - m_pSorted_luma[n - 1]));
 
1809
               if (min_error >= trial_solution.m_error)
 
1810
                  continue;
 
1811
            }
 
1812
 
 
1813
            memset(&m_temp_selectors[0], 0, n);
 
1814
 
 
1815
            for (uint c = 0; c < n; c++)
 
1816
               total_error += block_colors[0].squared_distance_rgb(pSrc_pixels[c]);
 
1817
         }
 
1818
         else if ((m_pSorted_luma[0] * 2) >= block_inten_midpoints[2])
 
1819
         {
 
1820
            if (m_pSorted_luma[0] > block_inten[3])
 
1821
            {
 
1822
               const uint min_error = (uint)(labs(m_pSorted_luma[0] - block_inten[3]));
 
1823
               if (min_error >= trial_solution.m_error)
 
1824
                  continue;
 
1825
            }
 
1826
 
 
1827
            memset(&m_temp_selectors[0], 3, n);
 
1828
 
 
1829
            for (uint c = 0; c < n; c++)
 
1830
               total_error += block_colors[3].squared_distance_rgb(pSrc_pixels[c]);
 
1831
         }
 
1832
         else
 
1833
         {
 
1834
            uint cur_selector = 0, c;
 
1835
            for (c = 0; c < n; c++)
 
1836
            {
 
1837
               const uint y = m_pSorted_luma[c];
 
1838
               while ((y * 2) >= block_inten_midpoints[cur_selector])
 
1839
                  if (++cur_selector > 2)
 
1840
                     goto done;
 
1841
               const uint sorted_pixel_index = m_pSorted_luma_indices[c];
 
1842
               m_temp_selectors[sorted_pixel_index] = static_cast<uint8>(cur_selector);
 
1843
               total_error += block_colors[cur_selector].squared_distance_rgb(pSrc_pixels[sorted_pixel_index]);
 
1844
            }
 
1845
done:
 
1846
            while (c < n)
 
1847
            {
 
1848
               const uint sorted_pixel_index = m_pSorted_luma_indices[c];
 
1849
               m_temp_selectors[sorted_pixel_index] = 3;
 
1850
               total_error += block_colors[3].squared_distance_rgb(pSrc_pixels[sorted_pixel_index]);
 
1851
               ++c;
 
1852
            }
 
1853
         }
 
1854
 
 
1855
         if (total_error < trial_solution.m_error)
 
1856
         {
 
1857
            trial_solution.m_error = total_error;
 
1858
            trial_solution.m_coords.m_inten_table = inten_table;
 
1859
            memcpy(trial_solution.m_selectors, m_temp_selectors, n);
 
1860
            trial_solution.m_valid = true;
 
1861
            if (!total_error)
 
1862
               break;
 
1863
         }
 
1864
      }
 
1865
      trial_solution.m_coords.m_unscaled_color = coords.m_unscaled_color;
 
1866
      trial_solution.m_coords.m_color4 = m_pParams->m_use_color4;
 
1867
      
 
1868
      bool success = false;
 
1869
      if (pBest_solution)
 
1870
      {
 
1871
         if (trial_solution.m_error < pBest_solution->m_error)
 
1872
         {
 
1873
            *pBest_solution = trial_solution;
 
1874
            success = true;
 
1875
         }
 
1876
      }
 
1877
 
 
1878
      return success;
 
1879
   }
 
1880
         
 
1881
   static uint etc1_decode_value(uint diff, uint inten, uint selector, uint packed_c)
 
1882
   {
 
1883
      //const uint limit = diff ? 32 : 16; limit;
 
1884
      //RG_ETC1_ASSERT((diff < 2) && (inten < 8) && (selector < 4) && (packed_c < limit));
 
1885
      int c;
 
1886
      if (diff)
 
1887
         c = (packed_c >> 2) | (packed_c << 3);
 
1888
      else 
 
1889
         c = packed_c | (packed_c << 4);
 
1890
      c += g_etc1_inten_tables[inten][selector];
 
1891
      c = rg_etc1::clamp<int>(c, 0, 255);
 
1892
      return c;
 
1893
   }
 
1894
 
 
1895
   static inline int mul_8bit(int a, int b) { int t = a*b + 128; return (t + (t >> 8)) >> 8; }
 
1896
 
 
1897
   void pack_etc1_block_init()
 
1898
   {
 
1899
      for (uint diff = 0; diff < 2; diff++)
 
1900
      {
 
1901
         const uint limit = diff ? 32 : 16;
 
1902
 
 
1903
         for (uint inten = 0; inten < 8; inten++)
 
1904
         {
 
1905
            for (uint selector = 0; selector < 4; selector++)
 
1906
            {
 
1907
               const uint inverse_table_index = diff + (inten << 1) + (selector << 4);
 
1908
               for (uint color = 0; color < 256; color++)
 
1909
               {
 
1910
                  uint best_error = cUINT32_MAX, best_packed_c = 0;
 
1911
                  for (uint packed_c = 0; packed_c < limit; packed_c++)
 
1912
                  {
 
1913
                     int v = etc1_decode_value(diff, inten, selector, packed_c);
 
1914
                     uint err = (uint)labs(v - static_cast<int>(color));
 
1915
                     if (err < best_error)
 
1916
                     {
 
1917
                        best_error = err;
 
1918
                        best_packed_c = packed_c;
 
1919
                        if (!best_error) 
 
1920
                           break;
 
1921
                     }
 
1922
                  }
 
1923
                  RG_ETC1_ASSERT(best_error <= 255);
 
1924
                  g_etc1_inverse_lookup[inverse_table_index][color] = static_cast<uint16>(best_packed_c | (best_error << 8));
 
1925
               }
 
1926
            }
 
1927
         }
 
1928
      }
 
1929
      
 
1930
      uint expand5[32];
 
1931
      for(int i = 0; i < 32; i++)
 
1932
         expand5[i] = (i << 3) | (i >> 2);
 
1933
 
 
1934
      for(int i = 0; i < 256 + 16; i++)
 
1935
      {
 
1936
         int v = clamp<int>(i - 8, 0, 255);
 
1937
         g_quant5_tab[i] = static_cast<uint8>(expand5[mul_8bit(v,31)]);
 
1938
      }
 
1939
   }
 
1940
 
 
1941
   // Packs solid color blocks efficiently using a set of small precomputed tables.
 
1942
   // For random 888 inputs, MSE results are better than Erricson's ETC1 packer in "slow" mode ~9.5% of the time, is slightly worse only ~.01% of the time, and is equal the rest of the time.
 
1943
   static uint64 pack_etc1_block_solid_color(etc1_block& block, const uint8* pColor, etc1_pack_params& pack_params)
 
1944
   {
 
1945
      // pack_params;
 
1946
      RG_ETC1_ASSERT(g_etc1_inverse_lookup[0][255]);
 
1947
            
 
1948
      static uint s_next_comp[4] = { 1, 2, 0, 1 };
 
1949
            
 
1950
      uint best_error = cUINT32_MAX, best_i = 0;
 
1951
      int best_x = 0, best_packed_c1 = 0, best_packed_c2 = 0;
 
1952
 
 
1953
      // For each possible 8-bit value, there is a precomputed list of diff/inten/selector configurations that allow that 8-bit value to be encoded with no error.
 
1954
      for (uint i = 0; i < 3; i++)
 
1955
      {
 
1956
         const uint c1 = pColor[s_next_comp[i]], c2 = pColor[s_next_comp[i + 1]];
 
1957
 
 
1958
         const int delta_range = 1;
 
1959
         for (int delta = -delta_range; delta <= delta_range; delta++)
 
1960
         {
 
1961
            const int c_plus_delta = rg_etc1::clamp<int>(pColor[i] + delta, 0, 255);
 
1962
 
 
1963
            const uint16* pTable;
 
1964
            if (!c_plus_delta)
 
1965
               pTable = g_color8_to_etc_block_config_0_255[0];
 
1966
            else if (c_plus_delta == 255)
 
1967
               pTable = g_color8_to_etc_block_config_0_255[1];
 
1968
            else
 
1969
               pTable = g_color8_to_etc_block_config_1_to_254[c_plus_delta - 1];
 
1970
 
 
1971
            do
 
1972
            {
 
1973
               const uint x = *pTable++;
 
1974
 
 
1975
#ifdef RG_ETC1_BUILD_DEBUG
 
1976
               const uint diff = x & 1;
 
1977
               const uint inten = (x >> 1) & 7;
 
1978
               const uint selector = (x >> 4) & 3;
 
1979
               const uint p0 = (x >> 8) & 255;
 
1980
               RG_ETC1_ASSERT(etc1_decode_value(diff, inten, selector, p0) == (uint)c_plus_delta);
 
1981
#endif
 
1982
 
 
1983
               const uint16* pInverse_table = g_etc1_inverse_lookup[x & 0xFF];
 
1984
               uint16 p1 = pInverse_table[c1];
 
1985
               uint16 p2 = pInverse_table[c2];
 
1986
               const uint trial_error = rg_etc1::square(c_plus_delta - pColor[i]) + rg_etc1::square(p1 >> 8) + rg_etc1::square(p2 >> 8);
 
1987
               if (trial_error < best_error)
 
1988
               {
 
1989
                  best_error = trial_error;
 
1990
                  best_x = x;
 
1991
                  best_packed_c1 = p1 & 0xFF;
 
1992
                  best_packed_c2 = p2 & 0xFF;
 
1993
                  best_i = i;
 
1994
                  if (!best_error)
 
1995
                     goto found_perfect_match;
 
1996
               }
 
1997
            } while (*pTable != 0xFFFF);
 
1998
         }
 
1999
      }
 
2000
found_perfect_match:
 
2001
 
 
2002
      const uint diff = best_x & 1;
 
2003
      const uint inten = (best_x >> 1) & 7;
 
2004
 
 
2005
      block.m_bytes[3] = static_cast<uint8>(((inten | (inten << 3)) << 2) | (diff << 1));
 
2006
                        
 
2007
      const uint etc1_selector = g_selector_index_to_etc1[(best_x >> 4) & 3];
 
2008
      *reinterpret_cast<uint16*>(&block.m_bytes[4]) = (etc1_selector & 2) ? 0xFFFF : 0;
 
2009
      *reinterpret_cast<uint16*>(&block.m_bytes[6]) = (etc1_selector & 1) ? 0xFFFF : 0;
 
2010
 
 
2011
      const uint best_packed_c0 = (best_x >> 8) & 255;
 
2012
      if (diff)
 
2013
      {
 
2014
         block.m_bytes[best_i] = static_cast<uint8>(best_packed_c0 << 3);
 
2015
         block.m_bytes[s_next_comp[best_i]] = static_cast<uint8>(best_packed_c1 << 3);
 
2016
         block.m_bytes[s_next_comp[best_i+1]] = static_cast<uint8>(best_packed_c2 << 3);
 
2017
      }
 
2018
      else
 
2019
      {
 
2020
         block.m_bytes[best_i] = static_cast<uint8>(best_packed_c0 | (best_packed_c0 << 4));
 
2021
         block.m_bytes[s_next_comp[best_i]] = static_cast<uint8>(best_packed_c1 | (best_packed_c1 << 4));
 
2022
         block.m_bytes[s_next_comp[best_i+1]] = static_cast<uint8>(best_packed_c2 | (best_packed_c2 << 4));
 
2023
      }
 
2024
 
 
2025
      return best_error;
 
2026
   }
 
2027
      
 
2028
   static uint pack_etc1_block_solid_color_constrained(
 
2029
      etc1_optimizer::results& results, 
 
2030
      uint num_colors, const uint8* pColor, 
 
2031
      etc1_pack_params& pack_params, 
 
2032
      bool use_diff,
 
2033
      const color_quad_u8* pBase_color5_unscaled)
 
2034
   {
 
2035
      RG_ETC1_ASSERT(g_etc1_inverse_lookup[0][255]);
 
2036
 
 
2037
      //pack_params;
 
2038
      static uint s_next_comp[4] = { 1, 2, 0, 1 };
 
2039
 
 
2040
      uint best_error = cUINT32_MAX, best_i = 0;
 
2041
      int best_x = 0, best_packed_c1 = 0, best_packed_c2 = 0;
 
2042
 
 
2043
      // For each possible 8-bit value, there is a precomputed list of diff/inten/selector configurations that allow that 8-bit value to be encoded with no error.
 
2044
      for (uint i = 0; i < 3; i++)
 
2045
      {
 
2046
         const uint c1 = pColor[s_next_comp[i]], c2 = pColor[s_next_comp[i + 1]];
 
2047
 
 
2048
         const int delta_range = 1;
 
2049
         for (int delta = -delta_range; delta <= delta_range; delta++)
 
2050
         {
 
2051
            const int c_plus_delta = rg_etc1::clamp<int>(pColor[i] + delta, 0, 255);
 
2052
 
 
2053
            const uint16* pTable;
 
2054
            if (!c_plus_delta)
 
2055
               pTable = g_color8_to_etc_block_config_0_255[0];
 
2056
            else if (c_plus_delta == 255)
 
2057
               pTable = g_color8_to_etc_block_config_0_255[1];
 
2058
            else
 
2059
               pTable = g_color8_to_etc_block_config_1_to_254[c_plus_delta - 1];
 
2060
 
 
2061
            do
 
2062
            {
 
2063
               const uint x = *pTable++;
 
2064
               const uint diff = x & 1;
 
2065
               if (static_cast<uint>(use_diff) != diff)
 
2066
               {
 
2067
                  if (*pTable == 0xFFFF)
 
2068
                     break;
 
2069
                  continue;
 
2070
               }
 
2071
 
 
2072
               if ((diff) && (pBase_color5_unscaled))
 
2073
               {
 
2074
                  const int p0 = (x >> 8) & 255;
 
2075
                  int delta = p0 - static_cast<int>(pBase_color5_unscaled->c[i]);
 
2076
                  if ((delta < cETC1ColorDeltaMin) || (delta > cETC1ColorDeltaMax))
 
2077
                  {
 
2078
                     if (*pTable == 0xFFFF)
 
2079
                        break;
 
2080
                     continue;
 
2081
                  }
 
2082
               }
 
2083
 
 
2084
#ifdef RG_ETC1_BUILD_DEBUG
 
2085
               {
 
2086
                  const uint inten = (x >> 1) & 7;
 
2087
                  const uint selector = (x >> 4) & 3;
 
2088
                  const uint p0 = (x >> 8) & 255;
 
2089
                  RG_ETC1_ASSERT(etc1_decode_value(diff, inten, selector, p0) == (uint)c_plus_delta);
 
2090
               }
 
2091
#endif
 
2092
 
 
2093
               const uint16* pInverse_table = g_etc1_inverse_lookup[x & 0xFF];
 
2094
               uint16 p1 = pInverse_table[c1];
 
2095
               uint16 p2 = pInverse_table[c2];
 
2096
 
 
2097
               if ((diff) && (pBase_color5_unscaled))
 
2098
               {
 
2099
                  int delta1 = (p1 & 0xFF) - static_cast<int>(pBase_color5_unscaled->c[s_next_comp[i]]);
 
2100
                  int delta2 = (p2 & 0xFF) - static_cast<int>(pBase_color5_unscaled->c[s_next_comp[i + 1]]);
 
2101
                  if ((delta1 < cETC1ColorDeltaMin) || (delta1 > cETC1ColorDeltaMax) || (delta2 < cETC1ColorDeltaMin) || (delta2 > cETC1ColorDeltaMax))
 
2102
                  {
 
2103
                     if (*pTable == 0xFFFF)
 
2104
                        break;
 
2105
                     continue;
 
2106
                  }
 
2107
               }
 
2108
 
 
2109
               const uint trial_error = rg_etc1::square(c_plus_delta - pColor[i]) + rg_etc1::square(p1 >> 8) + rg_etc1::square(p2 >> 8);
 
2110
               if (trial_error < best_error)
 
2111
               {
 
2112
                  best_error = trial_error;
 
2113
                  best_x = x;
 
2114
                  best_packed_c1 = p1 & 0xFF;
 
2115
                  best_packed_c2 = p2 & 0xFF;
 
2116
                  best_i = i;
 
2117
                  if (!best_error)
 
2118
                     goto found_perfect_match;
 
2119
               }
 
2120
            } while (*pTable != 0xFFFF);
 
2121
         }
 
2122
      }
 
2123
found_perfect_match:
 
2124
 
 
2125
      if (best_error == cUINT32_MAX)
 
2126
         return best_error;
 
2127
 
 
2128
      best_error *= num_colors;
 
2129
 
 
2130
      results.m_n = num_colors;
 
2131
      results.m_block_color4 = !(best_x & 1);
 
2132
      results.m_block_inten_table = (best_x >> 1) & 7;
 
2133
      memset(results.m_pSelectors, (best_x >> 4) & 3, num_colors);
 
2134
 
 
2135
      const uint best_packed_c0 = (best_x >> 8) & 255;
 
2136
      results.m_block_color_unscaled[best_i] = static_cast<uint8>(best_packed_c0);
 
2137
      results.m_block_color_unscaled[s_next_comp[best_i]] = static_cast<uint8>(best_packed_c1);
 
2138
      results.m_block_color_unscaled[s_next_comp[best_i + 1]] = static_cast<uint8>(best_packed_c2);
 
2139
      results.m_error = best_error;
 
2140
      
 
2141
      return best_error;
 
2142
   }
 
2143
 
 
2144
   // Function originally from RYG's public domain real-time DXT1 compressor, modified for 555.
 
2145
   static void dither_block_555(color_quad_u8* dest, const color_quad_u8* block)
 
2146
   {
 
2147
      int err[8],*ep1 = err,*ep2 = err+4;
 
2148
      uint8 *quant = g_quant5_tab+8;
 
2149
 
 
2150
      memset(dest, 0xFF, sizeof(color_quad_u8)*16);
 
2151
 
 
2152
      // process channels seperately
 
2153
      for(int ch=0;ch<3;ch++)
 
2154
      {
 
2155
         uint8* bp = (uint8*)block;
 
2156
         uint8* dp = (uint8*)dest;
 
2157
 
 
2158
         bp += ch; dp += ch;
 
2159
 
 
2160
         memset(err,0, sizeof(err));
 
2161
         for(int y = 0; y < 4; y++)
 
2162
         {
 
2163
            // pixel 0
 
2164
            dp[ 0] = quant[bp[ 0] + ((3*ep2[1] + 5*ep2[0]) >> 4)];
 
2165
            ep1[0] = bp[ 0] - dp[ 0];
 
2166
 
 
2167
            // pixel 1
 
2168
            dp[ 4] = quant[bp[ 4] + ((7*ep1[0] + 3*ep2[2] + 5*ep2[1] + ep2[0]) >> 4)];
 
2169
            ep1[1] = bp[ 4] - dp[ 4];
 
2170
 
 
2171
            // pixel 2
 
2172
            dp[ 8] = quant[bp[ 8] + ((7*ep1[1] + 3*ep2[3] + 5*ep2[2] + ep2[1]) >> 4)];
 
2173
            ep1[2] = bp[ 8] - dp[ 8];
 
2174
 
 
2175
            // pixel 3
 
2176
            dp[12] = quant[bp[12] + ((7*ep1[2] + 5*ep2[3] + ep2[2]) >> 4)];
 
2177
            ep1[3] = bp[12] - dp[12];
 
2178
 
 
2179
            // advance to next line
 
2180
            int* tmp = ep1; ep1 = ep2; ep2 = tmp;
 
2181
            bp += 16;
 
2182
            dp += 16;
 
2183
         }
 
2184
      }
 
2185
   }
 
2186
 
 
2187
   unsigned int pack_etc1_block(void* pETC1_block, const unsigned int* pSrc_pixels_rgba, etc1_pack_params& pack_params)
 
2188
   {
 
2189
      const color_quad_u8* pSrc_pixels = reinterpret_cast<const color_quad_u8*>(pSrc_pixels_rgba);
 
2190
      etc1_block& dst_block = *static_cast<etc1_block*>(pETC1_block);
 
2191
 
 
2192
#ifdef RG_ETC1_BUILD_DEBUG
 
2193
      // Ensure all alpha values are 0xFF.
 
2194
      for (uint i = 0; i < 16; i++)
 
2195
      {
 
2196
         RG_ETC1_ASSERT(pSrc_pixels[i].a == 255);
 
2197
      }
 
2198
#endif
 
2199
 
 
2200
      color_quad_u8 src_pixel0(pSrc_pixels[0]);
 
2201
 
 
2202
      // Check for solid block.
 
2203
      const uint32 first_pixel_u32 = pSrc_pixels->m_u32;
 
2204
      int r;
 
2205
      for (r = 15; r >= 1; --r)
 
2206
         if (pSrc_pixels[r].m_u32 != first_pixel_u32)
 
2207
            break;
 
2208
      if (!r)
 
2209
         return static_cast<unsigned int>(16 * pack_etc1_block_solid_color(dst_block, &pSrc_pixels[0].r, pack_params));
 
2210
      
 
2211
      color_quad_u8 dithered_pixels[16];
 
2212
      if (pack_params.m_dithering)
 
2213
      {
 
2214
         dither_block_555(dithered_pixels, pSrc_pixels);
 
2215
         pSrc_pixels = dithered_pixels;
 
2216
      }
 
2217
 
 
2218
      etc1_optimizer optimizer;
 
2219
 
 
2220
      uint64 best_error = cUINT64_MAX;
 
2221
      uint best_flip = false, best_use_color4 = false;
 
2222
      
 
2223
      uint8 best_selectors[2][8];
 
2224
      etc1_optimizer::results best_results[2];
 
2225
      for (uint i = 0; i < 2; i++)
 
2226
      {
 
2227
         best_results[i].m_n = 8;
 
2228
         best_results[i].m_pSelectors = best_selectors[i];
 
2229
      }
 
2230
      
 
2231
      uint8 selectors[3][8];
 
2232
      etc1_optimizer::results results[3];
 
2233
      
 
2234
      for (uint i = 0; i < 3; i++)
 
2235
      {
 
2236
         results[i].m_n = 8;
 
2237
         results[i].m_pSelectors = selectors[i];
 
2238
      }
 
2239
            
 
2240
      color_quad_u8 subblock_pixels[8];
 
2241
 
 
2242
      etc1_optimizer::params params(pack_params);
 
2243
      params.m_num_src_pixels = 8;
 
2244
      params.m_pSrc_pixels = subblock_pixels;
 
2245
 
 
2246
      for (uint flip = 0; flip < 2; flip++)
 
2247
      {
 
2248
         for (uint use_color4 = 0; use_color4 < 2; use_color4++)
 
2249
         {
 
2250
            uint64 trial_error = 0;
 
2251
 
 
2252
            uint subblock;
 
2253
            for (subblock = 0; subblock < 2; subblock++)
 
2254
            {
 
2255
               if (flip)
 
2256
                  memcpy(subblock_pixels, pSrc_pixels + subblock * 8, sizeof(color_quad_u8) * 8);
 
2257
               else
 
2258
               {
 
2259
                  const color_quad_u8* pSrc_col = pSrc_pixels + subblock * 2;
 
2260
                  subblock_pixels[0] = pSrc_col[0]; subblock_pixels[1] = pSrc_col[4]; subblock_pixels[2] = pSrc_col[8]; subblock_pixels[3] = pSrc_col[12];
 
2261
                  subblock_pixels[4] = pSrc_col[1]; subblock_pixels[5] = pSrc_col[5]; subblock_pixels[6] = pSrc_col[9]; subblock_pixels[7] = pSrc_col[13];
 
2262
               }
 
2263
 
 
2264
               results[2].m_error = cUINT64_MAX;
 
2265
               if ((params.m_quality >= cMediumQuality) && ((subblock) || (use_color4)))
 
2266
               {
 
2267
                  const uint32 subblock_pixel0_u32 = subblock_pixels[0].m_u32;
 
2268
                  for (r = 7; r >= 1; --r)
 
2269
                     if (subblock_pixels[r].m_u32 != subblock_pixel0_u32)
 
2270
                        break;
 
2271
                  if (!r)
 
2272
                  {
 
2273
                     pack_etc1_block_solid_color_constrained(results[2], 8, &subblock_pixels[0].r, pack_params, !use_color4, (subblock && !use_color4) ? &results[0].m_block_color_unscaled : NULL);
 
2274
                  }
 
2275
               }
 
2276
 
 
2277
               params.m_use_color4 = (use_color4 != 0);
 
2278
               params.m_constrain_against_base_color5 = false;
 
2279
 
 
2280
               if ((!use_color4) && (subblock))
 
2281
               {
 
2282
                  params.m_constrain_against_base_color5 = true;
 
2283
                  params.m_base_color5 = results[0].m_block_color_unscaled;
 
2284
               }
 
2285
                              
 
2286
               if (params.m_quality == cHighQuality)
 
2287
               {
 
2288
                  static const int s_scan_delta_0_to_4[] = { -4, -3, -2, -1, 0, 1, 2, 3, 4 };
 
2289
                  params.m_scan_delta_size = RG_ETC1_ARRAY_SIZE(s_scan_delta_0_to_4);
 
2290
                  params.m_pScan_deltas = s_scan_delta_0_to_4;
 
2291
               }
 
2292
               else if (params.m_quality == cMediumQuality)
 
2293
               {
 
2294
                  static const int s_scan_delta_0_to_1[] = { -1, 0, 1 };
 
2295
                  params.m_scan_delta_size = RG_ETC1_ARRAY_SIZE(s_scan_delta_0_to_1);
 
2296
                  params.m_pScan_deltas = s_scan_delta_0_to_1;
 
2297
               }
 
2298
               else
 
2299
               {
 
2300
                  static const int s_scan_delta_0[] = { 0 };
 
2301
                  params.m_scan_delta_size = RG_ETC1_ARRAY_SIZE(s_scan_delta_0);
 
2302
                  params.m_pScan_deltas = s_scan_delta_0;
 
2303
               }
 
2304
               
 
2305
               optimizer.init(params, results[subblock]);
 
2306
               if (!optimizer.compute())
 
2307
                  break;
 
2308
                              
 
2309
               if (params.m_quality >= cMediumQuality)
 
2310
               {
 
2311
                  // TODO: Fix fairly arbitrary/unrefined thresholds that control how far away to scan for potentially better solutions.
 
2312
                  const uint refinement_error_thresh0 = 3000;
 
2313
                  const uint refinement_error_thresh1 = 6000;
 
2314
                  if (results[subblock].m_error > refinement_error_thresh0)
 
2315
                  {
 
2316
                     if (params.m_quality == cMediumQuality)
 
2317
                     {
 
2318
                        static const int s_scan_delta_2_to_3[] = { -3, -2, 2, 3 };
 
2319
                        params.m_scan_delta_size = RG_ETC1_ARRAY_SIZE(s_scan_delta_2_to_3);
 
2320
                        params.m_pScan_deltas = s_scan_delta_2_to_3;
 
2321
                     }
 
2322
                     else
 
2323
                     {
 
2324
                        static const int s_scan_delta_5_to_5[] = { -5, 5 };
 
2325
                        static const int s_scan_delta_5_to_8[] = { -8, -7, -6, -5, 5, 6, 7, 8 };
 
2326
                        if (results[subblock].m_error > refinement_error_thresh1)
 
2327
                        {
 
2328
                           params.m_scan_delta_size = RG_ETC1_ARRAY_SIZE(s_scan_delta_5_to_8);
 
2329
                           params.m_pScan_deltas = s_scan_delta_5_to_8;
 
2330
                        }
 
2331
                        else
 
2332
                        {
 
2333
                           params.m_scan_delta_size = RG_ETC1_ARRAY_SIZE(s_scan_delta_5_to_5);
 
2334
                           params.m_pScan_deltas = s_scan_delta_5_to_5;
 
2335
                        }
 
2336
                     }
 
2337
 
 
2338
                     if (!optimizer.compute())
 
2339
                        break;
 
2340
                  }
 
2341
 
 
2342
                  if (results[2].m_error < results[subblock].m_error)
 
2343
                     results[subblock] = results[2];
 
2344
               }
 
2345
                            
 
2346
               trial_error += results[subblock].m_error;
 
2347
               if (trial_error >= best_error)
 
2348
                  break;
 
2349
            }
 
2350
 
 
2351
            if (subblock < 2)
 
2352
               continue;
 
2353
 
 
2354
            best_error = trial_error;
 
2355
            best_results[0] = results[0];
 
2356
            best_results[1] = results[1];
 
2357
            best_flip = flip;
 
2358
            best_use_color4 = use_color4;
 
2359
            
 
2360
         } // use_color4
 
2361
 
 
2362
      } // flip
 
2363
 
 
2364
      int dr = best_results[1].m_block_color_unscaled.r - best_results[0].m_block_color_unscaled.r;
 
2365
      int dg = best_results[1].m_block_color_unscaled.g - best_results[0].m_block_color_unscaled.g;
 
2366
      int db = best_results[1].m_block_color_unscaled.b - best_results[0].m_block_color_unscaled.b;
 
2367
      RG_ETC1_ASSERT(best_use_color4 || (rg_etc1::minimum(dr, dg, db) >= cETC1ColorDeltaMin) && (rg_etc1::maximum(dr, dg, db) <= cETC1ColorDeltaMax));
 
2368
           
 
2369
      if (best_use_color4)
 
2370
      {
 
2371
         dst_block.m_bytes[0] = static_cast<uint8>(best_results[1].m_block_color_unscaled.r | (best_results[0].m_block_color_unscaled.r << 4));
 
2372
         dst_block.m_bytes[1] = static_cast<uint8>(best_results[1].m_block_color_unscaled.g | (best_results[0].m_block_color_unscaled.g << 4));
 
2373
         dst_block.m_bytes[2] = static_cast<uint8>(best_results[1].m_block_color_unscaled.b | (best_results[0].m_block_color_unscaled.b << 4));
 
2374
      }
 
2375
      else
 
2376
      {
 
2377
         if (dr < 0) dr += 8; dst_block.m_bytes[0] = static_cast<uint8>((best_results[0].m_block_color_unscaled.r << 3) | dr);
 
2378
         if (dg < 0) dg += 8; dst_block.m_bytes[1] = static_cast<uint8>((best_results[0].m_block_color_unscaled.g << 3) | dg);
 
2379
         if (db < 0) db += 8; dst_block.m_bytes[2] = static_cast<uint8>((best_results[0].m_block_color_unscaled.b << 3) | db);
 
2380
      }
 
2381
      
 
2382
      dst_block.m_bytes[3] = static_cast<uint8>( (best_results[1].m_block_inten_table << 2) | (best_results[0].m_block_inten_table << 5) | ((~best_use_color4 & 1) << 1) | best_flip );
 
2383
      
 
2384
      uint selector0 = 0, selector1 = 0;
 
2385
      if (best_flip)
 
2386
      {
 
2387
         // flipped:
 
2388
         // { 0, 0 }, { 1, 0 }, { 2, 0 }, { 3, 0 },               
 
2389
         // { 0, 1 }, { 1, 1 }, { 2, 1 }, { 3, 1 } 
 
2390
         //
 
2391
         // { 0, 2 }, { 1, 2 }, { 2, 2 }, { 3, 2 },
 
2392
         // { 0, 3 }, { 1, 3 }, { 2, 3 }, { 3, 3 }
 
2393
         const uint8* pSelectors0 = best_results[0].m_pSelectors;
 
2394
         const uint8* pSelectors1 = best_results[1].m_pSelectors;
 
2395
         for (int x = 3; x >= 0; --x)
 
2396
         {
 
2397
            uint b;
 
2398
            b = g_selector_index_to_etc1[pSelectors1[4 + x]];
 
2399
            selector0 = (selector0 << 1) | (b & 1); selector1 = (selector1 << 1) | (b >> 1);
 
2400
 
 
2401
            b = g_selector_index_to_etc1[pSelectors1[x]];
 
2402
            selector0 = (selector0 << 1) | (b & 1); selector1 = (selector1 << 1) | (b >> 1);
 
2403
 
 
2404
            b = g_selector_index_to_etc1[pSelectors0[4 + x]];
 
2405
            selector0 = (selector0 << 1) | (b & 1); selector1 = (selector1 << 1) | (b >> 1);
 
2406
 
 
2407
            b = g_selector_index_to_etc1[pSelectors0[x]];
 
2408
            selector0 = (selector0 << 1) | (b & 1); selector1 = (selector1 << 1) | (b >> 1);
 
2409
         }
 
2410
      }
 
2411
      else
 
2412
      {
 
2413
         // non-flipped:
 
2414
         // { 0, 0 }, { 0, 1 }, { 0, 2 }, { 0, 3 },
 
2415
         // { 1, 0 }, { 1, 1 }, { 1, 2 }, { 1, 3 }
 
2416
         //
 
2417
         // { 2, 0 }, { 2, 1 }, { 2, 2 }, { 2, 3 },
 
2418
         // { 3, 0 }, { 3, 1 }, { 3, 2 }, { 3, 3 }
 
2419
         for (int subblock = 1; subblock >= 0; --subblock)
 
2420
         {
 
2421
            const uint8* pSelectors = best_results[subblock].m_pSelectors + 4;
 
2422
            for (uint i = 0; i < 2; i++)
 
2423
            {
 
2424
               uint b;
 
2425
               b = g_selector_index_to_etc1[pSelectors[3]];
 
2426
               selector0 = (selector0 << 1) | (b & 1); selector1 = (selector1 << 1) | (b >> 1);
 
2427
 
 
2428
               b = g_selector_index_to_etc1[pSelectors[2]];
 
2429
               selector0 = (selector0 << 1) | (b & 1); selector1 = (selector1 << 1) | (b >> 1);
 
2430
 
 
2431
               b = g_selector_index_to_etc1[pSelectors[1]];
 
2432
               selector0 = (selector0 << 1) | (b & 1); selector1 = (selector1 << 1) | (b >> 1);
 
2433
 
 
2434
               b = g_selector_index_to_etc1[pSelectors[0]];
 
2435
               selector0 = (selector0 << 1) | (b & 1);selector1 = (selector1 << 1) | (b >> 1);
 
2436
 
 
2437
               pSelectors -= 4;
 
2438
            }
 
2439
         }
 
2440
      }
 
2441
                  
 
2442
      dst_block.m_bytes[4] = static_cast<uint8>(selector1 >> 8); dst_block.m_bytes[5] = static_cast<uint8>(selector1 & 0xFF);
 
2443
      dst_block.m_bytes[6] = static_cast<uint8>(selector0 >> 8); dst_block.m_bytes[7] = static_cast<uint8>(selector0 & 0xFF);
 
2444
 
 
2445
      return static_cast<unsigned int>(best_error);
 
2446
   }
 
2447
 
 
2448
} // namespace rg_etc1