test_cpp_raw.cpp

/*
 * Copyright (c) 2015, Luca Fulchir<luker@fenrirproject.org>, All rights reserved.
 *
 * This file is part of "libRaptorQ".
 *
 * libRaptorQ is free software: you can redistribute it and/or modify
 * it under the terms of the GNU Lesser General Public License as
 * published by the Free Software Foundation, either version 3
 * of the License, or (at your option) any later version.
 *
 * libRaptorQ 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
 * and a copy of the GNU Lesser General Public License
 * along with libRaptorQ.  If not, see <http://www.gnu.org/licenses/>.
 */

#include "../src/RaptorQ/RaptorQ_v1.hpp"
#include <algorithm>
#include <fstream>
#include <iostream>
#include <limits>
#include <random>
#include <stdlib.h>
#include <vector>

// Demonstration of how to use the C++ interface
// it's pretty simple, we generate some input,
// then encode, drop some packets (source and repair)
// and finally decode everything.


// mysize is bytes.
template <typename in_enc_align, typename out_enc_align, typename out_dec_align>
bool decode (const uint32_t mysize, std::mt19937_64 &rnd, float drop_prob,
														const uint8_t overhead);
template <typename in_enc_align, typename out_enc_align, typename out_dec_align>
bool decode (const uint32_t mysize, std::mt19937_64 &rnd, float drop_prob,
														const uint8_t overhead)
{
	// define the alignment of the input and output data, for
	// decoder and encoder.
	// note that this is independent from the "mysize" argument,
	// which is always in bytes.
	// used as template arguments
	//typedef uint8_t			in_enc_align;
	//typedef uint16_t			out_enc_align;
	typedef out_enc_align	in_dec_align;
	//typedef uint32_t			out_dec_align;
	// NOTE:  out_enc_align is the same as in_dec_align so that we
	// can simulate data trnsmision just by passing along a vector, but
	// they do not need to be the same.

	// make sure the size is a multiple of its iterator. not much sense
	// otherwise.
	assert ((mysize % sizeof(in_enc_align)) == 0);

	std::vector<in_enc_align> myvec;

	// initialize vector with random data
	// fill remaining data (more than "mysize" bytes) with zeros
	std::uniform_int_distribution<uint8_t> distr(0, 0xFF);
	myvec.reserve (static_cast<size_t> (
				std::ceil(static_cast<float> (mysize) / sizeof(in_enc_align))));
	in_enc_align tmp = 0;
	uint8_t shift = 0;
	for (uint32_t i = 0; i < mysize; ++i) {
		tmp += static_cast<in_enc_align> (distr(rnd)) << shift * 8;
		//tmp += static_cast<in_enc_align> (i) << shift * 8;
		++shift;
		if (shift >= sizeof(in_enc_align)) {
			myvec.push_back (tmp);
			shift = 0;
			tmp = 0;
		}
	}
	if (shift != 0)
		myvec.push_back (tmp);

	// done initializing random data.



	// std::pair<symbol id (esi), symbol>
	std::vector<std::pair<uint32_t, std::vector<out_enc_align>>> encoded;

	// symbol size must be multiple of iterator.
	// symbol size max = uint16_t
	uint16_t max_symsize = static_cast<uint16_t>(
			std::min (static_cast<size_t> (myvec.size()),
				static_cast<size_t> (std::numeric_limits<uint16_t>::max()
													/ sizeof(in_enc_align))));
	std::uniform_int_distribution<uint16_t> sym_distr (1, max_symsize);
	const uint16_t in_aligned_symbol_size = sym_distr (rnd);
	uint16_t symbol_size = in_aligned_symbol_size * sizeof(in_enc_align);
	const uint16_t out_aligned_symbol_size = static_cast<uint16_t> (std::ceil (
					static_cast<float> (symbol_size) / sizeof (out_enc_align)));
	auto enc_it = myvec.begin();

	RaptorQ::Encoder<typename std::vector<in_enc_align>::iterator,
							typename std::vector<out_enc_align>::iterator> enc (
								enc_it, myvec.end(), in_aligned_symbol_size);
	std::cout << "Size: " << mysize << " symbols: " <<
								static_cast<int32_t>(enc.symbols()) <<
								" symbol size: " <<
								static_cast<int32_t>(enc.symbol_size()) << "\n";
	if (!enc.compute_sync()) {
		std::cout << "Enc-RaptorQ failure! really bad!\n";
		return false;
	}

	if (drop_prob > static_cast<float> (90.0))
		drop_prob = 90.0;	// this is still too high probably.


	{
		std::uniform_real_distribution<float> drop (0.0, 100.0);
		// start encoding
		int32_t repair = overhead;
		// Now get the source and repair symbols.
		// make sure that at the end we end with "block.symbols() + overhead"
		// symbols, so that decoding is possible
		for (auto sym_it = enc.begin_source(); sym_it != enc.end_source();
																	++sym_it) {
			float dropped = drop (rnd);
			if (dropped <= drop_prob) {
				// we dropped one source symbol, we need one more repair.
				++repair;
				continue;
			}
			// create a place where to save our source symbol
			std::vector<out_enc_align> source_sym;
			source_sym.reserve (out_aligned_symbol_size);
			source_sym.insert (source_sym.begin(), out_aligned_symbol_size, 0);
			auto it = source_sym.begin();
			// save the symbol
			auto written = (*sym_it) (it, source_sym.end());
			if (written != out_aligned_symbol_size) {
				std::cout << written << "-vs-" << out_aligned_symbol_size <<
									" Could not get the whole source symbol!\n";
				return false;
			}
			// finally add it to the encoded vector
			encoded.emplace_back ((*sym_it).id(), std::move(source_sym));
		}
		// now get (overhead + source_symbol_lost) repair symbols.
		std::cout << "Source Packet lost: " << repair - overhead << "\n";
		auto sym_it = enc.begin_repair();
		for (; repair >= 0 && sym_it != enc.end_repair (enc.max_repair());
																	++sym_it) {
			// repair symbols can be lost, too!
			float dropped = drop (rnd);
			if (dropped <= drop_prob) {
				continue;
			}
			--repair;
			// create a place where to save our source symbol
			std::vector<out_enc_align> repair_sym;
			repair_sym.reserve (out_aligned_symbol_size);
			repair_sym.insert (repair_sym.begin(), out_aligned_symbol_size, 0);
			auto it = repair_sym.begin();
			// save the repair symbol
			auto written = (*sym_it) (it, repair_sym.end());
			if (written != out_aligned_symbol_size) {
				std::cout << written << "-vs-" << out_aligned_symbol_size <<
									" Could not get the whole repair symbol!\n";
				return false;
			}
			// finally add it to the encoded vector
			encoded.emplace_back ((*sym_it).id(), std::move(repair_sym));
		}
		if (sym_it == enc.end_repair (enc.max_repair())) {
			// we dropped waaaay too many symbols! how much are you planning to
			// lose, again???
			std::cout << "Maybe losing " << drop_prob << "% is too much?\n";
			return false;
		}
	}
	// encoding done. now "encoded" is the vector with the trnasmitted data.
	// let's decode it

	using Decoder_type = RaptorQ::Decoder<
								typename std::vector<in_dec_align>::iterator,
								typename std::vector<out_dec_align>::iterator>;
	Decoder_type dec (mysize, symbol_size, Decoder_type::Report::COMPLETE);


	std::vector<out_dec_align> received;
	size_t out_size = static_cast<size_t> (
				std::ceil(static_cast<float>(mysize) / sizeof(out_dec_align)));
	received.reserve (out_size);
	// make sure that there's enough place in "received" to get the
	// whole decoded data.
	for (uint32_t i = 0; i < out_size; ++i)
		received.push_back (static_cast<out_dec_align> (0));

	for (size_t i = 0; i < encoded.size(); ++i) {
		auto it = encoded[i].second.begin();
		auto err = dec.add_symbol (it, encoded[i].second.end(),
															encoded[i].first);
		if (err != RaptorQ::Error::NONE && err != RaptorQ::Error::NOT_NEEDED) {
			std::cout << "error adding?\n";
			abort();
		}
	}

	auto res = dec.decode();
	if (res != Decoder_type::Decoder_Result::DECODED) {
		std::cout << "Couldn't decode.\n";
		return false;
	}

	auto re_it = received.begin();
	// decode all blocks
	// you can actually call ".decode(...)" as many times
	// as you want. It will only start decoding once
	// it has enough data.
	auto decoded = dec.decode_bytes (re_it, received.end(), 0, 0);

	if (decoded.first != mysize) {
		if (decoded.first == 0) {
			std::cout << "Couldn't decode, RaptorQ Algorithm failure. "
															"Can't Retry.\n";
			return true;
		} else {
			std::cout << "Partial Decoding? This should not have happened: " <<
									decoded.first  << " vs " << mysize << "\n";
		}
		return false;
	} else {
		std::cout << "Decoded: " << mysize << "\n";
	}
	// byte-wise check: did we actually decode everything the right way?
	uint8_t *in, *out;
	in = reinterpret_cast<uint8_t *> (myvec.data());
	out = reinterpret_cast<uint8_t *> (received.data());
	for (uint64_t i = 0; i < mysize; ++i) {
		if (in[i] != out[i]) {
			std::cout << "FAILED, but we thought otherwise! " << mysize << " - "
											<< drop_prob << " at " << i << "\n";
			return false;
		}
	}

	return true;
}

uint32_t rnd_size (std::mt19937_64 &rnd, uint8_t size);
uint32_t rnd_size (std::mt19937_64 &rnd, uint8_t size)
{
	std::uniform_int_distribution<uint32_t> distr(1, 10000);
	uint32_t ret;
	do {
		ret = distr (rnd);
	} while ((ret % size) != 0);
	return ret;
}

int main (void)
{
	// get a random number generator
	std::mt19937_64 rnd;
	std::ifstream rand("/dev/urandom");
	uint64_t seed = 0;
	rand.read (reinterpret_cast<char *> (&seed), sizeof(seed));
	rand.close ();
	rnd.seed (seed);


	// encode and decoder
	for (size_t i = 0; i < 1000; ++i) {
		std::cout << "08-08-08\n";
		bool ret = decode<uint8_t, uint8_t, uint8_t> (
								rnd_size (rnd, sizeof(uint8_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "08-08-16\n";
		ret = decode<uint8_t, uint8_t, uint16_t> (
								rnd_size (rnd, sizeof(uint8_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "08-08-32\n";
		ret = decode<uint8_t, uint8_t, uint32_t> (
								rnd_size (rnd, sizeof(uint8_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "08-16-08\n";
		ret = decode<uint8_t, uint16_t, uint8_t> (
								rnd_size (rnd, sizeof(uint8_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "08-16-16\n";
		ret = decode<uint8_t, uint16_t, uint16_t> (
								rnd_size (rnd, sizeof(uint8_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "08-16-32\n";
		ret = decode<uint8_t, uint16_t, uint32_t> (
								rnd_size (rnd, sizeof(uint8_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "08-32-08\n";
		ret = decode<uint8_t, uint32_t, uint8_t> (
								rnd_size (rnd, sizeof(uint8_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "08-32-16\n";
		ret = decode<uint8_t, uint32_t, uint16_t> (
								rnd_size (rnd, sizeof(uint8_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "08-32-32\n";
		ret = decode<uint8_t, uint32_t, uint32_t> (
								rnd_size (rnd, sizeof(uint8_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "16-08-08\n";
		ret = decode<uint16_t, uint8_t, uint8_t> (
								rnd_size (rnd, sizeof(uint16_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "16-08-16\n";
		ret = decode<uint16_t, uint8_t, uint16_t> (
								rnd_size (rnd, sizeof(uint16_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "16-08-32\n";
		ret = decode<uint16_t, uint8_t, uint32_t> (
								rnd_size (rnd, sizeof(uint16_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "16-16-08\n";
		ret = decode<uint16_t, uint16_t, uint8_t> (
								rnd_size (rnd, sizeof(uint16_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "16-16-16\n";
		ret = decode<uint16_t, uint16_t, uint16_t> (
								rnd_size (rnd, sizeof(uint16_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "16-16-32\n";
		ret = decode<uint16_t, uint16_t, uint32_t> (
								rnd_size (rnd, sizeof(uint16_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "16-32-08\n";
		ret = decode<uint16_t, uint32_t, uint8_t> (
								rnd_size (rnd, sizeof(uint16_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "16-32-16\n";
		ret = decode<uint16_t, uint32_t, uint16_t> (
								rnd_size (rnd, sizeof(uint16_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "16-32-32\n";
		ret = decode<uint16_t, uint32_t, uint32_t> (
								rnd_size (rnd, sizeof(uint16_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "32-08-08\n";
		ret = decode<uint32_t, uint8_t, uint8_t> (
								rnd_size (rnd, sizeof(uint32_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "32-08-16\n";
		ret = decode<uint32_t, uint8_t, uint16_t> (
								rnd_size (rnd, sizeof(uint32_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "32-08-32\n";
		ret = decode<uint32_t, uint8_t, uint32_t> (
								rnd_size (rnd, sizeof(uint32_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "32-16-08\n";
		ret = decode<uint32_t, uint16_t, uint8_t> (
								rnd_size (rnd, sizeof(uint32_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "32-16-16\n";
		ret = decode<uint32_t, uint16_t, uint16_t> (
								rnd_size (rnd, sizeof(uint32_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "32-16-32\n";
		ret = decode<uint32_t, uint16_t, uint32_t> (
								rnd_size (rnd, sizeof(uint32_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "32-32-08\n";
		ret = decode<uint32_t, uint32_t, uint8_t> (
								rnd_size (rnd, sizeof(uint32_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "32-32-16\n";
		ret = decode<uint32_t, uint32_t, uint16_t> (
								rnd_size (rnd, sizeof(uint32_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;
		std::cout << "32-32-32\n";
		ret = decode<uint32_t, uint32_t, uint32_t> (
								rnd_size (rnd, sizeof(uint32_t)), rnd, 20.0, 4);
		if (!ret)
			return -1;

	}
	std::cout << "All tests succesfull!\n";
	return 0;
}