Mercurial > hg > wm

comparison Meerwald/cmp_dugad_sig.c @ 8:f83ef905a63d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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

Repositories maintained by Peter Meerwald, pmeerw@pmeerw.net.