/* 
** Copyright 2004, Marcus Overhagen. All rights reserved.
** Distributed under the terms of the MIT license.
**
** This file generates the lookup tables used by image
** format conversion from planar YCbCr420 format to
** linear RGB32 format in file gfx_conv_c_lookup.cpp
*/

#include <strings.h>
#include <stdio.h>

double precision = 32768.0;
int shift = 15;

/*
        R = 1.164(Y - 16) + 1.596(Cr - 128)
        G = 1.164(Y - 16) - 0.813(Cr - 128) - 0.392(Cb - 128)
        B = 1.164(Y - 16) + 2.017(Cb - 128)
*/

int main()
{
        printf("const int32 Cb_Gtab[256] = {\n\t");
        for (int i = 0; i < 256; i++) {
                int     Cr_R = (- 128 + i) * 1.596 * precision + 0.5;
                int     Cr_G = (- 128 + i) * -0.813 * precision + 0.5;
                int     Cb_G = (- 128 + i) * -0.392 * precision + 0.5;
                int     Cb_B = (- 128 + i) * 2.017 * precision + 0.5;
                printf("%d, ", Cb_G);   
                if ((i % 10) == 9)
                        printf("\n\t");
         }
        printf("\n};\n\n");

        printf("const int32 Cb_Btab[256] = {\n\t");
        for (int i = 0; i < 256; i++) {
                int     Cr_R = (- 128 + i) * 1.596 * precision + 0.5;
                int     Cr_G = (- 128 + i) * -0.813 * precision + 0.5;
                int     Cb_G = (- 128 + i) * -0.392 * precision + 0.5;
                int     Cb_B = (- 128 + i) * 2.017 * precision + 0.5;
                printf("%d, ", Cb_B);   
                if ((i % 10) == 9)
                        printf("\n\t");
         }
        printf("\n};\n\n");

        printf("const int32 Cr_Rtab[256] = {\n\t");
        for (int i = 0; i < 256; i++) {
                int     Cr_R = (- 128 + i) * 1.596 * precision + 0.5;
                int     Cr_G = (- 128 + i) * -0.813 * precision + 0.5;
                int     Cb_G = (- 128 + i) * -0.392 * precision + 0.5;
                int     Cb_B = (- 128 + i) * 2.017 * precision + 0.5;
                printf("%d, ", Cr_R);   
                if ((i % 10) == 9)
                        printf("\n\t");
         }
        printf("\n};\n\n");

        printf("const int32 Cr_Gtab[256] = {\n\t");
        for (int i = 0; i < 256; i++) {
                int     Cr_R = (- 128 + i) * 1.596 * precision + 0.5;
                int     Cr_G = (- 128 + i) * -0.813 * precision + 0.5;
                int     Cb_G = (- 128 + i) * -0.392 * precision + 0.5;
                int     Cb_B = (- 128 + i) * 2.017 * precision + 0.5;
                printf("%d, ", Cr_G);   
                if ((i % 10) == 9)
                        printf("\n\t");
         }
        printf("\n};\n\n");

        printf("const int32 Ytab[256] = {\n\t");
        for (int i = 0; i < 256; i++) {
                int Y = (i - 16) * 1.164 * precision + 0.5;
                printf("%d, ", Y);      
                if ((i % 10) == 9)
                        printf("\n\t");
         }
        printf("\n};\n\n");

        int Ymin = (0 - 16) * 1.164 * precision + 0.5;
        int Ymax = (255 - 16) * 1.164 * precision + 0.5;

        int Cr_Rmin = (- 128 + 0) * 1.596 * precision + 0.5;
        int Cr_Rmax = (- 128 + 255) * 1.596 * precision + 0.5;
        int Cr_Gmin = (- 128 + 255) * -0.813 * precision + 0.5;
        int Cr_Gmax = (- 128 + 0) * -0.813 * precision + 0.5;
        int Cb_Gmin = (- 128 + 255) * -0.392 * precision + 0.5;
        int Cb_Gmax = (- 128 + 0) * -0.392 * precision + 0.5;
        int Cb_Bmin = (- 128 + 0) * 2.017 * precision + 0.5;
        int Cb_Bmax = (- 128 + 255) * 2.017 * precision + 0.5;

        printf("Cr_Rmin %d\n", Cr_Rmin);
        printf("Cr_Rmax %d\n", Cr_Rmax);
        printf("Cr_Gmin %d\n", Cr_Gmin);
        printf("Cr_Gmax %d\n", Cr_Gmax);
        printf("Cb_Gmin %d\n", Cb_Gmin);
        printf("Cb_Gmax %d\n", Cb_Gmax);
        printf("Cb_Bmin %d\n", Cb_Bmin);
        printf("Cb_Bmax %d\n", Cb_Bmax);
        
        int Rmax = (Ymax + Cr_Rmax) >> shift;
        int Rmin = (Ymin + Cr_Rmin) >> shift;
        int Gmax = (Ymax + Cr_Gmax + Cb_Gmax) >> shift;
        int Gmin = (Ymin + Cr_Gmin + Cb_Gmin) >> shift;
        int Bmax = (Ymax + Cb_Bmax) >> shift;
        int Bmin = (Ymin + Cb_Bmin) >> shift;
        
        printf("Rmax %d\n", Rmax);
        printf("Rmin %d\n", Rmin);
        printf("Gmax %d\n", Gmax);
        printf("Gmin %d\n", Gmin);
        printf("Bmax %d\n", Bmax);
        printf("Bmin %d\n", Bmin);
        
        int num_Rneg = (Rmin < 0) ? -Rmin : 0;
        int num_Rpos = (Rmax > 255) ? Rmax - 255 : 0;
        int num_Rent = num_Rneg + 256 + num_Rpos;

        printf("num_Rneg %d\n", num_Rneg);
        printf("num_Rpos %d\n", num_Rpos);
        printf("num_Rent %d\n", num_Rent);
        
        printf("const uint32 Rsat[%d] = {\n\t", num_Rent);
        for (int i = 0; i < num_Rneg; i++) {
                printf("0x000000, ");   
                if ((i % 10) == 9)
                        printf("\n\t");
        }
        printf("\n\t");
        for (int i = 0; i < 256; i++) {
                printf("0x%06x, ", i << 16);    
                if ((i % 10) == 9)
                        printf("\n\t");
        }
        printf("\n\t");
        for (int i = 0; i < num_Rpos; i++) {
                printf("0x%06x, ", 255 << 16);  
                if ((i % 10) == 9)
                        printf("\n\t");
        }
        printf("\n};\n");
        printf("const uint32 *Rtab = &Rsat[%d];\n\n", num_Rneg);
        
        
        
        int num_Gneg = (Gmin < 0) ? -Gmin : 0;
        int num_Gpos = (Gmax > 255) ? Gmax - 255 : 0;
        int num_Gent = num_Gneg + 256 + num_Gpos;

        printf("num_Gneg %d\n", num_Gneg);
        printf("num_Gpos %d\n", num_Gpos);
        printf("num_Gent %d\n", num_Gent);
        
        printf("const uint32 Gsat[%d] = {\n\t", num_Gent);
        for (int i = 0; i < num_Gneg; i++) {
                printf("0x000000, ");   
                if ((i % 10) == 9)
                        printf("\n\t");
        }
        printf("\n\t");
        for (int i = 0; i < 256; i++) {
                printf("0x%06x, ", i << 8);     
                if ((i % 10) == 9)
                        printf("\n\t");
        }
        printf("\n\t");
        for (int i = 0; i < num_Gpos; i++) {
                printf("0x%06x, ", 255 << 8);   
                if ((i % 10) == 9)
                        printf("\n\t");
        }
        printf("\n};\n");
        printf("const uint32 *Gtab = &Gsat[%d];\n\n", num_Gneg);
        
        
        int num_Bneg = (Bmin < 0) ? -Bmin : 0;
        int num_Bpos = (Bmax > 255) ? Bmax - 255 : 0;
        int num_Bent = num_Bneg + 256 + num_Bpos;

        printf("num_Bneg %d\n", num_Bneg);
        printf("num_Bpos %d\n", num_Bpos);
        printf("num_Bent %d\n", num_Bent);
        
        printf("const uint32 Bsat[%d] = {\n\t", num_Bent);
        for (int i = 0; i < num_Bneg; i++) {
                printf("0x000000, ");   
                if ((i % 10) == 9)
                        printf("\n\t");
        }
        printf("\n\t");
        for (int i = 0; i < 256; i++) {
                printf("0x%06x, ", i << 0);     
                if ((i % 10) == 9)
                        printf("\n\t");
        }
        printf("\n\t");
        for (int i = 0; i < num_Bpos; i++) {
                printf("0x%06x, ", 255 << 0);   
                if ((i % 10) == 9)
                        printf("\n\t");
        }
        printf("\n};\n");
        printf("const uint32 *Btab = &Bsat[%d];\n\n", num_Bneg);
        
        return 0;
}