1727 lines
49 KiB
C
1727 lines
49 KiB
C
/**************************************************************************
|
|
*
|
|
* Copyright (C) 1999-2008 Brian Paul All Rights Reserved.
|
|
* Copyright (c) 2008 VMware, Inc.
|
|
*
|
|
* Permission is hereby granted, free of charge, to any person obtaining a
|
|
* copy of this software and associated documentation files (the "Software"),
|
|
* to deal in the Software without restriction, including without limitation
|
|
* the rights to use, copy, modify, merge, publish, distribute, sublicense,
|
|
* and/or sell copies of the Software, and to permit persons to whom the
|
|
* Software is furnished to do so, subject to the following conditions:
|
|
*
|
|
* The above copyright notice and this permission notice shall be included
|
|
* in all copies or substantial portions of the Software.
|
|
*
|
|
* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
|
|
* OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
|
|
* FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
|
|
* THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR
|
|
* OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,
|
|
* ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR
|
|
* OTHER DEALINGS IN THE SOFTWARE.
|
|
*
|
|
**************************************************************************/
|
|
|
|
#include "util/format/u_format.h"
|
|
#include "util/format/u_format_fxt1.h"
|
|
#include "util/format/u_format_pack.h"
|
|
#include "util/format_srgb.h"
|
|
#include "util/u_math.h"
|
|
|
|
#define RCOMP 0
|
|
#define GCOMP 1
|
|
#define BCOMP 2
|
|
#define ACOMP 3
|
|
|
|
#define FXT1_BLOCK_SIZE 16
|
|
|
|
static void
|
|
fxt1_encode (uint32_t width, uint32_t height, int32_t comps,
|
|
const void *source, int32_t srcRowStride,
|
|
void *dest, int32_t destRowStride);
|
|
|
|
static void
|
|
fxt1_decode_1 (const void *texture, int32_t stride,
|
|
int32_t i, int32_t j, uint8_t *rgba);
|
|
|
|
/***************************************************************************\
|
|
* FXT1 encoder
|
|
*
|
|
* The encoder was built by reversing the decoder,
|
|
* and is vaguely based on Texus2 by 3dfx. Note that this code
|
|
* is merely a proof of concept, since it is highly UNoptimized;
|
|
* moreover, it is sub-optimal due to initial conditions passed
|
|
* to Lloyd's algorithm (the interpolation modes are even worse).
|
|
\***************************************************************************/
|
|
|
|
|
|
#define MAX_COMP 4 /* ever needed maximum number of components in texel */
|
|
#define MAX_VECT 4 /* ever needed maximum number of base vectors to find */
|
|
#define N_TEXELS 32 /* number of texels in a block (always 32) */
|
|
#define LL_N_REP 50 /* number of iterations in lloyd's vq */
|
|
#define LL_RMS_D 10 /* fault tolerance (maximum delta) */
|
|
#define LL_RMS_E 255 /* fault tolerance (maximum error) */
|
|
#define ALPHA_TS 2 /* alpha threshold: (255 - ALPHA_TS) deemed opaque */
|
|
static const uint32_t zero = 0;
|
|
#define ISTBLACK(v) (memcmp(&(v), &zero, sizeof(zero)) == 0)
|
|
|
|
/*
|
|
* Define a 64-bit unsigned integer type and macros
|
|
*/
|
|
#if 1
|
|
|
|
#define FX64_NATIVE 1
|
|
|
|
typedef uint64_t Fx64;
|
|
|
|
#define FX64_MOV32(a, b) a = b
|
|
#define FX64_OR32(a, b) a |= b
|
|
#define FX64_SHL(a, c) a <<= c
|
|
|
|
#else
|
|
|
|
#define FX64_NATIVE 0
|
|
|
|
typedef struct {
|
|
uint32_t lo, hi;
|
|
} Fx64;
|
|
|
|
#define FX64_MOV32(a, b) a.lo = b
|
|
#define FX64_OR32(a, b) a.lo |= b
|
|
|
|
#define FX64_SHL(a, c) \
|
|
do { \
|
|
if ((c) >= 32) { \
|
|
a.hi = a.lo << ((c) - 32); \
|
|
a.lo = 0; \
|
|
} else { \
|
|
a.hi = (a.hi << (c)) | (a.lo >> (32 - (c))); \
|
|
a.lo <<= (c); \
|
|
} \
|
|
} while (0)
|
|
|
|
#endif
|
|
|
|
|
|
#define F(i) (float)1 /* can be used to obtain an oblong metric: 0.30 / 0.59 / 0.11 */
|
|
#define SAFECDOT 1 /* for paranoids */
|
|
|
|
#define MAKEIVEC(NV, NC, IV, B, V0, V1) \
|
|
do { \
|
|
/* compute interpolation vector */ \
|
|
float d2 = 0.0F; \
|
|
float rd2; \
|
|
\
|
|
for (i = 0; i < NC; i++) { \
|
|
IV[i] = (V1[i] - V0[i]) * F(i); \
|
|
d2 += IV[i] * IV[i]; \
|
|
} \
|
|
rd2 = (float)NV / d2; \
|
|
B = 0; \
|
|
for (i = 0; i < NC; i++) { \
|
|
IV[i] *= F(i); \
|
|
B -= IV[i] * V0[i]; \
|
|
IV[i] *= rd2; \
|
|
} \
|
|
B = B * rd2 + 0.5f; \
|
|
} while (0)
|
|
|
|
#define CALCCDOT(TEXEL, NV, NC, IV, B, V)\
|
|
do { \
|
|
float dot = 0.0F; \
|
|
for (i = 0; i < NC; i++) { \
|
|
dot += V[i] * IV[i]; \
|
|
} \
|
|
TEXEL = (int32_t)(dot + B); \
|
|
if (SAFECDOT) { \
|
|
if (TEXEL < 0) { \
|
|
TEXEL = 0; \
|
|
} else if (TEXEL > NV) { \
|
|
TEXEL = NV; \
|
|
} \
|
|
} \
|
|
} while (0)
|
|
|
|
|
|
static int32_t
|
|
fxt1_bestcol (float vec[][MAX_COMP], int32_t nv,
|
|
uint8_t input[MAX_COMP], int32_t nc)
|
|
{
|
|
int32_t i, j, best = -1;
|
|
float err = 1e9; /* big enough */
|
|
|
|
for (j = 0; j < nv; j++) {
|
|
float e = 0.0F;
|
|
for (i = 0; i < nc; i++) {
|
|
e += (vec[j][i] - input[i]) * (vec[j][i] - input[i]);
|
|
}
|
|
if (e < err) {
|
|
err = e;
|
|
best = j;
|
|
}
|
|
}
|
|
|
|
return best;
|
|
}
|
|
|
|
|
|
static int32_t
|
|
fxt1_worst (float vec[MAX_COMP],
|
|
uint8_t input[N_TEXELS][MAX_COMP], int32_t nc, int32_t n)
|
|
{
|
|
int32_t i, k, worst = -1;
|
|
float err = -1.0F; /* small enough */
|
|
|
|
for (k = 0; k < n; k++) {
|
|
float e = 0.0F;
|
|
for (i = 0; i < nc; i++) {
|
|
e += (vec[i] - input[k][i]) * (vec[i] - input[k][i]);
|
|
}
|
|
if (e > err) {
|
|
err = e;
|
|
worst = k;
|
|
}
|
|
}
|
|
|
|
return worst;
|
|
}
|
|
|
|
|
|
static int32_t
|
|
fxt1_variance (uint8_t input[N_TEXELS / 2][MAX_COMP], int32_t nc)
|
|
{
|
|
const int n = N_TEXELS / 2;
|
|
int32_t i, k, best = 0;
|
|
int32_t sx, sx2;
|
|
double var, maxvar = -1; /* small enough */
|
|
double teenth = 1.0 / n;
|
|
|
|
for (i = 0; i < nc; i++) {
|
|
sx = sx2 = 0;
|
|
for (k = 0; k < n; k++) {
|
|
int32_t t = input[k][i];
|
|
sx += t;
|
|
sx2 += t * t;
|
|
}
|
|
var = sx2 * teenth - sx * sx * teenth * teenth;
|
|
if (maxvar < var) {
|
|
maxvar = var;
|
|
best = i;
|
|
}
|
|
}
|
|
|
|
return best;
|
|
}
|
|
|
|
|
|
static int32_t
|
|
fxt1_choose (float vec[][MAX_COMP], int32_t nv,
|
|
uint8_t input[N_TEXELS][MAX_COMP], int32_t nc, int32_t n)
|
|
{
|
|
#if 0
|
|
/* Choose colors from a grid.
|
|
*/
|
|
int32_t i, j;
|
|
|
|
for (j = 0; j < nv; j++) {
|
|
int32_t m = j * (n - 1) / (nv - 1);
|
|
for (i = 0; i < nc; i++) {
|
|
vec[j][i] = input[m][i];
|
|
}
|
|
}
|
|
#else
|
|
/* Our solution here is to find the darkest and brightest colors in
|
|
* the 8x4 tile and use those as the two representative colors.
|
|
* There are probably better algorithms to use (histogram-based).
|
|
*/
|
|
int32_t i, j, k;
|
|
int32_t minSum = 2000; /* big enough */
|
|
int32_t maxSum = -1; /* small enough */
|
|
int32_t minCol = 0; /* phoudoin: silent compiler! */
|
|
int32_t maxCol = 0; /* phoudoin: silent compiler! */
|
|
|
|
struct {
|
|
int32_t flag;
|
|
int32_t key;
|
|
int32_t freq;
|
|
int32_t idx;
|
|
} hist[N_TEXELS];
|
|
int32_t lenh = 0;
|
|
|
|
memset(hist, 0, sizeof(hist));
|
|
|
|
for (k = 0; k < n; k++) {
|
|
int32_t l;
|
|
int32_t key = 0;
|
|
int32_t sum = 0;
|
|
for (i = 0; i < nc; i++) {
|
|
key <<= 8;
|
|
key |= input[k][i];
|
|
sum += input[k][i];
|
|
}
|
|
for (l = 0; l < n; l++) {
|
|
if (!hist[l].flag) {
|
|
/* alloc new slot */
|
|
hist[l].flag = !0;
|
|
hist[l].key = key;
|
|
hist[l].freq = 1;
|
|
hist[l].idx = k;
|
|
lenh = l + 1;
|
|
break;
|
|
} else if (hist[l].key == key) {
|
|
hist[l].freq++;
|
|
break;
|
|
}
|
|
}
|
|
if (minSum > sum) {
|
|
minSum = sum;
|
|
minCol = k;
|
|
}
|
|
if (maxSum < sum) {
|
|
maxSum = sum;
|
|
maxCol = k;
|
|
}
|
|
}
|
|
|
|
if (lenh <= nv) {
|
|
for (j = 0; j < lenh; j++) {
|
|
for (i = 0; i < nc; i++) {
|
|
vec[j][i] = (float)input[hist[j].idx][i];
|
|
}
|
|
}
|
|
for (; j < nv; j++) {
|
|
for (i = 0; i < nc; i++) {
|
|
vec[j][i] = vec[0][i];
|
|
}
|
|
}
|
|
return 0;
|
|
}
|
|
|
|
for (j = 0; j < nv; j++) {
|
|
for (i = 0; i < nc; i++) {
|
|
vec[j][i] = ((nv - 1 - j) * input[minCol][i] + j * input[maxCol][i] + (nv - 1) / 2) / (float)(nv - 1);
|
|
}
|
|
}
|
|
#endif
|
|
|
|
return !0;
|
|
}
|
|
|
|
|
|
static int32_t
|
|
fxt1_lloyd (float vec[][MAX_COMP], int32_t nv,
|
|
uint8_t input[N_TEXELS][MAX_COMP], int32_t nc, int32_t n)
|
|
{
|
|
/* Use the generalized lloyd's algorithm for VQ:
|
|
* find 4 color vectors.
|
|
*
|
|
* for each sample color
|
|
* sort to nearest vector.
|
|
*
|
|
* replace each vector with the centroid of its matching colors.
|
|
*
|
|
* repeat until RMS doesn't improve.
|
|
*
|
|
* if a color vector has no samples, or becomes the same as another
|
|
* vector, replace it with the color which is farthest from a sample.
|
|
*
|
|
* vec[][MAX_COMP] initial vectors and resulting colors
|
|
* nv number of resulting colors required
|
|
* input[N_TEXELS][MAX_COMP] input texels
|
|
* nc number of components in input / vec
|
|
* n number of input samples
|
|
*/
|
|
|
|
int32_t sum[MAX_VECT][MAX_COMP]; /* used to accumulate closest texels */
|
|
int32_t cnt[MAX_VECT]; /* how many times a certain vector was chosen */
|
|
float error, lasterror = 1e9;
|
|
|
|
int32_t i, j, k, rep;
|
|
|
|
/* the quantizer */
|
|
for (rep = 0; rep < LL_N_REP; rep++) {
|
|
/* reset sums & counters */
|
|
for (j = 0; j < nv; j++) {
|
|
for (i = 0; i < nc; i++) {
|
|
sum[j][i] = 0;
|
|
}
|
|
cnt[j] = 0;
|
|
}
|
|
error = 0;
|
|
|
|
/* scan whole block */
|
|
for (k = 0; k < n; k++) {
|
|
#if 1
|
|
int32_t best = -1;
|
|
float err = 1e9; /* big enough */
|
|
/* determine best vector */
|
|
for (j = 0; j < nv; j++) {
|
|
float e = (vec[j][0] - input[k][0]) * (vec[j][0] - input[k][0]) +
|
|
(vec[j][1] - input[k][1]) * (vec[j][1] - input[k][1]) +
|
|
(vec[j][2] - input[k][2]) * (vec[j][2] - input[k][2]);
|
|
if (nc == 4) {
|
|
e += (vec[j][3] - input[k][3]) * (vec[j][3] - input[k][3]);
|
|
}
|
|
if (e < err) {
|
|
err = e;
|
|
best = j;
|
|
}
|
|
}
|
|
#else
|
|
int32_t best = fxt1_bestcol(vec, nv, input[k], nc, &err);
|
|
#endif
|
|
assert(best >= 0);
|
|
/* add in closest color */
|
|
for (i = 0; i < nc; i++) {
|
|
sum[best][i] += input[k][i];
|
|
}
|
|
/* mark this vector as used */
|
|
cnt[best]++;
|
|
/* accumulate error */
|
|
error += err;
|
|
}
|
|
|
|
/* check RMS */
|
|
if ((error < LL_RMS_E) ||
|
|
((error < lasterror) && ((lasterror - error) < LL_RMS_D))) {
|
|
return !0; /* good match */
|
|
}
|
|
lasterror = error;
|
|
|
|
/* move each vector to the barycenter of its closest colors */
|
|
for (j = 0; j < nv; j++) {
|
|
if (cnt[j]) {
|
|
float div = 1.0F / cnt[j];
|
|
for (i = 0; i < nc; i++) {
|
|
vec[j][i] = div * sum[j][i];
|
|
}
|
|
} else {
|
|
/* this vec has no samples or is identical with a previous vec */
|
|
int32_t worst = fxt1_worst(vec[j], input, nc, n);
|
|
for (i = 0; i < nc; i++) {
|
|
vec[j][i] = input[worst][i];
|
|
}
|
|
}
|
|
}
|
|
}
|
|
|
|
return 0; /* could not converge fast enough */
|
|
}
|
|
|
|
|
|
static void
|
|
fxt1_quantize_CHROMA (uint32_t *cc,
|
|
uint8_t input[N_TEXELS][MAX_COMP])
|
|
{
|
|
const int32_t n_vect = 4; /* 4 base vectors to find */
|
|
const int32_t n_comp = 3; /* 3 components: R, G, B */
|
|
float vec[MAX_VECT][MAX_COMP];
|
|
int32_t i, j, k;
|
|
Fx64 hi; /* high quadword */
|
|
uint32_t lohi, lolo; /* low quadword: hi dword, lo dword */
|
|
|
|
if (fxt1_choose(vec, n_vect, input, n_comp, N_TEXELS) != 0) {
|
|
fxt1_lloyd(vec, n_vect, input, n_comp, N_TEXELS);
|
|
}
|
|
|
|
FX64_MOV32(hi, 4); /* cc-chroma = "010" + unused bit */
|
|
for (j = n_vect - 1; j >= 0; j--) {
|
|
for (i = 0; i < n_comp; i++) {
|
|
/* add in colors */
|
|
FX64_SHL(hi, 5);
|
|
FX64_OR32(hi, (uint32_t)(vec[j][i] / 8.0F));
|
|
}
|
|
}
|
|
((Fx64 *)cc)[1] = hi;
|
|
|
|
lohi = lolo = 0;
|
|
/* right microtile */
|
|
for (k = N_TEXELS - 1; k >= N_TEXELS/2; k--) {
|
|
lohi <<= 2;
|
|
lohi |= fxt1_bestcol(vec, n_vect, input[k], n_comp);
|
|
}
|
|
/* left microtile */
|
|
for (; k >= 0; k--) {
|
|
lolo <<= 2;
|
|
lolo |= fxt1_bestcol(vec, n_vect, input[k], n_comp);
|
|
}
|
|
cc[1] = lohi;
|
|
cc[0] = lolo;
|
|
}
|
|
|
|
|
|
static void
|
|
fxt1_quantize_ALPHA0 (uint32_t *cc,
|
|
uint8_t input[N_TEXELS][MAX_COMP],
|
|
uint8_t reord[N_TEXELS][MAX_COMP], int32_t n)
|
|
{
|
|
const int32_t n_vect = 3; /* 3 base vectors to find */
|
|
const int32_t n_comp = 4; /* 4 components: R, G, B, A */
|
|
float vec[MAX_VECT][MAX_COMP];
|
|
int32_t i, j, k;
|
|
Fx64 hi; /* high quadword */
|
|
uint32_t lohi, lolo; /* low quadword: hi dword, lo dword */
|
|
|
|
/* the last vector indicates zero */
|
|
for (i = 0; i < n_comp; i++) {
|
|
vec[n_vect][i] = 0;
|
|
}
|
|
|
|
/* the first n texels in reord are guaranteed to be non-zero */
|
|
if (fxt1_choose(vec, n_vect, reord, n_comp, n) != 0) {
|
|
fxt1_lloyd(vec, n_vect, reord, n_comp, n);
|
|
}
|
|
|
|
FX64_MOV32(hi, 6); /* alpha = "011" + lerp = 0 */
|
|
for (j = n_vect - 1; j >= 0; j--) {
|
|
/* add in alphas */
|
|
FX64_SHL(hi, 5);
|
|
FX64_OR32(hi, (uint32_t)(vec[j][ACOMP] / 8.0F));
|
|
}
|
|
for (j = n_vect - 1; j >= 0; j--) {
|
|
for (i = 0; i < n_comp - 1; i++) {
|
|
/* add in colors */
|
|
FX64_SHL(hi, 5);
|
|
FX64_OR32(hi, (uint32_t)(vec[j][i] / 8.0F));
|
|
}
|
|
}
|
|
((Fx64 *)cc)[1] = hi;
|
|
|
|
lohi = lolo = 0;
|
|
/* right microtile */
|
|
for (k = N_TEXELS - 1; k >= N_TEXELS/2; k--) {
|
|
lohi <<= 2;
|
|
lohi |= fxt1_bestcol(vec, n_vect + 1, input[k], n_comp);
|
|
}
|
|
/* left microtile */
|
|
for (; k >= 0; k--) {
|
|
lolo <<= 2;
|
|
lolo |= fxt1_bestcol(vec, n_vect + 1, input[k], n_comp);
|
|
}
|
|
cc[1] = lohi;
|
|
cc[0] = lolo;
|
|
}
|
|
|
|
|
|
static void
|
|
fxt1_quantize_ALPHA1 (uint32_t *cc,
|
|
uint8_t input[N_TEXELS][MAX_COMP])
|
|
{
|
|
const int32_t n_vect = 3; /* highest vector number in each microtile */
|
|
const int32_t n_comp = 4; /* 4 components: R, G, B, A */
|
|
float vec[1 + 1 + 1][MAX_COMP]; /* 1.5 extrema for each sub-block */
|
|
float b, iv[MAX_COMP]; /* interpolation vector */
|
|
int32_t i, j, k;
|
|
Fx64 hi; /* high quadword */
|
|
uint32_t lohi, lolo; /* low quadword: hi dword, lo dword */
|
|
|
|
int32_t minSum;
|
|
int32_t maxSum;
|
|
int32_t minColL = 0, maxColL = 0;
|
|
int32_t minColR = 0, maxColR = 0;
|
|
int32_t sumL = 0, sumR = 0;
|
|
int32_t nn_comp;
|
|
/* Our solution here is to find the darkest and brightest colors in
|
|
* the 4x4 tile and use those as the two representative colors.
|
|
* There are probably better algorithms to use (histogram-based).
|
|
*/
|
|
nn_comp = n_comp;
|
|
while ((minColL == maxColL) && nn_comp) {
|
|
minSum = 2000; /* big enough */
|
|
maxSum = -1; /* small enough */
|
|
for (k = 0; k < N_TEXELS / 2; k++) {
|
|
int32_t sum = 0;
|
|
for (i = 0; i < nn_comp; i++) {
|
|
sum += input[k][i];
|
|
}
|
|
if (minSum > sum) {
|
|
minSum = sum;
|
|
minColL = k;
|
|
}
|
|
if (maxSum < sum) {
|
|
maxSum = sum;
|
|
maxColL = k;
|
|
}
|
|
sumL += sum;
|
|
}
|
|
|
|
nn_comp--;
|
|
}
|
|
|
|
nn_comp = n_comp;
|
|
while ((minColR == maxColR) && nn_comp) {
|
|
minSum = 2000; /* big enough */
|
|
maxSum = -1; /* small enough */
|
|
for (k = N_TEXELS / 2; k < N_TEXELS; k++) {
|
|
int32_t sum = 0;
|
|
for (i = 0; i < nn_comp; i++) {
|
|
sum += input[k][i];
|
|
}
|
|
if (minSum > sum) {
|
|
minSum = sum;
|
|
minColR = k;
|
|
}
|
|
if (maxSum < sum) {
|
|
maxSum = sum;
|
|
maxColR = k;
|
|
}
|
|
sumR += sum;
|
|
}
|
|
|
|
nn_comp--;
|
|
}
|
|
|
|
/* choose the common vector (yuck!) */
|
|
{
|
|
int32_t j1, j2;
|
|
int32_t v1 = 0, v2 = 0;
|
|
float err = 1e9; /* big enough */
|
|
float tv[2 * 2][MAX_COMP]; /* 2 extrema for each sub-block */
|
|
for (i = 0; i < n_comp; i++) {
|
|
tv[0][i] = input[minColL][i];
|
|
tv[1][i] = input[maxColL][i];
|
|
tv[2][i] = input[minColR][i];
|
|
tv[3][i] = input[maxColR][i];
|
|
}
|
|
for (j1 = 0; j1 < 2; j1++) {
|
|
for (j2 = 2; j2 < 4; j2++) {
|
|
float e = 0.0F;
|
|
for (i = 0; i < n_comp; i++) {
|
|
e += (tv[j1][i] - tv[j2][i]) * (tv[j1][i] - tv[j2][i]);
|
|
}
|
|
if (e < err) {
|
|
err = e;
|
|
v1 = j1;
|
|
v2 = j2;
|
|
}
|
|
}
|
|
}
|
|
for (i = 0; i < n_comp; i++) {
|
|
vec[0][i] = tv[1 - v1][i];
|
|
vec[1][i] = (tv[v1][i] * sumL + tv[v2][i] * sumR) / (sumL + sumR);
|
|
vec[2][i] = tv[5 - v2][i];
|
|
}
|
|
}
|
|
|
|
/* left microtile */
|
|
cc[0] = 0;
|
|
if (minColL != maxColL) {
|
|
/* compute interpolation vector */
|
|
MAKEIVEC(n_vect, n_comp, iv, b, vec[0], vec[1]);
|
|
|
|
/* add in texels */
|
|
lolo = 0;
|
|
for (k = N_TEXELS / 2 - 1; k >= 0; k--) {
|
|
int32_t texel;
|
|
/* interpolate color */
|
|
CALCCDOT(texel, n_vect, n_comp, iv, b, input[k]);
|
|
/* add in texel */
|
|
lolo <<= 2;
|
|
lolo |= texel;
|
|
}
|
|
|
|
cc[0] = lolo;
|
|
}
|
|
|
|
/* right microtile */
|
|
cc[1] = 0;
|
|
if (minColR != maxColR) {
|
|
/* compute interpolation vector */
|
|
MAKEIVEC(n_vect, n_comp, iv, b, vec[2], vec[1]);
|
|
|
|
/* add in texels */
|
|
lohi = 0;
|
|
for (k = N_TEXELS - 1; k >= N_TEXELS / 2; k--) {
|
|
int32_t texel;
|
|
/* interpolate color */
|
|
CALCCDOT(texel, n_vect, n_comp, iv, b, input[k]);
|
|
/* add in texel */
|
|
lohi <<= 2;
|
|
lohi |= texel;
|
|
}
|
|
|
|
cc[1] = lohi;
|
|
}
|
|
|
|
FX64_MOV32(hi, 7); /* alpha = "011" + lerp = 1 */
|
|
for (j = n_vect - 1; j >= 0; j--) {
|
|
/* add in alphas */
|
|
FX64_SHL(hi, 5);
|
|
FX64_OR32(hi, (uint32_t)(vec[j][ACOMP] / 8.0F));
|
|
}
|
|
for (j = n_vect - 1; j >= 0; j--) {
|
|
for (i = 0; i < n_comp - 1; i++) {
|
|
/* add in colors */
|
|
FX64_SHL(hi, 5);
|
|
FX64_OR32(hi, (uint32_t)(vec[j][i] / 8.0F));
|
|
}
|
|
}
|
|
((Fx64 *)cc)[1] = hi;
|
|
}
|
|
|
|
|
|
static void
|
|
fxt1_quantize_HI (uint32_t *cc,
|
|
uint8_t input[N_TEXELS][MAX_COMP],
|
|
uint8_t reord[N_TEXELS][MAX_COMP], int32_t n)
|
|
{
|
|
const int32_t n_vect = 6; /* highest vector number */
|
|
const int32_t n_comp = 3; /* 3 components: R, G, B */
|
|
float b = 0.0F; /* phoudoin: silent compiler! */
|
|
float iv[MAX_COMP]; /* interpolation vector */
|
|
int32_t i, k;
|
|
uint32_t hihi; /* high quadword: hi dword */
|
|
|
|
int32_t minSum = 2000; /* big enough */
|
|
int32_t maxSum = -1; /* small enough */
|
|
int32_t minCol = 0; /* phoudoin: silent compiler! */
|
|
int32_t maxCol = 0; /* phoudoin: silent compiler! */
|
|
|
|
/* Our solution here is to find the darkest and brightest colors in
|
|
* the 8x4 tile and use those as the two representative colors.
|
|
* There are probably better algorithms to use (histogram-based).
|
|
*/
|
|
for (k = 0; k < n; k++) {
|
|
int32_t sum = 0;
|
|
for (i = 0; i < n_comp; i++) {
|
|
sum += reord[k][i];
|
|
}
|
|
if (minSum > sum) {
|
|
minSum = sum;
|
|
minCol = k;
|
|
}
|
|
if (maxSum < sum) {
|
|
maxSum = sum;
|
|
maxCol = k;
|
|
}
|
|
}
|
|
|
|
hihi = 0; /* cc-hi = "00" */
|
|
for (i = 0; i < n_comp; i++) {
|
|
/* add in colors */
|
|
hihi <<= 5;
|
|
hihi |= reord[maxCol][i] >> 3;
|
|
}
|
|
for (i = 0; i < n_comp; i++) {
|
|
/* add in colors */
|
|
hihi <<= 5;
|
|
hihi |= reord[minCol][i] >> 3;
|
|
}
|
|
cc[3] = hihi;
|
|
cc[0] = cc[1] = cc[2] = 0;
|
|
|
|
/* compute interpolation vector */
|
|
if (minCol != maxCol) {
|
|
MAKEIVEC(n_vect, n_comp, iv, b, reord[minCol], reord[maxCol]);
|
|
}
|
|
|
|
/* add in texels */
|
|
for (k = N_TEXELS - 1; k >= 0; k--) {
|
|
int32_t t = k * 3;
|
|
uint32_t *kk = (uint32_t *)((char *)cc + t / 8);
|
|
int32_t texel = n_vect + 1; /* transparent black */
|
|
|
|
if (!ISTBLACK(input[k])) {
|
|
if (minCol != maxCol) {
|
|
/* interpolate color */
|
|
CALCCDOT(texel, n_vect, n_comp, iv, b, input[k]);
|
|
/* add in texel */
|
|
kk[0] |= texel << (t & 7);
|
|
}
|
|
} else {
|
|
/* add in texel */
|
|
kk[0] |= texel << (t & 7);
|
|
}
|
|
}
|
|
}
|
|
|
|
|
|
static void
|
|
fxt1_quantize_MIXED1 (uint32_t *cc,
|
|
uint8_t input[N_TEXELS][MAX_COMP])
|
|
{
|
|
const int32_t n_vect = 2; /* highest vector number in each microtile */
|
|
const int32_t n_comp = 3; /* 3 components: R, G, B */
|
|
uint8_t vec[2 * 2][MAX_COMP]; /* 2 extrema for each sub-block */
|
|
float b, iv[MAX_COMP]; /* interpolation vector */
|
|
int32_t i, j, k;
|
|
Fx64 hi; /* high quadword */
|
|
uint32_t lohi, lolo; /* low quadword: hi dword, lo dword */
|
|
|
|
int32_t minSum;
|
|
int32_t maxSum;
|
|
int32_t minColL = 0, maxColL = -1;
|
|
int32_t minColR = 0, maxColR = -1;
|
|
|
|
/* Our solution here is to find the darkest and brightest colors in
|
|
* the 4x4 tile and use those as the two representative colors.
|
|
* There are probably better algorithms to use (histogram-based).
|
|
*/
|
|
minSum = 2000; /* big enough */
|
|
maxSum = -1; /* small enough */
|
|
for (k = 0; k < N_TEXELS / 2; k++) {
|
|
if (!ISTBLACK(input[k])) {
|
|
int32_t sum = 0;
|
|
for (i = 0; i < n_comp; i++) {
|
|
sum += input[k][i];
|
|
}
|
|
if (minSum > sum) {
|
|
minSum = sum;
|
|
minColL = k;
|
|
}
|
|
if (maxSum < sum) {
|
|
maxSum = sum;
|
|
maxColL = k;
|
|
}
|
|
}
|
|
}
|
|
minSum = 2000; /* big enough */
|
|
maxSum = -1; /* small enough */
|
|
for (; k < N_TEXELS; k++) {
|
|
if (!ISTBLACK(input[k])) {
|
|
int32_t sum = 0;
|
|
for (i = 0; i < n_comp; i++) {
|
|
sum += input[k][i];
|
|
}
|
|
if (minSum > sum) {
|
|
minSum = sum;
|
|
minColR = k;
|
|
}
|
|
if (maxSum < sum) {
|
|
maxSum = sum;
|
|
maxColR = k;
|
|
}
|
|
}
|
|
}
|
|
|
|
/* left microtile */
|
|
if (maxColL == -1) {
|
|
/* all transparent black */
|
|
cc[0] = ~0u;
|
|
for (i = 0; i < n_comp; i++) {
|
|
vec[0][i] = 0;
|
|
vec[1][i] = 0;
|
|
}
|
|
} else {
|
|
cc[0] = 0;
|
|
for (i = 0; i < n_comp; i++) {
|
|
vec[0][i] = input[minColL][i];
|
|
vec[1][i] = input[maxColL][i];
|
|
}
|
|
if (minColL != maxColL) {
|
|
/* compute interpolation vector */
|
|
MAKEIVEC(n_vect, n_comp, iv, b, vec[0], vec[1]);
|
|
|
|
/* add in texels */
|
|
lolo = 0;
|
|
for (k = N_TEXELS / 2 - 1; k >= 0; k--) {
|
|
int32_t texel = n_vect + 1; /* transparent black */
|
|
if (!ISTBLACK(input[k])) {
|
|
/* interpolate color */
|
|
CALCCDOT(texel, n_vect, n_comp, iv, b, input[k]);
|
|
}
|
|
/* add in texel */
|
|
lolo <<= 2;
|
|
lolo |= texel;
|
|
}
|
|
cc[0] = lolo;
|
|
}
|
|
}
|
|
|
|
/* right microtile */
|
|
if (maxColR == -1) {
|
|
/* all transparent black */
|
|
cc[1] = ~0u;
|
|
for (i = 0; i < n_comp; i++) {
|
|
vec[2][i] = 0;
|
|
vec[3][i] = 0;
|
|
}
|
|
} else {
|
|
cc[1] = 0;
|
|
for (i = 0; i < n_comp; i++) {
|
|
vec[2][i] = input[minColR][i];
|
|
vec[3][i] = input[maxColR][i];
|
|
}
|
|
if (minColR != maxColR) {
|
|
/* compute interpolation vector */
|
|
MAKEIVEC(n_vect, n_comp, iv, b, vec[2], vec[3]);
|
|
|
|
/* add in texels */
|
|
lohi = 0;
|
|
for (k = N_TEXELS - 1; k >= N_TEXELS / 2; k--) {
|
|
int32_t texel = n_vect + 1; /* transparent black */
|
|
if (!ISTBLACK(input[k])) {
|
|
/* interpolate color */
|
|
CALCCDOT(texel, n_vect, n_comp, iv, b, input[k]);
|
|
}
|
|
/* add in texel */
|
|
lohi <<= 2;
|
|
lohi |= texel;
|
|
}
|
|
cc[1] = lohi;
|
|
}
|
|
}
|
|
|
|
FX64_MOV32(hi, 9 | (vec[3][GCOMP] & 4) | ((vec[1][GCOMP] >> 1) & 2)); /* chroma = "1" */
|
|
for (j = 2 * 2 - 1; j >= 0; j--) {
|
|
for (i = 0; i < n_comp; i++) {
|
|
/* add in colors */
|
|
FX64_SHL(hi, 5);
|
|
FX64_OR32(hi, vec[j][i] >> 3);
|
|
}
|
|
}
|
|
((Fx64 *)cc)[1] = hi;
|
|
}
|
|
|
|
|
|
static void
|
|
fxt1_quantize_MIXED0 (uint32_t *cc,
|
|
uint8_t input[N_TEXELS][MAX_COMP])
|
|
{
|
|
const int32_t n_vect = 3; /* highest vector number in each microtile */
|
|
const int32_t n_comp = 3; /* 3 components: R, G, B */
|
|
uint8_t vec[2 * 2][MAX_COMP]; /* 2 extrema for each sub-block */
|
|
float b, iv[MAX_COMP]; /* interpolation vector */
|
|
int32_t i, j, k;
|
|
Fx64 hi; /* high quadword */
|
|
uint32_t lohi, lolo; /* low quadword: hi dword, lo dword */
|
|
|
|
int32_t minColL = 0, maxColL = 0;
|
|
int32_t minColR = 0, maxColR = 0;
|
|
#if 0
|
|
int32_t minSum;
|
|
int32_t maxSum;
|
|
|
|
/* Our solution here is to find the darkest and brightest colors in
|
|
* the 4x4 tile and use those as the two representative colors.
|
|
* There are probably better algorithms to use (histogram-based).
|
|
*/
|
|
minSum = 2000; /* big enough */
|
|
maxSum = -1; /* small enough */
|
|
for (k = 0; k < N_TEXELS / 2; k++) {
|
|
int32_t sum = 0;
|
|
for (i = 0; i < n_comp; i++) {
|
|
sum += input[k][i];
|
|
}
|
|
if (minSum > sum) {
|
|
minSum = sum;
|
|
minColL = k;
|
|
}
|
|
if (maxSum < sum) {
|
|
maxSum = sum;
|
|
maxColL = k;
|
|
}
|
|
}
|
|
minSum = 2000; /* big enough */
|
|
maxSum = -1; /* small enough */
|
|
for (; k < N_TEXELS; k++) {
|
|
int32_t sum = 0;
|
|
for (i = 0; i < n_comp; i++) {
|
|
sum += input[k][i];
|
|
}
|
|
if (minSum > sum) {
|
|
minSum = sum;
|
|
minColR = k;
|
|
}
|
|
if (maxSum < sum) {
|
|
maxSum = sum;
|
|
maxColR = k;
|
|
}
|
|
}
|
|
#else
|
|
int32_t minVal;
|
|
int32_t maxVal;
|
|
int32_t maxVarL = fxt1_variance(input, n_comp);
|
|
int32_t maxVarR = fxt1_variance(&input[N_TEXELS / 2], n_comp);
|
|
|
|
/* Scan the channel with max variance for lo & hi
|
|
* and use those as the two representative colors.
|
|
*/
|
|
minVal = 2000; /* big enough */
|
|
maxVal = -1; /* small enough */
|
|
for (k = 0; k < N_TEXELS / 2; k++) {
|
|
int32_t t = input[k][maxVarL];
|
|
if (minVal > t) {
|
|
minVal = t;
|
|
minColL = k;
|
|
}
|
|
if (maxVal < t) {
|
|
maxVal = t;
|
|
maxColL = k;
|
|
}
|
|
}
|
|
minVal = 2000; /* big enough */
|
|
maxVal = -1; /* small enough */
|
|
for (; k < N_TEXELS; k++) {
|
|
int32_t t = input[k][maxVarR];
|
|
if (minVal > t) {
|
|
minVal = t;
|
|
minColR = k;
|
|
}
|
|
if (maxVal < t) {
|
|
maxVal = t;
|
|
maxColR = k;
|
|
}
|
|
}
|
|
#endif
|
|
|
|
/* left microtile */
|
|
cc[0] = 0;
|
|
for (i = 0; i < n_comp; i++) {
|
|
vec[0][i] = input[minColL][i];
|
|
vec[1][i] = input[maxColL][i];
|
|
}
|
|
if (minColL != maxColL) {
|
|
/* compute interpolation vector */
|
|
MAKEIVEC(n_vect, n_comp, iv, b, vec[0], vec[1]);
|
|
|
|
/* add in texels */
|
|
lolo = 0;
|
|
for (k = N_TEXELS / 2 - 1; k >= 0; k--) {
|
|
int32_t texel;
|
|
/* interpolate color */
|
|
CALCCDOT(texel, n_vect, n_comp, iv, b, input[k]);
|
|
/* add in texel */
|
|
lolo <<= 2;
|
|
lolo |= texel;
|
|
}
|
|
|
|
/* funky encoding for LSB of green */
|
|
if ((int32_t)((lolo >> 1) & 1) != (((vec[1][GCOMP] ^ vec[0][GCOMP]) >> 2) & 1)) {
|
|
for (i = 0; i < n_comp; i++) {
|
|
vec[1][i] = input[minColL][i];
|
|
vec[0][i] = input[maxColL][i];
|
|
}
|
|
lolo = ~lolo;
|
|
}
|
|
|
|
cc[0] = lolo;
|
|
}
|
|
|
|
/* right microtile */
|
|
cc[1] = 0;
|
|
for (i = 0; i < n_comp; i++) {
|
|
vec[2][i] = input[minColR][i];
|
|
vec[3][i] = input[maxColR][i];
|
|
}
|
|
if (minColR != maxColR) {
|
|
/* compute interpolation vector */
|
|
MAKEIVEC(n_vect, n_comp, iv, b, vec[2], vec[3]);
|
|
|
|
/* add in texels */
|
|
lohi = 0;
|
|
for (k = N_TEXELS - 1; k >= N_TEXELS / 2; k--) {
|
|
int32_t texel;
|
|
/* interpolate color */
|
|
CALCCDOT(texel, n_vect, n_comp, iv, b, input[k]);
|
|
/* add in texel */
|
|
lohi <<= 2;
|
|
lohi |= texel;
|
|
}
|
|
|
|
/* funky encoding for LSB of green */
|
|
if ((int32_t)((lohi >> 1) & 1) != (((vec[3][GCOMP] ^ vec[2][GCOMP]) >> 2) & 1)) {
|
|
for (i = 0; i < n_comp; i++) {
|
|
vec[3][i] = input[minColR][i];
|
|
vec[2][i] = input[maxColR][i];
|
|
}
|
|
lohi = ~lohi;
|
|
}
|
|
|
|
cc[1] = lohi;
|
|
}
|
|
|
|
FX64_MOV32(hi, 8 | (vec[3][GCOMP] & 4) | ((vec[1][GCOMP] >> 1) & 2)); /* chroma = "1" */
|
|
for (j = 2 * 2 - 1; j >= 0; j--) {
|
|
for (i = 0; i < n_comp; i++) {
|
|
/* add in colors */
|
|
FX64_SHL(hi, 5);
|
|
FX64_OR32(hi, vec[j][i] >> 3);
|
|
}
|
|
}
|
|
((Fx64 *)cc)[1] = hi;
|
|
}
|
|
|
|
|
|
static void
|
|
fxt1_quantize (uint32_t *cc, const uint8_t *lines[], int32_t comps)
|
|
{
|
|
int32_t trualpha;
|
|
uint8_t reord[N_TEXELS][MAX_COMP];
|
|
|
|
uint8_t input[N_TEXELS][MAX_COMP];
|
|
int32_t i, k, l;
|
|
|
|
if (comps == 3) {
|
|
/* make the whole block opaque */
|
|
memset(input, -1, sizeof(input));
|
|
}
|
|
|
|
/* 8 texels each line */
|
|
for (l = 0; l < 4; l++) {
|
|
for (k = 0; k < 4; k++) {
|
|
for (i = 0; i < comps; i++) {
|
|
input[k + l * 4][i] = *lines[l]++;
|
|
}
|
|
}
|
|
for (; k < 8; k++) {
|
|
for (i = 0; i < comps; i++) {
|
|
input[k + l * 4 + 12][i] = *lines[l]++;
|
|
}
|
|
}
|
|
}
|
|
|
|
/* block layout:
|
|
* 00, 01, 02, 03, 08, 09, 0a, 0b
|
|
* 10, 11, 12, 13, 18, 19, 1a, 1b
|
|
* 04, 05, 06, 07, 0c, 0d, 0e, 0f
|
|
* 14, 15, 16, 17, 1c, 1d, 1e, 1f
|
|
*/
|
|
|
|
/* [dBorca]
|
|
* stupidity flows forth from this
|
|
*/
|
|
l = N_TEXELS;
|
|
trualpha = 0;
|
|
if (comps == 4) {
|
|
/* skip all transparent black texels */
|
|
l = 0;
|
|
for (k = 0; k < N_TEXELS; k++) {
|
|
/* test all components against 0 */
|
|
if (!ISTBLACK(input[k])) {
|
|
/* texel is not transparent black */
|
|
memcpy(reord[l], input[k], 4);
|
|
if (reord[l][ACOMP] < (255 - ALPHA_TS)) {
|
|
/* non-opaque texel */
|
|
trualpha = !0;
|
|
}
|
|
l++;
|
|
}
|
|
}
|
|
}
|
|
|
|
#if 0
|
|
if (trualpha) {
|
|
fxt1_quantize_ALPHA0(cc, input, reord, l);
|
|
} else if (l == 0) {
|
|
cc[0] = cc[1] = cc[2] = -1;
|
|
cc[3] = 0;
|
|
} else if (l < N_TEXELS) {
|
|
fxt1_quantize_HI(cc, input, reord, l);
|
|
} else {
|
|
fxt1_quantize_CHROMA(cc, input);
|
|
}
|
|
(void)fxt1_quantize_ALPHA1;
|
|
(void)fxt1_quantize_MIXED1;
|
|
(void)fxt1_quantize_MIXED0;
|
|
#else
|
|
if (trualpha) {
|
|
fxt1_quantize_ALPHA1(cc, input);
|
|
} else if (l == 0) {
|
|
cc[0] = cc[1] = cc[2] = ~0u;
|
|
cc[3] = 0;
|
|
} else if (l < N_TEXELS) {
|
|
fxt1_quantize_MIXED1(cc, input);
|
|
} else {
|
|
fxt1_quantize_MIXED0(cc, input);
|
|
}
|
|
(void)fxt1_quantize_ALPHA0;
|
|
(void)fxt1_quantize_HI;
|
|
(void)fxt1_quantize_CHROMA;
|
|
#endif
|
|
}
|
|
|
|
|
|
|
|
/**
|
|
* Upscale an image by replication, not (typical) stretching.
|
|
* We use this when the image width or height is less than a
|
|
* certain size (4, 8) and we need to upscale an image.
|
|
*/
|
|
static void
|
|
upscale_teximage2d(int32_t inWidth, int32_t inHeight,
|
|
int32_t outWidth, int32_t outHeight,
|
|
int32_t comps, const uint8_t *src, int32_t srcRowStride,
|
|
uint8_t *dest )
|
|
{
|
|
int32_t i, j, k;
|
|
|
|
assert(outWidth >= inWidth);
|
|
assert(outHeight >= inHeight);
|
|
#if 0
|
|
assert(inWidth == 1 || inWidth == 2 || inHeight == 1 || inHeight == 2);
|
|
assert((outWidth & 3) == 0);
|
|
assert((outHeight & 3) == 0);
|
|
#endif
|
|
|
|
for (i = 0; i < outHeight; i++) {
|
|
const int32_t ii = i % inHeight;
|
|
for (j = 0; j < outWidth; j++) {
|
|
const int32_t jj = j % inWidth;
|
|
for (k = 0; k < comps; k++) {
|
|
dest[(i * outWidth + j) * comps + k]
|
|
= src[ii * srcRowStride + jj * comps + k];
|
|
}
|
|
}
|
|
}
|
|
}
|
|
|
|
|
|
static void
|
|
fxt1_encode (uint32_t width, uint32_t height, int32_t comps,
|
|
const void *source, int32_t srcRowStride,
|
|
void *dest, int32_t destRowStride)
|
|
{
|
|
uint32_t x, y;
|
|
const uint8_t *data;
|
|
uint32_t *encoded = (uint32_t *)dest;
|
|
void *newSource = NULL;
|
|
|
|
assert(comps == 3 || comps == 4);
|
|
|
|
/* Replicate image if width is not M8 or height is not M4 */
|
|
if ((width & 7) | (height & 3)) {
|
|
int32_t newWidth = (width + 7) & ~7;
|
|
int32_t newHeight = (height + 3) & ~3;
|
|
newSource = malloc(comps * newWidth * newHeight * sizeof(uint8_t));
|
|
if (!newSource)
|
|
return;
|
|
upscale_teximage2d(width, height, newWidth, newHeight,
|
|
comps, (const uint8_t *) source,
|
|
srcRowStride, (uint8_t *) newSource);
|
|
source = newSource;
|
|
width = newWidth;
|
|
height = newHeight;
|
|
srcRowStride = comps * newWidth;
|
|
}
|
|
|
|
data = (const uint8_t *) source;
|
|
destRowStride = (destRowStride - width * 2) / 4;
|
|
for (y = 0; y < height; y += 4) {
|
|
uint32_t offs = 0 + (y + 0) * srcRowStride;
|
|
for (x = 0; x < width; x += 8) {
|
|
const uint8_t *lines[4];
|
|
lines[0] = &data[offs];
|
|
lines[1] = lines[0] + srcRowStride;
|
|
lines[2] = lines[1] + srcRowStride;
|
|
lines[3] = lines[2] + srcRowStride;
|
|
offs += 8 * comps;
|
|
fxt1_quantize(encoded, lines, comps);
|
|
/* 128 bits per 8x4 block */
|
|
encoded += 4;
|
|
}
|
|
encoded += destRowStride;
|
|
}
|
|
|
|
free(newSource);
|
|
}
|
|
|
|
|
|
/***************************************************************************\
|
|
* FXT1 decoder
|
|
*
|
|
* The decoder is based on GL_3DFX_texture_compression_FXT1
|
|
* specification and serves as a concept for the encoder.
|
|
\***************************************************************************/
|
|
|
|
|
|
/* lookup table for scaling 5 bit colors up to 8 bits */
|
|
static const uint8_t _rgb_scale_5[] = {
|
|
0, 8, 16, 25, 33, 41, 49, 58,
|
|
66, 74, 82, 90, 99, 107, 115, 123,
|
|
132, 140, 148, 156, 165, 173, 181, 189,
|
|
197, 206, 214, 222, 230, 239, 247, 255
|
|
};
|
|
|
|
/* lookup table for scaling 6 bit colors up to 8 bits */
|
|
static const uint8_t _rgb_scale_6[] = {
|
|
0, 4, 8, 12, 16, 20, 24, 28,
|
|
32, 36, 40, 45, 49, 53, 57, 61,
|
|
65, 69, 73, 77, 81, 85, 89, 93,
|
|
97, 101, 105, 109, 113, 117, 121, 125,
|
|
130, 134, 138, 142, 146, 150, 154, 158,
|
|
162, 166, 170, 174, 178, 182, 186, 190,
|
|
194, 198, 202, 206, 210, 215, 219, 223,
|
|
227, 231, 235, 239, 243, 247, 251, 255
|
|
};
|
|
|
|
|
|
#define CC_SEL(cc, which) (((uint32_t *)(cc))[(which) / 32] >> ((which) & 31))
|
|
#define UP5(c) _rgb_scale_5[(c) & 31]
|
|
#define UP6(c, b) _rgb_scale_6[(((c) & 31) << 1) | ((b) & 1)]
|
|
#define LERP(n, t, c0, c1) (((n) - (t)) * (c0) + (t) * (c1) + (n) / 2) / (n)
|
|
|
|
|
|
static void
|
|
fxt1_decode_1HI (const uint8_t *code, int32_t t, uint8_t *rgba)
|
|
{
|
|
const uint32_t *cc;
|
|
|
|
t *= 3;
|
|
cc = (const uint32_t *)(code + t / 8);
|
|
t = (cc[0] >> (t & 7)) & 7;
|
|
|
|
if (t == 7) {
|
|
rgba[RCOMP] = rgba[GCOMP] = rgba[BCOMP] = rgba[ACOMP] = 0;
|
|
} else {
|
|
uint8_t r, g, b;
|
|
cc = (const uint32_t *)(code + 12);
|
|
if (t == 0) {
|
|
b = UP5(CC_SEL(cc, 0));
|
|
g = UP5(CC_SEL(cc, 5));
|
|
r = UP5(CC_SEL(cc, 10));
|
|
} else if (t == 6) {
|
|
b = UP5(CC_SEL(cc, 15));
|
|
g = UP5(CC_SEL(cc, 20));
|
|
r = UP5(CC_SEL(cc, 25));
|
|
} else {
|
|
b = LERP(6, t, UP5(CC_SEL(cc, 0)), UP5(CC_SEL(cc, 15)));
|
|
g = LERP(6, t, UP5(CC_SEL(cc, 5)), UP5(CC_SEL(cc, 20)));
|
|
r = LERP(6, t, UP5(CC_SEL(cc, 10)), UP5(CC_SEL(cc, 25)));
|
|
}
|
|
rgba[RCOMP] = r;
|
|
rgba[GCOMP] = g;
|
|
rgba[BCOMP] = b;
|
|
rgba[ACOMP] = 255;
|
|
}
|
|
}
|
|
|
|
|
|
static void
|
|
fxt1_decode_1CHROMA (const uint8_t *code, int32_t t, uint8_t *rgba)
|
|
{
|
|
const uint32_t *cc;
|
|
uint32_t kk;
|
|
|
|
cc = (const uint32_t *)code;
|
|
if (t & 16) {
|
|
cc++;
|
|
t &= 15;
|
|
}
|
|
t = (cc[0] >> (t * 2)) & 3;
|
|
|
|
t *= 15;
|
|
cc = (const uint32_t *)(code + 8 + t / 8);
|
|
kk = cc[0] >> (t & 7);
|
|
rgba[BCOMP] = UP5(kk);
|
|
rgba[GCOMP] = UP5(kk >> 5);
|
|
rgba[RCOMP] = UP5(kk >> 10);
|
|
rgba[ACOMP] = 255;
|
|
}
|
|
|
|
|
|
static void
|
|
fxt1_decode_1MIXED (const uint8_t *code, int32_t t, uint8_t *rgba)
|
|
{
|
|
const uint32_t *cc;
|
|
uint32_t col[2][3];
|
|
int32_t glsb, selb;
|
|
|
|
cc = (const uint32_t *)code;
|
|
if (t & 16) {
|
|
t &= 15;
|
|
t = (cc[1] >> (t * 2)) & 3;
|
|
/* col 2 */
|
|
col[0][BCOMP] = (*(const uint32_t *)(code + 11)) >> 6;
|
|
col[0][GCOMP] = CC_SEL(cc, 99);
|
|
col[0][RCOMP] = CC_SEL(cc, 104);
|
|
/* col 3 */
|
|
col[1][BCOMP] = CC_SEL(cc, 109);
|
|
col[1][GCOMP] = CC_SEL(cc, 114);
|
|
col[1][RCOMP] = CC_SEL(cc, 119);
|
|
glsb = CC_SEL(cc, 126);
|
|
selb = CC_SEL(cc, 33);
|
|
} else {
|
|
t = (cc[0] >> (t * 2)) & 3;
|
|
/* col 0 */
|
|
col[0][BCOMP] = CC_SEL(cc, 64);
|
|
col[0][GCOMP] = CC_SEL(cc, 69);
|
|
col[0][RCOMP] = CC_SEL(cc, 74);
|
|
/* col 1 */
|
|
col[1][BCOMP] = CC_SEL(cc, 79);
|
|
col[1][GCOMP] = CC_SEL(cc, 84);
|
|
col[1][RCOMP] = CC_SEL(cc, 89);
|
|
glsb = CC_SEL(cc, 125);
|
|
selb = CC_SEL(cc, 1);
|
|
}
|
|
|
|
if (CC_SEL(cc, 124) & 1) {
|
|
/* alpha[0] == 1 */
|
|
|
|
if (t == 3) {
|
|
/* zero */
|
|
rgba[RCOMP] = rgba[BCOMP] = rgba[GCOMP] = rgba[ACOMP] = 0;
|
|
} else {
|
|
uint8_t r, g, b;
|
|
if (t == 0) {
|
|
b = UP5(col[0][BCOMP]);
|
|
g = UP5(col[0][GCOMP]);
|
|
r = UP5(col[0][RCOMP]);
|
|
} else if (t == 2) {
|
|
b = UP5(col[1][BCOMP]);
|
|
g = UP6(col[1][GCOMP], glsb);
|
|
r = UP5(col[1][RCOMP]);
|
|
} else {
|
|
b = (UP5(col[0][BCOMP]) + UP5(col[1][BCOMP])) / 2;
|
|
g = (UP5(col[0][GCOMP]) + UP6(col[1][GCOMP], glsb)) / 2;
|
|
r = (UP5(col[0][RCOMP]) + UP5(col[1][RCOMP])) / 2;
|
|
}
|
|
rgba[RCOMP] = r;
|
|
rgba[GCOMP] = g;
|
|
rgba[BCOMP] = b;
|
|
rgba[ACOMP] = 255;
|
|
}
|
|
} else {
|
|
/* alpha[0] == 0 */
|
|
uint8_t r, g, b;
|
|
if (t == 0) {
|
|
b = UP5(col[0][BCOMP]);
|
|
g = UP6(col[0][GCOMP], glsb ^ selb);
|
|
r = UP5(col[0][RCOMP]);
|
|
} else if (t == 3) {
|
|
b = UP5(col[1][BCOMP]);
|
|
g = UP6(col[1][GCOMP], glsb);
|
|
r = UP5(col[1][RCOMP]);
|
|
} else {
|
|
b = LERP(3, t, UP5(col[0][BCOMP]), UP5(col[1][BCOMP]));
|
|
g = LERP(3, t, UP6(col[0][GCOMP], glsb ^ selb),
|
|
UP6(col[1][GCOMP], glsb));
|
|
r = LERP(3, t, UP5(col[0][RCOMP]), UP5(col[1][RCOMP]));
|
|
}
|
|
rgba[RCOMP] = r;
|
|
rgba[GCOMP] = g;
|
|
rgba[BCOMP] = b;
|
|
rgba[ACOMP] = 255;
|
|
}
|
|
}
|
|
|
|
|
|
static void
|
|
fxt1_decode_1ALPHA (const uint8_t *code, int32_t t, uint8_t *rgba)
|
|
{
|
|
const uint32_t *cc;
|
|
uint8_t r, g, b, a;
|
|
|
|
cc = (const uint32_t *)code;
|
|
if (CC_SEL(cc, 124) & 1) {
|
|
/* lerp == 1 */
|
|
uint32_t col0[4];
|
|
|
|
if (t & 16) {
|
|
t &= 15;
|
|
t = (cc[1] >> (t * 2)) & 3;
|
|
/* col 2 */
|
|
col0[BCOMP] = (*(const uint32_t *)(code + 11)) >> 6;
|
|
col0[GCOMP] = CC_SEL(cc, 99);
|
|
col0[RCOMP] = CC_SEL(cc, 104);
|
|
col0[ACOMP] = CC_SEL(cc, 119);
|
|
} else {
|
|
t = (cc[0] >> (t * 2)) & 3;
|
|
/* col 0 */
|
|
col0[BCOMP] = CC_SEL(cc, 64);
|
|
col0[GCOMP] = CC_SEL(cc, 69);
|
|
col0[RCOMP] = CC_SEL(cc, 74);
|
|
col0[ACOMP] = CC_SEL(cc, 109);
|
|
}
|
|
|
|
if (t == 0) {
|
|
b = UP5(col0[BCOMP]);
|
|
g = UP5(col0[GCOMP]);
|
|
r = UP5(col0[RCOMP]);
|
|
a = UP5(col0[ACOMP]);
|
|
} else if (t == 3) {
|
|
b = UP5(CC_SEL(cc, 79));
|
|
g = UP5(CC_SEL(cc, 84));
|
|
r = UP5(CC_SEL(cc, 89));
|
|
a = UP5(CC_SEL(cc, 114));
|
|
} else {
|
|
b = LERP(3, t, UP5(col0[BCOMP]), UP5(CC_SEL(cc, 79)));
|
|
g = LERP(3, t, UP5(col0[GCOMP]), UP5(CC_SEL(cc, 84)));
|
|
r = LERP(3, t, UP5(col0[RCOMP]), UP5(CC_SEL(cc, 89)));
|
|
a = LERP(3, t, UP5(col0[ACOMP]), UP5(CC_SEL(cc, 114)));
|
|
}
|
|
} else {
|
|
/* lerp == 0 */
|
|
|
|
if (t & 16) {
|
|
cc++;
|
|
t &= 15;
|
|
}
|
|
t = (cc[0] >> (t * 2)) & 3;
|
|
|
|
if (t == 3) {
|
|
/* zero */
|
|
r = g = b = a = 0;
|
|
} else {
|
|
uint32_t kk;
|
|
cc = (const uint32_t *)code;
|
|
a = UP5(cc[3] >> (t * 5 + 13));
|
|
t *= 15;
|
|
cc = (const uint32_t *)(code + 8 + t / 8);
|
|
kk = cc[0] >> (t & 7);
|
|
b = UP5(kk);
|
|
g = UP5(kk >> 5);
|
|
r = UP5(kk >> 10);
|
|
}
|
|
}
|
|
rgba[RCOMP] = r;
|
|
rgba[GCOMP] = g;
|
|
rgba[BCOMP] = b;
|
|
rgba[ACOMP] = a;
|
|
}
|
|
|
|
|
|
static void
|
|
fxt1_decode_1 (const void *texture, int32_t stride, /* in pixels */
|
|
int32_t i, int32_t j, uint8_t *rgba)
|
|
{
|
|
static void (*decode_1[]) (const uint8_t *, int32_t, uint8_t *) = {
|
|
fxt1_decode_1HI, /* cc-high = "00?" */
|
|
fxt1_decode_1HI, /* cc-high = "00?" */
|
|
fxt1_decode_1CHROMA, /* cc-chroma = "010" */
|
|
fxt1_decode_1ALPHA, /* alpha = "011" */
|
|
fxt1_decode_1MIXED, /* mixed = "1??" */
|
|
fxt1_decode_1MIXED, /* mixed = "1??" */
|
|
fxt1_decode_1MIXED, /* mixed = "1??" */
|
|
fxt1_decode_1MIXED /* mixed = "1??" */
|
|
};
|
|
|
|
const uint8_t *code = (const uint8_t *)texture +
|
|
((j / 4) * (stride / 8) + (i / 8)) * 16;
|
|
int32_t mode = CC_SEL(code, 125);
|
|
int32_t t = i & 7;
|
|
|
|
if (t & 4) {
|
|
t += 12;
|
|
}
|
|
t += (j & 3) * 4;
|
|
|
|
decode_1[mode](code, t, rgba);
|
|
}
|
|
|
|
/*
|
|
* Pixel fetch within a block.
|
|
*/
|
|
|
|
void
|
|
util_format_fxt1_rgb_fetch_rgba_8unorm(uint8_t *restrict dst, const uint8_t *restrict src, unsigned i, unsigned j)
|
|
{
|
|
fxt1_decode_1(src, 0, i, j, dst);
|
|
}
|
|
|
|
void
|
|
util_format_fxt1_rgba_fetch_rgba_8unorm(uint8_t *restrict dst, const uint8_t *restrict src, unsigned i, unsigned j)
|
|
{
|
|
fxt1_decode_1(src, 0, i, j, dst);
|
|
dst[3] = 0xff;
|
|
}
|
|
|
|
void
|
|
util_format_fxt1_rgb_fetch_rgba(void *restrict in_dst, const uint8_t *restrict src, unsigned i, unsigned j)
|
|
{
|
|
float *dst = in_dst;
|
|
uint8_t tmp[4];
|
|
fxt1_decode_1(src, 0, i, j, tmp);
|
|
dst[0] = ubyte_to_float(tmp[0]);
|
|
dst[1] = ubyte_to_float(tmp[1]);
|
|
dst[2] = ubyte_to_float(tmp[2]);
|
|
dst[3] = 1.0;
|
|
}
|
|
|
|
void
|
|
util_format_fxt1_rgba_fetch_rgba(void *restrict in_dst, const uint8_t *restrict src, unsigned i, unsigned j)
|
|
{
|
|
float *dst = in_dst;
|
|
uint8_t tmp[4];
|
|
fxt1_decode_1(src, 0, i, j, tmp);
|
|
dst[0] = ubyte_to_float(tmp[0]);
|
|
dst[1] = ubyte_to_float(tmp[1]);
|
|
dst[2] = ubyte_to_float(tmp[2]);
|
|
dst[3] = ubyte_to_float(tmp[3]);
|
|
}
|
|
|
|
/*
|
|
* Block decompression.
|
|
*/
|
|
|
|
static inline void
|
|
util_format_fxtn_rgb_unpack_rgba_8unorm(uint8_t *restrict dst_row, unsigned dst_stride,
|
|
const uint8_t *restrict src_row, unsigned src_stride,
|
|
unsigned width, unsigned height,
|
|
boolean rgba)
|
|
{
|
|
const unsigned bw = 8, bh = 4, comps = 4;
|
|
unsigned x, y, i, j;
|
|
for (y = 0; y < height; y += bh) {
|
|
const uint8_t *src = src_row;
|
|
for (x = 0; x < width; x += bw) {
|
|
for (j = 0; j < bh; ++j) {
|
|
for (i = 0; i < bw; ++i) {
|
|
uint8_t *dst = dst_row + (y + j) * dst_stride / sizeof(*dst_row) + (x + i) * comps;
|
|
fxt1_decode_1(src, 0, i, j, dst);
|
|
if (!rgba)
|
|
dst[3] = 0xff;
|
|
}
|
|
}
|
|
src += FXT1_BLOCK_SIZE;
|
|
}
|
|
src_row += src_stride;
|
|
}
|
|
}
|
|
|
|
void
|
|
util_format_fxt1_rgb_unpack_rgba_8unorm(uint8_t *restrict dst_row, unsigned dst_stride,
|
|
const uint8_t *restrict src_row, unsigned src_stride,
|
|
unsigned width, unsigned height)
|
|
{
|
|
util_format_fxtn_rgb_unpack_rgba_8unorm(dst_row, dst_stride,
|
|
src_row, src_stride,
|
|
width, height,
|
|
false);
|
|
}
|
|
|
|
void
|
|
util_format_fxt1_rgba_unpack_rgba_8unorm(uint8_t *restrict dst_row, unsigned dst_stride,
|
|
const uint8_t *restrict src_row, unsigned src_stride,
|
|
unsigned width, unsigned height)
|
|
{
|
|
util_format_fxtn_rgb_unpack_rgba_8unorm(dst_row, dst_stride,
|
|
src_row, src_stride,
|
|
width, height,
|
|
true);
|
|
}
|
|
|
|
static inline void
|
|
util_format_fxtn_rgb_unpack_rgba_float(float *dst_row, unsigned dst_stride,
|
|
const uint8_t *restrict src_row, unsigned src_stride,
|
|
unsigned width, unsigned height,
|
|
boolean rgba)
|
|
{
|
|
const unsigned bw = 8, bh = 4, comps = 4;
|
|
unsigned x, y, i, j;
|
|
for (y = 0; y < height; y += 4) {
|
|
const uint8_t *src = src_row;
|
|
for (x = 0; x < width; x += 8) {
|
|
for (j = 0; j < bh; ++j) {
|
|
for (i = 0; i < bw; ++i) {
|
|
float *dst = dst_row + (y + j)*dst_stride/sizeof(*dst_row) + (x + i) * comps;
|
|
uint8_t tmp[4];
|
|
fxt1_decode_1(src, 0, i, j, tmp);
|
|
dst[0] = ubyte_to_float(tmp[0]);
|
|
dst[1] = ubyte_to_float(tmp[1]);
|
|
dst[2] = ubyte_to_float(tmp[2]);
|
|
if (rgba)
|
|
dst[3] = ubyte_to_float(tmp[3]);
|
|
else
|
|
dst[3] = 1.0;
|
|
}
|
|
}
|
|
src += FXT1_BLOCK_SIZE;
|
|
}
|
|
src_row += src_stride;
|
|
}
|
|
}
|
|
|
|
void
|
|
util_format_fxt1_rgb_unpack_rgba_float(void *restrict dst_row, unsigned dst_stride,
|
|
const uint8_t *restrict src_row, unsigned src_stride,
|
|
unsigned width, unsigned height)
|
|
{
|
|
util_format_fxtn_rgb_unpack_rgba_float(dst_row, dst_stride,
|
|
src_row, src_stride,
|
|
width, height,
|
|
false);
|
|
}
|
|
|
|
void
|
|
util_format_fxt1_rgba_unpack_rgba_float(void *restrict dst_row, unsigned dst_stride,
|
|
const uint8_t *restrict src_row, unsigned src_stride,
|
|
unsigned width, unsigned height)
|
|
{
|
|
util_format_fxtn_rgb_unpack_rgba_float(dst_row, dst_stride,
|
|
src_row, src_stride,
|
|
width, height,
|
|
true);
|
|
}
|
|
|
|
/*
|
|
* Block compression.
|
|
*/
|
|
|
|
void
|
|
util_format_fxt1_rgb_pack_rgba_8unorm(uint8_t *restrict dst_row, unsigned dst_stride,
|
|
const uint8_t *restrict src, unsigned src_stride,
|
|
unsigned width, unsigned height)
|
|
{
|
|
/* The encoder for FXT1_RGB wants 24bpp packed rgb, so make a temporary to do that.
|
|
*/
|
|
int temp_stride = width * 3;
|
|
uint8_t *temp = malloc(height * temp_stride);
|
|
if (!temp)
|
|
return;
|
|
|
|
for (int y = 0; y < height; y++) {
|
|
for (int x = 0; x < width; x++) {
|
|
temp[y * temp_stride + x * 3 + 0] = src[x * 4 + 0];
|
|
temp[y * temp_stride + x * 3 + 1] = src[x * 4 + 1];
|
|
temp[y * temp_stride + x * 3 + 2] = src[x * 4 + 2];
|
|
}
|
|
src += src_stride;
|
|
}
|
|
|
|
fxt1_encode(width, height, 3, temp, temp_stride, dst_row, dst_stride);
|
|
|
|
free(temp);
|
|
}
|
|
|
|
void
|
|
util_format_fxt1_rgba_pack_rgba_8unorm(uint8_t *restrict dst_row, unsigned dst_stride,
|
|
const uint8_t *restrict src, unsigned src_stride,
|
|
unsigned width, unsigned height)
|
|
{
|
|
fxt1_encode(width, height, 4, src, src_stride, dst_row, dst_stride);
|
|
}
|
|
|
|
void
|
|
util_format_fxt1_rgb_pack_rgba_float(uint8_t *restrict dst_row, unsigned dst_stride,
|
|
const float *restrict src, unsigned src_stride,
|
|
unsigned width, unsigned height)
|
|
{
|
|
int temp_stride = width * 4;
|
|
uint8_t *temp = malloc(height * temp_stride);
|
|
if (!temp)
|
|
return;
|
|
|
|
util_format_r8g8b8a8_unorm_pack_rgba_float(temp, temp_stride,
|
|
src, src_stride,
|
|
width, height);
|
|
|
|
util_format_fxt1_rgb_pack_rgba_8unorm(dst_row, dst_stride,
|
|
temp, temp_stride,
|
|
width, height);
|
|
|
|
free(temp);
|
|
}
|
|
|
|
void
|
|
util_format_fxt1_rgba_pack_rgba_float(uint8_t *restrict dst_row, unsigned dst_stride,
|
|
const float *restrict src, unsigned src_stride,
|
|
unsigned width, unsigned height)
|
|
{
|
|
int temp_stride = width * 4;
|
|
uint8_t *temp = malloc(height * temp_stride);
|
|
if (!temp)
|
|
return;
|
|
|
|
util_format_r8g8b8a8_unorm_pack_rgba_float(temp, temp_stride,
|
|
src, src_stride,
|
|
width, height);
|
|
|
|
util_format_fxt1_rgba_pack_rgba_8unorm(dst_row, dst_stride,
|
|
temp, temp_stride,
|
|
width, height);
|
|
|
|
free(temp);
|
|
}
|