444.c
上传用户:dangjiwu
上传日期:2013-07-19
资源大小:42019k
文件大小:32k
源码类别:

Symbian

开发平台:

Visual C++

  1. /* ***** BEGIN LICENSE BLOCK *****
  2.  * Source last modified: $Id: 444.c,v 1.2.8.1 2004/07/09 02:00:07 hubbe Exp $
  3.  * 
  4.  * Portions Copyright (c) 1995-2004 RealNetworks, Inc. All Rights Reserved.
  5.  * 
  6.  * The contents of this file, and the files included with this file,
  7.  * are subject to the current version of the RealNetworks Public
  8.  * Source License (the "RPSL") available at
  9.  * http://www.helixcommunity.org/content/rpsl unless you have licensed
  10.  * the file under the current version of the RealNetworks Community
  11.  * Source License (the "RCSL") available at
  12.  * http://www.helixcommunity.org/content/rcsl, in which case the RCSL
  13.  * will apply. You may also obtain the license terms directly from
  14.  * RealNetworks.  You may not use this file except in compliance with
  15.  * the RPSL or, if you have a valid RCSL with RealNetworks applicable
  16.  * to this file, the RCSL.  Please see the applicable RPSL or RCSL for
  17.  * the rights, obligations and limitations governing use of the
  18.  * contents of the file.
  19.  * 
  20.  * Alternatively, the contents of this file may be used under the
  21.  * terms of the GNU General Public License Version 2 or later (the
  22.  * "GPL") in which case the provisions of the GPL are applicable
  23.  * instead of those above. If you wish to allow use of your version of
  24.  * this file only under the terms of the GPL, and not to allow others
  25.  * to use your version of this file under the terms of either the RPSL
  26.  * or RCSL, indicate your decision by deleting the provisions above
  27.  * and replace them with the notice and other provisions required by
  28.  * the GPL. If you do not delete the provisions above, a recipient may
  29.  * use your version of this file under the terms of any one of the
  30.  * RPSL, the RCSL or the GPL.
  31.  * 
  32.  * This file is part of the Helix DNA Technology. RealNetworks is the
  33.  * developer of the Original Code and owns the copyrights in the
  34.  * portions it created.
  35.  * 
  36.  * This file, and the files included with this file, is distributed
  37.  * and made available on an 'AS IS' basis, WITHOUT WARRANTY OF ANY
  38.  * KIND, EITHER EXPRESS OR IMPLIED, AND REALNETWORKS HEREBY DISCLAIMS
  39.  * ALL SUCH WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES
  40.  * OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE, QUIET
  41.  * ENJOYMENT OR NON-INFRINGEMENT.
  42.  * 
  43.  * Technology Compatibility Kit Test Suite(s) Location:
  44.  *    http://www.helixcommunity.org/content/tck
  45.  * 
  46.  * Contributor(s):
  47.  * 
  48.  * ***** END LICENSE BLOCK ***** */
  49. #include <stdlib.h>
  50. #include "nostatic/colorlib.h"
  51. #include "nostatic/yuv.h"
  52. /*
  53.  * Format-conversion routines.
  54.  * Use:
  55.  *  int XXXXtoYYYYEx (unsigned char *dest_ptr, int dest_width, int dest_height,
  56.  *      int dest_pitch, int dest_x, int dest_y, int dest_dx, int dest_dy,
  57.  *      unsigned char *src_ptr, int src_width, int src_height, int src_pitch,
  58.  *      int src_x, int src_y, int src_dx, int src_dy);
  59.  * Input:
  60.  *  dest_ptr - pointer to a destination buffer
  61.  *  dest_width, dest_height - width/height of the destination image (pixels)
  62.  *  dest_pitch - pitch of the dest. buffer (in bytes; <0 - if bottom up image)
  63.  *  dest_x, dest_y, dest_dx, dest_dy - destination rectangle (pixels)
  64.  *  src_ptr - pointer to an input image
  65.  *  src_width, src_height - width/height of the input image (pixels)
  66.  *  src_pitch - pitch of the source buffer (in bytes; <0 - if bottom up image)
  67.  *  src_x, src_y, src_dx, src_dy - source rectangle (pixels)
  68.  * Returns:
  69.  *  0 - if success; -1 if failure.
  70.  * Notes:
  71.  *  a) In all cases, pointers to the source and destination buffers must be
  72.  *     DWORD aligned, and both pitch parameters must be multiple of 4!!!
  73.  *  b) Converters that deal with YUV 4:2:2, or 4:2:0 formats may also require
  74.  *     rectangle parameters (x,y,dx,dy) to be multiple of 2. Failure to provide
  75.  *     aligned rectangles will result in partially converted image.
  76.  *  c) Currently only scale factors of 1:1 and 2:1 are supported; if the rates
  77.  *     dest_dx/src_dx & dest_dy/src_dy are neither 1, or 2, the converters
  78.  *     will fail.
  79.  */
  80. #define BYTESPERPIXEL 2
  81. /* Indexing tables */
  82. static const int next[3] = {1, 2, 0};                 /* (ibuf+1)%3 table */
  83. static const int next2[3] = {2, 0, 1};                /* (ibuf+2)%3 table */
  84. /* Half resolution converter.  */
  85. static void halfI420toRGB444 (unsigned char *d1,
  86. int dest_x, unsigned char *sy1, 
  87. unsigned char *sy2, unsigned char *su, 
  88. unsigned char *sv, int src_x, int dx, int isodd,
  89.       color_data_t* pcd)
  90. {
  91.   /* check if we have misaligned input/output: */
  92.   if ((src_x ^ dest_x) & 1) {
  93.     
  94.     ; /* not implemented yet */
  95.     
  96.   } else {
  97.     
  98.     /* aligned input/output: */
  99.     register int y11, y12;
  100.     register int rv1, guv1, bu1, rv2, guv2, bu2;
  101.     register int i;
  102.     
  103.     /* source width should always be multiples of 4 */
  104.     /*    assert((src_x & 3)==0); */
  105.     /* convert all integral 2x4 blocks: */
  106.     if (isodd)
  107.       for (i = 0; i < dx/4; i ++) {
  108. /* first 2x2 block */
  109. rv1 = pcd->rvtab[sv[0]];
  110. guv1 = pcd->gutab[su[0]] + pcd->gvtab[sv[0]];
  111. bu1 = pcd->butab[su[0]];
  112. y11 = (pcd->ytab[sy1[0]] + pcd->ytab[sy1[1]] + pcd->ytab[sy2[0]] + pcd->ytab[sy2[1]])>>2; /* avg */
  113. /* second 2x2 block */
  114. rv2  = pcd->rvtab[sv[1]];
  115. guv2 = pcd->gutab[su[1]] + pcd->gvtab[sv[1]];
  116. bu2  = pcd->butab[su[1]];
  117. y12 = (pcd->ytab[sy1[2]] + pcd->ytab[sy1[3]] + pcd->ytab[sy2[2]] + pcd->ytab[sy2[3]])>>2;
  118. /* output 4 bytes at a time */
  119. *(unsigned int *)(d1+0) =
  120.   (CLIP4[y11 + bu1  + DITH4L] << 0)  |
  121.   (CLIP4[y11 + guv1 + DITH4L] << 4)  |
  122.   (CLIP4[y11 + rv1  + DITH4L] << 8) |
  123.   (CLIP4[y12 + bu2  + DITH4H] << 16) |
  124.   (CLIP4[y12 + guv2 + DITH4H] << 20) |
  125.   (CLIP4[y12 + rv2  + DITH4H] << 24) ;
  126. /* next 4x2 block */
  127. sy1 += 4; sy2 += 4;
  128. su += 2; sv += 2;
  129. d1 += 2*BPP2;
  130.       }
  131.     else
  132.       for (i = 0; i < dx/4; i ++) {
  133. /* first 2x2 block */
  134. rv1 = pcd->rvtab[sv[0]];
  135. guv1 = pcd->gutab[su[0]] + pcd->gvtab[sv[0]];
  136. bu1 = pcd->butab[su[0]];
  137. y11 = (pcd->ytab[sy1[0]] + pcd->ytab[sy1[1]] + pcd->ytab[sy2[0]] + pcd->ytab[sy2[1]])>>2; /* avg */
  138. /* second 2x2 block */
  139. rv2  = pcd->rvtab[sv[1]];
  140. guv2 = pcd->gutab[su[1]] + pcd->gvtab[sv[1]];
  141. bu2  = pcd->butab[su[1]];
  142. y12 = (pcd->ytab[sy1[2]] + pcd->ytab[sy1[3]] + pcd->ytab[sy2[2]] + pcd->ytab[sy2[3]])>>2;
  143. /* output 4 bytes at a time */
  144. *(unsigned int *)(d1+0) =
  145.   (CLIP4[y11 + bu1  + DITH4H] << 0)  |
  146.   (CLIP4[y11 + guv1 + DITH4H] << 4)  |
  147.   (CLIP4[y11 + rv1  + DITH4H] << 8) |
  148.   (CLIP4[y12 + bu2  + DITH4L] << 16) |
  149.   (CLIP4[y12 + guv2 + DITH4L] << 20) |
  150.   (CLIP4[y12 + rv2  + DITH4L] << 24) ;
  151. /* next 4x2 block */
  152. sy1 += 4; sy2 += 4;
  153. su += 2; sv += 2;
  154. d1 += 2*BPP2;
  155.       }
  156.   }
  157. }
  158. static void dblineI420toRGB444 (unsigned char *d1, unsigned char *d2, 
  159. int dest_x, unsigned char *sy1, 
  160. unsigned char *sy2, unsigned char *su, 
  161. unsigned char *sv, int src_x, int dx,
  162. color_data_t* pcd)
  163. {
  164.     /* check if we have misaligned input/output: */
  165.     if ((src_x ^ dest_x) & 1) {
  166.         ; /* not implemented yet */
  167.     } else {
  168.         /* aligned input/output: */
  169.         register int y11, y12, y21, y22;
  170.         register int rv, guv, bu;
  171.         register int i;
  172.         /* convert first 2x1 block: */
  173.         if (dest_x & 1) {
  174.             rv  = pcd->rvtab[sv[0]];
  175.             guv = pcd->gutab[su[0]] + pcd->gvtab[sv[0]];
  176.             bu  = pcd->butab[su[0]];
  177.             y11 = pcd->ytab[sy1[0]];
  178.             y21 = pcd->ytab[sy2[0]];
  179.             /* output 2 bytes at a time */
  180.             *(unsigned short *)(d1+0) =
  181.                 (CLIP4[y11 + bu  + DITH4L] << 0)  |
  182.                 (CLIP4[y11 + guv + DITH4L] << 4)  |
  183.                 (CLIP4[y11 + rv  + DITH4L] << 8) ;
  184.             *(unsigned short *)(d2+0) =
  185.                 (CLIP4[y21 + bu  + DITH4H] << 0)  |
  186.                 (CLIP4[y21 + guv + DITH4H] << 4)  |
  187.                 (CLIP4[y21 + rv  + DITH4H] << 8) ;
  188.             sy1 += 1; sy2 += 1;
  189.             su += 1; sv += 1;
  190.             d1 += BPP2; d2 += BPP2;
  191.             dest_x ++; dx --;
  192.         }
  193.         /* convert first 2x2 block: */
  194.         if (dest_x & 2) {
  195.             rv  = pcd->rvtab[sv[0]];
  196.             guv = pcd->gutab[su[0]] + pcd->gvtab[sv[0]];
  197.             bu  = pcd->butab[su[0]];
  198.             y11 = pcd->ytab[sy1[0]];
  199.             y12 = pcd->ytab[sy1[1]];
  200.             /* output 4 bytes at a time */
  201.             *(unsigned int *)(d1+0) =
  202.                 (CLIP4[y11 + bu  + DITH4L] << 0)  |
  203.                 (CLIP4[y11 + guv + DITH4L] << 4)  |
  204.                 (CLIP4[y11 + rv  + DITH4L] << 8) |
  205.                 (CLIP4[y12 + bu  + DITH4H] << 16) |
  206.                 (CLIP4[y12 + guv + DITH4H] << 20) |
  207.                 (CLIP4[y12 + rv  + DITH4H] << 24) ;
  208.             y21 = pcd->ytab[sy2[0]];
  209.             y22 = pcd->ytab[sy2[1]];
  210.             *(unsigned int *)(d2+0) =
  211.                 (CLIP4[y21 + bu  + DITH4H] << 0)  |
  212.                 (CLIP4[y21 + guv + DITH4H] << 4)  |
  213.                 (CLIP4[y21 + rv  + DITH4H] << 8) |
  214.                 (CLIP4[y22 + bu  + DITH4L] << 16) |
  215.                 (CLIP4[y22 + guv + DITH4L] << 20) |
  216.                 (CLIP4[y22 + rv  + DITH4L] << 24) ;
  217.             sy1 += 2; sy2 += 2;
  218.             su += 1; sv += 1;
  219.             d1 += 2*BPP2; d2 += 2*BPP2;
  220.             dest_x += 2; dx -= 2;
  221.         }
  222.         /* convert all integral 2x4 blocks: */
  223.         for (i = 0; i < dx/4; i ++) {
  224.             /* first 2x2 block */
  225.             rv = pcd->rvtab[sv[0]];
  226.             guv = pcd->gutab[su[0]] + pcd->gvtab[sv[0]];
  227.             bu = pcd->butab[su[0]];
  228.             y11 = pcd->ytab[sy1[0]];
  229.             y12 = pcd->ytab[sy1[1]];
  230.             /* output 4 bytes at a time */
  231.             *(unsigned int *)(d1+0) =
  232.                 (CLIP4[y11 + bu  + DITH4L] << 0)  |
  233.                 (CLIP4[y11 + guv + DITH4L] << 4)  |
  234.                 (CLIP4[y11 + rv  + DITH4L] << 8) |
  235.                 (CLIP4[y12 + bu  + DITH4H] << 16) |
  236.                 (CLIP4[y12 + guv + DITH4H] << 20) |
  237.                 (CLIP4[y12 + rv  + DITH4H] << 24) ;
  238.             y21 = pcd->ytab[sy2[0]];
  239.             y22 = pcd->ytab[sy2[1]];
  240.             *(unsigned int *)(d2+0) =
  241.                 (CLIP4[y21 + bu  + DITH4H] << 0)  |
  242.                 (CLIP4[y21 + guv + DITH4H] << 4)  |
  243.                 (CLIP4[y21 + rv  + DITH4H] << 8) |
  244.                 (CLIP4[y22 + bu  + DITH4L] << 16) |
  245.                 (CLIP4[y22 + guv + DITH4L] << 20) |
  246.                 (CLIP4[y22 + rv  + DITH4L] << 24) ;
  247.             /* second 2x2 block */
  248.             rv  = pcd->rvtab[sv[1]];
  249.             guv = pcd->gutab[su[1]] + pcd->gvtab[sv[1]];
  250.             bu  = pcd->butab[su[1]];
  251.             y11 = pcd->ytab[sy1[2]];
  252.             y12 = pcd->ytab[sy1[3]];
  253.             *(unsigned int *)(d1+2*BPP2) =
  254.                 (CLIP4[y11 + bu  + DITH4L] << 0)  |
  255.                 (CLIP4[y11 + guv + DITH4L] << 4)  |
  256.                 (CLIP4[y11 + rv  + DITH4L] << 8) |
  257.                 (CLIP4[y12 + bu  + DITH4H] << 16) |
  258.                 (CLIP4[y12 + guv + DITH4H] << 20) |
  259.                 (CLIP4[y12 + rv  + DITH4H] << 24) ;
  260.             y21 = pcd->ytab[sy2[2]];
  261.             y22 = pcd->ytab[sy2[3]];
  262.             *(unsigned int *)(d2+2*BPP2) =
  263.                 (CLIP4[y21 + bu  + DITH4H] << 0)  |
  264.                 (CLIP4[y21 + guv + DITH4H] << 4)  |
  265.                 (CLIP4[y21 + rv  + DITH4H] << 8) |
  266.                 (CLIP4[y22 + bu  + DITH4L] << 16) |
  267.                 (CLIP4[y22 + guv + DITH4L] << 20) |
  268.                 (CLIP4[y22 + rv  + DITH4L] << 24) ;
  269.             /* next 4x2 block */
  270.             sy1 += 4; sy2 += 4;
  271.             su += 2; sv += 2;
  272.             d1 += 4*BPP2; d2 += 4*BPP2;
  273.         }
  274.         /* convert last 2x2 block: */
  275.         if (dx & 2) {
  276.             rv = pcd->rvtab[sv[0]];
  277.             guv = pcd->gutab[su[0]] + pcd->gvtab[sv[0]];
  278.             bu = pcd->butab[su[0]];
  279.             y11 = pcd->ytab[sy1[0]];
  280.             y12 = pcd->ytab[sy1[1]];
  281.             /* output 4 bytes at a time */
  282.             *(unsigned int *)(d1+0) =
  283.                 (CLIP4[y11 + bu  + DITH4L] << 0)  |
  284.                 (CLIP4[y11 + guv + DITH4L] << 4)  |
  285.                 (CLIP4[y11 + rv  + DITH4L] << 8) |
  286.                 (CLIP4[y12 + bu  + DITH4H] << 16) |
  287.                 (CLIP4[y12 + guv + DITH4H] << 20) |
  288.                 (CLIP4[y12 + rv  + DITH4H] << 24) ;
  289.             y21 = pcd->ytab[sy2[0]];
  290.             y22 = pcd->ytab[sy2[1]];
  291.             *(unsigned int *)(d2+0) =
  292.                 (CLIP4[y21 + bu  + DITH4H] << 0)  |
  293.                 (CLIP4[y21 + guv + DITH4H] << 4)  |
  294.                 (CLIP4[y21 + rv  + DITH4H] << 8) |
  295.                 (CLIP4[y22 + bu  + DITH4L] << 16) |
  296.                 (CLIP4[y22 + guv + DITH4L] << 20) |
  297.                 (CLIP4[y22 + rv  + DITH4L] << 24) ;
  298.             sy1 += 2; sy2 += 2;
  299.             su += 1; sv += 1;
  300.             d1 += 2*BPP2; d2 += 2*BPP2;
  301.             dx -= 2;
  302.         }
  303.         /* convert last 2x1 block: */
  304.         if (dx & 1) {
  305.             rv  = pcd->rvtab[sv[0]];
  306.             guv = pcd->gutab[su[0]] + pcd->gvtab[sv[0]];
  307.             bu  = pcd->butab[su[0]];
  308.             y11 = pcd->ytab[sy1[0]];
  309.             y21 = pcd->ytab[sy2[0]];
  310.             /* output 2 bytes at a time */
  311.             *(unsigned short *)(d1+0) =
  312.                 (CLIP4[y11 + bu  + DITH4L] << 0)  |
  313.                 (CLIP4[y11 + guv + DITH4L] << 4)  |
  314.                 (CLIP4[y11 + rv  + DITH4L] << 8) ;
  315.             *(unsigned short *)(d2+0) =
  316.                 (CLIP4[y21 + bu  + DITH4H] << 0)  |
  317.                 (CLIP4[y21 + guv + DITH4H] << 4)  |
  318.                 (CLIP4[y21 + rv  + DITH4H] << 8) ;
  319.         }
  320.     }
  321. }
  322. #define DITH1X4L 0x01002004
  323. #define DITH1X4H 0x0300600C
  324. #define DITH2X4L 0x02004008
  325. #define DITH2X4H 0x0600C018
  326. #define DITH4X4L 0x04008010
  327. #define DITH4X4H 0x0C018030
  328. /* RGB444 version: */
  329. static void dblineXRGBtoRGB444x2 (unsigned char *d1, unsigned char *d2, int dest_x,
  330.     unsigned char *s1, unsigned char *s2, int src_x, int dx,
  331. color_data_t* pcd)
  332. {
  333.     register unsigned int a, b, c, d, e, f;
  334.     register unsigned int w, x, y, z;
  335.     /* convert first 2x1 block: */
  336.     if (dest_x & 1) {
  337.         /* Input pels  =>   Output pels
  338.          *  a b             w  x
  339.          *  c d             w' x'
  340.          */
  341.         /* top line */
  342.         a = *(unsigned int *)s1;
  343.         b = *(unsigned int *)(s1+BPP4);
  344.         w = a     + DITH1X4L;
  345.         x = a + b + DITH2X4H;
  346.         /* pack and store */
  347.         *(unsigned int *)d1 =
  348.             ((w & 0x000000f0) >> 4)  |
  349.             ((w & 0x00078000) >> 11) |
  350.             ((w & 0x3c000000) >> 18) |
  351.             ((x & 0x000001e0) << 11) |
  352.             ((x & 0x000f0000) << 4)  |
  353.             ((x & 0x78000000) >> 3);
  354.         /* bottom line */
  355.         c = *(unsigned int *)s2;
  356.         d = *(unsigned int *)(s2+BPP4);
  357.         w = a + c         + DITH2X4H;
  358.         x = a + b + c + d + DITH4X4L;
  359.         /* pack and store */
  360.         *(unsigned int *)d2 =
  361.             ((w & 0x000001e0) >> 5)  |
  362.             ((w & 0x000f0000) >> 12) |
  363.             ((w & 0x78000000) >> 19) |
  364.             ((x & 0x000003c0) << 10) |
  365.             ((x & 0x001e0000) << 3)  |
  366.             ((x & 0xf0000000) >> 4);
  367.         /* bump pointers to next block */
  368.         s1 += BPP4; s2 += BPP4;
  369.         d1 += 2*BPP2; d2 += 2*BPP2;
  370.         dx -= 1;
  371.     }
  372.     /* process all integral 2x2 blocks: */
  373.     while (dx > 2) {    /* we need to leave at least one block */
  374.         /*
  375.          * Input stored as 00 RRRRRRRR 000 GGGGGGGG 000 BBBBBBBB
  376.          *
  377.          * Input pels       Output pels
  378.          *  a b e           w  x  y  z
  379.          *  c d f           w' x' y' z'
  380.          */
  381.         /* top line */
  382.         a = *(unsigned int *)s1;
  383.         b = *(unsigned int *)(s1+BPP4);
  384.         e = *(unsigned int *)(s1+2*BPP4);
  385.         w = a     + DITH1X4L;
  386.         x = a + b + DITH2X4H;
  387.         y = b     + DITH1X4L;
  388.         z = b + e + DITH2X4H;
  389.         /* pack and store */
  390.         *(unsigned int *)d1 =
  391.             ((w & 0x000000f0) >> 4)  |
  392.             ((w & 0x00078000) >> 11) |
  393.             ((w & 0x3c000000) >> 18) |
  394.             ((x & 0x000001e0) << 11) |
  395.             ((x & 0x000f0000) << 4)  |
  396.             ((x & 0x78000000) >> 3);
  397.         *(unsigned int *)(d1+2*BPP2) =
  398.             ((y & 0x000000f0) >> 4)  |
  399.             ((y & 0x00078000) >> 11) |
  400.             ((y & 0x3c000000) >> 18) |
  401.             ((z & 0x000001e0) << 11) |
  402.             ((z & 0x000f0000) << 4)  |
  403.             ((z & 0x78000000) >> 3);
  404.         /* bottom line */
  405.         c = *(unsigned int *)s2;
  406.         d = *(unsigned int *)(s2+BPP4);
  407.         f = *(unsigned int *)(s2+2*BPP4);
  408.         w = a + c         + DITH2X4H;
  409.         x = a + b + c + d + DITH4X4L;
  410.         y = b + d         + DITH2X4H;
  411.         z = b + e + d + f + DITH4X4L;
  412.         /* pack and store */
  413.         *(unsigned int *)d2 =
  414.             ((w & 0x000001e0) >> 5)  |
  415.             ((w & 0x000f0000) >> 12) |
  416.             ((w & 0x78000000) >> 19) |
  417.             ((x & 0x000003c0) << 10) |
  418.             ((x & 0x001e0000) << 3)  |
  419.             ((x & 0xf0000000) >> 4);
  420.         *(unsigned int *)(d2+2*BPP2) =
  421.             ((y & 0x000001e0) >> 5)  |
  422.             ((y & 0x000f0000) >> 12) |
  423.             ((y & 0x78000000) >> 19) |
  424.             ((z & 0x000003c0) << 10) |
  425.             ((z & 0x001e0000) << 3)  |
  426.             ((z & 0xf0000000) >> 4);
  427.         /* next 2x2 input block */
  428.         s1 += 2*BPP4; s2 += 2*BPP4;
  429.         d1 += 4*BPP2; d2 += 4*BPP2;
  430.         dx -= 2;
  431.     }
  432.     /* check if this is the last 2x2 block: */
  433.     if (dx > 1) {
  434.         /*
  435.          * For last 4 output pels, repeat final input pel
  436.          * for offscreen input.  Equivalent to pixel-doubling the
  437.          * last output pel.
  438.          */
  439.         /* top line */
  440.         a = *(unsigned int *)s1;
  441.         b = *(unsigned int *)(s1+BPP4);
  442.         e = b;      /* repeat last input pel */
  443.         w = a     + DITH1X4L;
  444.         x = a + b + DITH2X4H;
  445.         y = b     + DITH1X4L;
  446.         z = b + e + DITH2X4H;
  447.         /* pack and store */
  448.         *(unsigned int *)d1 =
  449.             ((w & 0x000000f0) >> 4)  |
  450.             ((w & 0x00078000) >> 11) |
  451.             ((w & 0x3c000000) >> 18) |
  452.             ((x & 0x000001e0) << 11) |
  453.             ((x & 0x000f0000) << 4)  |
  454.             ((x & 0x78000000) >> 3);
  455.         *(unsigned int *)(d1+2*BPP2) =
  456.             ((y & 0x000000f0) >> 4)  |
  457.             ((y & 0x00078000) >> 11) |
  458.             ((y & 0x3c000000) >> 18) |
  459.             ((z & 0x000001e0) << 11) |
  460.             ((z & 0x000f0000) << 4)  |
  461.             ((z & 0x78000000) >> 3);
  462.         /* bottom line */
  463.         c = *(unsigned int *)s2;
  464.         d = *(unsigned int *)(s2+BPP4);
  465.         f = d;      /* repeat last input pel */
  466.         w = a + c         + DITH2X4H;
  467.         x = a + b + c + d + DITH4X4L;
  468.         y = b + d         + DITH2X4H;
  469.         z = b + e + d + f + DITH4X4L;
  470.         *(unsigned int *)d2 =
  471.             ((w & 0x000001e0) >> 5)  |
  472.             ((w & 0x000f0000) >> 12) |
  473.             ((w & 0x78000000) >> 19) |
  474.             ((x & 0x000003c0) << 10) |
  475.             ((x & 0x001e0000) << 3)  |
  476.             ((x & 0xf0000000) >> 4);
  477.         *(unsigned int *)(d2+2*BPP2) =
  478.             ((y & 0x000001e0) >> 5)  |
  479.             ((y & 0x000f0000) >> 12) |
  480.             ((y & 0x78000000) >> 19) |
  481.             ((z & 0x000003c0) << 10) |
  482.             ((z & 0x001e0000) << 3)  |
  483.             ((z & 0xf0000000) >> 4);
  484.     } else {
  485.         /* last 2x1 block: */
  486.         /* Input pels  =>   Output pels
  487.          *  a b             w  x
  488.          *  c d             w' x'
  489.          */
  490.         /* top line */
  491.         a = *(unsigned int *)s1;
  492.         b = a;      /* repeat last input pel */
  493.         w = a     + DITH1X4L;
  494.         x = a + b + DITH2X4H;
  495.         /* pack and store */
  496.         *(unsigned int *)d1 =
  497.             ((w & 0x000000f0) >> 4)  |
  498.             ((w & 0x00078000) >> 11) |
  499.             ((w & 0x3c000000) >> 18) |
  500.             ((x & 0x000001e0) << 11) |
  501.             ((x & 0x000f0000) << 4)  |
  502.             ((x & 0x78000000) >> 3);
  503.         /* bottom line */
  504.         c = *(unsigned int *)s2;
  505.         d = c;      /* repeat last input pel */
  506.         w = a + c         + DITH2X4H;
  507.         x = a + b + c + d + DITH4X4L;
  508.         /* pack and store */
  509.         *(unsigned int *)d2 =
  510.             ((w & 0x000001e0) >> 5)  |
  511.             ((w & 0x000f0000) >> 12) |
  512.             ((w & 0x78000000) >> 19) |
  513.             ((x & 0x000003c0) << 10) |
  514.             ((x & 0x001e0000) << 3)  |
  515.             ((x & 0xf0000000) >> 4);
  516.     }
  517. }
  518. /*
  519.  * Convert two YUV lines into RGB linebufs.
  520.  * Produces two RGB lines per call.
  521.  * Output in padded RGB format, needed for SIMD interpolation.
  522.  */
  523. static void dblineI420toXRGB (unsigned char *d1, unsigned char *d2, int dest_x,
  524.     unsigned char *sy1, unsigned char *sy2, unsigned char *su, unsigned char *sv,
  525.       int src_x, int dx,
  526.       color_data_t* pcd)
  527. {
  528.     register int y1, y2, rv, guv, bu;
  529.     register int i;
  530.     /* convert first 2x1 block: */
  531.     if (src_x & 1) {
  532.         bu  = pcd->butab[su[0]];
  533.         guv = pcd->gutab[su[0]] + pcd->gvtab[sv[0]];
  534.         rv  = pcd->rvtab[sv[0]];
  535.         y1  = pcd->ytab[sy1[0]];
  536.         y2  = pcd->ytab[sy2[0]];
  537.         /* first line 00rrrrrr.rr000ggg.ggggg000.bbbbbbbb: */
  538.         *(unsigned int *)(d1+0) =
  539.             (CLIP8[y1 + bu])        |
  540.             (CLIP8[y1 + guv] << 11) |
  541.             (CLIP8[y1 + rv]  << 22) ;
  542.         /* second line 00rrrrrr.rr000ggg.ggggg000.bbbbbbbb: */
  543.         *(unsigned int *)(d2+0) =
  544.             (CLIP8[y2 + bu])        |
  545.             (CLIP8[y2 + guv] << 11) |
  546.             (CLIP8[y2 + rv]  << 22) ;
  547.         sy1 += 1; sy2 += 1;
  548.         su += 1; sv += 1;
  549.         d1 += BPP4; d2 += BPP4;
  550.         dx --;
  551.     }
  552.     /* convert all integral 2x2 blocks: */
  553.     for (i = 0; i < dx/2; i ++) {
  554.         bu = pcd->butab[su[0]];
  555.         guv = pcd->gutab[su[0]] + pcd->gvtab[sv[0]];
  556.         rv = pcd->rvtab[sv[0]];
  557.         y1 = pcd->ytab[sy1[0]];
  558.         y2 = pcd->ytab[sy2[0]];
  559.         /* first line 00rrrrrr.rr000ggg.ggggg000.bbbbbbbb: */
  560.         *(unsigned int *)(d1+0) =
  561.             (CLIP8[y1 + bu])        |
  562.             (CLIP8[y1 + guv] << 11) |
  563.             (CLIP8[y1 + rv]  << 22) ;
  564.         /* second line 00rrrrrr.rr000ggg.ggggg000.bbbbbbbb: */
  565.         *(unsigned int *)(d2+0) =
  566.             (CLIP8[y2 + bu])        |
  567.             (CLIP8[y2 + guv] << 11) |
  568.             (CLIP8[y2 + rv]  << 22) ;
  569.         /* 2nd row: */
  570.         y1 = pcd->ytab[sy1[1]];
  571.         y2 = pcd->ytab[sy2[1]];
  572.         /* first line 00rrrrrr.rr000ggg.ggggg000.bbbbbbbb: */
  573.         *(unsigned int *)(d1+BPP4) =
  574.             (CLIP8[y1 + bu])        |
  575.             (CLIP8[y1 + guv] << 11) |
  576.             (CLIP8[y1 + rv]  << 22) ;
  577.         /* second line 00rrrrrr.rr000ggg.ggggg000.bbbbbbbb: */
  578.         *(unsigned int *)(d2+BPP4) =
  579.             (CLIP8[y2 + bu])        |
  580.             (CLIP8[y2 + guv] << 11) |
  581.             (CLIP8[y2 + rv]  << 22) ;
  582.         /* next 2x2 block */
  583.         sy1 += 2; sy2 += 2;
  584.         su += 1; sv += 1;
  585.         d1 += 2*BPP4; d2 += 2*BPP4;
  586.     }
  587.     /* convert last 2x1 block: */
  588.     if (dx & 1) {
  589.         bu = pcd->butab[su[0]];
  590.         guv = pcd->gutab[su[0]] + pcd->gvtab[sv[0]];
  591.         rv = pcd->rvtab[sv[0]];
  592.         y1 = pcd->ytab[sy1[0]];
  593.         y2 = pcd->ytab[sy2[0]];
  594.         /* first line 00rrrrrr.rr000ggg.ggggg000.bbbbbbbb: */
  595.         *(unsigned int *)(d1+0) =
  596.             (CLIP8[y1 + bu])        |
  597.             (CLIP8[y1 + guv] << 11) |
  598.             (CLIP8[y1 + rv]  << 22) ;
  599.         /* second line BGR0 */
  600.         *(unsigned int *)(d2+0) =
  601.             (CLIP8[y2 + bu])        |
  602.             (CLIP8[y2 + guv] << 11) |
  603.             (CLIP8[y2 + rv]  << 22) ;
  604.     }
  605. }
  606. /*
  607.  * I420->RGB444 converter:
  608.  */
  609. int I420toRGB444 (unsigned char *dest_ptr, int dest_width, int dest_height,
  610.     int dest_pitch, int dest_x, int dest_y, int dest_dx, int dest_dy,
  611.     unsigned char *src_ptr, int src_width, int src_height, int src_pitch,
  612.     int src_x, int src_y, int src_dx, int src_dy,
  613.   color_data_t* pcd)
  614. {
  615.     /* scale factors: */
  616.     int scale_x, scale_y;
  617.     /* check arguments: */
  618.     if (!chk_args (dest_ptr, dest_width, dest_height, dest_pitch,
  619.    dest_x, dest_y, dest_dx, dest_dy, src_ptr, src_width, 
  620.    src_height, src_pitch, src_x, src_y, src_dx, src_dy, 
  621.    &scale_x, &scale_y))
  622.       return -1;
  623.     /* check if bottop-up images: */
  624.     if (dest_pitch < 0) dest_ptr -= (dest_height-1) * dest_pitch;
  625.     if (src_pitch <= 0) return -1;          /* not supported */
  626.     /* check if 1:1 scale: */
  627.     if (scale_x == 1 && scale_y == 1) {
  628.         /* color converter to use: */
  629.         void (*dbline) (unsigned char *d1, unsigned char *d2, int dest_x,
  630.             unsigned char *sy1, unsigned char *sy2, unsigned char *su, unsigned char *sv,
  631.             int src_x, int dx, color_data_t* pcd) = 
  632.             dblineI420toRGB444;
  633.         /* local variables: */
  634.         unsigned char *sy1, *sy2, *su, *sv, *d1, *d2;
  635.         register int j;
  636.         /* get pointers: */
  637.         sy1 = src_ptr + (src_x + src_y * src_pitch);        /* luma offset */
  638.         sy2 = sy1 + src_pitch;
  639.         su  = src_ptr + src_height * src_pitch + (src_x/2 + src_y/2 * src_pitch);
  640.         sv  = su + src_height * src_pitch / 4;
  641.         d1  = dest_ptr + dest_x * 2 + dest_y * dest_pitch; /* RGB offset */
  642.         d2  = d1 + dest_pitch;
  643.         /* convert a top line: */
  644.         if (src_y & 1) {                        /* chroma shift */
  645.             /* this is a bit inefficient, but who cares??? */
  646.             (* dbline) (d1, d1, dest_x, sy1, sy1, su, sv, src_x, src_dx, pcd);
  647.             sy1 += src_pitch;   sy2 += src_pitch;
  648.             su  += src_pitch/2; sv  += src_pitch/2;
  649.             d1  += dest_pitch;  d2  += dest_pitch;
  650.             dest_dy --;
  651.         }
  652.         /* convert aligned portion of the image: */
  653.         for (j = 0; j < dest_dy/2; j ++) {
  654.             /* convert two lines a time: */
  655.             (* dbline) (d1, d2, dest_x, sy1, sy2, su, sv, src_x, src_dx, pcd);
  656.             sy1 += src_pitch*2; sy2 += src_pitch*2;
  657.             su  += src_pitch/2; sv  += src_pitch/2;
  658.             d1  += dest_pitch*2; d2 += dest_pitch*2;
  659.         }
  660.         /* convert bottom line (if dest_dy is odd): */
  661.         if (dest_dy & 1) {
  662.             /* convert one line only: */
  663.             (* dbline) (d1, d1, dest_x, sy1, sy1, su, sv, src_x, src_dx, pcd);
  664.         }
  665.         return 0;
  666.     }
  667.     /* check if 2:1 scale: */
  668.     if (scale_x == 2 && scale_y == 2) {
  669.         /* color converter & interpolator to use: */
  670.         void (*cvt) (unsigned char *d1, unsigned char *d2, int dest_x,
  671.             unsigned char *sy1, unsigned char *sy2, unsigned char *su, unsigned char *sv,
  672.             int src_x, int dx, color_data_t* pcd) = dblineI420toXRGB;
  673.         void (*x2) (unsigned char *d1, unsigned char *d2, int dest_x,
  674.             unsigned char *s1, unsigned char *s2, int src_x, int dx, color_data_t* pcd) = dblineXRGBtoRGB444x2;
  675.         /* local variables: */
  676.         unsigned char *sy1, *sy2, *su, *sv, *d1, *d2;
  677.         register int dy = src_dy;
  678.         /* line buffers (we want them to be as compact as possible): */
  679.         int ibuf = 0;                           /* circular buffer index */
  680.         unsigned char * buf[3];                 /* actual pointers  */
  681.         buf[0] = pcd->linebuf;
  682.         buf[1] = pcd->linebuf + src_dx * BPP4;
  683.         buf[2] = pcd->linebuf + 2 * src_dx * BPP4;
  684.         /* get pointers: */
  685.         sy1 = src_ptr + (src_x + src_y * src_pitch);        /* luma offset */
  686.         sy2 = sy1 + src_pitch;
  687.         su  = src_ptr + src_height * src_pitch + (src_x/2 + src_y/2 * src_pitch);
  688.         sv  = su + src_height * src_pitch / 4;
  689.         d1  = dest_ptr + dest_x * BYTESPERPIXEL + dest_y * dest_pitch; /* RGB offset */
  690.         d2  = d1 + dest_pitch;
  691.         /* check if we have misaligned top line: */
  692.         if (src_y & 1) {
  693.             /* convert an odd first line: */
  694.             (*cvt) (buf[ibuf], buf[ibuf], 0, sy1, sy1, su, sv, src_x, src_dx, pcd);
  695.             sy1 += src_pitch;   sy2 += src_pitch;
  696.             su  += src_pitch/2; sv  += src_pitch/2;
  697.             dy --;
  698.         } else {
  699.             /* otherwise, convert first two lines: */
  700.             (*cvt) (buf[next[ibuf]], buf[next2[ibuf]], 0, sy1, sy2, su, sv, src_x, src_dx, pcd);
  701.             sy1 += src_pitch*2; sy2 += src_pitch*2;
  702.             su  += src_pitch/2; sv  += src_pitch/2;
  703.             ibuf = next[ibuf];      /* skip first interpolation: */
  704.             (*x2) (d1, d2, dest_x, buf[ibuf], buf[next[ibuf]], 0, src_dx, pcd);
  705.             d1  += dest_pitch*2; d2  += dest_pitch*2;
  706.             ibuf = next[ibuf];
  707.             dy -= 2;
  708.         }
  709.         /*
  710.          * Convert & interpolate the main portion of image:
  711.          *
  712.          *  source:      temp.store:                destination:
  713.          *
  714.          *               buf[ibuf] -------    /--> d1
  715.          *                                  x2  --> d2
  716.          *  s1 --   /-> buf[next[ibuf]] -<    /--> d1'=d1+2*pitch
  717.          *        cvt                       x2 ---> d2'=d2+2*pitch
  718.          *  s2 --/   -> buf[next2[ibuf]] /
  719.          */
  720.         while (dy >= 2) {
  721.             /* convert two lines into XRGB buffers: */
  722.             (*cvt) (buf[next[ibuf]], buf[next2[ibuf]], 0, sy1, sy2, su, sv, src_x, src_dx, pcd);
  723.             sy1 += src_pitch*2; sy2 += src_pitch*2;
  724.             su  += src_pitch/2; sv  += src_pitch/2;
  725.             /* interpolate first line: */
  726.             (*x2) (d1, d2, dest_x, buf[ibuf], buf[next[ibuf]], 0, src_dx, pcd);
  727.             d1  += dest_pitch*2; d2  += dest_pitch*2;
  728.             ibuf = next[ibuf];
  729.             /* interpolate second one: */
  730.             (*x2) (d1, d2, dest_x, buf[ibuf], buf[next[ibuf]], 0, src_dx, pcd);
  731.             d1  += dest_pitch*2; d2  += dest_pitch*2;
  732.             ibuf = next[ibuf];
  733.             dy -= 2;
  734.         }
  735.         /* check the # of remaining rows: */
  736.         if (dy & 1) {
  737.             /* convert the last odd line: */
  738.             (*cvt) (buf[next[ibuf]], buf[next[ibuf]], 0, sy1, sy1, su, sv, src_x, src_dx, pcd);
  739.             /* interpolate first line: */
  740.             (*x2) (d1, d2, dest_x, buf[ibuf], buf[next[ibuf]], 0, src_dx, pcd);
  741.             d1  += dest_pitch*2; d2  += dest_pitch*2;
  742.             ibuf = next[ibuf];
  743.             /* replicate the last line: */
  744.             (*x2) (d1, d2, dest_x, buf[ibuf], buf[ibuf], 0, src_dx, pcd);
  745.         } else {
  746.             /* replicate the last line: */
  747.             (*x2) (d1, d2, dest_x, buf[ibuf], buf[ibuf], 0, src_dx, pcd);
  748.         }
  749.         return 0;
  750.     }
  751.     /* check for 1/2 X scaling */
  752.     if (scale_x == 0 && scale_y == 0) {
  753.         /* local variables: */
  754.         unsigned char *sy1, *sy2, *su, *sv, *d1, *d2;
  755.         register int j;
  756.         /* get pointers: */
  757.         sy1 = src_ptr + (src_x + src_y * src_pitch);        /* luma offset */
  758.         sy2 = sy1 + src_pitch;
  759.         su  = src_ptr + src_height * src_pitch + (src_x/2 + src_y/2 * src_pitch);
  760.         sv  = su + src_height * src_pitch / 4;
  761.         d1  = dest_ptr + dest_x * 2 + dest_y * dest_pitch; /* RGB offset */
  762.         d2  = d1 + dest_pitch;
  763.         /* convert aligned portion of the image: */
  764.         for (j = 0; j < dest_dy>>1; j ++) {
  765.             /* convert two lines a time: */
  766.             halfI420toRGB444 (d1, dest_x, sy1, sy2, su, sv, src_x, src_dx, 0, pcd);
  767.             sy1 += src_pitch*2; sy2 += src_pitch*2;
  768.             su  += src_pitch/2; sv  += src_pitch/2;
  769.             d1  += dest_pitch;
  770.             halfI420toRGB444 (d1, dest_x, sy1, sy2, su, sv, src_x, src_dx, 1, pcd);
  771.             sy1 += src_pitch*2; sy2 += src_pitch*2;
  772.             su  += src_pitch/2; sv  += src_pitch/2;
  773.             d1  += dest_pitch;
  774.         }
  775.         return 0;
  776.     }
  777.     /* conversion is not supported */
  778.     return -1;
  779. }
  780. /*
  781.  * InitI420toRGB444 initializes the minimum set of clip tables necessary to perform
  782.  * color conversion to 12 bit RGB.  The 8-bit clip table is necessary for 2x
  783.  * interpolation.
  784.  *
  785.  * This routine also calls SetSrcI420Colors, which initializes color conversion tables
  786.  * and adds color balance.
  787.  */
  788. void InitI420toRGB444(float brightness, float contrast, float saturation, float hue,
  789.       color_data_t* pcd)
  790. {
  791.     InitColorData(pcd);
  792.     SetSrcI420Colors (brightness, contrast, saturation, hue, pcd);
  793. }
  794. void UninitI420toRGB444( ){
  795.     /* noop */
  796. }