kopia lustrzana https://github.com/Guenael/rtlsdr-wsprd
Guenael, VA2GKA
GitHub
commit
3d456a562d
|
|
@ -324,7 +324,7 @@ void postSpots(uint32_t n_results) {
|
||||||
// "Table 'wsprnet_db.activity' doesn't exist" reported on web site...
|
// "Table 'wsprnet_db.activity' doesn't exist" reported on web site...
|
||||||
// Anyone has doc about this?
|
// Anyone has doc about this?
|
||||||
if (n_results == 0) {
|
if (n_results == 0) {
|
||||||
snprintf(url, sizeof(url) - 1, "http://wsprnet.org/post?function=wsprstat&rcall=%s&rgrid=%s&rqrg=%.6f&tpct=%.2f&tqrg=%.6f&dbm=%d&version=rtlsdr-wsprd_v0.4.2&mode=2",
|
snprintf(url, sizeof(url) - 1, "http://wsprnet.org/post?function=wsprstat&rcall=%s&rgrid=%s&rqrg=%.6f&tpct=%.2f&tqrg=%.6f&dbm=%d&version=rtlsdr-050&mode=2",
|
||||||
dec_options.rcall,
|
dec_options.rcall,
|
||||||
dec_options.rloc,
|
dec_options.rloc,
|
||||||
rx_options.realfreq / 1e6,
|
rx_options.realfreq / 1e6,
|
||||||
|
|
@ -347,7 +347,7 @@ void postSpots(uint32_t n_results) {
|
||||||
}
|
}
|
||||||
|
|
||||||
for (uint32_t i = 0; i < n_results; i++) {
|
for (uint32_t i = 0; i < n_results; i++) {
|
||||||
snprintf(url, sizeof(url) - 1, "http://wsprnet.org/post?function=wspr&rcall=%s&rgrid=%s&rqrg=%.6f&date=%02d%02d%02d&time=%02d%02d&sig=%.0f&dt=%.1f&tqrg=%.6f&tcall=%s&tgrid=%s&dbm=%s&version=rtlsdr-wsprd_v0.4.2&mode=2",
|
snprintf(url, sizeof(url) - 1, "http://wsprnet.org/post?function=wspr&rcall=%s&rgrid=%s&rqrg=%.6f&date=%02d%02d%02d&time=%02d%02d&sig=%.0f&dt=%.1f&tqrg=%.6f&tcall=%s&tgrid=%s&dbm=%s&version=rtlsdr-050&mode=2",
|
||||||
dec_options.rcall,
|
dec_options.rcall,
|
||||||
dec_options.rloc,
|
dec_options.rloc,
|
||||||
dec_results[i].freq,
|
dec_results[i].freq,
|
||||||
|
|
@ -612,10 +612,10 @@ int32_t decoderSelfTest() {
|
||||||
unsigned char symbols[162];
|
unsigned char symbols[162];
|
||||||
char message[] = "K1JT FN20QI 20";
|
char message[] = "K1JT FN20QI 20";
|
||||||
char hashtab[32768*13] = {0};
|
char hashtab[32768*13] = {0};
|
||||||
//char loctab[32768*5] = {0}; // EVAL: code update from wsprd
|
char loctab[32768*5] = {0}; // EVAL: code update from wsprd
|
||||||
|
|
||||||
// Compute sympbols from the message
|
// Compute sympbols from the message
|
||||||
get_wspr_channel_symbols(message, hashtab, symbols);
|
get_wspr_channel_symbols(message, hashtab, loctab, symbols);
|
||||||
|
|
||||||
float f0 = 50.0;
|
float f0 = 50.0;
|
||||||
float t0 = 2.0; // Caution!! Possible buffer overflow with the index calculation (no user input here!)
|
float t0 = 2.0; // Caution!! Possible buffer overflow with the index calculation (no user input here!)
|
||||||
|
|
@ -721,28 +721,53 @@ int main(int argc, char **argv) {
|
||||||
case 'f': // Frequency
|
case 'f': // Frequency
|
||||||
if (!strcasecmp(optarg, "LF")) {
|
if (!strcasecmp(optarg, "LF")) {
|
||||||
rx_options.dialfreq = 136000;
|
rx_options.dialfreq = 136000;
|
||||||
|
// Implicit direct sampling for HF bands & lower
|
||||||
|
if (!rx_options.directsampling)
|
||||||
|
rx_options.directsampling = 2;
|
||||||
} else if (!strcasecmp(optarg, "MF")) {
|
} else if (!strcasecmp(optarg, "MF")) {
|
||||||
rx_options.dialfreq = 474200;
|
rx_options.dialfreq = 474200;
|
||||||
|
if (!rx_options.directsampling)
|
||||||
|
rx_options.directsampling = 2;
|
||||||
} else if (!strcasecmp(optarg, "160m")) {
|
} else if (!strcasecmp(optarg, "160m")) {
|
||||||
rx_options.dialfreq = 1836600;
|
rx_options.dialfreq = 1836600;
|
||||||
|
if (!rx_options.directsampling)
|
||||||
|
rx_options.directsampling = 2;
|
||||||
} else if (!strcasecmp(optarg, "80m")) {
|
} else if (!strcasecmp(optarg, "80m")) {
|
||||||
rx_options.dialfreq = 3592600;
|
rx_options.dialfreq = 3592600;
|
||||||
|
if (!rx_options.directsampling)
|
||||||
|
rx_options.directsampling = 2;
|
||||||
} else if (!strcasecmp(optarg, "60m")) {
|
} else if (!strcasecmp(optarg, "60m")) {
|
||||||
rx_options.dialfreq = 5287200;
|
rx_options.dialfreq = 5287200;
|
||||||
|
if (!rx_options.directsampling)
|
||||||
|
rx_options.directsampling = 2;
|
||||||
} else if (!strcasecmp(optarg, "40m")) {
|
} else if (!strcasecmp(optarg, "40m")) {
|
||||||
rx_options.dialfreq = 7038600;
|
rx_options.dialfreq = 7038600;
|
||||||
|
if (!rx_options.directsampling)
|
||||||
|
rx_options.directsampling = 2;
|
||||||
} else if (!strcasecmp(optarg, "30m")) {
|
} else if (!strcasecmp(optarg, "30m")) {
|
||||||
rx_options.dialfreq = 10138700;
|
rx_options.dialfreq = 10138700;
|
||||||
|
if (!rx_options.directsampling)
|
||||||
|
rx_options.directsampling = 2;
|
||||||
} else if (!strcasecmp(optarg, "20m")) {
|
} else if (!strcasecmp(optarg, "20m")) {
|
||||||
rx_options.dialfreq = 14095600;
|
rx_options.dialfreq = 14095600;
|
||||||
|
if (!rx_options.directsampling)
|
||||||
|
rx_options.directsampling = 2;
|
||||||
} else if (!strcasecmp(optarg, "17m")) {
|
} else if (!strcasecmp(optarg, "17m")) {
|
||||||
rx_options.dialfreq = 18104600;
|
rx_options.dialfreq = 18104600;
|
||||||
|
if (!rx_options.directsampling)
|
||||||
|
rx_options.directsampling = 2;
|
||||||
} else if (!strcasecmp(optarg, "15m")) {
|
} else if (!strcasecmp(optarg, "15m")) {
|
||||||
rx_options.dialfreq = 21094600;
|
rx_options.dialfreq = 21094600;
|
||||||
|
if (!rx_options.directsampling)
|
||||||
|
rx_options.directsampling = 2;
|
||||||
} else if (!strcasecmp(optarg, "12m")) {
|
} else if (!strcasecmp(optarg, "12m")) {
|
||||||
rx_options.dialfreq = 24924600;
|
rx_options.dialfreq = 24924600;
|
||||||
|
if (!rx_options.directsampling)
|
||||||
|
rx_options.directsampling = 2;
|
||||||
} else if (!strcasecmp(optarg, "10m")) {
|
} else if (!strcasecmp(optarg, "10m")) {
|
||||||
rx_options.dialfreq = 28124600;
|
rx_options.dialfreq = 28124600;
|
||||||
|
if (!rx_options.directsampling)
|
||||||
|
rx_options.directsampling = 2;
|
||||||
} else if (!strcasecmp(optarg, "6m")) {
|
} else if (!strcasecmp(optarg, "6m")) {
|
||||||
rx_options.dialfreq = 50293000;
|
rx_options.dialfreq = 50293000;
|
||||||
} else if (!strcasecmp(optarg, "4m")) {
|
} else if (!strcasecmp(optarg, "4m")) {
|
||||||
|
|
@ -822,18 +847,18 @@ int main(int argc, char **argv) {
|
||||||
if (rx_options.selftest == true) {
|
if (rx_options.selftest == true) {
|
||||||
if (decoderSelfTest()) {
|
if (decoderSelfTest()) {
|
||||||
fprintf(stdout, "Self-test SUCCESS!\n");
|
fprintf(stdout, "Self-test SUCCESS!\n");
|
||||||
exit(0);
|
return EXIT_SUCCESS;
|
||||||
}
|
}
|
||||||
else {
|
else {
|
||||||
fprintf(stderr, "Self-test FAILED!\n");
|
fprintf(stderr, "Self-test FAILED!\n");
|
||||||
exit(1);
|
return EXIT_FAILURE;
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
if (rx_options.readfile == true) {
|
if (rx_options.readfile == true) {
|
||||||
fprintf(stdout, "Reading IQ file: %s\n", rx_options.filename);
|
fprintf(stdout, "Reading IQ file: %s\n", rx_options.filename);
|
||||||
decodeRecordedFile(rx_options.filename);
|
decodeRecordedFile(rx_options.filename);
|
||||||
exit(0);
|
return EXIT_SUCCESS;
|
||||||
}
|
}
|
||||||
|
|
||||||
if (rx_options.writefile == true) {
|
if (rx_options.writefile == true) {
|
||||||
|
|
@ -843,19 +868,19 @@ int main(int argc, char **argv) {
|
||||||
if (rx_options.dialfreq == 0) {
|
if (rx_options.dialfreq == 0) {
|
||||||
fprintf(stderr, "Please specify a dial frequency.\n");
|
fprintf(stderr, "Please specify a dial frequency.\n");
|
||||||
fprintf(stderr, " --help for usage...\n");
|
fprintf(stderr, " --help for usage...\n");
|
||||||
exit(1);
|
return EXIT_FAILURE;
|
||||||
}
|
}
|
||||||
|
|
||||||
if (dec_options.rcall[0] == 0) {
|
if (dec_options.rcall[0] == 0) {
|
||||||
fprintf(stderr, "Please specify your callsign.\n");
|
fprintf(stderr, "Please specify your callsign.\n");
|
||||||
fprintf(stderr, " --help for usage...\n");
|
fprintf(stderr, " --help for usage...\n");
|
||||||
exit(1);
|
return EXIT_FAILURE;
|
||||||
}
|
}
|
||||||
|
|
||||||
if (dec_options.rloc[0] == 0) {
|
if (dec_options.rloc[0] == 0) {
|
||||||
fprintf(stderr, "Please specify your locator.\n");
|
fprintf(stderr, "Please specify your locator.\n");
|
||||||
fprintf(stderr, " --help for usage...\n");
|
fprintf(stderr, " --help for usage...\n");
|
||||||
exit(1);
|
return EXIT_FAILURE;
|
||||||
}
|
}
|
||||||
|
|
||||||
/* Calcule shift offset */
|
/* Calcule shift offset */
|
||||||
|
|
@ -961,7 +986,7 @@ int main(int argc, char **argv) {
|
||||||
struct tm *gtm = gmtime(&rawtime);
|
struct tm *gtm = gmtime(&rawtime);
|
||||||
|
|
||||||
/* Print used parameter */
|
/* Print used parameter */
|
||||||
printf("\nStarting rtlsdr-wsprd (%04d-%02d-%02d, %02d:%02dz) -- Version 0.4.2\n",
|
printf("\nStarting rtlsdr-wsprd (%04d-%02d-%02d, %02d:%02dz) -- Version 0.5.0\n",
|
||||||
gtm->tm_year + 1900, gtm->tm_mon + 1, gtm->tm_mday, gtm->tm_hour, gtm->tm_min);
|
gtm->tm_year + 1900, gtm->tm_mon + 1, gtm->tm_mday, gtm->tm_hour, gtm->tm_min);
|
||||||
printf(" Callsign : %s\n", dec_options.rcall);
|
printf(" Callsign : %s\n", dec_options.rcall);
|
||||||
printf(" Locator : %s\n", dec_options.rloc);
|
printf(" Locator : %s\n", dec_options.rloc);
|
||||||
|
|
|
||||||
11
wsprd/fano.c
11
wsprd/fano.c
|
|
@ -166,8 +166,7 @@ int fano(unsigned int *metric, // Final path metric (returned value)
|
||||||
*/
|
*/
|
||||||
while (ngamma >= t + delta) t += delta;
|
while (ngamma >= t + delta) t += delta;
|
||||||
}
|
}
|
||||||
// Move forward
|
np[1].gamma = ngamma; // Move forward
|
||||||
np[1].gamma = ngamma;
|
|
||||||
np[1].encstate = np->encstate << 1;
|
np[1].encstate = np->encstate << 1;
|
||||||
if (++np == (lastnode + 1)) {
|
if (++np == (lastnode + 1)) {
|
||||||
break; // Done!
|
break; // Done!
|
||||||
|
|
@ -194,14 +193,11 @@ int fano(unsigned int *metric, // Final path metric (returned value)
|
||||||
np->encstate++; // Set low bit
|
np->encstate++; // Set low bit
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
np->i = 0;
|
np->i = 0; // Start with best branch
|
||||||
// Start with best branch
|
|
||||||
continue;
|
continue;
|
||||||
}
|
}
|
||||||
|
|
||||||
// Threshold violated, can't go forward
|
// Threshold violated, can't go forward
|
||||||
for (;;) {
|
for (;;) { // Look backward
|
||||||
// Look backward
|
|
||||||
if (np == nodes || np[-1].gamma < t) {
|
if (np == nodes || np[-1].gamma < t) {
|
||||||
/* Can't back up either.
|
/* Can't back up either.
|
||||||
* Relax threshold and and look
|
* Relax threshold and and look
|
||||||
|
|
@ -214,7 +210,6 @@ int fano(unsigned int *metric, // Final path metric (returned value)
|
||||||
}
|
}
|
||||||
break;
|
break;
|
||||||
}
|
}
|
||||||
|
|
||||||
// Back up
|
// Back up
|
||||||
if (--np < tail && np->i != 1) {
|
if (--np < tail && np->i != 1) {
|
||||||
np->i++; // Search next best branch
|
np->i++; // Search next best branch
|
||||||
|
|
|
||||||
|
|
@ -4,114 +4,136 @@
|
||||||
* should be normalized to have rms amplitude equal to "symbol_scale".
|
* should be normalized to have rms amplitude equal to "symbol_scale".
|
||||||
********************************************************************************/
|
********************************************************************************/
|
||||||
|
|
||||||
#pragma GCC diagnostic push
|
//float symbol_scale[5]={42.6, 53.3, 72.7, 100.2, 125.4};
|
||||||
#pragma GCC diagnostic ignored "-Wmissing-braces"
|
float metric_tables[5][256]={
|
||||||
|
{0.9782, 0.9695, 0.9689, 0.9669, 0.9666, 0.9653, 0.9638, 0.9618, 0.9599, 0.9601,
|
||||||
// float symbol_scale[4]={42.6, 53.3, 72.7, 100.2};
|
0.9592, 0.9570, 0.9556, 0.9540, 0.9525, 0.9527, 0.9486, 0.9477, 0.9450, 0.9436,
|
||||||
static const float metric_tables[4][256] = {
|
0.9424, 0.9400, 0.9381, 0.9360, 0.9340, 0.9316, 0.9301, 0.9272, 0.9254, 0.9224,
|
||||||
0.9782, 0.9695, 0.9689, 0.9669, 0.9666, 0.9653, 0.9638, 0.9618, 0.9599, 0.9601,
|
0.9196, 0.9171, 0.9154, 0.9123, 0.9076, 0.9061, 0.9030, 0.9000, 0.8965, 0.8934,
|
||||||
0.9592, 0.9570, 0.9556, 0.9540, 0.9525, 0.9527, 0.9486, 0.9477, 0.9450, 0.9436,
|
0.8903, 0.8874, 0.8834, 0.8792, 0.8760, 0.8726, 0.8685, 0.8639, 0.8599, 0.8550,
|
||||||
0.9424, 0.9400, 0.9381, 0.9360, 0.9340, 0.9316, 0.9301, 0.9272, 0.9254, 0.9224,
|
0.8504, 0.8459, 0.8422, 0.8364, 0.8320, 0.8262, 0.8215, 0.8159, 0.8111, 0.8052,
|
||||||
0.9196, 0.9171, 0.9154, 0.9123, 0.9076, 0.9061, 0.9030, 0.9000, 0.8965, 0.8934,
|
0.7996, 0.7932, 0.7878, 0.7812, 0.7745, 0.7685, 0.7616, 0.7550, 0.7479, 0.7405,
|
||||||
0.8903, 0.8874, 0.8834, 0.8792, 0.8760, 0.8726, 0.8685, 0.8639, 0.8599, 0.8550,
|
0.7336, 0.7255, 0.7184, 0.7102, 0.7016, 0.6946, 0.6860, 0.6769, 0.6687, 0.6598,
|
||||||
0.8504, 0.8459, 0.8422, 0.8364, 0.8320, 0.8262, 0.8215, 0.8159, 0.8111, 0.8052,
|
0.6503, 0.6416, 0.6325, 0.6219, 0.6122, 0.6016, 0.5920, 0.5818, 0.5711, 0.5606,
|
||||||
0.7996, 0.7932, 0.7878, 0.7812, 0.7745, 0.7685, 0.7616, 0.7550, 0.7479, 0.7405,
|
0.5487, 0.5374, 0.5266, 0.5142, 0.5020, 0.4908, 0.4784, 0.4663, 0.4532, 0.4405,
|
||||||
0.7336, 0.7255, 0.7184, 0.7102, 0.7016, 0.6946, 0.6860, 0.6769, 0.6687, 0.6598,
|
0.4271, 0.4144, 0.4006, 0.3865, 0.3731, 0.3594, 0.3455, 0.3304, 0.3158, 0.3009,
|
||||||
0.6503, 0.6416, 0.6325, 0.6219, 0.6122, 0.6016, 0.5920, 0.5818, 0.5711, 0.5606,
|
0.2858, 0.2708, 0.2560, 0.2399, 0.2233, 0.2074, 0.1919, 0.1756, 0.1590, 0.1427,
|
||||||
0.5487, 0.5374, 0.5266, 0.5142, 0.5020, 0.4908, 0.4784, 0.4663, 0.4532, 0.4405,
|
0.1251, 0.1074, 0.0905, 0.0722, 0.0550, 0.0381, 0.0183, 0.0000, -0.0185, -0.0391,
|
||||||
0.4271, 0.4144, 0.4006, 0.3865, 0.3731, 0.3594, 0.3455, 0.3304, 0.3158, 0.3009,
|
-0.0571, -0.0760, -0.0966, -0.1160, -0.1370, -0.1584, -0.1787, -0.1999, -0.2214, -0.2423,
|
||||||
0.2858, 0.2708, 0.2560, 0.2399, 0.2233, 0.2074, 0.1919, 0.1756, 0.1590, 0.1427,
|
-0.2643, -0.2879, -0.3114, -0.3336, -0.3568, -0.3806, -0.4050, -0.4293, -0.4552, -0.4798,
|
||||||
0.1251, 0.1074, 0.0905, 0.0722, 0.0550, 0.0381, 0.0183, 0.0000, -0.0185, -0.0391,
|
-0.5046, -0.5296, -0.5564, -0.5836, -0.6093, -0.6372, -0.6645, -0.6933, -0.7208, -0.7495,
|
||||||
-0.0571, -0.0760, -0.0966, -0.1160, -0.1370, -0.1584, -0.1787, -0.1999, -0.2214, -0.2423,
|
-0.7763, -0.8065, -0.8378, -0.8660, -0.8964, -0.9293, -0.9592, -0.9907, -1.0214, -1.0509,
|
||||||
-0.2643, -0.2879, -0.3114, -0.3336, -0.3568, -0.3806, -0.4050, -0.4293, -0.4552, -0.4798,
|
-1.0850, -1.1168, -1.1528, -1.1847, -1.2157, -1.2511, -1.2850, -1.3174, -1.3540, -1.3900,
|
||||||
-0.5046, -0.5296, -0.5564, -0.5836, -0.6093, -0.6372, -0.6645, -0.6933, -0.7208, -0.7495,
|
-1.4201, -1.4580, -1.4956, -1.5292, -1.5683, -1.6030, -1.6411, -1.6789, -1.7147, -1.7539,
|
||||||
-0.7763, -0.8065, -0.8378, -0.8660, -0.8964, -0.9293, -0.9592, -0.9907, -1.0214, -1.0509,
|
-1.7887, -1.8289, -1.8699, -1.9043, -1.9469, -1.9849, -2.0267, -2.0610, -2.1028, -2.1391,
|
||||||
-1.0850, -1.1168, -1.1528, -1.1847, -1.2157, -1.2511, -1.2850, -1.3174, -1.3540, -1.3900,
|
-2.1855, -2.2215, -2.2712, -2.3033, -2.3440, -2.3870, -2.4342, -2.4738, -2.5209, -2.5646,
|
||||||
-1.4201, -1.4580, -1.4956, -1.5292, -1.5683, -1.6030, -1.6411, -1.6789, -1.7147, -1.7539,
|
-2.6016, -2.6385, -2.6868, -2.7356, -2.7723, -2.8111, -2.8524, -2.9009, -2.9428, -2.9879,
|
||||||
-1.7887, -1.8289, -1.8699, -1.9043, -1.9469, -1.9849, -2.0267, -2.0610, -2.1028, -2.1391,
|
-3.0103, -3.0832, -3.1340, -3.1628, -3.2049, -3.2557, -3.3101, -3.3453, -3.4025, -3.4317,
|
||||||
-2.1855, -2.2215, -2.2712, -2.3033, -2.3440, -2.3870, -2.4342, -2.4738, -2.5209, -2.5646,
|
-3.4828, -3.5270, -3.5745, -3.6181, -3.6765, -3.7044, -3.7410, -3.8118, -3.8368, -3.9549,
|
||||||
-2.6016, -2.6385, -2.6868, -2.7356, -2.7723, -2.8111, -2.8524, -2.9009, -2.9428, -2.9879,
|
-3.9488, -3.9941, -4.0428, -4.0892, -4.1648, -4.1965, -4.1892, -4.2565, -4.3356, -4.3948,
|
||||||
-3.0103, -3.0832, -3.1340, -3.1628, -3.2049, -3.2557, -3.3101, -3.3453, -3.4025, -3.4317,
|
-4.4481, -4.4607, -4.5533, -4.5809, -4.5927, -5.1047},
|
||||||
-3.4828, -3.5270, -3.5745, -3.6181, -3.6765, -3.7044, -3.7410, -3.8118, -3.8368, -3.9549,
|
{0.9978, 0.9962, 0.9961, 0.9959, 0.9958, 0.9954, 0.9949, 0.9950, 0.9947, 0.9942,
|
||||||
-3.9488, -3.9941, -4.0428, -4.0892, -4.1648, -4.1965, -4.1892, -4.2565, -4.3356, -4.3948,
|
0.9940, 0.9939, 0.9933, 0.9931, 0.9928, 0.9924, 0.9921, 0.9916, 0.9911, 0.9909,
|
||||||
-4.4481, -4.4607, -4.5533, -4.5809, -4.5927, -5.1047,
|
0.9903, 0.9900, 0.9892, 0.9887, 0.9883, 0.9877, 0.9869, 0.9863, 0.9857, 0.9848,
|
||||||
0.9978, 0.9962, 0.9961, 0.9959, 0.9958, 0.9954, 0.9949, 0.9950, 0.9947, 0.9942,
|
0.9842, 0.9835, 0.9825, 0.9817, 0.9808, 0.9799, 0.9791, 0.9777, 0.9767, 0.9757,
|
||||||
0.9940, 0.9939, 0.9933, 0.9931, 0.9928, 0.9924, 0.9921, 0.9916, 0.9911, 0.9909,
|
0.9744, 0.9729, 0.9716, 0.9704, 0.9690, 0.9674, 0.9656, 0.9641, 0.9625, 0.9609,
|
||||||
0.9903, 0.9900, 0.9892, 0.9887, 0.9883, 0.9877, 0.9869, 0.9863, 0.9857, 0.9848,
|
0.9587, 0.9567, 0.9548, 0.9524, 0.9501, 0.9478, 0.9453, 0.9426, 0.9398, 0.9371,
|
||||||
0.9842, 0.9835, 0.9825, 0.9817, 0.9808, 0.9799, 0.9791, 0.9777, 0.9767, 0.9757,
|
0.9339, 0.9311, 0.9277, 0.9242, 0.9206, 0.9168, 0.9131, 0.9087, 0.9043, 0.8999,
|
||||||
0.9744, 0.9729, 0.9716, 0.9704, 0.9690, 0.9674, 0.9656, 0.9641, 0.9625, 0.9609,
|
0.8953, 0.8907, 0.8857, 0.8803, 0.8747, 0.8690, 0.8632, 0.8572, 0.8507, 0.8439,
|
||||||
0.9587, 0.9567, 0.9548, 0.9524, 0.9501, 0.9478, 0.9453, 0.9426, 0.9398, 0.9371,
|
0.8368, 0.8295, 0.8217, 0.8138, 0.8058, 0.7972, 0.7883, 0.7784, 0.7694, 0.7597,
|
||||||
0.9339, 0.9311, 0.9277, 0.9242, 0.9206, 0.9168, 0.9131, 0.9087, 0.9043, 0.8999,
|
0.7489, 0.7378, 0.7269, 0.7152, 0.7030, 0.6911, 0.6782, 0.6643, 0.6506, 0.6371,
|
||||||
0.8953, 0.8907, 0.8857, 0.8803, 0.8747, 0.8690, 0.8632, 0.8572, 0.8507, 0.8439,
|
0.6211, 0.6054, 0.5897, 0.5740, 0.5565, 0.5393, 0.5214, 0.5027, 0.4838, 0.4643,
|
||||||
0.8368, 0.8295, 0.8217, 0.8138, 0.8058, 0.7972, 0.7883, 0.7784, 0.7694, 0.7597,
|
0.4436, 0.4225, 0.4004, 0.3787, 0.3562, 0.3324, 0.3089, 0.2839, 0.2584, 0.2321,
|
||||||
0.7489, 0.7378, 0.7269, 0.7152, 0.7030, 0.6911, 0.6782, 0.6643, 0.6506, 0.6371,
|
0.2047, 0.1784, 0.1499, 0.1213, 0.0915, 0.0628, 0.0314, 0.0000, -0.0321, -0.0657,
|
||||||
0.6211, 0.6054, 0.5897, 0.5740, 0.5565, 0.5393, 0.5214, 0.5027, 0.4838, 0.4643,
|
-0.0977, -0.1324, -0.1673, -0.2036, -0.2387, -0.2768, -0.3150, -0.3538, -0.3936, -0.4327,
|
||||||
0.4436, 0.4225, 0.4004, 0.3787, 0.3562, 0.3324, 0.3089, 0.2839, 0.2584, 0.2321,
|
-0.4739, -0.5148, -0.5561, -0.6000, -0.6438, -0.6889, -0.7331, -0.7781, -0.8247, -0.8712,
|
||||||
0.2047, 0.1784, 0.1499, 0.1213, 0.0915, 0.0628, 0.0314, 0.0000, -0.0321, -0.0657,
|
-0.9177, -0.9677, -1.0142, -1.0631, -1.1143, -1.1686, -1.2169, -1.2680, -1.3223, -1.3752,
|
||||||
-0.0977, -0.1324, -0.1673, -0.2036, -0.2387, -0.2768, -0.3150, -0.3538, -0.3936, -0.4327,
|
-1.4261, -1.4806, -1.5356, -1.5890, -1.6462, -1.7041, -1.7591, -1.8124, -1.8735, -1.9311,
|
||||||
-0.4739, -0.5148, -0.5561, -0.6000, -0.6438, -0.6889, -0.7331, -0.7781, -0.8247, -0.8712,
|
-1.9891, -2.0459, -2.1048, -2.1653, -2.2248, -2.2855, -2.3466, -2.4079, -2.4668, -2.5263,
|
||||||
-0.9177, -0.9677, -1.0142, -1.0631, -1.1143, -1.1686, -1.2169, -1.2680, -1.3223, -1.3752,
|
-2.5876, -2.6507, -2.7142, -2.7761, -2.8366, -2.8995, -2.9620, -3.0279, -3.0973, -3.1576,
|
||||||
-1.4261, -1.4806, -1.5356, -1.5890, -1.6462, -1.7041, -1.7591, -1.8124, -1.8735, -1.9311,
|
-3.2238, -3.2890, -3.3554, -3.4215, -3.4805, -3.5518, -3.6133, -3.6812, -3.7473, -3.8140,
|
||||||
-1.9891, -2.0459, -2.1048, -2.1653, -2.2248, -2.2855, -2.3466, -2.4079, -2.4668, -2.5263,
|
-3.8781, -3.9450, -4.0184, -4.0794, -4.1478, -4.2241, -4.2853, -4.3473, -4.4062, -4.4839,
|
||||||
-2.5876, -2.6507, -2.7142, -2.7761, -2.8366, -2.8995, -2.9620, -3.0279, -3.0973, -3.1576,
|
-4.5539, -4.6202, -4.6794, -4.7478, -4.8309, -4.9048, -4.9669, -5.0294, -5.1194, -5.1732,
|
||||||
-3.2238, -3.2890, -3.3554, -3.4215, -3.4805, -3.5518, -3.6133, -3.6812, -3.7473, -3.8140,
|
-5.2378, -5.3094, -5.3742, -5.4573, -5.5190, -5.5728, -5.6637, -5.7259, -5.7843, -5.8854,
|
||||||
-3.8781, -3.9450, -4.0184, -4.0794, -4.1478, -4.2241, -4.2853, -4.3473, -4.4062, -4.4839,
|
-5.9553, -6.0054, -6.0656, -6.1707, -6.2241, -6.3139, -6.3393, -6.4356, -6.5153, -6.5758,
|
||||||
-4.5539, -4.6202, -4.6794, -4.7478, -4.8309, -4.9048, -4.9669, -5.0294, -5.1194, -5.1732,
|
-6.6506, -6.7193, -6.7542, -6.8942, -6.9219, -6.9605, -7.1013, -7.1895, -7.1549, -7.2799,
|
||||||
-5.2378, -5.3094, -5.3742, -5.4573, -5.5190, -5.5728, -5.6637, -5.7259, -5.7843, -5.8854,
|
-7.4119, -7.4608, -7.5256, -7.5879, -7.7598, -8.4120},
|
||||||
-5.9553, -6.0054, -6.0656, -6.1707, -6.2241, -6.3139, -6.3393, -6.4356, -6.5153, -6.5758,
|
{0.9999, 0.9998, 0.9998, 0.9998, 0.9998, 0.9998, 0.9997, 0.9997, 0.9997, 0.9997,
|
||||||
-6.6506, -6.7193, -6.7542, -6.8942, -6.9219, -6.9605, -7.1013, -7.1895, -7.1549, -7.2799,
|
0.9997, 0.9996, 0.9996, 0.9996, 0.9995, 0.9995, 0.9994, 0.9994, 0.9994, 0.9993,
|
||||||
-7.4119, -7.4608, -7.5256, -7.5879, -7.7598, -8.4120,
|
0.9993, 0.9992, 0.9991, 0.9991, 0.9990, 0.9989, 0.9988, 0.9988, 0.9988, 0.9986,
|
||||||
0.9999, 0.9998, 0.9998, 0.9998, 0.9998, 0.9998, 0.9997, 0.9997, 0.9997, 0.9997,
|
0.9985, 0.9984, 0.9983, 0.9982, 0.9980, 0.9979, 0.9977, 0.9976, 0.9974, 0.9971,
|
||||||
0.9997, 0.9996, 0.9996, 0.9996, 0.9995, 0.9995, 0.9994, 0.9994, 0.9994, 0.9993,
|
0.9969, 0.9968, 0.9965, 0.9962, 0.9960, 0.9957, 0.9953, 0.9950, 0.9947, 0.9941,
|
||||||
0.9993, 0.9992, 0.9991, 0.9991, 0.9990, 0.9989, 0.9988, 0.9988, 0.9988, 0.9986,
|
0.9937, 0.9933, 0.9928, 0.9922, 0.9917, 0.9911, 0.9904, 0.9897, 0.9890, 0.9882,
|
||||||
0.9985, 0.9984, 0.9983, 0.9982, 0.9980, 0.9979, 0.9977, 0.9976, 0.9974, 0.9971,
|
0.9874, 0.9863, 0.9855, 0.9843, 0.9832, 0.9819, 0.9806, 0.9792, 0.9777, 0.9760,
|
||||||
0.9969, 0.9968, 0.9965, 0.9962, 0.9960, 0.9957, 0.9953, 0.9950, 0.9947, 0.9941,
|
0.9743, 0.9724, 0.9704, 0.9683, 0.9659, 0.9634, 0.9609, 0.9581, 0.9550, 0.9516,
|
||||||
0.9937, 0.9933, 0.9928, 0.9922, 0.9917, 0.9911, 0.9904, 0.9897, 0.9890, 0.9882,
|
0.9481, 0.9446, 0.9406, 0.9363, 0.9317, 0.9270, 0.9218, 0.9160, 0.9103, 0.9038,
|
||||||
0.9874, 0.9863, 0.9855, 0.9843, 0.9832, 0.9819, 0.9806, 0.9792, 0.9777, 0.9760,
|
0.8972, 0.8898, 0.8822, 0.8739, 0.8647, 0.8554, 0.8457, 0.8357, 0.8231, 0.8115,
|
||||||
0.9743, 0.9724, 0.9704, 0.9683, 0.9659, 0.9634, 0.9609, 0.9581, 0.9550, 0.9516,
|
0.7984, 0.7854, 0.7704, 0.7556, 0.7391, 0.7210, 0.7038, 0.6840, 0.6633, 0.6408,
|
||||||
0.9481, 0.9446, 0.9406, 0.9363, 0.9317, 0.9270, 0.9218, 0.9160, 0.9103, 0.9038,
|
0.6174, 0.5939, 0.5678, 0.5410, 0.5137, 0.4836, 0.4524, 0.4193, 0.3850, 0.3482,
|
||||||
0.8972, 0.8898, 0.8822, 0.8739, 0.8647, 0.8554, 0.8457, 0.8357, 0.8231, 0.8115,
|
0.3132, 0.2733, 0.2315, 0.1891, 0.1435, 0.0980, 0.0493, 0.0000, -0.0510, -0.1052,
|
||||||
0.7984, 0.7854, 0.7704, 0.7556, 0.7391, 0.7210, 0.7038, 0.6840, 0.6633, 0.6408,
|
-0.1593, -0.2177, -0.2759, -0.3374, -0.4005, -0.4599, -0.5266, -0.5935, -0.6626, -0.7328,
|
||||||
0.6174, 0.5939, 0.5678, 0.5410, 0.5137, 0.4836, 0.4524, 0.4193, 0.3850, 0.3482,
|
-0.8051, -0.8757, -0.9498, -1.0271, -1.1019, -1.1816, -1.2642, -1.3459, -1.4295, -1.5077,
|
||||||
0.3132, 0.2733, 0.2315, 0.1891, 0.1435, 0.0980, 0.0493, 0.0000, -0.0510, -0.1052,
|
-1.5958, -1.6818, -1.7647, -1.8548, -1.9387, -2.0295, -2.1152, -2.2154, -2.3011, -2.3904,
|
||||||
-0.1593, -0.2177, -0.2759, -0.3374, -0.4005, -0.4599, -0.5266, -0.5935, -0.6626, -0.7328,
|
-2.4820, -2.5786, -2.6730, -2.7652, -2.8616, -2.9546, -3.0526, -3.1445, -3.2445, -3.3416,
|
||||||
-0.8051, -0.8757, -0.9498, -1.0271, -1.1019, -1.1816, -1.2642, -1.3459, -1.4295, -1.5077,
|
-3.4357, -3.5325, -3.6324, -3.7313, -3.8225, -3.9209, -4.0248, -4.1278, -4.2261, -4.3193,
|
||||||
-1.5958, -1.6818, -1.7647, -1.8548, -1.9387, -2.0295, -2.1152, -2.2154, -2.3011, -2.3904,
|
-4.4220, -4.5262, -4.6214, -4.7242, -4.8234, -4.9245, -5.0298, -5.1250, -5.2232, -5.3267,
|
||||||
-2.4820, -2.5786, -2.6730, -2.7652, -2.8616, -2.9546, -3.0526, -3.1445, -3.2445, -3.3416,
|
-5.4332, -5.5342, -5.6431, -5.7270, -5.8401, -5.9350, -6.0407, -6.1418, -6.2363, -6.3384,
|
||||||
-3.4357, -3.5325, -3.6324, -3.7313, -3.8225, -3.9209, -4.0248, -4.1278, -4.2261, -4.3193,
|
-6.4536, -6.5429, -6.6582, -6.7433, -6.8438, -6.9478, -7.0789, -7.1894, -7.2714, -7.3815,
|
||||||
-4.4220, -4.5262, -4.6214, -4.7242, -4.8234, -4.9245, -5.0298, -5.1250, -5.2232, -5.3267,
|
-7.4810, -7.5575, -7.6852, -7.8071, -7.8580, -7.9724, -8.1000, -8.2207, -8.2867, -8.4017,
|
||||||
-5.4332, -5.5342, -5.6431, -5.7270, -5.8401, -5.9350, -6.0407, -6.1418, -6.2363, -6.3384,
|
-8.5287, -8.6347, -8.7082, -8.8319, -8.9448, -9.0355, -9.1885, -9.2095, -9.2863, -9.4186,
|
||||||
-6.4536, -6.5429, -6.6582, -6.7433, -6.8438, -6.9478, -7.0789, -7.1894, -7.2714, -7.3815,
|
-9.5064, -9.6386, -9.7207, -9.8286, -9.9453, -10.0701, -10.1735, -10.3001, -10.2858, -10.5427,
|
||||||
-7.4810, -7.5575, -7.6852, -7.8071, -7.8580, -7.9724, -8.1000, -8.2207, -8.2867, -8.4017,
|
-10.5982, -10.7361, -10.7042, -10.9212, -11.0097, -11.0469, -11.1155, -11.2812, -11.3472, -11.4988,
|
||||||
-8.5287, -8.6347, -8.7082, -8.8319, -8.9448, -9.0355, -9.1885, -9.2095, -9.2863, -9.4186,
|
-11.5327, -11.6692, -11.9376, -11.8606, -12.1372, -13.2539},
|
||||||
-9.5064, -9.6386, -9.7207, -9.8286, -9.9453, -10.0701, -10.1735, -10.3001, -10.2858, -10.5427,
|
{1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000,
|
||||||
-10.5982, -10.7361, -10.7042, -10.9212, -11.0097, -11.0469, -11.1155, -11.2812, -11.3472, -11.4988,
|
1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000,
|
||||||
-11.5327, -11.6692, -11.9376, -11.8606, -12.1372, -13.2539,
|
1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000,
|
||||||
1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000,
|
0.9999, 0.9999, 0.9999, 0.9999, 0.9999, 0.9999, 0.9999, 0.9999, 0.9999, 0.9999,
|
||||||
1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000,
|
0.9999, 0.9998, 0.9998, 0.9998, 0.9998, 0.9997, 0.9997, 0.9997, 0.9997, 0.9996,
|
||||||
1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000,
|
0.9996, 0.9995, 0.9995, 0.9994, 0.9994, 0.9993, 0.9992, 0.9991, 0.9991, 0.9989,
|
||||||
0.9999, 0.9999, 0.9999, 0.9999, 0.9999, 0.9999, 0.9999, 0.9999, 0.9999, 0.9999,
|
0.9988, 0.9986, 0.9985, 0.9983, 0.9981, 0.9980, 0.9977, 0.9974, 0.9971, 0.9968,
|
||||||
0.9999, 0.9998, 0.9998, 0.9998, 0.9998, 0.9997, 0.9997, 0.9997, 0.9997, 0.9996,
|
0.9965, 0.9962, 0.9956, 0.9950, 0.9948, 0.9941, 0.9933, 0.9926, 0.9919, 0.9910,
|
||||||
0.9996, 0.9995, 0.9995, 0.9994, 0.9994, 0.9993, 0.9992, 0.9991, 0.9991, 0.9989,
|
0.9899, 0.9889, 0.9877, 0.9863, 0.9845, 0.9829, 0.9811, 0.9791, 0.9769, 0.9741,
|
||||||
0.9988, 0.9986, 0.9985, 0.9983, 0.9981, 0.9980, 0.9977, 0.9974, 0.9971, 0.9968,
|
0.9716, 0.9684, 0.9645, 0.9611, 0.9563, 0.9519, 0.9463, 0.9406, 0.9344, 0.9272,
|
||||||
0.9965, 0.9962, 0.9956, 0.9950, 0.9948, 0.9941, 0.9933, 0.9926, 0.9919, 0.9910,
|
0.9197, 0.9107, 0.9016, 0.8903, 0.8791, 0.8653, 0.8523, 0.8357, 0.8179, 0.7988,
|
||||||
0.9899, 0.9889, 0.9877, 0.9863, 0.9845, 0.9829, 0.9811, 0.9791, 0.9769, 0.9741,
|
0.7779, 0.7562, 0.7318, 0.7024, 0.6753, 0.6435, 0.6089, 0.5700, 0.5296, 0.4860,
|
||||||
0.9716, 0.9684, 0.9645, 0.9611, 0.9563, 0.9519, 0.9463, 0.9406, 0.9344, 0.9272,
|
0.4366, 0.3855, 0.3301, 0.2735, 0.2114, 0.1443, 0.0682, 0.0000, -0.0715, -0.1604,
|
||||||
0.9197, 0.9107, 0.9016, 0.8903, 0.8791, 0.8653, 0.8523, 0.8357, 0.8179, 0.7988,
|
-0.2478, -0.3377, -0.4287, -0.5277, -0.6291, -0.7384, -0.8457, -0.9559, -1.0742, -1.1913,
|
||||||
0.7779, 0.7562, 0.7318, 0.7024, 0.6753, 0.6435, 0.6089, 0.5700, 0.5296, 0.4860,
|
-1.3110, -1.4238, -1.5594, -1.6854, -1.8093, -1.9414, -2.0763, -2.2160, -2.3611, -2.4876,
|
||||||
0.4366, 0.3855, 0.3301, 0.2735, 0.2114, 0.1443, 0.0682, 0.0000, -0.0715, -0.1604,
|
-2.6374, -2.7710, -2.9225, -3.0591, -3.2077, -3.3452, -3.4916, -3.6316, -3.7735, -3.9296,
|
||||||
-0.2478, -0.3377, -0.4287, -0.5277, -0.6291, -0.7384, -0.8457, -0.9559, -1.0742, -1.1913,
|
-4.0682, -4.2334, -4.3607, -4.5270, -4.6807, -4.8108, -4.9753, -5.1212, -5.2631, -5.4042,
|
||||||
-1.3110, -1.4238, -1.5594, -1.6854, -1.8093, -1.9414, -2.0763, -2.2160, -2.3611, -2.4876,
|
-5.5510, -5.7227, -5.8794, -6.0244, -6.1677, -6.3271, -6.4862, -6.6130, -6.7449, -6.9250,
|
||||||
-2.6374, -2.7710, -2.9225, -3.0591, -3.2077, -3.3452, -3.4916, -3.6316, -3.7735, -3.9296,
|
-7.1232, -7.1736, -7.3628, -7.5596, -7.6906, -7.8129, -7.9817, -8.1440, -8.3016, -8.4797,
|
||||||
-4.0682, -4.2334, -4.3607, -4.5270, -4.6807, -4.8108, -4.9753, -5.1212, -5.2631, -5.4042,
|
-8.5734, -8.7692, -8.9198, -9.0610, -9.1746, -9.3536, -9.5939, -9.6957, -9.8475, -9.9639,
|
||||||
-5.5510, -5.7227, -5.8794, -6.0244, -6.1677, -6.3271, -6.4862, -6.6130, -6.7449, -6.9250,
|
-10.1730, -10.2427, -10.4573, -10.5413, -10.7303, -10.9339, -11.0215, -11.2047, -11.2894, -11.4572,
|
||||||
-7.1232, -7.1736, -7.3628, -7.5596, -7.6906, -7.8129, -7.9817, -8.1440, -8.3016, -8.4797,
|
-11.6256, -11.7794, -11.8801, -12.1717, -12.2354, -12.3686, -12.6195, -12.6527, -12.8247, -12.9560,
|
||||||
-8.5734, -8.7692, -8.9198, -9.0610, -9.1746, -9.3536, -9.5939, -9.6957, -9.8475, -9.9639,
|
-13.3265, -13.1667, -13.4274, -13.6064, -13.5515, -13.9501, -13.9926, -14.4049, -14.1653, -14.4348,
|
||||||
-10.1730, -10.2427, -10.4573, -10.5413, -10.7303, -10.9339, -11.0215, -11.2047, -11.2894, -11.4572,
|
-14.7983, -14.7807, -15.2349, -15.3536, -15.3026, -15.2739, -15.7170, -16.2161, -15.9185, -15.9490,
|
||||||
-11.6256, -11.7794, -11.8801, -12.1717, -12.2354, -12.3686, -12.6195, -12.6527, -12.8247, -12.9560,
|
-16.6258, -16.5568, -16.4318, -16.7999, -16.4101, -17.6393, -17.7643, -17.2644, -17.5973, -17.0403,
|
||||||
-13.3265, -13.1667, -13.4274, -13.6064, -13.5515, -13.9501, -13.9926, -14.4049, -14.1653, -14.4348,
|
-17.7039, -18.0073, -18.1840, -18.3848, -18.6286, -20.7063},
|
||||||
-14.7983, -14.7807, -15.2349, -15.3536, -15.3026, -15.2739, -15.7170, -16.2161, -15.9185, -15.9490,
|
{1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000,
|
||||||
-16.6258, -16.5568, -16.4318, -16.7999, -16.4101, -17.6393, -17.7643, -17.2644, -17.5973, -17.0403,
|
1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000,
|
||||||
-17.7039, -18.0073, -18.1840, -18.3848, -18.6286, -20.7063};
|
1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000, 1.0000,
|
||||||
|
0.9999, 0.9999, 0.9999, 0.9999, 0.9999, 0.9999, 0.9999, 0.9999, 0.9999, 0.9999,
|
||||||
#pragma GCC diagnostic pop
|
0.9999, 0.9998, 0.9998, 0.9998, 0.9998, 0.9997, 0.9997, 0.9997, 0.9997, 0.9996,
|
||||||
|
0.9996, 0.9995, 0.9995, 0.9994, 0.9994, 0.9993, 0.9992, 0.9991, 0.9991, 0.9989,
|
||||||
|
0.9988, 0.9986, 0.9985, 0.9983, 0.9981, 0.9980, 0.9977, 0.9974, 0.9971, 0.9968,
|
||||||
|
0.9965, 0.9962, 0.9956, 0.9950, 0.9948, 0.9941, 0.9933, 0.9926, 0.9919, 0.9910,
|
||||||
|
0.9899, 0.9889, 0.9877, 0.9863, 0.9845, 0.9829, 0.9811, 0.9791, 0.9769, 0.9741,
|
||||||
|
0.9716, 0.9684, 0.9645, 0.9611, 0.9563, 0.9519, 0.9463, 0.9406, 0.9344, 0.9272,
|
||||||
|
0.9197, 0.9107, 0.9016, 0.8903, 0.8791, 0.8653, 0.8523, 0.8357, 0.8179, 0.7988,
|
||||||
|
0.7779, 0.7562, 0.7318, 0.7024, 0.6753, 0.6435, 0.6089, 0.5700, 0.5296, 0.4860,
|
||||||
|
0.4366, 0.3855, 0.3301, 0.2735, 0.2114, 0.1443, 0.0682, 0.0000, -0.0715, -0.1604,
|
||||||
|
-0.2478, -0.3377, -0.4287, -0.5277, -0.6291, -0.7384, -0.8457, -0.9559, -1.0742, -1.1913,
|
||||||
|
-1.3110, -1.4238, -1.5594, -1.6854, -1.8093, -1.9414, -2.0763, -2.2160, -2.3611, -2.4876,
|
||||||
|
-2.6374, -2.7710, -2.9225, -3.0591, -3.2077, -3.3452, -3.4916, -3.6316, -3.7735, -3.9296,
|
||||||
|
-4.0682, -4.2334, -4.3607, -4.5270, -4.6807, -4.8108, -4.9753, -5.1212, -5.2631, -5.4042,
|
||||||
|
-5.5510, -5.7227, -5.8794, -6.0244, -6.1677, -6.3271, -6.4862, -6.6130, -6.7449, -6.9250,
|
||||||
|
-7.1232, -7.1736, -7.3628, -7.5596, -7.6906, -7.8129, -7.9817, -8.1440, -8.3016, -8.4797,
|
||||||
|
-8.5734, -8.7692, -8.9198, -9.0610, -9.1746, -9.3536, -9.5939, -9.6957, -9.8475, -9.9639,
|
||||||
|
-10.1730, -10.2427, -10.4573, -11.7794, -11.8801, -12.1717, -12.2354, -12.3686, -12.6195, -12.6527,
|
||||||
|
-12.8247, -12.9560, -13.3265, -13.1667, -13.4274, -13.6064, -13.5515, -13.9501, -13.9926, -14.4049,
|
||||||
|
-14.1653, -14.4348, -14.7983, -14.7807, -15.2349, -15.3536, -15.3026, -15.2739, -15.7170, -16.2161,
|
||||||
|
-15.9185, -15.9490, -16.6258, -16.5568, -16.4318, -16.7999, -16.4101, -17.6393, -17.7643, -17.2644,
|
||||||
|
-17.5973, -17.0403, -17.7039, -18.0073, -18.1840, -18.3848, -18.6286, -20.7063, 1.43370769e-019,
|
||||||
|
2.64031087e-006, 6.6908396e+031, 1.77537994e+028, 2.79322819e+020, 1.94326e-019,
|
||||||
|
0.00019371575, 2.80722121e-041}};
|
||||||
|
|
|
||||||
|
|
@ -449,10 +449,3 @@ uint32_t nhash(const void *key, size_t length, uint32_t initval) {
|
||||||
|
|
||||||
return c;
|
return c;
|
||||||
}
|
}
|
||||||
|
|
||||||
/*
|
|
||||||
* Fortran argument compatible wrapper
|
|
||||||
*/
|
|
||||||
uint32_t nhash_(const void *key, size_t const *length, uint32_t const *initval) {
|
|
||||||
return nhash(key, *length, *initval);
|
|
||||||
}
|
|
||||||
|
|
|
||||||
|
|
@ -1,4 +1,3 @@
|
||||||
#pragma once
|
#pragma once
|
||||||
|
|
||||||
uint32_t nhash(const void* key, size_t length, uint32_t initval);
|
uint32_t nhash(const void* key, size_t length, uint32_t initval);
|
||||||
uint32_t nhash_(void const* key, size_t const* length, uint32_t const* initval);
|
|
||||||
|
|
|
||||||
807
wsprd/wsprd.c
807
wsprd/wsprd.c
Plik diff jest za duży
Load Diff
|
|
@ -35,37 +35,71 @@
|
||||||
|
|
||||||
/* Option & config of decoder (Shared with the wsprd code) */
|
/* Option & config of decoder (Shared with the wsprd code) */
|
||||||
struct decoder_options {
|
struct decoder_options {
|
||||||
uint32_t freq; // Dial frequency
|
int freq; // Dial frequency
|
||||||
char rcall[13]; // Callsign of the RX station
|
char rcall[13]; // Callsign of the RX station
|
||||||
char rloc[7]; // Locator of the RX station
|
char rloc[7]; // Locator of the RX station
|
||||||
uint32_t quickmode; // Decoder option & tweak
|
int quickmode; // Decoder option & tweak
|
||||||
uint32_t usehashtable; // ''
|
int usehashtable; // ''
|
||||||
uint32_t npasses; // ''
|
int npasses; // ''
|
||||||
uint32_t subtraction; // ''
|
int subtraction; // ''
|
||||||
|
};
|
||||||
|
|
||||||
|
struct cand {
|
||||||
|
float freq;
|
||||||
|
float snr;
|
||||||
|
int shift;
|
||||||
|
float drift;
|
||||||
|
float sync;
|
||||||
};
|
};
|
||||||
|
|
||||||
struct decoder_results {
|
struct decoder_results {
|
||||||
double freq;
|
double freq;
|
||||||
float sync;
|
float sync;
|
||||||
float snr;
|
float snr;
|
||||||
float dt;
|
float dt;
|
||||||
float drift;
|
float drift;
|
||||||
int32_t jitter;
|
int jitter;
|
||||||
char message[23];
|
char message[23];
|
||||||
char call[13];
|
char call[13];
|
||||||
char loc[7];
|
char loc[7];
|
||||||
char pwr[3];
|
char pwr[3];
|
||||||
uint32_t cycles;
|
int cycles;
|
||||||
};
|
};
|
||||||
|
|
||||||
void sync_and_demodulate(float *id, float *qd, long np,
|
void sync_and_demodulate(float *id,
|
||||||
uint8_t *symbols, float *freq1, float fstep,
|
float *qd,
|
||||||
int32_t *shift1, int32_t lagmin, int32_t lagmax, int32_t lagstep,
|
long np,
|
||||||
float *drift1, int32_t symfac, float *sync, int32_t mode);
|
unsigned char *symbols,
|
||||||
void subtract_signal(float *id, float *qd, long np,
|
float *freq,
|
||||||
float f0, int32_t shift0, float drift0, uint8_t *channel_symbols);
|
int ifmin,
|
||||||
void subtract_signal2(float *id, float *qd, long np,
|
int ifmax,
|
||||||
float f0, int32_t shift0, float drift0, uint8_t *channel_symbols);
|
float fstep,
|
||||||
int32_t wspr_decode(float *idat, float *qdat, uint32_t npoints,
|
int *shift,
|
||||||
struct decoder_options options, struct decoder_results *decodes,
|
int lagmin,
|
||||||
int32_t *n_results);
|
int lagmax,
|
||||||
|
int lagstep,
|
||||||
|
float *drift,
|
||||||
|
int symfac,
|
||||||
|
float *sync,
|
||||||
|
int mode);
|
||||||
|
void subtract_signal(float *id,
|
||||||
|
float *qd,
|
||||||
|
long np,
|
||||||
|
float f0,
|
||||||
|
int shift,
|
||||||
|
float drift,
|
||||||
|
unsigned char *channel_symbols);
|
||||||
|
void subtract_signal2(float *id,
|
||||||
|
float *qd,
|
||||||
|
long np,
|
||||||
|
float f0,
|
||||||
|
int shift,
|
||||||
|
float drift,
|
||||||
|
unsigned char *channel_symbols);
|
||||||
|
int wspr_decode(float *idat,
|
||||||
|
float *qdat,
|
||||||
|
int samples,
|
||||||
|
struct decoder_options options,
|
||||||
|
struct decoder_results *decodes,
|
||||||
|
int *n_results);
|
||||||
|
|
||||||
|
|
|
||||||
|
|
@ -152,7 +152,7 @@ int unpackgrid(int32_t ngrid, char *grid) {
|
||||||
}
|
}
|
||||||
|
|
||||||
int unpackpfx(int32_t nprefix, char *call) {
|
int unpackpfx(int32_t nprefix, char *call) {
|
||||||
char nc, pfx[4] = "", tmpcall[7] = "";
|
char nc, pfx[4] = {'\0'}, tmpcall[7];
|
||||||
int i;
|
int i;
|
||||||
int32_t n;
|
int32_t n;
|
||||||
|
|
||||||
|
|
@ -172,7 +172,8 @@ int unpackpfx(int32_t nprefix, char *call) {
|
||||||
n = n / 37;
|
n = n / 37;
|
||||||
}
|
}
|
||||||
|
|
||||||
strcpy(call, pfx);
|
char *p = strrchr(pfx, ' ');
|
||||||
|
strcpy(call, p ? p + 1 : pfx);
|
||||||
strncat(call, "/", 1);
|
strncat(call, "/", 1);
|
||||||
strncat(call, tmpcall, strlen(tmpcall));
|
strncat(call, tmpcall, strlen(tmpcall));
|
||||||
|
|
||||||
|
|
@ -234,9 +235,9 @@ int floatcomp(const void *elem1, const void *elem2) {
|
||||||
return *(const float *)elem1 > *(const float *)elem2;
|
return *(const float *)elem1 > *(const float *)elem2;
|
||||||
}
|
}
|
||||||
|
|
||||||
int unpk_(signed char *message, char *hashtab, char *call_loc_pow, char *call, char *loc, char *pwr, char *callsign) {
|
int unpk_(signed char *message, char *hashtab, char *loctab, char *call_loc_pow, char *call, char *loc, char *pwr, char *callsign) {
|
||||||
int n1, n2, n3, ndbm, ihash, nadd, noprint = 0;
|
int n1, n2, n3, ndbm, ihash, nadd, noprint = 0;
|
||||||
char grid[5], grid6[7], cdbm[3];
|
char grid[5], grid6[7], cdbm[4];
|
||||||
|
|
||||||
unpack50(message, &n1, &n2);
|
unpack50(message, &n1, &n2);
|
||||||
if (!unpackcall(n1, callsign))
|
if (!unpackcall(n1, callsign))
|
||||||
|
|
@ -250,13 +251,10 @@ int unpk_(signed char *message, char *hashtab, char *call_loc_pow, char *call, c
|
||||||
/*
|
/*
|
||||||
Based on the value of ntype, decide whether this is a Type 1, 2, or
|
Based on the value of ntype, decide whether this is a Type 1, 2, or
|
||||||
3 message.
|
3 message.
|
||||||
|
|
||||||
* Type 1: 6 digit call, grid, power - ntype is positive and is a member
|
* Type 1: 6 digit call, grid, power - ntype is positive and is a member
|
||||||
of the set {0,3,7,10,13,17,20...60}
|
of the set {0,3,7,10,13,17,20...60}
|
||||||
|
|
||||||
* Type 2: extended callsign, power - ntype is positive but not
|
* Type 2: extended callsign, power - ntype is positive but not
|
||||||
a member of the set of allowed powers
|
a member of the set of allowed powers
|
||||||
|
|
||||||
* Type 3: hash, 6 digit grid, power - ntype is negative.
|
* Type 3: hash, 6 digit grid, power - ntype is negative.
|
||||||
*/
|
*/
|
||||||
|
|
||||||
|
|
@ -274,26 +272,23 @@ int unpk_(signed char *message, char *hashtab, char *call_loc_pow, char *call, c
|
||||||
strncat(call_loc_pow, "\0", 1);
|
strncat(call_loc_pow, "\0", 1);
|
||||||
ihash = nhash(callsign, strlen(callsign), (uint32_t)146);
|
ihash = nhash(callsign, strlen(callsign), (uint32_t)146);
|
||||||
strcpy(hashtab + ihash * 13, callsign);
|
strcpy(hashtab + ihash * 13, callsign);
|
||||||
|
strcpy(loctab + ihash * 5, grid);
|
||||||
|
|
||||||
memset(call, 0, strlen(callsign) + 1);
|
memset(call, 0, sizeof(char) * strlen(callsign) + 1);
|
||||||
memset(loc, 0, strlen(grid) + 1);
|
memset(loc, 0, sizeof(char) * strlen(grid) + 1);
|
||||||
memset(pwr, 0, 2 + 1);
|
memset(pwr, 0, sizeof(char) * 2 + 1);
|
||||||
strncat(call, callsign, strlen(callsign));
|
strncat(call, callsign, strlen(callsign));
|
||||||
strncat(call, "\0", 1);
|
strncat(call, "\0", 1);
|
||||||
strncat(loc, grid, strlen(grid));
|
strncat(loc, grid, strlen(grid));
|
||||||
strncat(loc, "\0", 1);
|
strncat(loc, "\0", 1);
|
||||||
strncat(pwr, cdbm, 2);
|
strncat(pwr, cdbm, 2);
|
||||||
strncat(pwr, "\0", 1);
|
strncat(pwr, "\0", 1);
|
||||||
|
|
||||||
} else {
|
} else {
|
||||||
nadd = nu;
|
nadd = nu;
|
||||||
if (nu > 3)
|
if (nu > 3) nadd = nu - 3;
|
||||||
nadd = nu - 3;
|
if (nu > 7) nadd = nu - 7;
|
||||||
if (nu > 7)
|
|
||||||
nadd = nu - 7;
|
|
||||||
n3 = n2 / 128 + 32768 * (nadd - 1);
|
n3 = n2 / 128 + 32768 * (nadd - 1);
|
||||||
if (!unpackpfx(n3, callsign))
|
if (!unpackpfx(n3, callsign)) return 1;
|
||||||
return 1;
|
|
||||||
ndbm = ntype - nadd;
|
ndbm = ntype - nadd;
|
||||||
memset(call_loc_pow, 0, sizeof(char) * 23);
|
memset(call_loc_pow, 0, sizeof(char) * 23);
|
||||||
sprintf(cdbm, "%2d", ndbm);
|
sprintf(cdbm, "%2d", ndbm);
|
||||||
|
|
@ -305,33 +300,23 @@ int unpk_(signed char *message, char *hashtab, char *call_loc_pow, char *call, c
|
||||||
if (nu == 0 || nu == 3 || nu == 7) { // make sure power is OK
|
if (nu == 0 || nu == 3 || nu == 7) { // make sure power is OK
|
||||||
ihash = nhash(callsign, strlen(callsign), (uint32_t)146);
|
ihash = nhash(callsign, strlen(callsign), (uint32_t)146);
|
||||||
strcpy(hashtab + ihash * 13, callsign);
|
strcpy(hashtab + ihash * 13, callsign);
|
||||||
} else {
|
} else
|
||||||
noprint = 1;
|
noprint = 1;
|
||||||
}
|
|
||||||
|
|
||||||
memset(call, 0, strlen(callsign) + 1);
|
|
||||||
memset(loc, 0, 1);
|
|
||||||
memset(pwr, 0, 2 + 1);
|
|
||||||
strncat(call, callsign, strlen(callsign));
|
|
||||||
strncat(call, "\0", 1);
|
|
||||||
strncat(loc, "\0", 1);
|
|
||||||
strncat(pwr, cdbm, 2);
|
|
||||||
strncat(pwr, "\0", 1);
|
|
||||||
}
|
}
|
||||||
} else if (ntype < 0) {
|
} else if (ntype < 0) {
|
||||||
ndbm = -(ntype + 1);
|
ndbm = -(ntype + 1);
|
||||||
memset(grid6, 0, sizeof(char) * 7);
|
memset(grid6, 0, sizeof(char) * 7);
|
||||||
strncat(grid6, callsign + 5, 1);
|
// size_t len=strlen(callsign);
|
||||||
strncat(grid6, callsign, 5);
|
size_t len = 6;
|
||||||
|
strncat(grid6, callsign + len - 1, 1);
|
||||||
|
strncat(grid6, callsign, len - 1);
|
||||||
int nu = ndbm % 10;
|
int nu = ndbm % 10;
|
||||||
if ((nu == 0 || nu == 3 || nu == 7) &&
|
if ((nu != 0 && nu != 3 && nu != 7) ||
|
||||||
(isalpha(grid6[0]) && isalpha(grid6[1]) && isdigit(grid6[2]) && isdigit(grid6[3]))) {
|
!isalpha(grid6[0]) || !isalpha(grid6[1]) ||
|
||||||
|
!isdigit(grid6[2]) || !isdigit(grid6[3])) {
|
||||||
// not testing 4'th and 5'th chars because of this case: <PA0SKT/2> JO33 40
|
// not testing 4'th and 5'th chars because of this case: <PA0SKT/2> JO33 40
|
||||||
// grid is only 4 chars even though this is a hashed callsign...
|
// grid is only 4 chars even though this is a hashed callsign...
|
||||||
// isalpha(grid6[4]) && isalpha(grid6[5]) ) ) {
|
// isalpha(grid6[4]) && isalpha(grid6[5]) ) ) {
|
||||||
ihash = nhash(callsign, strlen(callsign), (uint32_t)146);
|
|
||||||
strcpy(hashtab + ihash * 13, callsign);
|
|
||||||
} else {
|
|
||||||
noprint = 1;
|
noprint = 1;
|
||||||
}
|
}
|
||||||
|
|
||||||
|
|
@ -351,9 +336,9 @@ int unpk_(signed char *message, char *hashtab, char *call_loc_pow, char *call, c
|
||||||
strncat(call_loc_pow, cdbm, 2);
|
strncat(call_loc_pow, cdbm, 2);
|
||||||
strncat(call_loc_pow, "\0", 1);
|
strncat(call_loc_pow, "\0", 1);
|
||||||
|
|
||||||
memset(call, 0, strlen(callsign) + 1);
|
memset(call, 0, sizeof(char) * strlen(callsign) + 1);
|
||||||
memset(loc, 0, strlen(grid6) + 1);
|
memset(loc, 0, sizeof(char) * strlen(grid6) + 1);
|
||||||
memset(pwr, 0, 2 + 1);
|
memset(pwr, 0, sizeof(char) * 2 + 1);
|
||||||
strncat(call, callsign, strlen(callsign));
|
strncat(call, callsign, strlen(callsign));
|
||||||
strncat(call, "\0", 1);
|
strncat(call, "\0", 1);
|
||||||
strncat(loc, grid6, strlen(grid6));
|
strncat(loc, grid6, strlen(grid6));
|
||||||
|
|
@ -362,7 +347,7 @@ int unpk_(signed char *message, char *hashtab, char *call_loc_pow, char *call, c
|
||||||
strncat(pwr, "\0", 1);
|
strncat(pwr, "\0", 1);
|
||||||
|
|
||||||
// I don't know what to do with these... They show up as "A000AA" grids.
|
// I don't know what to do with these... They show up as "A000AA" grids.
|
||||||
if (ntype == -64)
|
if (ntype == -64)
|
||||||
noprint = 1;
|
noprint = 1;
|
||||||
}
|
}
|
||||||
return noprint;
|
return noprint;
|
||||||
|
|
|
||||||
|
|
@ -39,4 +39,4 @@ void deinterleave(unsigned char *sym);
|
||||||
int doublecomp(const void *elem1, const void *elem2);
|
int doublecomp(const void *elem1, const void *elem2);
|
||||||
int floatcomp(const void *elem1, const void *elem2);
|
int floatcomp(const void *elem1, const void *elem2);
|
||||||
|
|
||||||
int unpk_(signed char *message, char *hashtab, char *call_loc_pow, char *call, char *loc, char *pwr, char *callsign);
|
int unpk_(signed char *message, char *hashtab, char *loctab, char *call_loc_pow, char *call, char *loc, char *pwr, char *callsign);
|
||||||
|
|
|
||||||
|
|
@ -38,7 +38,7 @@ char get_callsign_character_code(char ch) {
|
||||||
return -1;
|
return -1;
|
||||||
}
|
}
|
||||||
|
|
||||||
long unsigned int pack_grid4_power(char *grid4, int power) {
|
long unsigned int pack_grid4_power(char const *grid4, int power) {
|
||||||
long unsigned int m;
|
long unsigned int m;
|
||||||
|
|
||||||
m = (179 - 10 * grid4[0] - grid4[2]) * 180 + 10 * grid4[1] + grid4[3];
|
m = (179 - 10 * grid4[0] - grid4[2]) * 180 + 10 * grid4[1] + grid4[3];
|
||||||
|
|
@ -46,31 +46,23 @@ long unsigned int pack_grid4_power(char *grid4, int power) {
|
||||||
return m;
|
return m;
|
||||||
}
|
}
|
||||||
|
|
||||||
long unsigned int pack_call(char *callsign) {
|
long unsigned int pack_call(char const *callsign) {
|
||||||
int i;
|
unsigned int i;
|
||||||
long unsigned int n;
|
long unsigned int n;
|
||||||
char call6[6];
|
char call6[6];
|
||||||
memset(call6, 32, sizeof(char) * 6);
|
memset(call6, ' ', sizeof(call6));
|
||||||
// callsign is 6 characters in length. Exactly.
|
// callsign is 6 characters in length. Exactly.
|
||||||
int call_len = strlen(callsign);
|
size_t call_len = strlen(callsign);
|
||||||
if (call_len > 6) {
|
if (call_len > 6) {
|
||||||
return 0;
|
return 0;
|
||||||
}
|
}
|
||||||
if (isdigit(*(callsign + 2))) {
|
if (isdigit(callsign[2])) {
|
||||||
for (i = 0; i < 6; i++) {
|
for (i = 0; i < call_len; i++) {
|
||||||
if (callsign[i] == 0) {
|
call6[i] = callsign[i];
|
||||||
call6[i] = 32;
|
|
||||||
} else {
|
|
||||||
call6[i] = *(callsign + i);
|
|
||||||
}
|
|
||||||
}
|
}
|
||||||
} else if (isdigit(*(callsign + 1))) {
|
} else if (isdigit(callsign[1])) {
|
||||||
for (i = 0; i < 6; i++) {
|
for (i = 1; i < call_len + 1; i++) {
|
||||||
if (i == 0 || callsign[i - 1] == 0) {
|
call6[i] = callsign[i - 1];
|
||||||
call6[i] = 32;
|
|
||||||
} else {
|
|
||||||
call6[i] = *(callsign + i - 1);
|
|
||||||
}
|
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
for (i = 0; i < 6; i++) {
|
for (i = 0; i < 6; i++) {
|
||||||
|
|
@ -86,17 +78,16 @@ long unsigned int pack_call(char *callsign) {
|
||||||
}
|
}
|
||||||
|
|
||||||
void pack_prefix(char *callsign, int32_t *n, int32_t *m, int32_t *nadd) {
|
void pack_prefix(char *callsign, int32_t *n, int32_t *m, int32_t *nadd) {
|
||||||
int i;
|
size_t i;
|
||||||
char *call6;
|
char *call6 = calloc(7, sizeof(char));
|
||||||
call6 = malloc(sizeof(char) * 6);
|
size_t i1 = strcspn(callsign, "/");
|
||||||
memset(call6, 32, sizeof(char) * 6);
|
|
||||||
int i1 = strcspn(callsign, "/");
|
|
||||||
|
|
||||||
if (callsign[i1 + 2] == 0) {
|
if (callsign[i1 + 2] == 0) {
|
||||||
// single char suffix
|
// single char suffix
|
||||||
for (i = 0; i < i1; i++) {
|
for (i = 0; i < i1; i++) {
|
||||||
call6[i] = callsign[i];
|
call6[i] = callsign[i];
|
||||||
}
|
}
|
||||||
|
call6[i] = '\0';
|
||||||
*n = pack_call(call6);
|
*n = pack_call(call6);
|
||||||
*nadd = 1;
|
*nadd = 1;
|
||||||
int nc = callsign[i1 + 1];
|
int nc = callsign[i1 + 1];
|
||||||
|
|
@ -118,11 +109,10 @@ void pack_prefix(char *callsign, int32_t *n, int32_t *m, int32_t *nadd) {
|
||||||
*m = 10 * (callsign[i1 + 1] - 48) + (callsign[i1 + 2] - 48);
|
*m = 10 * (callsign[i1 + 1] - 48) + (callsign[i1 + 2] - 48);
|
||||||
*m = 60000 + 26 + *m;
|
*m = 60000 + 26 + *m;
|
||||||
} else {
|
} else {
|
||||||
char *pfx = strtok(callsign, "/");
|
char const *pfx = strtok(callsign, "/");
|
||||||
*n = pack_call(strtok(NULL, " "));
|
char const *call = strtok(NULL, " ");
|
||||||
//call6 = strtok(NULL, " ");
|
*n = pack_call(call);
|
||||||
//*n = pack_call(call6);
|
size_t plen = strlen(pfx);
|
||||||
int plen = strlen(pfx);
|
|
||||||
if (plen == 1) {
|
if (plen == 1) {
|
||||||
*m = 36;
|
*m = 36;
|
||||||
*m = 37 * (*m) + 36;
|
*m = 37 * (*m) + 36;
|
||||||
|
|
@ -148,6 +138,7 @@ void pack_prefix(char *callsign, int32_t *n, int32_t *m, int32_t *nadd) {
|
||||||
*nadd = 1;
|
*nadd = 1;
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
free(call6);
|
||||||
}
|
}
|
||||||
|
|
||||||
void interleave(unsigned char *sym) {
|
void interleave(unsigned char *sym) {
|
||||||
|
|
@ -169,10 +160,11 @@ void interleave(unsigned char *sym) {
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
int get_wspr_channel_symbols(char *rawmessage, char *hashtab, unsigned char *symbols) {
|
int get_wspr_channel_symbols(char *rawmessage, char *hashtab, char *loctab, unsigned char *symbols) {
|
||||||
int m = 0, n = 0, ntype = 0;
|
int m = 0, ntype = 0;
|
||||||
|
long unsigned int n = 0;
|
||||||
int i, j, ihash;
|
int i, j, ihash;
|
||||||
unsigned char pr3[162] = {
|
unsigned char pr3vector[162] = {
|
||||||
1, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0,
|
1, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 1, 0,
|
||||||
0, 1, 0, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1,
|
0, 1, 0, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 1,
|
||||||
0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 1, 1, 0, 1, 0, 0, 0, 1,
|
0, 0, 0, 0, 0, 0, 1, 0, 1, 1, 0, 0, 1, 1, 0, 1, 0, 0, 0, 1,
|
||||||
|
|
@ -181,7 +173,8 @@ int get_wspr_channel_symbols(char *rawmessage, char *hashtab, unsigned char *sym
|
||||||
0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 1, 1, 0, 1, 1, 0, 0, 1, 1,
|
0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 1, 1, 0, 1, 1, 0, 0, 1, 1,
|
||||||
0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 1,
|
0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 1,
|
||||||
0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 1, 1, 0, 0, 0, 1, 1, 0,
|
0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 1, 1, 0, 0, 0, 1, 1, 0,
|
||||||
0, 0};
|
0, 0 };
|
||||||
|
|
||||||
int nu[10] = {0, -1, 1, 0, -1, 2, 1, 0, -1, 1};
|
int nu[10] = {0, -1, 1, 0, -1, 2, 1, 0, -1, 1};
|
||||||
char *callsign, *grid, *powstr;
|
char *callsign, *grid, *powstr;
|
||||||
char grid4[5], message[23];
|
char grid4[5], message[23];
|
||||||
|
|
@ -193,16 +186,16 @@ int get_wspr_channel_symbols(char *rawmessage, char *hashtab, unsigned char *sym
|
||||||
i++;
|
i++;
|
||||||
}
|
}
|
||||||
|
|
||||||
int i1 = strcspn(message, " ");
|
size_t i1 = strcspn(message, " ");
|
||||||
int i2 = strcspn(message, "/");
|
size_t i2 = strcspn(message, "/");
|
||||||
int i3 = strcspn(message, "<");
|
size_t i3 = strcspn(message, "<");
|
||||||
int i4 = strcspn(message, ">");
|
size_t i4 = strcspn(message, ">");
|
||||||
int mlen = strlen(message);
|
size_t mlen = strlen(message);
|
||||||
|
|
||||||
// Use the presence and/or absence of "<" and "/" to decide what
|
// Use the presence and/or absence of "<" and "/" to decide what
|
||||||
// type of message. No sanity checks! Beware!
|
// type of message. No sanity checks! Beware!
|
||||||
|
|
||||||
if ((i1 > 3) & (i1 < 7) & (i2 == mlen) & (i3 == mlen)) {
|
if (i1 > 3 && i1 < 7 && i2 == mlen && i3 == mlen) {
|
||||||
// Type 1 message: K9AN EN50 33
|
// Type 1 message: K9AN EN50 33
|
||||||
// xxnxxxx xxnn nn
|
// xxnxxxx xxnn nn
|
||||||
callsign = strtok(message, " ");
|
callsign = strtok(message, " ");
|
||||||
|
|
@ -232,8 +225,8 @@ int get_wspr_channel_symbols(char *rawmessage, char *hashtab, unsigned char *sym
|
||||||
ihash = nhash(callsign, strlen(callsign), (uint32_t)146);
|
ihash = nhash(callsign, strlen(callsign), (uint32_t)146);
|
||||||
m = 128 * ihash + ntype + 64;
|
m = 128 * ihash + ntype + 64;
|
||||||
|
|
||||||
char grid6[6];
|
char grid6[7];
|
||||||
memset(grid6, 32, sizeof(char) * 6);
|
memset(grid6, 0, sizeof(char) * 7);
|
||||||
j = strlen(grid);
|
j = strlen(grid);
|
||||||
for (i = 0; i < j - 1; i++) {
|
for (i = 0; i < j - 1; i++) {
|
||||||
grid6[i] = grid[i + 1];
|
grid6[i] = grid[i + 1];
|
||||||
|
|
@ -243,7 +236,7 @@ int get_wspr_channel_symbols(char *rawmessage, char *hashtab, unsigned char *sym
|
||||||
} else if (i2 < mlen) { // just looks for a right slash
|
} else if (i2 < mlen) { // just looks for a right slash
|
||||||
// Type 2: PJ4/K1ABC 37
|
// Type 2: PJ4/K1ABC 37
|
||||||
callsign = strtok(message, " ");
|
callsign = strtok(message, " ");
|
||||||
if (strlen(callsign) < i2) return 0; // guards against pathological case
|
if (i2 == 0 || i2 > strlen(callsign)) return 0; // guards against pathological case
|
||||||
powstr = strtok(NULL, " ");
|
powstr = strtok(NULL, " ");
|
||||||
int power = atoi(powstr);
|
int power = atoi(powstr);
|
||||||
if (power < 0) power = 0;
|
if (power < 0) power = 0;
|
||||||
|
|
@ -283,8 +276,6 @@ int get_wspr_channel_symbols(char *rawmessage, char *hashtab, unsigned char *sym
|
||||||
// make sure that the 11-byte data vector is unpackable
|
// make sure that the 11-byte data vector is unpackable
|
||||||
// unpack it with the routine that the decoder will use and display
|
// unpack it with the routine that the decoder will use and display
|
||||||
// the result. let the operator decide whether it worked.
|
// the result. let the operator decide whether it worked.
|
||||||
// char hashtab[32768][13];
|
|
||||||
// memset(hashtab,0,sizeof(char)*32768*13);
|
|
||||||
|
|
||||||
char *check_call_loc_pow, *check_callsign, *call, *loc, *pwr;
|
char *check_call_loc_pow, *check_callsign, *call, *loc, *pwr;
|
||||||
check_call_loc_pow = malloc(sizeof(char) * 23);
|
check_call_loc_pow = malloc(sizeof(char) * 23);
|
||||||
|
|
@ -294,21 +285,23 @@ int get_wspr_channel_symbols(char *rawmessage, char *hashtab, unsigned char *sym
|
||||||
pwr = malloc(sizeof(char) * 3);
|
pwr = malloc(sizeof(char) * 3);
|
||||||
signed char check_data[11];
|
signed char check_data[11];
|
||||||
memcpy(check_data, data, sizeof(char) * 11);
|
memcpy(check_data, data, sizeof(char) * 11);
|
||||||
unpk_(check_data, hashtab, check_call_loc_pow, call, loc, pwr, check_callsign);
|
|
||||||
|
unpk_(check_data, hashtab, loctab, check_call_loc_pow, call, loc, pwr, check_callsign);
|
||||||
// printf("Will decode as: %s\n",check_call_loc_pow);
|
// printf("Will decode as: %s\n",check_call_loc_pow);
|
||||||
|
|
||||||
unsigned int nbytes = 11; // The message with tail is packed into 11 bytes.
|
unsigned int nbytes = 11; // The message with tail is packed into almost 11 bytes.
|
||||||
unsigned int nencoded = (nbytes * 2 * 8); // This is how much encode() writes
|
unsigned char channelbits[nbytes * 8 * 2]; /* 162 rounded up */
|
||||||
unsigned char channelbits[nencoded];
|
memset(channelbits, 0, sizeof(char) * nbytes * 8 * 2);
|
||||||
memset(channelbits, 0, sizeof(char) * nencoded);
|
|
||||||
|
|
||||||
encode(channelbits, data, nbytes);
|
encode(channelbits, data, nbytes);
|
||||||
|
|
||||||
interleave(channelbits);
|
interleave(channelbits);
|
||||||
|
|
||||||
for (i = 0; i < 162; i++) {
|
for (i = 0; i < 162; i++) {
|
||||||
symbols[i] = 2 * channelbits[i] + pr3[i];
|
symbols[i] = 2 * channelbits[i] + pr3vector[i];
|
||||||
}
|
}
|
||||||
|
free(check_call_loc_pow);
|
||||||
|
free(check_callsign);
|
||||||
return 1;
|
return 1;
|
||||||
}
|
}
|
||||||
|
|
||||||
|
|
|
||||||
|
|
@ -2,8 +2,8 @@
|
||||||
|
|
||||||
char get_locator_character_code(char ch);
|
char get_locator_character_code(char ch);
|
||||||
char get_callsign_character_code(char ch);
|
char get_callsign_character_code(char ch);
|
||||||
long unsigned int pack_grid4_power(char *grid4, int power);
|
long unsigned int pack_grid4_power(char const *grid4, int power);
|
||||||
long unsigned int pack_call(char *callsign);
|
long unsigned int pack_call(char const *callsign);
|
||||||
void pack_prefix(char *callsign, int32_t *n, int32_t *m, int32_t *nadd);
|
void pack_prefix(char *callsign, int32_t *n, int32_t *m, int32_t *nadd);
|
||||||
void interleave(unsigned char *sym);
|
void interleave(unsigned char *sym);
|
||||||
int get_wspr_channel_symbols(char *message, char *hashtab, unsigned char *symbols);
|
int get_wspr_channel_symbols(char *rawmessage, char *hashtab, char *loctab, unsigned char *symbols);
|
||||||
|
|
|
||||||
Ładowanie…
Reference in New Issue