Fix line endlings

1 files changed, 769 insertions, 769 deletions

diff --git a/indra/llmessage/patch_dct.cpp b/indra/llmessage/patch_dct.cpp
index 10b3e5dd77..728fe84537 100644
--- a/indra/llmessage/patch_dct.cpp
+++ b/indra/llmessage/patch_dct.cpp
@@ -1,769 +1,769 @@
-/**

- * @file patch_dct.cpp

- * @brief DCT patch.

- *

- * $LicenseInfo:firstyear=2000&license=viewerlgpl$

- * Second Life Viewer Source Code

- * Copyright (C) 2010, Linden Research, Inc.

- *

- * This library is free software; you can redistribute it and/or

- * modify it under the terms of the GNU Lesser General Public

- * License as published by the Free Software Foundation;

- * version 2.1 of the License only.

- *

- * This library is distributed in the hope that it will be useful,

- * but WITHOUT ANY WARRANTY; without even the implied warranty of

- * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU

- * Lesser General Public License for more details.

- *

- * You should have received a copy of the GNU Lesser General Public

- * License along with this library; if not, write to the Free Software

- * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA

- *

- * Linden Research, Inc., 945 Battery Street, San Francisco, CA  94111  USA

- * $/LicenseInfo$

- */

-

-#include "linden_common.h"

-

-#include "llmath.h"

-//#include "vmath.h"

-#include "v3math.h"

-#include "patch_dct.h"

-

-typedef struct s_patch_compress_global_data

-{

-    S32 patch_size;

-    S32 patch_stride;

-    U32 charptr;

-    S32 layer_type;

-} PCGD;

-

-PCGD    gPatchCompressGlobalData;

-

-void reset_patch_compressor(void)

-{

-    PCGD *pcp = &gPatchCompressGlobalData;

-

-    pcp->charptr = 0;

-}

-

-S32 gCurrentSize = 0;

-

-F32 gPatchQuantizeTable[LARGE_PATCH_SIZE*LARGE_PATCH_SIZE];

-

-void build_patch_quantize_table(S32 size)

-{

-    S32 i, j;

-    for (j = 0; j < size; j++)

-    {

-        for (i = 0; i < size; i++)

-        {

-            gPatchQuantizeTable[j*size + i] = 1.f/(1.f + 2.f*(i+j));

-        }

-    }

-}

-

-F32 gPatchCosines[LARGE_PATCH_SIZE*LARGE_PATCH_SIZE];

-

-void setup_patch_cosines(S32 size)

-{

-    S32 n, u;

-    F32 oosob = F_PI*0.5f/size;

-

-    for (u = 0; u < size; u++)

-    {

-        for (n = 0; n < size; n++)

-        {

-            gPatchCosines[u*size+n] = cosf((2.f*n+1.f)*u*oosob);

-        }

-    }

-}

-

-S32 gCopyMatrix[LARGE_PATCH_SIZE*LARGE_PATCH_SIZE];

-

-void build_copy_matrix(S32 size)

-{

-    S32 i, j, count;

-    bool    b_diag = false;

-    bool    b_right = true;

-

-    i = 0;

-    j = 0;

-    count = 0;

-

-    while (  (i < size)

-           &&(j < size))

-    {

-        gCopyMatrix[j*size + i] = count;

-

-        count++;

-

-        if (!b_diag)

-        {

-            if (b_right)

-            {

-                if (i < size - 1)

-                    i++;

-                else

-                    j++;

-                b_right = false;

-                b_diag = true;

-            }

-            else

-            {

-                if (j < size - 1)

-                    j++;

-                else

-                    i++;

-                b_right = true;

-                b_diag = true;

-            }

-        }

-        else

-        {

-            if (b_right)

-            {

-                i++;

-                j--;

-                if (  (i == size - 1)

-                    ||(j == 0))

-                {

-                    b_diag = false;

-                }

-            }

-            else

-            {

-                i--;

-                j++;

-                if (  (i == 0)

-                    ||(j == size - 1))

-                {

-                    b_diag = false;

-                }

-            }

-        }

-    }

-}

-

-

-void init_patch_compressor(S32 patch_size, S32 patch_stride, S32 layer_type)

-{

-    PCGD *pcp = &gPatchCompressGlobalData;

-

-    pcp->charptr = 0;

-

-    pcp->patch_size = patch_size;

-    pcp->patch_stride = patch_stride;

-    pcp->layer_type = layer_type;

-

-    if (patch_size != gCurrentSize)

-    {

-        gCurrentSize = patch_size;

-        build_patch_quantize_table(patch_size);

-        setup_patch_cosines(patch_size);

-        build_copy_matrix(patch_size);

-    }

-}

-

-void prescan_patch(F32 *patch, LLPatchHeader *php, F32 &zmax, F32 &zmin)

-{

-    S32     i, j;

-    PCGD    *pcp = &gPatchCompressGlobalData;

-    S32     stride = pcp->patch_stride;

-    S32     size = pcp->patch_size;

-    S32     jstride;

-

-    zmax = -99999999.f;

-    zmin = 99999999.f;

-

-    for (j = 0; j < size; j++)

-    {

-        jstride = j*stride;

-        for (i = 0; i < size; i++)

-        {

-            if (*(patch + jstride + i) > zmax)

-            {

-                zmax = *(patch + jstride + i);

-            }

-            if (*(patch + jstride + i) < zmin)

-            {

-                zmin = *(patch + jstride + i);

-            }

-        }

-    }

-

-    php->dc_offset = zmin;

-    php->range = (U16) ((zmax - zmin) + 1.f);

-}

-

-void dct_line(F32 *linein, F32 *lineout, S32 line)

-{

-    S32 u;

-    F32 total;

-    F32 *pcp = gPatchCosines;

-    S32 line_size = line*NORMAL_PATCH_SIZE;

-

-#ifdef _PATCH_SIZE_16_AND_32_ONLY

-    F32 *tlinein, *tpcp;

-

-    tlinein = linein + line_size;

-

-    total = *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein);

-

-    *(lineout + line_size) = OO_SQRT2*total;

-

-    for (u = 1; u < NORMAL_PATCH_SIZE; u++)

-    {

-        tlinein = linein + line_size;

-        tpcp = pcp + (u<<4);

-

-        total = *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein)*(*tpcp);

-

-        *(lineout + line_size + u) = total;

-    }

-#else

-    S32 n;

-    S32 size = gPatchCompressGlobalData.patch_size;

-    total = 0.f;

-    for (n = 0; n < size; n++)

-    {

-        total += linein[line_size + n];

-    }

-    lineout[line_size] = OO_SQRT2*total;

-

-    for (u = 1; u < size; u++)

-    {

-        total = 0.f;

-        for (n = 0; n < size; n++)

-        {

-            total += linein[line_size + n]*pcp[u*size+n];

-        }

-        lineout[line_size + u] = total;

-    }

-#endif

-}

-

-void dct_line_large(F32 *linein, F32 *lineout, S32 line)

-{

-    S32 u;

-    F32 total;

-    F32 *pcp = gPatchCosines;

-    S32 line_size = line*LARGE_PATCH_SIZE;

-

-    F32 *tlinein, *tpcp;

-

-    tlinein = linein + line_size;

-

-    total = *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein++);

-    total += *(tlinein);

-

-    *(lineout + line_size) = OO_SQRT2*total;

-

-    for (u = 1; u < LARGE_PATCH_SIZE; u++)

-    {

-        tlinein = linein + line_size;

-        tpcp = pcp + (u<<5);

-

-        total = *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein++)*(*(tpcp++));

-        total += *(tlinein)*(*tpcp);

-

-        *(lineout + line_size + u) = total;

-    }

-}

-

-inline void dct_column(F32 *linein, S32 *lineout, S32 column)

-{

-    S32 u;

-    F32 total;

-    F32 oosob = 2.f/16.f;

-    F32 *pcp = gPatchCosines;

-    S32 *copy_matrix = gCopyMatrix;

-    F32 *qt = gPatchQuantizeTable;

-

-#ifdef _PATCH_SIZE_16_AND_32_ONLY

-    F32 *tlinein, *tpcp;

-    S32 sizeu;

-

-    tlinein = linein + column;

-

-    total = *(tlinein);

-    total += *(tlinein += NORMAL_PATCH_SIZE);

-    total += *(tlinein += NORMAL_PATCH_SIZE);

-    total += *(tlinein += NORMAL_PATCH_SIZE);

-

-    total += *(tlinein += NORMAL_PATCH_SIZE);

-    total += *(tlinein += NORMAL_PATCH_SIZE);

-    total += *(tlinein += NORMAL_PATCH_SIZE);

-    total += *(tlinein += NORMAL_PATCH_SIZE);

-

-    total += *(tlinein += NORMAL_PATCH_SIZE);

-    total += *(tlinein += NORMAL_PATCH_SIZE);

-    total += *(tlinein += NORMAL_PATCH_SIZE);

-    total += *(tlinein += NORMAL_PATCH_SIZE);

-

-    total += *(tlinein += NORMAL_PATCH_SIZE);

-    total += *(tlinein += NORMAL_PATCH_SIZE);

-    total += *(tlinein += NORMAL_PATCH_SIZE);

-    total += *(tlinein += NORMAL_PATCH_SIZE);

-

-    *(lineout + *(copy_matrix + column)) = (S32)(OO_SQRT2*total*oosob*(*(qt + column)));

-

-    for (u = 1; u < NORMAL_PATCH_SIZE; u++)

-    {

-        tlinein = linein + column;

-        tpcp = pcp + (u<<4);

-

-        total = *(tlinein)*(*(tpcp++));

-        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));

-

-        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));

-

-        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));

-

-        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp));

-

-        sizeu = NORMAL_PATCH_SIZE*u + column;

-

-        *(lineout + *(copy_matrix + sizeu)) = (S32)(total*oosob*(*(qt+sizeu)));

-    }

-#else

-    S32 size = gPatchCompressGlobalData.patch_size;

-    F32 oosob = 2.f/size;

-    S32 n;

-    total = 0.f;

-    for (n = 0; n < size; n++)

-    {

-        total += linein[size*n + column];

-    }

-    lineout[copy_matrix[column]] = OO_SQRT2*total*oosob*qt[column];

-

-    for (u = 1; u < size; u++)

-    {

-        total = 0.f;

-        for (n = 0; n < size; n++)

-        {

-            total += linein[size*n + column]*pcp[u*size+n];

-        }

-        lineout[copy_matrix[size*u + column]] = total*oosob*qt[size*u + column];

-    }

-#endif

-}

-

-inline void dct_column_large(F32 *linein, S32 *lineout, S32 column)

-{

-    S32 u;

-    F32 total;

-    F32 oosob = 2.f/32.f;

-    F32 *pcp = gPatchCosines;

-    S32 *copy_matrix = gCopyMatrix;

-    F32 *qt = gPatchQuantizeTable;

-

-    F32 *tlinein, *tpcp;

-    S32 sizeu;

-

-    tlinein = linein + column;

-

-    total = *(tlinein);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-    total += *(tlinein += LARGE_PATCH_SIZE);

-

-    *(lineout + *(copy_matrix + column)) = (S32)(OO_SQRT2*total*oosob*(*(qt + column)));

-

-    for (u = 1; u < LARGE_PATCH_SIZE; u++)

-    {

-        tlinein = linein + column;

-        tpcp = pcp + (u<<5);

-

-        total = *(tlinein)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));

-        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp));

-

-        sizeu = LARGE_PATCH_SIZE*u + column;

-

-        *(lineout + *(copy_matrix + sizeu)) = (S32)(total*oosob*(*(qt+sizeu)));

-    }

-}

-

-inline void dct_patch(F32 *block, S32 *cpatch)

-{

-    F32 temp[NORMAL_PATCH_SIZE*NORMAL_PATCH_SIZE];

-

-#ifdef _PATCH_SIZE_16_AND_32_ONLY

-    dct_line(block, temp, 0);

-    dct_line(block, temp, 1);

-    dct_line(block, temp, 2);

-    dct_line(block, temp, 3);

-

-    dct_line(block, temp, 4);

-    dct_line(block, temp, 5);

-    dct_line(block, temp, 6);

-    dct_line(block, temp, 7);

-

-    dct_line(block, temp, 8);

-    dct_line(block, temp, 9);

-    dct_line(block, temp, 10);

-    dct_line(block, temp, 11);

-

-    dct_line(block, temp, 12);

-    dct_line(block, temp, 13);

-    dct_line(block, temp, 14);

-    dct_line(block, temp, 15);

-

-    dct_column(temp, cpatch, 0);

-    dct_column(temp, cpatch, 1);

-    dct_column(temp, cpatch, 2);

-    dct_column(temp, cpatch, 3);

-

-    dct_column(temp, cpatch, 4);

-    dct_column(temp, cpatch, 5);

-    dct_column(temp, cpatch, 6);

-    dct_column(temp, cpatch, 7);

-

-    dct_column(temp, cpatch, 8);

-    dct_column(temp, cpatch, 9);

-    dct_column(temp, cpatch, 10);

-    dct_column(temp, cpatch, 11);

-

-    dct_column(temp, cpatch, 12);

-    dct_column(temp, cpatch, 13);

-    dct_column(temp, cpatch, 14);

-    dct_column(temp, cpatch, 15);

-#else

-    S32 i;

-    S32 size = gPatchCompressGlobalData.patch_size;

-    for (i = 0; i < size; i++)

-    {

-        dct_line(block, temp, i);

-    }

-    for (i = 0; i < size; i++)

-    {

-        dct_column(temp, cpatch, i);

-    }

-#endif

-}

-

-inline void dct_patch_large(F32 *block, S32 *cpatch)

-{

-    F32 temp[LARGE_PATCH_SIZE*LARGE_PATCH_SIZE];

-

-    dct_line_large(block, temp, 0);

-    dct_line_large(block, temp, 1);

-    dct_line_large(block, temp, 2);

-    dct_line_large(block, temp, 3);

-

-    dct_line_large(block, temp, 4);

-    dct_line_large(block, temp, 5);

-    dct_line_large(block, temp, 6);

-    dct_line_large(block, temp, 7);

-

-    dct_line_large(block, temp, 8);

-    dct_line_large(block, temp, 9);

-    dct_line_large(block, temp, 10);

-    dct_line_large(block, temp, 11);

-

-    dct_line_large(block, temp, 12);

-    dct_line_large(block, temp, 13);

-    dct_line_large(block, temp, 14);

-    dct_line_large(block, temp, 15);

-

-    dct_line_large(block, temp, 16);

-    dct_line_large(block, temp, 17);

-    dct_line_large(block, temp, 18);

-    dct_line_large(block, temp, 19);

-

-    dct_line_large(block, temp, 20);

-    dct_line_large(block, temp, 21);

-    dct_line_large(block, temp, 22);

-    dct_line_large(block, temp, 23);

-

-    dct_line_large(block, temp, 24);

-    dct_line_large(block, temp, 25);

-    dct_line_large(block, temp, 26);

-    dct_line_large(block, temp, 27);

-

-    dct_line_large(block, temp, 28);

-    dct_line_large(block, temp, 29);

-    dct_line_large(block, temp, 30);

-    dct_line_large(block, temp, 31);

-

-    dct_column_large(temp, cpatch, 0);

-    dct_column_large(temp, cpatch, 1);

-    dct_column_large(temp, cpatch, 2);

-    dct_column_large(temp, cpatch, 3);

-

-    dct_column_large(temp, cpatch, 4);

-    dct_column_large(temp, cpatch, 5);

-    dct_column_large(temp, cpatch, 6);

-    dct_column_large(temp, cpatch, 7);

-

-    dct_column_large(temp, cpatch, 8);

-    dct_column_large(temp, cpatch, 9);

-    dct_column_large(temp, cpatch, 10);

-    dct_column_large(temp, cpatch, 11);

-

-    dct_column_large(temp, cpatch, 12);

-    dct_column_large(temp, cpatch, 13);

-    dct_column_large(temp, cpatch, 14);

-    dct_column_large(temp, cpatch, 15);

-

-    dct_column_large(temp, cpatch, 16);

-    dct_column_large(temp, cpatch, 17);

-    dct_column_large(temp, cpatch, 18);

-    dct_column_large(temp, cpatch, 19);

-

-    dct_column_large(temp, cpatch, 20);

-    dct_column_large(temp, cpatch, 21);

-    dct_column_large(temp, cpatch, 22);

-    dct_column_large(temp, cpatch, 23);

-

-    dct_column_large(temp, cpatch, 24);

-    dct_column_large(temp, cpatch, 25);

-    dct_column_large(temp, cpatch, 26);

-    dct_column_large(temp, cpatch, 27);

-

-    dct_column_large(temp, cpatch, 28);

-    dct_column_large(temp, cpatch, 29);

-    dct_column_large(temp, cpatch, 30);

-    dct_column_large(temp, cpatch, 31);

-}

-

-void compress_patch(F32 *patch, S32 *cpatch, LLPatchHeader *php, S32 prequant)

-{

-    S32     i, j;

-    PCGD    *pcp = &gPatchCompressGlobalData;

-    S32     stride = pcp->patch_stride;

-    S32     size = pcp->patch_size;

-    F32     block[LARGE_PATCH_SIZE*LARGE_PATCH_SIZE], *tblock;

-    F32     *tpatch;

-

-    S32     wordsize = prequant;

-    F32     oozrange = 1.f/php->range;

-

-    F32     dc = php->dc_offset;

-

-    S32     range = (1<<prequant);

-    F32     premult = oozrange*range;

-//  F32     sub = (F32)(1<<(prequant - 1));

-    F32     sub = (F32)(1<<(prequant - 1)) + dc*premult;

-

-    php->quant_wbits = wordsize - 2;

-    php->quant_wbits |= (prequant - 2)<<4;

-

-    for (j = 0; j < size; j++)

-    {

-        tblock = block + j*size;

-        tpatch = patch + j*stride;

-        for (i = 0; i < size; i++)

-        {

-//          block[j*size + i] = (patch[j*stride + i] - dc)*premult - sub;

-            *(tblock++) = *(tpatch++)*premult - sub;

-        }

-    }

-

-    if (size == 16)

-        dct_patch(block, cpatch);

-    else

-        dct_patch_large(block, cpatch);

-}

-

-void get_patch_group_header(LLGroupHeader *gopp)

-{

-    PCGD    *pcp = &gPatchCompressGlobalData;

-    gopp->stride = pcp->patch_stride;

-    gopp->patch_size = pcp->patch_size;

-    gopp->layer_type = pcp->layer_type;

-}

+/**
+ * @file patch_dct.cpp
+ * @brief DCT patch.
+ *
+ * $LicenseInfo:firstyear=2000&license=viewerlgpl$
+ * Second Life Viewer Source Code
+ * Copyright (C) 2010, Linden Research, Inc.
+ *
+ * This library is free software; you can redistribute it and/or
+ * modify it under the terms of the GNU Lesser General Public
+ * License as published by the Free Software Foundation;
+ * version 2.1 of the License only.
+ *
+ * This library is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+ * Lesser General Public License for more details.
+ *
+ * You should have received a copy of the GNU Lesser General Public
+ * License along with this library; if not, write to the Free Software
+ * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
+ *
+ * Linden Research, Inc., 945 Battery Street, San Francisco, CA  94111  USA
+ * $/LicenseInfo$
+ */
+
+#include "linden_common.h"
+
+#include "llmath.h"
+//#include "vmath.h"
+#include "v3math.h"
+#include "patch_dct.h"
+
+typedef struct s_patch_compress_global_data
+{
+    S32 patch_size;
+    S32 patch_stride;
+    U32 charptr;
+    S32 layer_type;
+} PCGD;
+
+PCGD    gPatchCompressGlobalData;
+
+void reset_patch_compressor(void)
+{
+    PCGD *pcp = &gPatchCompressGlobalData;
+
+    pcp->charptr = 0;
+}
+
+S32 gCurrentSize = 0;
+
+F32 gPatchQuantizeTable[LARGE_PATCH_SIZE*LARGE_PATCH_SIZE];
+
+void build_patch_quantize_table(S32 size)
+{
+    S32 i, j;
+    for (j = 0; j < size; j++)
+    {
+        for (i = 0; i < size; i++)
+        {
+            gPatchQuantizeTable[j*size + i] = 1.f/(1.f + 2.f*(i+j));
+        }
+    }
+}
+
+F32 gPatchCosines[LARGE_PATCH_SIZE*LARGE_PATCH_SIZE];
+
+void setup_patch_cosines(S32 size)
+{
+    S32 n, u;
+    F32 oosob = F_PI*0.5f/size;
+
+    for (u = 0; u < size; u++)
+    {
+        for (n = 0; n < size; n++)
+        {
+            gPatchCosines[u*size+n] = cosf((2.f*n+1.f)*u*oosob);
+        }
+    }
+}
+
+S32 gCopyMatrix[LARGE_PATCH_SIZE*LARGE_PATCH_SIZE];
+
+void build_copy_matrix(S32 size)
+{
+    S32 i, j, count;
+    bool    b_diag = false;
+    bool    b_right = true;
+
+    i = 0;
+    j = 0;
+    count = 0;
+
+    while (  (i < size)
+           &&(j < size))
+    {
+        gCopyMatrix[j*size + i] = count;
+
+        count++;
+
+        if (!b_diag)
+        {
+            if (b_right)
+            {
+                if (i < size - 1)
+                    i++;
+                else
+                    j++;
+                b_right = false;
+                b_diag = true;
+            }
+            else
+            {
+                if (j < size - 1)
+                    j++;
+                else
+                    i++;
+                b_right = true;
+                b_diag = true;
+            }
+        }
+        else
+        {
+            if (b_right)
+            {
+                i++;
+                j--;
+                if (  (i == size - 1)
+                    ||(j == 0))
+                {
+                    b_diag = false;
+                }
+            }
+            else
+            {
+                i--;
+                j++;
+                if (  (i == 0)
+                    ||(j == size - 1))
+                {
+                    b_diag = false;
+                }
+            }
+        }
+    }
+}
+
+
+void init_patch_compressor(S32 patch_size, S32 patch_stride, S32 layer_type)
+{
+    PCGD *pcp = &gPatchCompressGlobalData;
+
+    pcp->charptr = 0;
+
+    pcp->patch_size = patch_size;
+    pcp->patch_stride = patch_stride;
+    pcp->layer_type = layer_type;
+
+    if (patch_size != gCurrentSize)
+    {
+        gCurrentSize = patch_size;
+        build_patch_quantize_table(patch_size);
+        setup_patch_cosines(patch_size);
+        build_copy_matrix(patch_size);
+    }
+}
+
+void prescan_patch(F32 *patch, LLPatchHeader *php, F32 &zmax, F32 &zmin)
+{
+    S32     i, j;
+    PCGD    *pcp = &gPatchCompressGlobalData;
+    S32     stride = pcp->patch_stride;
+    S32     size = pcp->patch_size;
+    S32     jstride;
+
+    zmax = -99999999.f;
+    zmin = 99999999.f;
+
+    for (j = 0; j < size; j++)
+    {
+        jstride = j*stride;
+        for (i = 0; i < size; i++)
+        {
+            if (*(patch + jstride + i) > zmax)
+            {
+                zmax = *(patch + jstride + i);
+            }
+            if (*(patch + jstride + i) < zmin)
+            {
+                zmin = *(patch + jstride + i);
+            }
+        }
+    }
+
+    php->dc_offset = zmin;
+    php->range = (U16) ((zmax - zmin) + 1.f);
+}
+
+void dct_line(F32 *linein, F32 *lineout, S32 line)
+{
+    S32 u;
+    F32 total;
+    F32 *pcp = gPatchCosines;
+    S32 line_size = line*NORMAL_PATCH_SIZE;
+
+#ifdef _PATCH_SIZE_16_AND_32_ONLY
+    F32 *tlinein, *tpcp;
+
+    tlinein = linein + line_size;
+
+    total = *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein);
+
+    *(lineout + line_size) = OO_SQRT2*total;
+
+    for (u = 1; u < NORMAL_PATCH_SIZE; u++)
+    {
+        tlinein = linein + line_size;
+        tpcp = pcp + (u<<4);
+
+        total = *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein)*(*tpcp);
+
+        *(lineout + line_size + u) = total;
+    }
+#else
+    S32 n;
+    S32 size = gPatchCompressGlobalData.patch_size;
+    total = 0.f;
+    for (n = 0; n < size; n++)
+    {
+        total += linein[line_size + n];
+    }
+    lineout[line_size] = OO_SQRT2*total;
+
+    for (u = 1; u < size; u++)
+    {
+        total = 0.f;
+        for (n = 0; n < size; n++)
+        {
+            total += linein[line_size + n]*pcp[u*size+n];
+        }
+        lineout[line_size + u] = total;
+    }
+#endif
+}
+
+void dct_line_large(F32 *linein, F32 *lineout, S32 line)
+{
+    S32 u;
+    F32 total;
+    F32 *pcp = gPatchCosines;
+    S32 line_size = line*LARGE_PATCH_SIZE;
+
+    F32 *tlinein, *tpcp;
+
+    tlinein = linein + line_size;
+
+    total = *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein++);
+    total += *(tlinein);
+
+    *(lineout + line_size) = OO_SQRT2*total;
+
+    for (u = 1; u < LARGE_PATCH_SIZE; u++)
+    {
+        tlinein = linein + line_size;
+        tpcp = pcp + (u<<5);
+
+        total = *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein++)*(*(tpcp++));
+        total += *(tlinein)*(*tpcp);
+
+        *(lineout + line_size + u) = total;
+    }
+}
+
+inline void dct_column(F32 *linein, S32 *lineout, S32 column)
+{
+    S32 u;
+    F32 total;
+    F32 oosob = 2.f/16.f;
+    F32 *pcp = gPatchCosines;
+    S32 *copy_matrix = gCopyMatrix;
+    F32 *qt = gPatchQuantizeTable;
+
+#ifdef _PATCH_SIZE_16_AND_32_ONLY
+    F32 *tlinein, *tpcp;
+    S32 sizeu;
+
+    tlinein = linein + column;
+
+    total = *(tlinein);
+    total += *(tlinein += NORMAL_PATCH_SIZE);
+    total += *(tlinein += NORMAL_PATCH_SIZE);
+    total += *(tlinein += NORMAL_PATCH_SIZE);
+
+    total += *(tlinein += NORMAL_PATCH_SIZE);
+    total += *(tlinein += NORMAL_PATCH_SIZE);
+    total += *(tlinein += NORMAL_PATCH_SIZE);
+    total += *(tlinein += NORMAL_PATCH_SIZE);
+
+    total += *(tlinein += NORMAL_PATCH_SIZE);
+    total += *(tlinein += NORMAL_PATCH_SIZE);
+    total += *(tlinein += NORMAL_PATCH_SIZE);
+    total += *(tlinein += NORMAL_PATCH_SIZE);
+
+    total += *(tlinein += NORMAL_PATCH_SIZE);
+    total += *(tlinein += NORMAL_PATCH_SIZE);
+    total += *(tlinein += NORMAL_PATCH_SIZE);
+    total += *(tlinein += NORMAL_PATCH_SIZE);
+
+    *(lineout + *(copy_matrix + column)) = (S32)(OO_SQRT2*total*oosob*(*(qt + column)));
+
+    for (u = 1; u < NORMAL_PATCH_SIZE; u++)
+    {
+        tlinein = linein + column;
+        tpcp = pcp + (u<<4);
+
+        total = *(tlinein)*(*(tpcp++));
+        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));
+
+        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));
+
+        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));
+
+        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += NORMAL_PATCH_SIZE)*(*(tpcp));
+
+        sizeu = NORMAL_PATCH_SIZE*u + column;
+
+        *(lineout + *(copy_matrix + sizeu)) = (S32)(total*oosob*(*(qt+sizeu)));
+    }
+#else
+    S32 size = gPatchCompressGlobalData.patch_size;
+    F32 oosob = 2.f/size;
+    S32 n;
+    total = 0.f;
+    for (n = 0; n < size; n++)
+    {
+        total += linein[size*n + column];
+    }
+    lineout[copy_matrix[column]] = OO_SQRT2*total*oosob*qt[column];
+
+    for (u = 1; u < size; u++)
+    {
+        total = 0.f;
+        for (n = 0; n < size; n++)
+        {
+            total += linein[size*n + column]*pcp[u*size+n];
+        }
+        lineout[copy_matrix[size*u + column]] = total*oosob*qt[size*u + column];
+    }
+#endif
+}
+
+inline void dct_column_large(F32 *linein, S32 *lineout, S32 column)
+{
+    S32 u;
+    F32 total;
+    F32 oosob = 2.f/32.f;
+    F32 *pcp = gPatchCosines;
+    S32 *copy_matrix = gCopyMatrix;
+    F32 *qt = gPatchQuantizeTable;
+
+    F32 *tlinein, *tpcp;
+    S32 sizeu;
+
+    tlinein = linein + column;
+
+    total = *(tlinein);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+    total += *(tlinein += LARGE_PATCH_SIZE);
+
+    *(lineout + *(copy_matrix + column)) = (S32)(OO_SQRT2*total*oosob*(*(qt + column)));
+
+    for (u = 1; u < LARGE_PATCH_SIZE; u++)
+    {
+        tlinein = linein + column;
+        tpcp = pcp + (u<<5);
+
+        total = *(tlinein)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp++));
+        total += *(tlinein += LARGE_PATCH_SIZE)*(*(tpcp));
+
+        sizeu = LARGE_PATCH_SIZE*u + column;
+
+        *(lineout + *(copy_matrix + sizeu)) = (S32)(total*oosob*(*(qt+sizeu)));
+    }
+}
+
+inline void dct_patch(F32 *block, S32 *cpatch)
+{
+    F32 temp[NORMAL_PATCH_SIZE*NORMAL_PATCH_SIZE];
+
+#ifdef _PATCH_SIZE_16_AND_32_ONLY
+    dct_line(block, temp, 0);
+    dct_line(block, temp, 1);
+    dct_line(block, temp, 2);
+    dct_line(block, temp, 3);
+
+    dct_line(block, temp, 4);
+    dct_line(block, temp, 5);
+    dct_line(block, temp, 6);
+    dct_line(block, temp, 7);
+
+    dct_line(block, temp, 8);
+    dct_line(block, temp, 9);
+    dct_line(block, temp, 10);
+    dct_line(block, temp, 11);
+
+    dct_line(block, temp, 12);
+    dct_line(block, temp, 13);
+    dct_line(block, temp, 14);
+    dct_line(block, temp, 15);
+
+    dct_column(temp, cpatch, 0);
+    dct_column(temp, cpatch, 1);
+    dct_column(temp, cpatch, 2);
+    dct_column(temp, cpatch, 3);
+
+    dct_column(temp, cpatch, 4);
+    dct_column(temp, cpatch, 5);
+    dct_column(temp, cpatch, 6);
+    dct_column(temp, cpatch, 7);
+
+    dct_column(temp, cpatch, 8);
+    dct_column(temp, cpatch, 9);
+    dct_column(temp, cpatch, 10);
+    dct_column(temp, cpatch, 11);
+
+    dct_column(temp, cpatch, 12);
+    dct_column(temp, cpatch, 13);
+    dct_column(temp, cpatch, 14);
+    dct_column(temp, cpatch, 15);
+#else
+    S32 i;
+    S32 size = gPatchCompressGlobalData.patch_size;
+    for (i = 0; i < size; i++)
+    {
+        dct_line(block, temp, i);
+    }
+    for (i = 0; i < size; i++)
+    {
+        dct_column(temp, cpatch, i);
+    }
+#endif
+}
+
+inline void dct_patch_large(F32 *block, S32 *cpatch)
+{
+    F32 temp[LARGE_PATCH_SIZE*LARGE_PATCH_SIZE];
+
+    dct_line_large(block, temp, 0);
+    dct_line_large(block, temp, 1);
+    dct_line_large(block, temp, 2);
+    dct_line_large(block, temp, 3);
+
+    dct_line_large(block, temp, 4);
+    dct_line_large(block, temp, 5);
+    dct_line_large(block, temp, 6);
+    dct_line_large(block, temp, 7);
+
+    dct_line_large(block, temp, 8);
+    dct_line_large(block, temp, 9);
+    dct_line_large(block, temp, 10);
+    dct_line_large(block, temp, 11);
+
+    dct_line_large(block, temp, 12);
+    dct_line_large(block, temp, 13);
+    dct_line_large(block, temp, 14);
+    dct_line_large(block, temp, 15);
+
+    dct_line_large(block, temp, 16);
+    dct_line_large(block, temp, 17);
+    dct_line_large(block, temp, 18);
+    dct_line_large(block, temp, 19);
+
+    dct_line_large(block, temp, 20);
+    dct_line_large(block, temp, 21);
+    dct_line_large(block, temp, 22);
+    dct_line_large(block, temp, 23);
+
+    dct_line_large(block, temp, 24);
+    dct_line_large(block, temp, 25);
+    dct_line_large(block, temp, 26);
+    dct_line_large(block, temp, 27);
+
+    dct_line_large(block, temp, 28);
+    dct_line_large(block, temp, 29);
+    dct_line_large(block, temp, 30);
+    dct_line_large(block, temp, 31);
+
+    dct_column_large(temp, cpatch, 0);
+    dct_column_large(temp, cpatch, 1);
+    dct_column_large(temp, cpatch, 2);
+    dct_column_large(temp, cpatch, 3);
+
+    dct_column_large(temp, cpatch, 4);
+    dct_column_large(temp, cpatch, 5);
+    dct_column_large(temp, cpatch, 6);
+    dct_column_large(temp, cpatch, 7);
+
+    dct_column_large(temp, cpatch, 8);
+    dct_column_large(temp, cpatch, 9);
+    dct_column_large(temp, cpatch, 10);
+    dct_column_large(temp, cpatch, 11);
+
+    dct_column_large(temp, cpatch, 12);
+    dct_column_large(temp, cpatch, 13);
+    dct_column_large(temp, cpatch, 14);
+    dct_column_large(temp, cpatch, 15);
+
+    dct_column_large(temp, cpatch, 16);
+    dct_column_large(temp, cpatch, 17);
+    dct_column_large(temp, cpatch, 18);
+    dct_column_large(temp, cpatch, 19);
+
+    dct_column_large(temp, cpatch, 20);
+    dct_column_large(temp, cpatch, 21);
+    dct_column_large(temp, cpatch, 22);
+    dct_column_large(temp, cpatch, 23);
+
+    dct_column_large(temp, cpatch, 24);
+    dct_column_large(temp, cpatch, 25);
+    dct_column_large(temp, cpatch, 26);
+    dct_column_large(temp, cpatch, 27);
+
+    dct_column_large(temp, cpatch, 28);
+    dct_column_large(temp, cpatch, 29);
+    dct_column_large(temp, cpatch, 30);
+    dct_column_large(temp, cpatch, 31);
+}
+
+void compress_patch(F32 *patch, S32 *cpatch, LLPatchHeader *php, S32 prequant)
+{
+    S32     i, j;
+    PCGD    *pcp = &gPatchCompressGlobalData;
+    S32     stride = pcp->patch_stride;
+    S32     size = pcp->patch_size;
+    F32     block[LARGE_PATCH_SIZE*LARGE_PATCH_SIZE], *tblock;
+    F32     *tpatch;
+
+    S32     wordsize = prequant;
+    F32     oozrange = 1.f/php->range;
+
+    F32     dc = php->dc_offset;
+
+    S32     range = (1<<prequant);
+    F32     premult = oozrange*range;
+//  F32     sub = (F32)(1<<(prequant - 1));
+    F32     sub = (F32)(1<<(prequant - 1)) + dc*premult;
+
+    php->quant_wbits = wordsize - 2;
+    php->quant_wbits |= (prequant - 2)<<4;
+
+    for (j = 0; j < size; j++)
+    {
+        tblock = block + j*size;
+        tpatch = patch + j*stride;
+        for (i = 0; i < size; i++)
+        {
+//          block[j*size + i] = (patch[j*stride + i] - dc)*premult - sub;
+            *(tblock++) = *(tpatch++)*premult - sub;
+        }
+    }
+
+    if (size == 16)
+        dct_patch(block, cpatch);
+    else
+        dct_patch_large(block, cpatch);
+}
+
+void get_patch_group_header(LLGroupHeader *gopp)
+{
+    PCGD    *pcp = &gPatchCompressGlobalData;
+    gopp->stride = pcp->patch_stride;
+    gopp->patch_size = pcp->patch_size;
+    gopp->layer_type = pcp->layer_type;
+}