src/common/tmprivat.h

/* RCSid: $Id$ */
/*
 * Private header file for tone mapping routines.
 */

/* ====================================================================
 * The Radiance Software License, Version 1.0
 *
 * Copyright (c) 1990 - 2002 The Regents of the University of California,
 * through Lawrence Berkeley National Laboratory.   All rights reserved.
 *
 * Redistribution and use in source and binary forms, with or without
 * modification, are permitted provided that the following conditions
 * are met:
 *
 * 1. Redistributions of source code must retain the above copyright
 *         notice, this list of conditions and the following disclaimer.
 *
 * 2. Redistributions in binary form must reproduce the above copyright
 *       notice, this list of conditions and the following disclaimer in
 *       the documentation and/or other materials provided with the
 *       distribution.
 *
 * 3. The end-user documentation included with the redistribution,
 *           if any, must include the following acknowledgment:
 *             "This product includes Radiance software
 *                 (http://radsite.lbl.gov/)
 *                 developed by the Lawrence Berkeley National Laboratory
 *               (http://www.lbl.gov/)."
 *       Alternately, this acknowledgment may appear in the software itself,
 *       if and wherever such third-party acknowledgments normally appear.
 *
 * 4. The names "Radiance," "Lawrence Berkeley National Laboratory"
 *       and "The Regents of the University of California" must
 *       not be used to endorse or promote products derived from this
 *       software without prior written permission. For written
 *       permission, please contact [email protected].
 *
 * 5. Products derived from this software may not be called "Radiance",
 *       nor may "Radiance" appear in their name, without prior written
 *       permission of Lawrence Berkeley National Laboratory.
 *
 * THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESSED OR IMPLIED
 * WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
 * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
 * DISCLAIMED.   IN NO EVENT SHALL Lawrence Berkeley National Laboratory OR
 * ITS CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
 * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
 * LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF
 * USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
 * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
 * OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT
 * OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
 * SUCH DAMAGE.
 * ====================================================================
 *
 * This software consists of voluntary contributions made by many
 * individuals on behalf of Lawrence Berkeley National Laboratory.   For more
 * information on Lawrence Berkeley National Laboratory, please see
 * <http://www.lbl.gov/>.
 */

#ifndef MEM_PTR
#define MEM_PTR         void *
#endif
#include        "color.h"
#include        "tonemap.h"

#ifdef __cplusplus
extern "C" {
#endif

                                /* required constants */
#ifndef M_LN2
#define M_LN2           0.69314718055994530942
#endif
#ifndef M_LN10
#define M_LN10          2.30258509299404568402
#endif
                                /* minimum values and defaults */
#define MINGAM          0.75
#define DEFGAM          2.2
#define MINLDDYN        2.
#define DEFLDDYN        32.
#define MINLDMAX        1.
#define DEFLDMAX        100.

#define BRT2SCALE(l2)   (int)(M_LN2*TM_BRTSCALE*(l2) + ((l2)>0 ? .5 : -.5))

#define HISTEP          8               /* steps in BRTSCALE for each bin */

#define MINBRT          (-16*TM_BRTSCALE)       /* minimum usable brightness */
#define MINLUM          (1.125352e-7)           /* tmLuminance(MINBRT) */

#define LMESLOWER       (5.62e-3)               /* lower mesopic limit */
#define LMESUPPER       (5.62)                  /* upper mesopic limit */
#if     (TM_BRTSCALE==128)
#define BMESLOWER       (-663)                  /* encoded LMESLOWER */
#define BMESUPPER       (221)                   /* encoded LMESUPPER */
#else
#define BMESLOWER       ((int)(-5.18*TM_BRTSCALE-.5))
#define BMESUPPER       ((int)(1.73*TM_BRTSCALE+.5))
#endif

#ifndef int4
#define int4            int                     /* assume 32-bit integers */
#endif
                                                /* approximate scotopic lum. */
#define SCO_rf          0.062
#define SCO_gf          0.608
#define SCO_bf          0.330
#define scotlum(c)      (SCO_rf*(c)[RED] + SCO_gf*(c)[GRN] + SCO_bf*(c)[BLU])
#define normscot(c)     ( (     (int4)(SCO_rf*256.+.5)*(c)[RED] + \
                                (int4)(SCO_gf*256.+.5)*(c)[GRN] + \
                                (int4)(SCO_bf*256.+.5)*(c)[BLU] ) >> 8 )

#ifndef malloc
MEM_PTR         malloc();
MEM_PTR         calloc();
#endif
extern int      tmErrorReturn();

                                                /* lookup for mesopic scaling */
extern BYTE     tmMesofact[BMESUPPER-BMESLOWER];

extern void     tmMkMesofact();                 /* build tmMesofact */

#define returnErr(code) return(tmErrorReturn(funcName,code))
#define returnOK        return(TM_E_OK)

#define FEQ(a,b)        ((a) < (b)+1e-5 && (b) < (a)+1e-5)

#define PRIMEQ(p1,p2)   (FEQ((p1)[0][0],(p2)[0][0])&&FEQ((p1)[0][1],(p2)[0][1])\
                        &&FEQ((p1)[1][0],(p2)[1][0])&&FEQ((p1)[1][1],(p2)[1][1])\
                        &&FEQ((p1)[2][0],(p2)[2][0])&&FEQ((p1)[2][1],(p2)[2][1])\
                        &&FEQ((p1)[3][0],(p2)[3][0])&&FEQ((p1)[3][1],(p2)[3][1]))

#ifdef __cplusplus
}
#endif
Revision:	3.6
Committed:	Sat Feb 22 02:07:22 2003 UTC (22 years, 8 months ago) by greg
Content type:	text/plain
Branch:	MAIN
Changes since 3.5:	+77 -12 lines
Log Message:	Changes and check-in for 3.5 release Includes new source files and modifications not recorded for many years See ray/doc/notes/ReleaseNotes for notes between 3.1 and 3.5 release
#	User	Rev	Content
1	greg	3.6	/* RCSid: $Id$ */
2	greg	3.1	/*
3			* Private header file for tone mapping routines.
4			*/
5
6	greg	3.6	/* ====================================================================
7			* The Radiance Software License, Version 1.0
8			*
9			* Copyright (c) 1990 - 2002 The Regents of the University of California,
10			* through Lawrence Berkeley National Laboratory. All rights reserved.
11			*
12			* Redistribution and use in source and binary forms, with or without
13			* modification, are permitted provided that the following conditions
14			* are met:
15			*
16			* 1. Redistributions of source code must retain the above copyright
17			* notice, this list of conditions and the following disclaimer.
18			*
19			* 2. Redistributions in binary form must reproduce the above copyright
20			* notice, this list of conditions and the following disclaimer in
21			* the documentation and/or other materials provided with the
22			* distribution.
23			*
24			* 3. The end-user documentation included with the redistribution,
25			* if any, must include the following acknowledgment:
26			* "This product includes Radiance software
27			* (http://radsite.lbl.gov/)
28			* developed by the Lawrence Berkeley National Laboratory
29			* (http://www.lbl.gov/)."
30			* Alternately, this acknowledgment may appear in the software itself,
31			* if and wherever such third-party acknowledgments normally appear.
32			*
33			* 4. The names "Radiance," "Lawrence Berkeley National Laboratory"
34			* and "The Regents of the University of California" must
35			* not be used to endorse or promote products derived from this
36			* software without prior written permission. For written
37			* permission, please contact [email protected].
38			*
39			* 5. Products derived from this software may not be called "Radiance",
40			* nor may "Radiance" appear in their name, without prior written
41			* permission of Lawrence Berkeley National Laboratory.
42			*
43			* THIS SOFTWARE IS PROVIDED ``AS IS'' AND ANY EXPRESSED OR IMPLIED
44			* WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
45			* OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
46			* DISCLAIMED. IN NO EVENT SHALL Lawrence Berkeley National Laboratory OR
47			* ITS CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
48			* SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
49			* LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF
50			* USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
51			* ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
52			* OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT
53			* OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
54			* SUCH DAMAGE.
55			* ====================================================================
56			*
57			* This software consists of voluntary contributions made by many
58			* individuals on behalf of Lawrence Berkeley National Laboratory. For more
59			* information on Lawrence Berkeley National Laboratory, please see
60			* <http://www.lbl.gov/>.
61			*/
62
63	greg	3.4	#ifndef MEM_PTR
64	greg	3.6	#define MEM_PTR void *
65	greg	3.4	#endif
66	gwlarson	3.5	#include "color.h"
67	greg	3.1	#include "tonemap.h"
68	greg	3.2
69	greg	3.6	#ifdef __cplusplus
70			extern "C" {
71			#endif
72
73	greg	3.2	/* required constants */
74			#ifndef M_LN2
75			#define M_LN2 0.69314718055994530942
76			#endif
77			#ifndef M_LN10
78			#define M_LN10 2.30258509299404568402
79			#endif
80	greg	3.1	/* minimum values and defaults */
81			#define MINGAM 0.75
82			#define DEFGAM 2.2
83			#define MINLDDYN 2.
84			#define DEFLDDYN 32.
85			#define MINLDMAX 1.
86			#define DEFLDMAX 100.
87
88	greg	3.6	#define BRT2SCALE(l2) (int)(M_LN2TM_BRTSCALE(l2) + ((l2)>0 ? .5 : -.5))
89	greg	3.1
90			#define HISTEP 8 /* steps in BRTSCALE for each bin */
91
92			#define MINBRT (-16TM_BRTSCALE) / minimum usable brightness */
93			#define MINLUM (1.125352e-7) /* tmLuminance(MINBRT) */
94
95			#define LMESLOWER (5.62e-3) /* lower mesopic limit */
96	greg	3.3	#define LMESUPPER (5.62) /* upper mesopic limit */
97			#if (TM_BRTSCALE==128)
98	greg	3.2	#define BMESLOWER (-663) /* encoded LMESLOWER */
99			#define BMESUPPER (221) /* encoded LMESUPPER */
100	greg	3.3	#else
101			#define BMESLOWER ((int)(-5.18*TM_BRTSCALE-.5))
102			#define BMESUPPER ((int)(1.73*TM_BRTSCALE+.5))
103			#endif
104	greg	3.1
105			#ifndef int4
106			#define int4 int /* assume 32-bit integers */
107			#endif
108			/* approximate scotopic lum. */
109			#define SCO_rf 0.062
110			#define SCO_gf 0.608
111			#define SCO_bf 0.330
112			#define scotlum(c) (SCO_rf(c)[RED] + SCO_gf(c)[GRN] + SCO_bf*(c)[BLU])
113			#define normscot(c) ( ( (int4)(SCO_rf256.+.5)(c)[RED] + \
114			(int4)(SCO_gf256.+.5)(c)[GRN] + \
115			(int4)(SCO_bf256.+.5)(c)[BLU] ) >> 8 )
116
117			#ifndef malloc
118	greg	3.4	MEM_PTR malloc();
119			MEM_PTR calloc();
120	greg	3.1	#endif
121			extern int tmErrorReturn();
122
123	greg	3.6	/* lookup for mesopic scaling */
124			extern BYTE tmMesofact[BMESUPPER-BMESLOWER];
125
126			extern void tmMkMesofact(); /* build tmMesofact */
127
128	greg	3.1	#define returnErr(code) return(tmErrorReturn(funcName,code))
129			#define returnOK return(TM_E_OK)
130
131			#define FEQ(a,b) ((a) < (b)+1e-5 && (b) < (a)+1e-5)
132
133	greg	3.6	#define PRIMEQ(p1,p2) (FEQ((p1)[0][0],(p2)[0][0])&&FEQ((p1)[0][1],(p2)[0][1])\
134			&&FEQ((p1)[1][0],(p2)[1][0])&&FEQ((p1)[1][1],(p2)[1][1])\
135			&&FEQ((p1)[2][0],(p2)[2][0])&&FEQ((p1)[2][1],(p2)[2][1])\
136			&&FEQ((p1)[3][0],(p2)[3][0])&&FEQ((p1)[3][1],(p2)[3][1]))
137
138			#ifdef __cplusplus
139			}
140			#endif