diff options
Diffstat (limited to 'tests/coding/coding_test.c')
| -rw-r--r-- | tests/coding/coding_test.c | 501 |
1 files changed, 501 insertions, 0 deletions
diff --git a/tests/coding/coding_test.c b/tests/coding/coding_test.c new file mode 100644 index 00000000..9a00f0d3 --- /dev/null +++ b/tests/coding/coding_test.c @@ -0,0 +1,501 @@ +/* + * (C) 2013 by Andreas Eversberg <jolly@eversberg.eu> + * (C) 2015 by Alexander Chemeris <Alexander.Chemeris@fairwaves.co> + * + * All Rights Reserved + * + * This program is free software; you can redistribute it and/or modify + * it under the terms of the GNU General Public License as published by + * the Free Software Foundation; either version 2 of the License, or + * (at your option) any later version. + * + * This program is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the + * GNU General Public License for more details. + * + * You should have received a copy of the GNU General Public License along + * with this program; if not, write to the Free Software Foundation, Inc., + * 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA. + */ + +#include <stdio.h> +#include <errno.h> +#include <unistd.h> +#include <string.h> +#include <stdlib.h> + +#include <osmocom/core/bits.h> +#include <osmocom/core/utils.h> + +#include <osmocom/coding/gsm0503_coding.h> + +#define ASSERT_TRUE(rc) \ + if (!(rc)) { \ + printf("Assert failed in %s:%d.\n", \ + __FILE__, __LINE__); \ + abort(); \ + } + +/* set condition to 1, to show debugging */ +#define printd if (0) printf + +static int ubits2sbits(ubit_t *ubits, sbit_t *sbits, int count) +{ + int i; + + for (i = 0; i < count; i++) { + if (*ubits == 0x23) { + ubits++; + sbits++; + continue; + } + + if ((*ubits++) & 1) { + *sbits++ = -127; + } else { + *sbits++ = 127; + } + } + + return count; +} + +static void test_xcch(uint8_t *l2) +{ + uint8_t result[23]; + ubit_t bursts_u[116 * 4]; + sbit_t bursts_s[116 * 4]; + int n_errors, n_bits_total; + + /* Encode L2 message */ + printd("Encoding: %s\n", osmo_hexdump(l2, 23)); + gsm0503_xcch_encode(bursts_u, l2); + + /* Prepare soft-bits */ + ubits2sbits(bursts_u, bursts_s, 116 * 4); + + printd("U-Bits:\n"); + printd("%s %02x %02x ", osmo_hexdump(bursts_u, 57), + bursts_u[57], bursts_u[58]); + printd("%s\n", osmo_hexdump(bursts_u + 59, 57)); + printd("%s %02x %02x ", osmo_hexdump(bursts_u + 116, 57), + bursts_u[57 + 116], bursts_u[58 + 116]); + printd("%s\n", osmo_hexdump(bursts_u + 59 + 116, 57)); + printd("%s %02x %02x ", osmo_hexdump(bursts_u + 232, 57), + bursts_u[57 + 232], bursts_u[58 + 232]); + printd("%s\n", osmo_hexdump(bursts_u + 59 + 232, 57)); + printd("%s %02x %02x ", osmo_hexdump(bursts_u + 348, 57), + bursts_u[57 + 348], bursts_u[58 + 348]); + printd("%s\n", osmo_hexdump(bursts_u + 59 + 348, 57)); + + printd("S-Bits:\n"); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s, 57), + (uint8_t)bursts_s[57], (uint8_t)bursts_s[58]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59, 57)); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s + 116, 57), + (uint8_t)bursts_s[57 + 116], (uint8_t)bursts_s[58 + 116]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59 + 116, 57)); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s + 232, 57), + (uint8_t)bursts_s[57 + 232], (uint8_t)bursts_s[58 + 232]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59 + 232, 57)); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s + 348, 57), + (uint8_t)bursts_s[57 + 348], (uint8_t)bursts_s[58 + 348]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59 + 348, 57)); + + /* Destroy some bits */ + memset(bursts_s, 0, 30); + memset(bursts_s + 116, 0, 30); + + /* Decode, correcting errors */ + gsm0503_xcch_decode(result, bursts_s, &n_errors, &n_bits_total); + printd("Decoded: %s\n", osmo_hexdump(result, 23)); + printf("xcch_decode: n_errors=%d n_bits_total=%d ber=%.2f\n", + n_errors, n_bits_total, (float) n_errors / n_bits_total); + + ASSERT_TRUE(n_bits_total == 456); + ASSERT_TRUE(!memcmp(l2, result, 23)); + + printd("\n"); +} + +static void test_rach(uint8_t bsic, uint8_t ra) +{ + uint8_t result; + ubit_t bursts_u[36]; + sbit_t bursts_s[36]; + + /* Encode L2 message */ + printd("Encoding: %02x\n", ra); + gsm0503_rach_encode(bursts_u, &ra, bsic); + + /* Prepare soft-bits */ + ubits2sbits(bursts_u, bursts_s, 36); + + printd("U-Bits:\n"); + printd("%s\n", osmo_hexdump(bursts_u, 36)); + + printd("S-Bits:\n"); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s, 36)); + + /* Destroy some bits */ + memset(bursts_s + 6, 0, 8); + + /* Decode, correcting errors */ + gsm0503_rach_decode(&result, bursts_s, bsic); + printd("Decoded: %02x\n", result); + + ASSERT_TRUE(ra == result); + + printd("\n"); +} + +static void test_sch(uint8_t *info) +{ + uint8_t result[4]; + ubit_t bursts_u[78]; + sbit_t bursts_s[78]; + + /* Zero bits 25 and above */ + info[3] &= 1; + result[3] = 0; + + /* Encode L2 message */ + printd("Encoding: %s\n", osmo_hexdump(info, 4)); + gsm0503_sch_encode(bursts_u, info); + + /* Prepare soft-bits */ + ubits2sbits(bursts_u, bursts_s, 78); + + printd("U-Bits:\n"); + printd("%s\n", osmo_hexdump(bursts_u, 78)); + + printd("S-Bits:\n"); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s, 78)); + + /* Destroy some bits */ + memset(bursts_s + 6, 0, 10); + + /* Decode, correcting errors */ + gsm0503_sch_decode(result, bursts_s); + printd("Decoded: %s\n", osmo_hexdump(result, 4)); + + ASSERT_TRUE(!memcmp(info, result, 4)); + + printd("\n"); +} + +static void test_fr(uint8_t *speech, int len) +{ + uint8_t result[33]; + ubit_t bursts_u[116 * 8]; + sbit_t bursts_s[116 * 8]; + int n_errors, n_bits_total; + int rc; + + memset(bursts_u, 0x23, sizeof(bursts_u)); + memset(bursts_s, 0, sizeof(bursts_s)); + + /* Encode L2 message */ + printd("Encoding: %s\n", osmo_hexdump(speech, len)); + gsm0503_tch_fr_encode(bursts_u, speech, len, 1); + + /* Prepare soft-bits */ + ubits2sbits(bursts_u, bursts_s, 116 * 8); + + printd("U-Bits:\n"); + printd("%s %02x %02x ", osmo_hexdump(bursts_u, 57), + bursts_u[57], bursts_u[58]); + printd("%s\n", osmo_hexdump(bursts_u + 59, 57)); + printd("%s %02x %02x ", osmo_hexdump(bursts_u + 116, 57), + bursts_u[57 + 116], bursts_u[58 + 116]); + printd("%s\n", osmo_hexdump(bursts_u + 59 + 116, 57)); + printd("%s %02x %02x ", osmo_hexdump(bursts_u + 232, 57), + bursts_u[57 + 232], bursts_u[58 + 232]); + printd("%s\n", osmo_hexdump(bursts_u + 59 + 232, 57)); + printd("%s %02x %02x ", osmo_hexdump(bursts_u + 348, 57), + bursts_u[57 + 348], bursts_u[58 + 348]); + printd("%s\n", osmo_hexdump(bursts_u + 59 + 348, 57)); + printd("%s %02x %02x ", osmo_hexdump(bursts_u + 464, 57), + bursts_u[57 + 464], bursts_u[58 + 464]); + printd("%s\n", osmo_hexdump(bursts_u + 59 + 464, 57)); + printd("%s %02x %02x ", osmo_hexdump(bursts_u + 580, 57), + bursts_u[57 + 580], bursts_u[58 + 580]); + printd("%s\n", osmo_hexdump(bursts_u + 59 + 580, 57)); + printd("%s %02x %02x ", osmo_hexdump(bursts_u + 696, 57), + bursts_u[57 + 696], bursts_u[58 + 696]); + printd("%s\n", osmo_hexdump(bursts_u + 59 + 696, 57)); + printd("%s %02x %02x ", osmo_hexdump(bursts_u + 812, 57), + bursts_u[57 + 812], bursts_u[58 + 812]); + printd("%s\n", osmo_hexdump(bursts_u + 59 + 812, 57)); + + printd("S-Bits:\n"); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s, 57), + (uint8_t)bursts_s[57], (uint8_t)bursts_s[58]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59, 57)); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s + 116, 57), + (uint8_t)bursts_s[57 + 116], (uint8_t)bursts_s[58 + 116]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59 + 116, 57)); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s + 232, 57), + (uint8_t)bursts_s[57 + 232], (uint8_t)bursts_s[58 + 232]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59 + 232, 57)); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s + 348, 57), + (uint8_t)bursts_s[57 + 348], (uint8_t)bursts_s[58 + 348]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59 + 348, 57)); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s + 464, 57), + (uint8_t)bursts_s[57 + 464], (uint8_t)bursts_s[58 + 464]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59 + 464, 57)); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s + 580, 57), + (uint8_t)bursts_s[57 + 580], (uint8_t)bursts_s[58 + 580]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59 + 580, 57)); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s + 696, 57), + (uint8_t)bursts_s[57 + 696], (uint8_t)bursts_s[58 + 696]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59 + 696, 57)); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s + 812, 57), + (uint8_t)bursts_s[57 + 812], (uint8_t)bursts_s[58 + 812]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59 + 812, 57)); + + /* Destroy some bits */ + memset(bursts_s + 6, 0, 20); + + /* Decode, correcting errors */ + rc = gsm0503_tch_fr_decode(result, bursts_s, 1, len == 31, + &n_errors, &n_bits_total); + printd("Decoded: %s\n", osmo_hexdump(result, len)); + printf("tch_fr_decode: n_errors=%d n_bits_total=%d ber=%.2f\n", + n_errors, n_bits_total, (float)n_errors/n_bits_total); + + ASSERT_TRUE(rc == len); + ASSERT_TRUE(!memcmp(speech, result, len)); + + printd("\n"); +} + +static void test_hr(uint8_t *speech, int len) +{ + uint8_t result[23]; + ubit_t bursts_u[116 * 6]; + sbit_t bursts_s[116 * 6]; + int n_errors, n_bits_total; + int rc; + + memset(bursts_u, 0x23, sizeof(bursts_u)); + memset(bursts_s, 0, sizeof(bursts_s)); + + /* Encode L2 message */ + printd("Encoding: %s\n", osmo_hexdump(speech, len)); + gsm0503_tch_hr_encode(bursts_u, speech, len); + + /* Prepare soft-bits */ + ubits2sbits(bursts_u, bursts_s, 116 * 6); + + printd("U-Bits:\n"); + printd("%s %02x %02x ", osmo_hexdump(bursts_u, 57), + bursts_u[57], bursts_u[58]); + printd("%s\n", osmo_hexdump(bursts_u + 59, 57)); + printd("%s %02x %02x ", osmo_hexdump(bursts_u + 116, 57), + bursts_u[57 + 116], bursts_u[58 + 116]); + printd("%s\n", osmo_hexdump(bursts_u + 59 + 116, 57)); + printd("%s %02x %02x ", osmo_hexdump(bursts_u + 232, 57), + bursts_u[57 + 232], bursts_u[58 + 232]); + printd("%s\n", osmo_hexdump(bursts_u + 59 + 232, 57)); + printd("%s %02x %02x ", osmo_hexdump(bursts_u + 348, 57), + bursts_u[57 + 348], bursts_u[58 + 348]); + printd("%s\n", osmo_hexdump(bursts_u + 59 + 348, 57)); + printd("%s %02x %02x ", osmo_hexdump(bursts_u + 464, 57), + bursts_u[57 + 464], bursts_u[58 + 464]); + printd("%s\n", osmo_hexdump(bursts_u + 59 + 464, 57)); + printd("%s %02x %02x ", osmo_hexdump(bursts_u + 580, 57), + bursts_u[57 + 580], bursts_u[58 + 580]); + printd("%s\n", osmo_hexdump(bursts_u + 59 + 580, 57)); + + printd("S-Bits:\n"); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s, 57), + (uint8_t)bursts_s[57], (uint8_t)bursts_s[58]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59, 57)); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s + 116, 57), + (uint8_t)bursts_s[57 + 116], (uint8_t)bursts_s[58 + 116]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59 + 116, 57)); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s + 232, 57), + (uint8_t)bursts_s[57 + 232], (uint8_t)bursts_s[58 + 232]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59 + 232, 57)); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s + 348, 57), + (uint8_t)bursts_s[57 + 348], (uint8_t)bursts_s[58 + 348]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59 + 348, 57)); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s + 464, 57), + (uint8_t)bursts_s[57 + 464], (uint8_t)bursts_s[58 + 464]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59 + 464, 57)); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s + 580, 57), + (uint8_t)bursts_s[57 + 580], (uint8_t)bursts_s[58 + 580]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59 + 580, 57)); + + /* Destroy some bits */ + memset(bursts_s + 6, 0, 20); + + /* Decode, correcting errors */ + rc = gsm0503_tch_hr_decode(result, bursts_s, 0, + &n_errors, &n_bits_total); + printd("Decoded: %s\n", osmo_hexdump(result, len)); + printf("tch_hr_decode: n_errors=%d n_bits_total=%d ber=%.2f\n", + n_errors, n_bits_total, (float)n_errors/n_bits_total); + + ASSERT_TRUE(rc == len); + ASSERT_TRUE(!memcmp(speech, result, len)); + + printd("\n"); +} + +static void test_pdtch(uint8_t *l2, int len) +{ + uint8_t result[len]; + ubit_t bursts_u[116 * 4]; + sbit_t bursts_s[116 * 4]; + int n_errors, n_bits_total; + int rc; + + /* Zero the not coded tail bits */ + switch (len) { + case 34: + case 54: + l2[len - 1] &= 0x7f; + result[len - 1] &= 0x7f; + break; + case 40: + l2[len - 1] &= 0x07; + result[len - 1] &= 0x07; + break; + } + + /* Encode L2 message */ + printd("Encoding: %s\n", osmo_hexdump(l2, len)); + gsm0503_pdtch_encode(bursts_u, l2, len); + + /* Prepare soft-bits */ + ubits2sbits(bursts_u, bursts_s, 116 * 4); + + printd("U-Bits:\n"); + printd("%s %02x %02x ", osmo_hexdump(bursts_u, 57), + bursts_u[57], bursts_u[58]); + printd("%s\n", osmo_hexdump(bursts_u + 59, 57)); + printd("%s %02x %02x ", osmo_hexdump(bursts_u + 116, 57), + bursts_u[57 + 116], bursts_u[58 + 116]); + printd("%s\n", osmo_hexdump(bursts_u + 59 + 116, 57)); + printd("%s %02x %02x ", osmo_hexdump(bursts_u + 232, 57), + bursts_u[57 + 232], bursts_u[58 + 232]); + printd("%s\n", osmo_hexdump(bursts_u + 59 + 232, 57)); + printd("%s %02x %02x ", osmo_hexdump(bursts_u + 348, 57), + bursts_u[57 + 348], bursts_u[58 + 348]); + printd("%s\n", osmo_hexdump(bursts_u + 59 + 348, 57)); + + printd("S-Bits:\n"); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s, 57), + (uint8_t)bursts_s[57], (uint8_t)bursts_s[58]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59, 57)); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s + 116, 57), + (uint8_t)bursts_s[57 + 116], (uint8_t)bursts_s[58 + 116]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59 + 116, 57)); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s + 232, 57), + (uint8_t)bursts_s[57 + 232], (uint8_t)bursts_s[58 + 232]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59 + 232, 57)); + printd("%s %02x %02x ", osmo_hexdump((uint8_t *)bursts_s + 348, 57), + (uint8_t)bursts_s[57 + 348], (uint8_t)bursts_s[58 + 348]); + printd("%s\n", osmo_hexdump((uint8_t *)bursts_s + 59 + 348, 57)); + + /* Decode */ + rc = gsm0503_pdtch_decode(result, bursts_s, NULL, + &n_errors, &n_bits_total); + printd("Decoded: %s\n", osmo_hexdump(result, len)); + printf("pdtch_decode: n_errors=%d n_bits_total=%d ber=%.2f\n", + n_errors, n_bits_total, (float)n_errors/n_bits_total); + + ASSERT_TRUE(rc == len); + ASSERT_TRUE(!memcmp(l2, result, len)); + + printd("\n"); +} + +uint8_t test_l2[][23] = { + /* Dummy frame */ + { 0x03, 0x03, 0x01, 0x00, 0x00, 0x00, 0x00, 0x00, + 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, + 0x00, 0x00, 0x00, 0x00, 0x00, 0x00, 0x00 }, + /* Random frame */ + { 0xa3, 0xaf, 0x5f, 0xc6, 0x36, 0x43, 0x44, 0xab, + 0xd9, 0x6d, 0x7d, 0x62, 0x24, 0xc9, 0xd2, 0x92, + 0xfa, 0x27, 0x5d, 0x71, 0x7a, 0x59, 0xa8 }, + /* jolly frame */ + { 0x01, 0x02, 0x03, 0x04, 0x05, 0x06, 0x07, 0x08, + 0x09, 0x0a, 0x0b, 0x0c, 0x0d, 0x0e, 0x0f, 0x10, + 0x11, 0x12, 0x13, 0x14, 0x15, 0x16, 0x17 }, +}; + +uint8_t test_macblock[][54] = { + /* Random frame */ + { 0xa3, 0xaf, 0x5f, 0xc6, 0x36, 0x43, 0x44, 0xab, + 0xd9, 0x6d, 0x7d, 0x62, 0x24, 0xc9, 0xd2, 0x92, + 0xfa, 0x27, 0x5d, 0x71, 0x7a, 0x59, 0xa8, 0x42, + 0xa3, 0xaf, 0x5f, 0xc6, 0x36, 0x43, 0x44, 0xab, + 0xa3, 0xaf, 0x5f, 0xc6, 0x36, 0x43, 0x44, 0xab, + 0xd9, 0x6d, 0x7d, 0x62, 0x24, 0xc9, 0xd2, 0x92, + 0xfa, 0x27, 0x5d, 0x71, 0x7a, 0xa8 }, + /* jolly frame */ + { 0x01, 0x02, 0x03, 0x04, 0x05, 0x06, 0x07, 0x08, + 0x09, 0x0a, 0x0b, 0x0c, 0x0d, 0x0e, 0x0f, 0x10, + 0x11, 0x12, 0x13, 0x14, 0x15, 0x16, 0x17 }, +}; + +uint8_t test_speech_fr[33]; +uint8_t test_speech_efr[31]; +uint8_t test_speech_hr[15]; + +int main(int argc, char **argv) +{ + int i, len_l2, len_mb; + + len_l2 = sizeof(test_l2) / sizeof(test_l2[0]); + len_mb = sizeof(test_macblock) / sizeof(test_macblock[0]); + + for (i = 0; i < len_l2; i++) + test_xcch(test_l2[i]); + + for (i = 0; i < 256; i++) { + test_rach(0x3f, i); + test_rach(0x00, i); + test_rach(0x1a, i); + } + + for (i = 0; i < len_l2; i++) + test_sch(test_l2[i]); + + for (i = 0; i < sizeof(test_speech_fr); i++) + test_speech_fr[i] = i; + test_speech_fr[0] = 0xd0; + test_fr(test_speech_fr, sizeof(test_speech_fr)); + + for (i = 0; i < sizeof(test_speech_efr); i++) + test_speech_efr[i] = i; + test_speech_efr[0] = 0xc0; + test_fr(test_speech_efr, sizeof(test_speech_efr)); + + for (i = 0; i < len_l2; i++) + test_fr(test_l2[i], sizeof(test_l2[0])); + + for (i = 0; i < sizeof(test_speech_hr); i++) + test_speech_hr[i] = i * 17; + test_speech_hr[0] = 0x00; + test_hr(test_speech_hr, sizeof(test_speech_hr)); + + for (i = 0; i < len_l2; i++) + test_hr(test_l2[i], sizeof(test_l2[0])); + + for (i = 0; i < len_mb; i++) { + test_pdtch(test_macblock[i], 23); + test_pdtch(test_macblock[i], 34); + test_pdtch(test_macblock[i], 40); + test_pdtch(test_macblock[i], 54); + } + + printf("Success\n"); + + return 0; +} |