1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
|
Index: netpbm-10.29/editor/pnmaddnoise.c
===================================================================
--- netpbm-10.29.orig/editor/pnmaddnoise.c
+++ netpbm-10.29/editor/pnmaddnoise.c
@@ -252,14 +252,14 @@ laplacian_noise(FILE * const ifp,
int rows, cols, format, row;
double rtemp, gtemp, btemp;
double u, u1;
- double INFINITY;
+ double infinity;
double const EPSILON = 1.0e-5;
int const forceplain = 0;
pnm_readpnminit( ifp, &cols, &rows, &maxval, &format );
xelrow = pnm_allocrow( cols );
- INFINITY = (double) maxval;
+ infinity = (double) maxval;
pnm_writepnminit( stdout, cols, rows, maxval, format, forceplain );
@@ -272,12 +272,12 @@ laplacian_noise(FILE * const ifp,
u = ( rand( ) & RANDOM_MASK ) / arand;
if ( u <= 0.5 ) {
if ( u <= EPSILON )
- rtemp = PPM_GETR( *xP ) - INFINITY;
+ rtemp = PPM_GETR( *xP ) - infinity;
else
rtemp = PPM_GETR( *xP ) + lsigma * log( 2.0 * u );
} else {
if ( ( u1 = 1.0 - u ) <= 0.5 * EPSILON )
- rtemp = PPM_GETR( *xP ) + INFINITY;
+ rtemp = PPM_GETR( *xP ) + infinity;
else
rtemp = PPM_GETR( *xP ) - lsigma * log( 2.0 * u1 );
}
@@ -285,12 +285,12 @@ laplacian_noise(FILE * const ifp,
u = ( rand( ) & RANDOM_MASK ) / arand;
if ( u <= 0.5 ) {
if ( u <= EPSILON )
- gtemp = PPM_GETG( *xP ) - INFINITY;
+ gtemp = PPM_GETG( *xP ) - infinity;
else
gtemp = PPM_GETG( *xP ) + lsigma * log( 2.0 * u );
} else {
if ( ( u1 = 1.0 - u ) <= 0.5 * EPSILON )
- gtemp = PPM_GETG( *xP ) + INFINITY;
+ gtemp = PPM_GETG( *xP ) + infinity;
else
gtemp = PPM_GETG( *xP ) - lsigma * log( 2.0 * u1 );
}
@@ -298,12 +298,12 @@ laplacian_noise(FILE * const ifp,
u = ( rand( ) & RANDOM_MASK ) / arand;
if ( u <= 0.5 ) {
if ( u <= EPSILON )
- btemp = PPM_GETB( *xP ) - INFINITY;
+ btemp = PPM_GETB( *xP ) - infinity;
else
btemp = PPM_GETB( *xP ) + lsigma * log( 2.0 * u );
} else {
if ( ( u1 = 1.0 - u ) <= 0.5 * EPSILON )
- btemp = PPM_GETB( *xP ) + INFINITY;
+ btemp = PPM_GETB( *xP ) + infinity;
else
btemp = PPM_GETB( *xP ) - lsigma * log( 2.0 * u1 );
}
@@ -318,12 +318,12 @@ laplacian_noise(FILE * const ifp,
u = ( rand( ) & RANDOM_MASK ) / arand;
if ( u <= 0.5 ) {
if ( u <= EPSILON )
- gtemp = PNM_GET1( *xP ) - INFINITY;
+ gtemp = PNM_GET1( *xP ) - infinity;
else
gtemp = PNM_GET1( *xP ) + lsigma * log( 2.0 * u );
} else {
if ( ( u1 = 1.0 - u ) <= 0.5 * EPSILON )
- gtemp = PNM_GET1( *xP ) + INFINITY;
+ gtemp = PNM_GET1( *xP ) + infinity;
else
gtemp = PNM_GET1( *xP ) - lsigma * log( 2.0 * u1 );
}
@@ -356,14 +356,14 @@ multiplicative_gaussian_noise(FILE * con
int rows, cols, format, row;
double rtemp, gtemp, btemp;
double uniform, rayleigh, gauss;
- double INFINITY;
+ double infinity;
double const EPSILON = 1.0e-5;
int const forceplain = 0;
pnm_readpnminit( ifp, &cols, &rows, &maxval, &format );
xelrow = pnm_allocrow( cols );
- INFINITY = (double) maxval;
+ infinity = (double) maxval;
pnm_writepnminit( stdout, cols, rows, maxval, format, forceplain );
@@ -375,7 +375,7 @@ multiplicative_gaussian_noise(FILE * con
case PPM_TYPE:
uniform = ( rand( ) & RANDOM_MASK ) / arand;
if ( uniform <= EPSILON )
- rayleigh = INFINITY;
+ rayleigh = infinity;
else
rayleigh = sqrt( -2.0 * log( uniform ) );
uniform = ( rand( ) & RANDOM_MASK ) / arand;
@@ -384,7 +384,7 @@ multiplicative_gaussian_noise(FILE * con
uniform = ( rand( ) & RANDOM_MASK ) / arand;
if ( uniform <= EPSILON )
- rayleigh = INFINITY;
+ rayleigh = infinity;
else
rayleigh = sqrt( -2.0 * log( uniform ) );
uniform = ( rand( ) & RANDOM_MASK ) / arand;
@@ -393,7 +393,7 @@ multiplicative_gaussian_noise(FILE * con
uniform = ( rand( ) & RANDOM_MASK ) / arand;
if ( uniform <= EPSILON )
- rayleigh = INFINITY;
+ rayleigh = infinity;
else
rayleigh = sqrt( -2.0 * log( uniform ) );
uniform = ( rand( ) & RANDOM_MASK ) / arand;
@@ -409,7 +409,7 @@ multiplicative_gaussian_noise(FILE * con
default:
uniform = ( rand( ) & RANDOM_MASK ) / arand;
if ( uniform <= EPSILON )
- rayleigh = INFINITY;
+ rayleigh = infinity;
else
rayleigh = sqrt( -2.0 * log( uniform ) );
uniform = ( rand( ) & RANDOM_MASK ) / arand;
|
GitWeb