4246f2ecdb0c43d408d94b76c84c065f78a0c90e
[goodguy/history.git] / cinelerra-5.1 / guicast / clip.h
1
2 /*
3  * CINELERRA
4  * Copyright (C) 2008 Adam Williams <broadcast at earthling dot net>
5  *
6  * This program is free software; you can redistribute it and/or modify
7  * it under the terms of the GNU General Public License as published by
8  * the Free Software Foundation; either version 2 of the License, or
9  * (at your option) any later version.
10  *
11  * This program is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14  * GNU General Public License for more details.
15  *
16  * You should have received a copy of the GNU General Public License
17  * along with this program; if not, write to the Free Software
18  * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
19  *
20  */
21
22 #ifndef __CLIP_H__
23 #define __CLIP_H__
24
25 // Math macros
26 #undef SQR
27 #define SQR(x) ((x) * (x))
28 #undef CLIP
29 #define CLIP(x, y, z) ((x) < (y) ? (y) : ((x) > (z) ? (z) : (x)))
30 #undef RECLIP
31 #define RECLIP(x, y, z) ((x) = ((x) < (y) ? (y) : ((x) > (z) ? (z) : (x))))
32 #undef CLAMP
33 #define CLAMP(x, y, z) ((x) = ((x) < (y) ? (y) : ((x) > (z) ? (z) : (x))))
34 #undef MAX
35 #define MAX(x, y) ((x) > (y) ? (x) : (y))
36 #undef MIN
37 #define MIN(x, y) ((x) < (y) ? (x) : (y))
38 #undef EQUIV
39 #define EQUIV(x, y) (fabs((x) - (y)) < 0.001)
40 #undef DISTANCE
41 #define DISTANCE(x1, y1, x2, y2) \
42 (sqrt(((x2) - (x1)) * ((x2) - (x1)) + ((y2) - (y1)) * ((y2) - (y1))))
43 #define TO_RAD(x) ((x) * 2 * M_PI / 360)
44 #define TO_DEG(x) ((x) * 360 / 2 / M_PI)
45
46 static inline int bclip(int &iv, int imn, int imx) {
47         return iv < imn ? imn : iv > imx ? imx : iv;
48 }
49 static inline float bclip(float &fv, float fmn, float fmx) {
50         return fv < fmn ? fmn : fv > fmx ? fmx : fv;
51 }
52 static inline double bclip(double &dv, double dmn, double dmx) {
53         return dv < dmn ? dmn : dv > dmx ? dmx : dv;
54 }
55 static inline void bclamp(int &iv, int imn, int imx) {
56         if( iv < imn ) iv = imn; else if( iv > imx ) iv = imx;
57 }
58 static inline void bclamp(float &fv, float fmn, float fmx) {
59         if( fv < fmn ) fv = fmn; else if( fv > fmx ) fv = fmx;
60 }
61 static inline void bclamp(double &dv, double dmn, double dmx) {
62         if( dv < dmn ) dv = dmn; else if( dv > dmx ) dv = dmx;
63 }
64 static inline void bc_rgb2yuv(int r, int g, int b, int &y, int &u, int &v, int max=255)
65 { //bt601, jpeg, unclipped
66         double mx = max, rr = r/mx, gg = g/mx, bb = b/mx;
67         y = (int)(( 0.29900*rr + 0.58700*gg + 0.11400*bb)       * mx + 0.5);
68         u = (int)((-0.16874*rr - 0.33126*gg + 0.50000*bb + 0.5) * mx + 0.5);
69         v = (int)(( 0.50000*rr - 0.41869*gg - 0.08131*bb + 0.5) * mx + 0.5);
70 }
71 static inline void bc_yuv2rgb(int y, int u, int v, int &r, int &g, int &b, int max=255)
72 {
73         int ofs = (max + 1) / 2;
74         double mx = max, yy = y/mx, uu = (u-ofs)/mx, vv = (v-ofs)/mx;
75         r = (int)((yy + 1.40200*vv)              * mx + 0.5);
76         g = (int)((yy - 0.34414*uu - 0.71414*vv) * mx + 0.5);
77         b = (int)((yy + 1.77200*uu)              * mx + 0.5);
78 }
79
80 #endif