Mercurial > hg > wm
comparison Meerwald/cmp_dugad_sig.c @ 8:f83ef905a63d
fixing many warnings
increase capacity for coordinates in bruyn
fix some uninit. variables
author | Peter Meerwald <pmeerw@cosy.sbg.ac.at> |
---|---|
date | Tue, 22 Apr 2008 13:36:05 +0200 |
parents | be303a3f5ea8 |
children |
comparison
equal
deleted
inserted
replaced
7:2b350281f8b0 | 8:f83ef905a63d |
---|---|
13 | 13 |
14 int main(int argc, char *argv[]) { | 14 int main(int argc, char *argv[]) { |
15 | 15 |
16 FILE *in = stdin; | 16 FILE *in = stdin; |
17 FILE *out = stdout; | 17 FILE *out = stdout; |
18 FILE *sig = NULL; | |
19 | 18 |
20 char output_name[MAXPATHLEN] = "(stdout)"; | 19 char output_name[MAXPATHLEN] = "(stdout)"; |
21 char input_name[MAXPATHLEN] = "(stdin)"; | 20 char input_name[MAXPATHLEN] = "(stdin)"; |
22 | 21 |
23 int c, i, n, ok; | 22 int c, i, n, ok; |
26 double diff; | 25 double diff; |
27 char line[32]; | 26 char line[32]; |
28 | 27 |
29 int correlation_only = 0; | 28 int correlation_only = 0; |
30 int verbose = 0; | 29 int verbose = 0; |
31 double corr; | |
32 | 30 |
33 progname = argv[0]; | 31 progname = argv[0]; |
34 | 32 |
35 while ((c = getopt(argc, argv, "h?Co:v:s:")) != EOF) { | 33 while ((c = getopt(argc, argv, "h?Co:v:s:")) != EOF) { |
36 switch (c) { | 34 switch (c) { |
66 if (argc > 1) { | 64 if (argc > 1) { |
67 usage(); | 65 usage(); |
68 exit(1); | 66 exit(1); |
69 } | 67 } |
70 | 68 |
71 if (argc == 1 && *argv[0] != '-') | 69 if (argc == 1 && *argv[0] != '-') { |
72 if ((in = fopen(argv[0], "r")) == NULL) { | 70 if ((in = fopen(argv[0], "r")) == NULL) { |
73 fprintf(stderr, "%s: unable to open input file %s\n", progname, argv[0]); | 71 fprintf(stderr, "%s: unable to open input file %s\n", progname, argv[0]); |
74 exit(1); | 72 exit(1); |
75 } | 73 } |
76 else | 74 else |
77 strcpy(input_name, argv[0]); | 75 strcpy(input_name, argv[0]); |
78 | 76 } |
77 | |
79 fgets(line, sizeof(line), in); | 78 fgets(line, sizeof(line), in); |
80 if (strspn(line, "DGWM") < 4) { | 79 if (strspn(line, "DGWM") < 4) { |
81 fprintf(stderr, "%s: watermark file %s invalid\n", progname, input_name); | 80 fprintf(stderr, "%s: watermark file %s invalid\n", progname, input_name); |
82 exit(1); | 81 exit(1); |
83 } | 82 } |
92 int m; | 91 int m; |
93 double z, v; | 92 double z, v; |
94 | 93 |
95 // HL subband | 94 // HL subband |
96 fscanf(in, "%d %lf %lf\n", &m, &z, &v); | 95 fscanf(in, "%d %lf %lf\n", &m, &z, &v); |
97 if (verbose && !correlation_only) | 96 if (verbose && !correlation_only) { |
98 if (m) | 97 if (m) |
99 fprintf(out, "%f %f\n", z / (double) m, (v * alpha) / (double) (1.0 * m)); | 98 fprintf(out, "%f %f\n", z / (double) m, (v * alpha) / (double) (1.0 * m)); |
100 else | 99 else |
101 fprintf(out, "0.0 0.0\n"); | 100 fprintf(out, "0.0 0.0\n"); |
101 } | |
102 if (m) { | 102 if (m) { |
103 ok += (z > v * alpha / (double) 1.0) ? 1 : 0; | 103 ok += (z > v * alpha / (double) 1.0) ? 1 : 0; |
104 diff += ((z - v * alpha) / (double) (1.0 * m)); | 104 diff += ((z - v * alpha) / (double) (1.0 * m)); |
105 } | 105 } |
106 else | 106 else |
107 n--; | 107 n--; |
108 | 108 |
109 // LH subband | 109 // LH subband |
110 fscanf(in, "%d %lf %lf\n", &m, &z, &v); | 110 fscanf(in, "%d %lf %lf\n", &m, &z, &v); |
111 if (verbose && !correlation_only) | 111 if (verbose && !correlation_only) { |
112 if (m) | 112 if (m) |
113 fprintf(out, "%f %f\n", z / (double) m, (v * alpha) / (double) (1.0 * m)); | 113 fprintf(out, "%f %f\n", z / (double) m, (v * alpha) / (double) (1.0 * m)); |
114 else | 114 else |
115 fprintf(out, "0.0 0.0\n"); | 115 fprintf(out, "0.0 0.0\n"); |
116 } | |
116 if (m) { | 117 if (m) { |
117 ok += (z > v * alpha / (double) 1.0) ? 1 : 0; | 118 ok += (z > v * alpha / (double) 1.0) ? 1 : 0; |
118 diff += ((z - v * alpha) / (double) (1.0 * m)); | 119 diff += ((z - v * alpha) / (double) (1.0 * m)); |
119 } | 120 } |
120 else | 121 else |
121 n--; | 122 n--; |
122 | 123 |
123 // HH subband | 124 // HH subband |
124 fscanf(in, "%d %lf %lf\n", &m, &z, &v); | 125 fscanf(in, "%d %lf %lf\n", &m, &z, &v); |
125 if (verbose && !correlation_only) | 126 if (verbose && !correlation_only) { |
126 if (m) | 127 if (m) |
127 fprintf(out, "%f %f\n", z / (double) m, (v * alpha) / (double) (1.0 * m)); | 128 fprintf(out, "%f %f\n", z / (double) m, (v * alpha) / (double) (1.0 * m)); |
128 else | 129 else |
129 fprintf(out, "0.0 0.0\n"); | 130 fprintf(out, "0.0 0.0\n"); |
130 | 131 } |
132 | |
131 if (m) { | 133 if (m) { |
132 ok += (z > v * alpha / (double) 1.0) ? 1 : 0; | 134 ok += (z > v * alpha / (double) 1.0) ? 1 : 0; |
133 diff += ((z - v * alpha) / (double) (1.0 * m)); | 135 diff += ((z - v * alpha) / (double) (1.0 * m)); |
134 } | 136 } |
135 else | 137 else |