/* Copyright (C) 2013 Matthias P. Braendli http://mpb.li EDI output, Protection, Fragmentation and Transport. (PFT) Are supported: Reed-Solomon and Fragmentation This implements part of PFT as defined ETSI TS 102 821. */ /* This file is part of CRC-DabMux. CRC-DabMux 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 3 of the License, or (at your option) any later version. CRC-DabMux 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 CRC-DabMux. If not, see . */ #include "config.h" #include #include #include #include #include #include "PFT.h" #include "ReedSolomon.h" using namespace std; typedef vector Chunk; RSPacket PFT::Protect(AFPacket af_packet) { RSPacket rs_packet; // number of chunks is ceil(afpacketsize / m_k) if (af_packet.size() % m_k == 0) { m_num_chunks = af_packet.size() / m_k; } else { m_num_chunks = af_packet.size() / m_k + 1; } const size_t zero_pad = m_num_chunks * m_k - af_packet.size(); // add zero padding to last chunk for (size_t i = 0; i < zero_pad; i++) { af_packet.push_back(0); } for (size_t i = 1; i < af_packet.size(); i+= m_k) { Chunk c(m_k + ParityBytes); // copy m_k bytes into new chunk memcpy(&c.front(), &af_packet[i], m_k); // calculate RS for chunk m_encoder.encode(&c.front(), c.size()); // append new chunk to the RS Packet rs_packet.insert(rs_packet.end(), c.begin(), c.end()); } return rs_packet; } vector< vector > PFT::ProtectAndFragment(AFPacket af_packet) { RSPacket rs_packet = Protect(af_packet); const size_t max_payload_size = ( m_num_chunks * ParityBytes ) / (m_m + 1); const size_t fragment_size = m_num_chunks * (m_k + ParityBytes) / max_payload_size; const size_t num_fragments = m_num_chunks * (m_k + ParityBytes) / fragment_size; vector< vector > fragments(num_fragments); for (size_t i = 0; i < num_fragments; i++) { fragments[i].resize(fragment_size); for (size_t j = 0; j < fragment_size; j++) { fragments[i][j] = rs_packet[j*num_fragments + i]; } } return fragments; }