/*
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;
}