lut.c
上传用户:weiyuanprp
上传日期:2020-05-20
资源大小:1169k
文件大小:3k
源码类别:

传真(Fax)编程

开发平台:

C/C++

  1. /* $Id: lut.c,v 1.1.1.1 2005/11/11 21:32:03 faxguy Exp $ */
  2. /*
  3.  * Copyright (c) 1990-1996 Sam Leffler
  4.  * Copyright (c) 1991-1996 Silicon Graphics, Inc.
  5.  * HylaFAX is a trademark of Silicon Graphics
  6.  *
  7.  * Permission to use, copy, modify, distribute, and sell this software and 
  8.  * its documentation for any purpose is hereby granted without fee, provided
  9.  * that (i) the above copyright notices and this permission notice appear in
  10.  * all copies of the software and related documentation, and (ii) the names of
  11.  * Sam Leffler and Silicon Graphics may not be used in any advertising or
  12.  * publicity relating to the software without the specific, prior written
  13.  * permission of Sam Leffler and Silicon Graphics.
  14.  * 
  15.  * THE SOFTWARE IS PROVIDED "AS-IS" AND WITHOUT WARRANTY OF ANY KIND, 
  16.  * EXPRESS, IMPLIED OR OTHERWISE, INCLUDING WITHOUT LIMITATION, ANY 
  17.  * WARRANTY OF MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE.  
  18.  * 
  19.  * IN NO EVENT SHALL SAM LEFFLER OR SILICON GRAPHICS BE LIABLE FOR
  20.  * ANY SPECIAL, INCIDENTAL, INDIRECT OR CONSEQUENTIAL DAMAGES OF ANY KIND,
  21.  * OR ANY DAMAGES WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR PROFITS,
  22.  * WHETHER OR NOT ADVISED OF THE POSSIBILITY OF DAMAGE, AND ON ANY THEORY OF 
  23.  * LIABILITY, ARISING OUT OF OR IN CONNECTION WITH THE USE OR PERFORMANCE 
  24.  * OF THIS SOFTWARE.
  25.  */
  26. #include "lut.h"
  27. float frand();
  28. lut *makelut(func,insteps,outsteps,stoclut)
  29. float (*func)();
  30. int insteps, outsteps, stoclut;
  31. {
  32.     lut *l;
  33.     int i;
  34.     float low, high, inspace;
  35.     l = (lut *)malloc(sizeof(lut));
  36.     l->insteps = insteps;
  37.     l->outsteps = outsteps;
  38.     l->stoclut = stoclut;
  39.     if(stoclut) {
  40.         l->flow = (float *)malloc(insteps*sizeof(float));
  41.         l->fhigh = (float *)malloc(insteps*sizeof(float));
  42.         inspace = insteps-1.0;
  43.         for(i=0; i<insteps; i++) {
  44.             low = (i-0.5)/inspace;
  45.             high = (i+0.5)/inspace;
  46.             if(low<0.0)
  47.                 low = 0.0;
  48.             if(high>1.0)
  49.                 high = 1.0;
  50.             l->flow[i] = (func)(low);
  51.             l->fhigh[i] = (func)(high);
  52.         }
  53.     } else { 
  54.         l->stab = (unsigned short *)malloc(insteps*sizeof(unsigned short));
  55.         inspace = insteps-1.0;
  56.         for(i=0; i<insteps; i++) 
  57.             l->stab[i] = (l->outsteps-1)*(func)(i/inspace)+0.5;
  58.     }
  59.     return l;
  60. }
  61. applylut(l,sptr,n)
  62. lut *l;
  63. unsigned short *sptr;
  64. int n;
  65. {
  66.     float delta, val;
  67.     float *fhigh, *flow;
  68.     unsigned short *stab;
  69.     float outspace;
  70.     int ival;
  71.     if(l->stoclut) {
  72.         fhigh = l->fhigh;
  73.         flow = l->flow;
  74.         outspace = l->outsteps-1;
  75.         while(n--) {
  76.             delta = fhigh[*sptr]-flow[*sptr];
  77.             val = outspace*(flow[*sptr]+frand()*delta);
  78.             ival = val;
  79.             if((val-ival)<frand()) 
  80.                 *sptr++ = ival;
  81.             else
  82.                 *sptr++ = ival+1;
  83.         }
  84.     } else {
  85.         stab = l->stab;
  86.         while(n--) {
  87.             *sptr = stab[*sptr];
  88.             sptr++;
  89.         }
  90.     }
  91. }
  92. void freelut(lut* l)
  93. {
  94.     if(l->stoclut) {
  95.         free((char*)l->flow);
  96.         free((char*)l->fhigh);
  97.     } else { 
  98.         free((char*)l->stab);
  99.     }
  100.     free((char*)l);
  101. }