Next Article in Journal / Special Issue
Design and Analysis of Systematic Batched Network Codes
Previous Article in Journal
FedDroidMeter: A Privacy Risk Evaluator for FL-Based Android Malware Classification Systems
Previous Article in Special Issue
Network Coding Approaches for Distributed Computation over Lossy Wireless Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

BAR: Blockwise Adaptive Recoding for Batched Network Coding †

1
Department of Information Engineering, The Chinese University of Hong Kong, Shatin, New Territories, Hong Kong, China
2
Department of Electronic and Computer Engineering, The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong, China
3
School of Science and Engineering, The Chinese University of Hong Kong, Shenzhen, Shenzhen 518172, China
4
Department of Computer Science, School of Computing, National University of Singapore, Singapore 119391, Singapore
5
Independent Researcher, Hong Kong, China
*
Authors to whom correspondence should be addressed.
This paper is an extended version of our papers published in Yin, H.H.F.; Yang, S.; Zhou, Q.; Yung, L.M.L. Adaptive Recoding for BATS Codes. In Proceedings of the 2016 IEEE International Symposium on Information Theory (ISIT), Barcelona, Spain, 10–15 July 2016; pp. 2349–2353; and Yin, H.H.F.; Ng, K.H. Impact of Packet Loss Rate Estimation on Blockwise Adaptive Recoding for Batched Network Coding. In Proceedings of the 2021 IEEE International Symposium on Information Theory (ISIT), Melbourne, VIC, Australia, 12–20 July 2021; pp. 1415–1420.
Entropy 2023, 25(7), 1054; https://doi.org/10.3390/e25071054
Submission received: 31 May 2023 / Revised: 23 June 2023 / Accepted: 28 June 2023 / Published: 13 July 2023
(This article belongs to the Special Issue Information Theory and Network Coding II)

Abstract

:
Multi-hop networks have become popular network topologies in various emerging Internet of Things (IoT) applications. Batched network coding (BNC) is a solution to reliable communications in such networks with packet loss. By grouping packets into small batches and restricting recoding to the packets belonging to the same batch; BNC has much smaller computational and storage requirements at intermediate nodes compared with direct application of random linear network coding. In this paper, we discuss a practical recoding scheme called blockwise adaptive recoding (BAR) which learns the latest channel knowledge from short observations so that BAR can adapt to fluctuations in channel conditions. Due to the low computational power of remote IoT devices, we focus on investigating practical concerns such as how to implement efficient BAR algorithms. We also design and investigate feedback schemes for BAR under imperfect feedback systems. Our numerical evaluations show that BAR has significant throughput gain for small batch sizes compared with existing baseline recoding schemes. More importantly, this gain is insensitive to inaccurate channel knowledge. This encouraging result suggests that BAR is suitable to be used in practice as the exact channel model and its parameters could be unknown and subject to changes from time to time.

1. Introduction

Noise, interference and congestion are common causes of packet loss in network communications. Usually, a packet has to travel through multiple hops before it can arrive at the destination node. Traditionally, the intermediate nodes apply the store-and-forward strategy. In order to maintain a reliable communication, retransmission is a common practice. A feedback mechanism is applied so that a network node can acknowledge that a packet is lost. However, due to the delay and bandwidth consumption of the feedback packets, retransmission schemes come with a cost of degraded system performance.
In the era of the Internet of things (IoT), there is a diversity of devices and network topologies. Embedded devices or microcomputers have been heavily deployed due to their mobility, lightweight design, and low power consumption. Multi-hop wireless networks have become a common network topology, highlighting the issues in reliable transmission as wireless links are more vulnerable to packet loss. The packet loss at each link accumulates and the chance of successfully receiving a packet at the destination drops exponentially. Fountain codes, such as in [1,2,3], can recover the lost packets without the need for feedback because of their ratelessness property. However, the throughput still degrades quickly if there is packet loss at each network link unless link-by-link feedback and retransmission are adopted.

1.1. Network Coding-Based Approaches

Random linear network coding (RLNC) [4,5], a simple realization of network coding [6,7,8], can achieve the capacity of multi-hop networks with packet loss even without the need for feedback [9,10]. Unfortunately, a direct application of RLNC induces an enormous overhead for the coefficient vectors, as well as high computational and storage costs in network coding operations at intermediate nodes, where the intermediate nodes are usually routers or embedded devices with low computational power and storage space.
Generation-based RLNC was proposed in [11] to resolve these issues. The input packets of the file are partitioned into multiple subsets called the generations, and RLNC is applied to each generation independently. This approach, however, cannot achieve an optimal theoretical rate. Practical concerns and solutions have been further investigated to improve this type of RLNC, such as decoding delay and complexity [12,13,14,15,16,17,18,19,20], packet size [21,22,23,24,25], and coefficient overhead [26,27,28].
Instead of partitioning into disjoint subsets, overlapped subsets were investigated in [29,30,31,32]. To further reduce the computational costs, the use of RLNC was restricted to small subsets of coded packets generated from the input packets in [33,34,35,36,37,38]. Example codes used for generating the coded packets include LDPC and fountain codes. This combination of coding theory and network coding is another variant of RLNC called batched network coding (BNC). BATS codes [38,39], a class of BNC, have a close-to-optimal achievable rate where the achievable rate is upper bounded by the expectation of the rank distribution of the batch transfer matrices that model the end-to-end network operations (packet erasures, network coding operations, etc.) on the batches [40]. This hints that the network coding operations, also known as recoding, have an impact on the throughput of BNC.

1.2. Challenges of Recoding in Practice

Baseline recoding is the simplest recoding scheme which generates the same number of recoded packets for every batch. Due to its simple and deterministic structure, baseline recoding appears in many BNC designs and analyses, such as [41,42,43,44,45,46]. However, the throughput of baseline recoding is not optimal with finite batch sizes [47]. The idea of adaptive recoding, aiming to outperform baseline recoding by generating different numbers of recoded packets for different batches, was proposed in [47] without truly optimizing the numbers. Two adaptive recoding optimization models for independent packet loss channels were then formulated independently in [48,49]. A unified adaptive recoding framework was proposed in [50], subsuming both optimization models and supporting other channel models under certain conditions.
Although adaptive recoding can be applied distributively with local network information, it is a challenge to obtain accurate local information when deploying adaptive recoding in real-world scenarios. Adaptive recoding requires two pieces of information: information distribution remaining in the received batches and the channel condition of the outgoing link.
The first piece of information may change over time if the channel condition of the incoming link varies. One reason for the variation is that the link quality can be affected by interference from users of other networks around the network node. We proposed a simple way to adapt to this variation in [49], grouping a few batches into a block and observing the distribution of received batches in this block. This approach was later called blockwise adaptive recoding (BAR) in [51,52].
The second piece of information may also vary from time to time. In some scenarios, such as deep-space [53,54,55] and underwater communications [56,57,58], feedback can be expensive or is not available; therefore, a feedbackless network is preferred. Without feedback, we cannot update our knowledge on the channel condition of the outgoing link. Although we may assume an unchanged channel condition and measure information such as the packet loss rate of the channel beforehand, this measurement, however, can be inaccurate due to observational errors or precision limits.

1.3. Contributions

In this paper, we focus on the practical design of applying BAR in real-world applications. Specifically, we answer the following questions in this paper:
  • How does the block size affect the throughput?
  • Is BAR sensitive to an inaccurate channel condition?
  • How can one calculate the components of BAR and solve the optimization efficiently?
  • How can one make use of link-by-link feedback if it is available?
The first question is related to the trade-off between throughput and latency: a larger block induces a longer delay but gives a higher throughput. We show by numerical evaluations that a small block size can already give a significant throughput gain compared with baseline recoding.
For the second question, we demonstrate that BAR performs very well with an independent packet loss model on channels with dependent packet loss. We also show that BAR is insensitive to an inaccurate packet loss rate. This is an encouraging result as this suggests that it is feasible to apply BAR in real-world applications.
The third question is important in practice as BAR is suppose to run at network nodes, usually routers or IoT devices with limited computational power, but also they may need to handle a huge amount of network traffic. Furthermore, by updating the knowledge of the incoming link from a short observation, we need to recalculate the components of BAR and solve the optimization problem again. In light of this, we want to reduce the number of computations to improve the reaction time and reduce the stress of congestion. We answer this question by proposing an on-demand dynamic programming approach to build the components and some implementation techniques to speed up the algorithm for BAR.
Lastly, for the fourth question, we consider both a perfect feedback system (e.g., the feedback passes through a side-channel with no packet loss) and a lossy feedback system (e.g., the feedback uses the reverse direction of the lossy channel for data transmission). We investigate a few ways to estimate the packet loss rate and show that the throughput can be further boosted by using feedback. Furthermore, a rough estimation is sufficient to catch up the variation in the channel condition. In other words, unless there is another application which requires a more accurate estimation on the packet loss rate, we may consider using an estimation with low computational cost, e.g., the maximum likelihood estimator.

1.4. Paper Organization and Nomenclature

The paper is organized as follows. We first formulate BAR in Section 2. Then, we discuss the implementation techniques for solving BAR efficiently and evaluate the throughput of different block sizes in Section 3. In Section 4, we demonstrate that BAR is insensitive to inaccurate channel models and investigate the use of feedback mechanisms. Lastly, we conclude the paper in Section 5.
Some specific terminology and notations appear frequently throughout the paper. We summarize some of the important terminology and frequently used notations in Table 1 and Table 2, respectively.

2. Blockwise Adaptive Recoding

In this section, we briefly introduce BNC and then formulate BAR.

2.1. Network Model

As some intermediate nodes may be hardware-implemented routers or not easily reachable for an upgrade, it is not required to deploy a BNC recoder at every intermediate node. The nodes that do not deploy a recoder are transparent to the BNC network as no network coding operations are performed. In the following text, we only consider intermediate nodes that have deployed BNC recoders.
It is not practical to assume every intermediate node knows the information of the whole network; thus, a distributed scheme that only requires local information is desirable. For example, the statistics of the incoming batches, the channel condition of the outgoing link, etc. In a general network, there may be more than one possible outgoing link to reach the destination. We can assign one recoder or one management unit for each outgoing link at an intermediate node [59,60]. In this way, we need a constraint to limit the number of recoded packets of certain batches sent via the outgoing links. The details are discussed in Section 2.4. In other words, we consider each route from the source to the destination separately as a line network.
Line networks are the fundamental building blocks of a general network. Conversely, a recoding scheme for line networks can be extended to general unicast networks and certain multicast networks [38,48]. A line network is a sequence of network nodes where the network links only exist between two neighbouring nodes. An example of a line network is illustrated in Figure 1. In this paper, we only consider line networks in our numerical evaluations.

2.2. Batched Network Coding

Suppose we want to send a file from a source node to a destination node through a multi-hop network. The file is divided into multiple input packets, where each packet is regarded as a vector over a fixed finite field. A BNC has three main components: the encoder, the recoder and the decoder.
The encoder of a BNC is applied at the source node to generate batches from the input packets, where each batch consists of a small number of coded packets. Recently, a reinforcement learning approach to optimize the generation of batches was proposed in [61]. Nevertheless, batches are commonly generated using the traditional approach as follows. To generate a batch, the encoder samples a predefined degree distribution to obtain a degree, where the degree is the number of input packets that constitute the batch. Depending on the application, there are various ways to formulate the degree distribution [62,63,64,65]. According to the degree, a set of packets is chosen randomly from the input packets. The size of the input packets may be obtained via certain optimizations, such as in [66], to minimize the overhead. Each packet in the batch is formed by taking random linear combinations on the chosen set of packets. The encoder generates M packets per batch, where M is known as the batch size.
Each packet in a batch has a coefficient vector attached to it. Two packets in a batch are defined as linearly independent of each other if and only if their coefficient vectors are linearly independent from each other. Immediately after a batch is generated, the packets within it are assigned as linearly independent from each other. This is accomplished by suitably choosing the initial coefficient vectors [59,67].
A recoder is applied at each intermediate node, performing network coding operations on the received batches to generate recoded packets. This procedure is known as recoding. Some packets of a batch may be lost when they pass through a network link. Each recoded packet of a batch is formed by taking a random linear combination of the received packets in a given batch. The number of recoded packets depends on the recoding scheme. For example, baseline recoding generates the same number of recoded packets for every batch. Optionally, we can also apply a recoder at the source node so that we can have more than M packets per batch at the beginning. After recoding, the recoded packets are sent to the next network node.
At the destination node, a decoder is applied to recover the input packets. Depending on the specific BNC, we can use different decoding algorithms, such as Gaussian elimination, belief propagation and inactivation [68,69].

2.3. Expected Rank Functions

The rank of a batch at a network node is defined by the number of linearly independent packets remaining in the batch, a measure of the amount of information carried by the batch. Adaptive recoding aims to maximize the sum of the expected value of the rank distribution of each batch arriving at the next network node. For simplicity, we called this expected value the expected rank.
For batch b, we denote its rank by r b and the number of recoded packets to be generated by t b . The expectation of r b at the next network node, denoted as E ( r b , t b ) , is known as the expected rank function. We have
E ( r , t ) = i = 0 t Pr ( X t = i ) j = 0 min { i , r } j ζ j i , r ,
where X t is the random variable of the number of packets of a batch received by the next network node when we send t packets for this batch at the current node, and ζ j i , r is the probability that a batch of rank r at the current node with i received packets at the next network node has rank j at the next network node. The exact formulation of ζ j i , r can be found in [38], which is ζ j i , r = ζ j i ζ j r ζ j j q ( i j ) ( r j ) , where q is the field size for the linear algebra operations and ζ j m = k = 0 j 1 ( 1 q m + k ) . It is convenient to use q = 2 8 in practice as each symbol in this field can be represented by 1 byte. For a sufficiently large field size, say q = 2 8 , ζ j i , r is very close to 1 if j = min { i , r } , and is very close to 0 otherwise. That is, we can approximate ζ j i , r by δ j , min { i , r } where δ · , · is the Kronecker delta. This approximation has also been used in the literature, see, e.g., [45,70,71,72,73,74,75,76].
Besides generating all recoded packets by taking random linear combinations, systematic recoding [39,47,59,67], which concerns received packets as recoded packets, can be applied to save computational time. Systematic recoding can achieve a nearly indistinguishable performance compared with methods which generate all recoded packets by taking random linear combinations [39]. Therefore, we can also use (1) to approximate the expected rank functions for systematic recoding accurately.
For the independent packet loss model with packet loss rate p, we have X t Binom ( t , 1 p ) , a binomial distribution. If p = 1 , then a store-and-forward technique can guarantee the maximal expected rank. If p = 0 , then no matter how many packets we transmit, the next network node must receive no packets. Thus, we assume 0 < p < 1 in this paper. It is easy to prove that the results in this paper are also valid for p = 0 or 1 when we define 0 0 : = 1 , which is a convention in combinatorics such that Binom ( t , 0 ) and Binom ( t , 1 ) are well-defined with correct interpretation. In the remaining text, we assume ζ j i , r = δ j , min { i , r } . That is, for the independent packet loss model, we have
E indep ( r , t ) = i = 0 t t i ( 1 p ) i p t i min { i , r } .
A demonstration of the accuracy of the approximation ζ j i , r δ j , min { i , r } can be found in Appendix A.
We also consider the expected rank functions for burst packet loss channels modelled by Gilbert–Elliott (GE) models [77,78], where the GE model was also used in other BNC literature such as [52,55,70]. A GE model is a two-state Markov chain, as illustrated in Figure 2. In each state, there is an independent event to decide whether a packet is lost or not. We define f ( s , i , t ) : = Pr ( S t = s , X t = i ) , where S t is the random variable of the state of the GE model after sending t packets of a batch. By exploiting the structure of the GE model, computation of f can be performed by dynamic programming. Then, we have
E GE ( r , t ) = i = 0 t ( f ( G , i , t ) + f ( B , i , t ) ) min { i , r } .
It is easy to see that it would take more steps to compute (3) than (2). Therefore, a natural question to ask is that for burst packet loss channels, is the throughput gap small between adaptive recoding with (2) and (3)? We demonstrate in Section 4.2 that the gap is small so we use (2) any time a nice throughput is received. Therefore, in our investigation we mainly focus on (2).
In the rest of this paper, we refer to E ( r , t ) as E indep ( r , t ) unless otherwise specified. From [50], we know that when the loss pattern follows a stationary stochastic process, the expected rank function E ( r , t ) is a non-negative, monotonically increasing concave function with respect to t, which is valid for arbitrary field sizes. Further, E ( r , 0 ) = 0 for all r. However, we need to calculate the values of E ( r , t ) or its supergradients to apply adaptive recoding in practice. To cope with this issue, we first investigate the recursive formula for E ( r , t ) .
We define the probability mass function of the binomial distribution Binom ( t , 1 p ) by
B p ( t , i ) = t i ( 1 p ) i p t i if 0 i t , 0 otherwise .
For integers r 0 and t 1 , we define
β p ( t , r ) = 1 if t r 1 , i = 0 r 1 t i ( 1 p ) i p t i = i = 0 r 1 B p ( t , i ) otherwise .
When t 0 , the function β p ( t , r ) is the partial sum of the probability masses of a binomial distribution Binom ( t , 1 p ) . The case where t = 1 is used in the approximation scheme in Section 3 and is discussed in that section.
The regularized incomplete beta function, defined as I x ( a , b ) : = 0 x t a 1 ( 1 t ) b 1 d t 0 1 t a 1 ( 1 t ) b 1 d t ([79], Equation 8.17.2), can be used to express the partial sum of the probability masses of a binomial distribution. When t r > 0 , we can apply ([79], Equation 8.17.5) and obtain
β p ( t , r ) = i = 0 r 1 t i ( 1 p ) i p t i = I p ( t r + 1 , r ) .
There are different implementations of I p ( · , · ) available for different languages. For example, the GNU Scientific Library [80] for C and C++, or the built-in function betainc in MATLAB. However, most available implementations consider non-negative real parameters and calculate different queries independently. This consideration is too general for our application, as we only need to query the integral points efficiently. In other words, this formula may be sufficient for prototyping or simulation, but it is not efficient enough for real-time deployment on devices with limited computational power. Nevertheless, this formula is useful for proving the following properties:
Lemma 1.
Assume 0 < p < 1 . Let Λ be an index set.
(a) 
B p ( t + 1 , i ) = ( 1 p ) B p ( t , i 1 ) + p B p ( t , i ) for i = 0 , 1 , , t ;
(b) 
β p ( t + 1 , r ) β p ( t , r ) where the equality holds if and only if t + 1 < r or t r = 0 ;
(c) 
β p ( t , r ) β p ( t + 1 , r + 1 ) where the equality holds if and only if t < r ;
(d) 
β p ( t , r + 1 ) β p ( t , r ) where the equality holds if and only if t < r ;
(e) 
1 max b Λ β p ( t b , r b ) β p ( t a + s , r a ) for all a Λ and any non-negative integer s;
(f) 
0 min b Λ β p ( t b , r b ) β p ( t a s , r a ) for all a Λ and any non-negative integer s such that t a s 1 .
Proof. 
See Appendix B.    □
With the notation of β p ( t , r ) , we can now write the recursive formula for E ( r , t ) .
Lemma 2.
E ( r , t + 1 ) = E ( r , t ) + ( 1 p ) β p ( t , r ) , where t and r are non-negative integers.
Proof. 
Let Y i be independent and identically distributed Bernoulli random variables, where Pr ( Y i = 1 ) = 1 p for all i. When Y i = 1 , the i-th packet is received by the next hop.
When we transmit one more packet at the current node, Y t + 1 indicates whether this packet is received by the next network node or not. If Y t + 1 = 0 , i.e., the packet is lost, then the expected rank will not change. If Y t + 1 = 1 , then the packet is linearly independent from all the already received packets at the next network node if the number of received packets at the next network node is less than r. That is, the rank of this batch at the next network node increases by 1 if i = 1 t Y i < r . Therefore, the increment of E ( r , t ) is Pr ( Y t + 1 = 1 , i = 1 t Y i < r ) . Note that i = 1 t Y i Binom ( t , 1 p ) . As Y i are all independent and identically distributed, we have Pr Y t + 1 = 1 , i = 1 t Y i < r = ( 1 p ) β p ( t , r ) .    □
The formula shown in Lemma 2 can be interpreted as a newly received packet that is linearly independent of all the already received packets with a probability tends to 1 unless the rank has already reached r. This can also be interpreted as ζ j i , r = δ j , min { i , r } with a probability tends to 1. The above lemma can be rewritten in a more useful form as stated below.
Lemma 3.
Let t and r be non-negative integers.
(a) 
E ( r , t + 1 ) = E ( r , t ) + ( 1 p ) if t < r ;
(b) 
E ( r , t ) = ( 1 p ) j = 0 t 1 β p ( j , r ) = ( 1 p ) min { r , t } + j = r t 1 β p ( j , r ) .
Proof. 
See Appendix C.    □

2.4. Blockwise Adaptive Recoding

The idea of adaptive recoding was presented in [47], and then independently formulated in [48,49]. The former formulation imposes an artificial upper bound on the number of recoded packets and then applies a probabilistic approach to avoid integer programming. The latter investigates the properties of the integer programming problem and proposed efficient algorithms to directly tackle this problem. These two formulations were unified in [50] as a general recoding framework for BNC. This framework requires the distribution of the ranks of the incoming batches, also called the incoming rank distribution. This distribution, however, is not known in advance, and can continually change due to environmental factors. A rank distribution inference approach was proposed in [81], but the long solving time hinders its application in real-time scenarios.
A more direct way to obtain up-to-date statistics is to use the ranks of the few latest batches, a trade-off between a latency of a few batches and the throughput of the whole transmission. This approach was proposed in [49], and later called BAR in [51,52]. In other words, BAR is a recoding scheme which groups batches into blocks and jointly optimizes the number of recoded packets for each batch in the block.
We first describe the adaptive recoding framework and its relation to BAR. We fix an intermediate network node. Let ( h 0 , h 1 , , h M ) be the incoming rank distribution, t r the number of recoded packets to be sent for a batch of rank r, and t avg the average number of recoded packets to be sent per batch. The value of t r is a non-negative real number that is interpreted as follows. Let ϵ = t r t r be the fractional part of t r . There is an ϵ chance to transmit t r + 1 recoded packets, and a 1 ϵ chance to transmit t r packets. That is, the fraction is the probability of transmitting one more packet. Similarly, E ( r , t r ) is defined as the linear interpolation by ( 1 ϵ ) E ( r , t r ) + ϵ E ( r , t r + 1 ) . The framework maximizes the expected rank of the batches at the next node, which is the optimization problem:
max t r 0 , r { 0 , 1 , , M } r = 0 M h r E ( r , t r ) s . t . r = 0 M h r t r = t avg .
For BAR, the incoming rank distribution is obtained from the recently received few batches. Let a block be a set of batches. We assume that the blocks at a network node are mutually disjoint. Suppose the node receives a block L . For each batch b L , let r b and t b be the rank of b and the number of recoded packets to be generated for b, respectively. Let t max L = t avg / | L | be the number of recoded packets to be transmitted for the block L . The batches of the same rank are considered individually with the notations r b and t b , and the total number of packets to be transmitted for a block is finite; therefore, we assume t b for each b L is a non-negative integer, and t max L is a positive integer. By dividing both the objective and the constraint of the framework by | L | , we obtain the simplest formulation of BAR:
max t b { 0 , 1 , 2 , } , b L b L E ( r b , t b ) s . t . b L t b = t max L .
To support scenarios with multiple outgoing links for the same batch, e.g., load balancing, we may impose an upper bound on the number of recoded packets per batch. Let t max b be a non-negative integer that represents the maximum number of recoded packets allowed to be transmitted for the batch b. This value may depend on the rank of b at the node. Subsequently, we can formulate the following optimization problem based on (8):
max t b { 0 , 1 , 2 , } , b L b L E ( r b , t b ) s . t . b L t b = t max L t b t max b , b L .
Note that we must have b L t max b t max L . In the case where this inequality does not hold, we can include more batches in the block to resolve this issue. When t max b is sufficiently large for all b L , (9) degenerates into (8).
The above optimization only depends on the local knowledge at the node. The batch rank r b can be known from the coefficient vectors of the received packets of batch b. As a remark, the value of t max L can affect the stability of the packet buffer. For a general network transmission scenario with multiple transmission sessions, the value of t max L can be determined by optimizing the utility of certain local network transmissions [82,83].
Though we do not discuss such optimizations in this paper, we consider solving BAR with a general value of t max L .
On the other hand, note that the solution to (9) may not be unique. We only need to obtain one solution for recoding purpose. In general, (9) is a non-linear integer programming problem. A linear programming variant of (9) can be formulated by using a technique in [81]. However, such a formulation has a huge amount of constraints and requires the values of E ( r b , t ) for all b L and all possible t to be calculated beforehand. We defer the discussion of this formulation to Appendix H.
A network node will keep receiving packets until it has received enough batches to form a block L . A packet buffer is used to store the received packets. Then, the node solves (9) to obtain the number of recoded packets for each batch in the block, i.e., { t b } b L . The node then generates and transmits t b -recoded packets for every batch b L . At the same time, the network node continually receives new packets. After all the recoded packets for the block L are transmitted, the node drops the block from its packet buffer and then repeats the procedure by considering another block.
We do not specify the transmission order of the packets. Within the same block, the ordering of packets can be shuffled to combat burst loss, e.g., [43,44,52]. Such shuffling can reduce the burst error length experienced by each batch so that the packet loss events are more “independent” from each other. On the other hand, we do not specify the rate control mechanism, as it should be separated as another module in the system. This can be reflected in BAR by choosing suitable expected rank functions, e.g., modifying the parameters in the GE model. BAR is only responsible for deciding the number of recoded packets per batch.
The size of a block depends on its application. For example, if an interleaver is applied to L batches, we can group the L batches as a block. When | L | = 1 , the only solution is t b = t max L , which degenerates into baseline recoding. Therefore, we need to use a block size of at least 2 in order to utilize the throughput enhancement of BAR. Intuitively, it is better to optimize (9) with a larger block size. However, the block size is related to the transmission latency as well as the computational and storage burdens at the network nodes. Note that we cannot conclude the exact rank of each batch in a block until the previous network node finishes sending all the packets of this block. That is, we need to wait for the previous network node to send the packets of all the batches in a block until we can solve the optimization problem. Numerical evaluations in Section 3.5 show that | L | = 2 already has obvious advantage over | L | = 1 , and it may not be necessary to use a block size larger than eight.

3. Implementation Techniques for Blockwise Adaptive Recoding

In this paper, we focus on the implementation and performance of BAR. Due to the non-linear integer programming structure of (9), we need to make use of certain properties of the model in order to solve it efficiently. The authors of [49] proposed greedy algorithms to solve (9), which were then generalized in [50] to solve (7). The greedy algorithms in [50] have an potential issue when certain probability masses in the incoming rank distribution are too small, as they may take too many iterations to find a feasible solution. The number of iterations is in the order of r = 0 M t r , depending on the solution to (7). That is, we cannot establish a bound on the time complexity as the incoming rank distribution can be arbitrary.
For BAR, we do not have this issue because the number of recoded packets in a block, t max L , is fixed.
In this section, we first discuss the greedy algorithm to solve (9) in Section 3.1. Then, we propose an approximation scheme in Section 3.2, and discuss its application to speed up the solver for practical implementations in Section 3.3. The algorithms in Section 3.1 and Section 3.3 are similar to that in [50], but they are modified to optimize BAR. Note that the algorithms in [50] are generalized from [49], so the correctness of the aforementioned modified algorithms is inherited directly from the generalized proofs in [50]. For the approximation scheme in Section 3.2, which did not appear in [50], a more detailed discussion is provided in this section.
The algorithms in this section frequently query and compare the values of ( 1 p ) β p ( t , r ) for different t { 1 , 0 , 1 , , t max L } and r { 0 , 1 , 2 , , M } . We suppose a lookup table is constructed so that the queries can be performed in O ( 1 ) time. The table is reusable if the packet loss rate of the outgoing link is unchanged. We only consider the subset { 1 , 0 , 1 , , t max L } × { 0 , 1 , 2 , , M } of the β p domain because
  • the maximum rank of a batch is M;
  • any t b cannot exceed t max L as b L t b = t max L .
The case t = 1 will be used by our approximation scheme so we keep it in the lookup table. We can build the table on-demand by dynamic programming, discussed in Section 3.4.

3.1. Greedy Algorithm

We first discuss the case t max L b L min { r b , t max b } . This condition means that the value of t max L is too small such that the node has just enough or even not enough time to forward the linearly independent packets received. It is trivial that every { t b } b L satisfying 0 t b min { r b , t max b } and b L t b = t max L is a solution to (9), because every such recoded packet gains 1 p to the expected rank by Lemma 3(a),where this gain is maximal according to the definition of β p ( t , r ) .
For t max L > b L min { r b , t max b } , we can initialize t b by min { r b , t max b } for every b L as every such recoded packet gains the maximal value 1 p to the expected rank. After this, the algorithm chooses the batch that can gain the most expected rank by sending one more recoded packet, and assigns one more recoded packet to it. The correctness is due to the concavity of the expected rank functions.
The above initialization reduces most iterations in the algorithm, as in practice, the difference between the number of recoded packets and the rank of the batch is not huge. Algorithm 1 is the improved greedy algorithm. Unlike the version in [50], the complexity of Algorithm 1 does not depend on the solution.
Algorithm 1: Solver for BAR.
Entropy 25 01054 i001
Theorem 1.
Algorithm 1 can be ran in O ( | L | + max { 0 , t max L b L min { r b , t max b } } log | L | )  time.
Proof. 
There are totally max { 0 , t max L b L min { r b , t max b } } iterations in the while loop. The query of μ = arg max b L T β p ( t b , r b ) can be implemented by using a binary heap. The initialization of the heap, i.e., heapify, takes O ( | L | ) time, which can be performed outside the loop. Each query in the loop takes O ( 1 ) time. The update from β p ( t μ , r μ ) into β p ( t μ + 1 , r μ ) , if t μ < t max μ 1 , takes O ( log | L | ) time. For t μ = t max μ 1 , we may remove the entry from the heap, taking the same time complexity as the update above. As b L min { r b , t max b } t max L by assumption, the algorithm will not query an empty heap. Therefore, the overall time complexity is O ( | L | + max { 0 , t max L b L min { r b , t max b } } log | L | ) .    □
In the algorithm, we assume that a lookup table for β p ( t , r ) is pre-computed. The table can be reused unless there is an update on the outgoing channel condition. Nevertheless, we will discuss an efficient way to construct the lookup table in Section 3.4, and the insignificance of the measurement or prediction errors of the loss probability of the outgoing channel in Section 4.
As the query arg max b L T β p ( t b , r b ) is run repeatedly and an update is performed after every query, we can use a binary heap as described in the proof in real implementation. Note that by Lemma 1(b), β p ( t μ , r μ ) β p ( t μ + 1 , r μ ) , so the update is a decrease key operation in a max-heap. In other words, a Fibonacci heap [84] cannot benefit from this operation here.

3.2. Equal Opportunity Approximation Scheme

Algorithm 1 increases t b step by step. From a geometric perspective, the algorithm finds a path from the interior of a compact convex polytope that models the feasible solutions to the facet H : b L t b = t max L . If we have a method to move a non-optimal feasible point on H towards an optimal point, together with a fast and accurate approximation to (8) or (9), then we can combine them to solve (9) faster than using Algorithm 1 directly. This idea is illustrated in Figure 3. A generalized tuning scheme can be found in [50] based on the algorithm in [49]. However, there is no approximation scheme proposed in [50].
We first give an approximation scheme in this subsection. The approximation is based on an observation of the solution for (8) that does not impose an upper boundary on t b : A batch of higher rank should have more recoded packets transmitted than a batch of lower rank. Unless most t max b are too small, the approximation for (8) is also a good approximation for (9).
Theorem 2.
Let L be a block where | L | 2 . If { t b } b L solves (8) and t b r b for all b L , then t m < t n for all m , n L such that r m < r n .
Proof. 
See Appendix D.    □
As we cannot generate any linearly independent packets for a batch of rank 0, we have E ( 0 , · ) = 0 . Therefore, we can exclude batches of rank 0 from L before we start the approximation. We define L = { b L : r b > 0 } L . When t max L > b L r b , we have t b r b for all b L . An easy way to obtain an approximation is to assign { t b } b L following the guidelines given in Theorem 2 by:
  • t b = 0 for all b L \ L ;
  • t b = r b + for all b L .
where = ( t max L b L r b ) / | L | . In the case where is not an integer, we can round it up for batches with higher ranks and round it down for those with lower ranks.
The above rules allocate the unassigned packets to batches equally after r b packets have been assigned to each batch b. Thus, we call this approach the equal opportunity approximation scheme. The steps of this scheme are summarized in Algorithm 2.
Note that we do not need to know the packet loss rate p to apply this approximation. That is, if we do not know the value of p, we can still apply this approximation to outperform baseline recoding.
Theorem 3.
Algorithm 2 approximates (8) in O ( | L | ) time. If t max L b L r b , then the algorithm solves (8).
Algorithm 2: Equal opportunity approximation scheme.
Entropy 25 01054 i002
Proof. 
It is easy to see that Algorithm 2 outputs { t b } b L which satisfies b L t b = t max L . That is, the output is a feasible solution of (8). Note that | L | | L | , so the assignments and the branches before the last for loop take O ( | L | ) time in total. The variable L after the first foreach loop equals | L | . Adding one to the number of recoded packets for r = mod L batches with the highest ranks can be performed in O ( | L | ) time. Therefore, the overall running time is O ( | L | ) .
If L = , i.e., the whole block is lost, then any feasible { t b } b L is a solution, and the optimal objective value is 0. If t max L b L r b , then the algorithm terminates with an output satisfying t b r b for all b L , which is an optimal solution.    □
For the step that adds 1 to the number of recoded packets for r = mod L batches with the highest ranks in the algorithm, the worst linear time case can be achieved by using introselect [85] (which is quickselect [86] with a random pivot, but changes to use the median of medians [87] pivot strategy when the complexity grows). We use the selection algorithm to find the r-th largest element, making use of its intermediate steps. During an iteration, one of the following three cases will occur. If the algorithm decides to search a part larger than the pivot, then the discarded part does not contain the largest r elements. If a part smaller than the pivot is selected, then the discarded part is part of the largest r elements. If the pivot is exactly the r-th largest element, then the part larger than the pivot together with the pivot are part of the largest r elements.
In practice, the batch size M is small. We can search these r batches with the highest ranks in O ( | L | + M ) time using a counting technique as an efficient alternative. The technique is to use part of the counting algorithm [88]. We first compute a histogram of the number of times each rank occurs, taking O ( M ) time for initialization and O ( | L | ) time to scan the block. Then, we can scan and count the frequencies of the histogram from the highest rank, and eliminate the part where the count exceeds mod | L | . This takes O ( M ) time. Lastly, we scan the ranks of the batches again in O ( | L | ) time. If included in the modified histogram, we add 1 to the corresponding t b and minus 1 to the corresponding frequency in the histogram.
Algorithm 2 is a ( 1 p ) -approximation algorithm, although the relative performance guarantee factor 1 p is not tight in general. However, this suggests that the smaller the packet loss rate p, the more accurate the output the algorithm gives. We defer this discussion to Appendix E.

3.3. Speed-Up via Approximation

In this subsection, we discuss the implementation that corrects an approximate solution to an optimal solution for (9). Algorithm 3 is a greedy algorithm that uses any feasible solution of (8) as a starting point. The foreach loop removes the exceeding recoded packets, assigning the released slot to another batch following the iterations in Algorithm 1. This can be regarded as mimicking a replacement of β p ( · , · ) with the smallest possible value for the batch b that violates the constraint t b t max b . After this, the intermediate solution is a feasible solution to (9). Then, the last loop finds an increase to the objective by reassigning some slots among the batches.
Algorithm 3: Solver for BAR via approximation.
Entropy 25 01054 i003
Note that the algorithm may query β p ( t a 1 , r a ) for a L . If t a = 0 , then it has access to the value β p ( 1 , r a ) . Recall that we defined β p ( 1 , · ) = 1 , the upper bound of β p ( · , · ) by (10). Therefore, these values act as barriers to prevent outputting a negative number of recoded packets.
Theorem 4.
Let { t b } b L be an approximate solution of (8) computed in O ( T approx ) time. Algorithm 3 can be run in O ( T approx + | L | + b L | t b * t b | log | L | ) time.
Proof. 
The assignments before the foreach loop takes O ( T approx + | L | ) time. There are a total of b L | t b * t b | / 2 iterations in the loops. The queries for the minimum and maximum values can be implemented using a min-heap and a max-heap, respectively. Similar to Algorithm 1, we can use binary heaps, taking O ( | L | ) initialization time, O ( 1 ) query time, and O ( log | L | ) update time. Each iteration contains at most two heap queries and four heap updates. The update of the set T can be performed implicitly by setting β p ( t a , r a ) to 0 during the heap updates for a T . The overall time complexity is then O ( T approx + | L | + b L | t b * t b | log | L | ) . □
In the last while loop, we need to query the minimum of β p ( t a 1 , r a ) and the maximum of β p ( t b , r b ) . It is clear that we need to decrease the key β p ( t b , r b ) to β p ( t b + 1 , r b ) in the max-heap, and increase the key β p ( t a 1 , r a ) to β p ( t a 2 , r a ) in the min-heap. However, we can omit the updates for batches a and b in the max-heap and min-heap, respectively, i.e., reduce from four heap updates to two heap updates. We defer this discussion to Appendix G.

3.4. Construction of the Lookup Table

In the above algorithms, we assume that we have a lookup table for the function β p ( · , · ) so that we can query its values quickly. In this subsection, we propose an on-demand approach to construct a lookup table by dynamic programming.
Due to the fact that i = 0 t B p ( t , i ) = 1 , we have
0 β p ( t , r ) 1 .
Furthermore, it is easy to see that
β p ( t , r ) = 0 if and only if r = 0 and t 0 ;
β p ( t , r ) = 1 if and only if t r 1 .
A tabular form of β p is illustrated in Figure 4 after introducing the boundaries 0 and 1 s.
Being a dynamic programming approach, we need the following recursive relations:
B p ( t + 1 , r ) = ( 1 p ) B p ( t , r 1 ) + p B p ( t , r ) for 0 r t ;
B p ( r , r ) = ( 1 p ) B p ( r 1 , r 1 ) for r > 0 ;
β p ( t , r ) = β p ( t , r 1 ) + B p ( t , r 1 ) for 1 < r t + 1 ,
where (13) is stated in Lemma 1(a); and (14) and (15) are by the definitions of B p ( t , r ) and β p ( t , r ) , respectively. The boundary conditions are B p ( 0 , 0 ) = 1 , B p ( i , 1 ) = 0 , and β p ( i , 1 ) = B p ( i , 0 ) for i = 0 , 1 , . The table can be built in-place in two stages. The first stage fills in B p ( y , x 1 ) at the ( y , x ) position of the table. The second stage finishes the table by using (15). Figure 5 illustrates the two stages where the arrows represent the recursive relations (13)–(15). As β p ( 0 , 1 ) = β p ( 1 , 2 ) = = 1 , the corresponding entries can be substituted in directly.
We can compute the values in the table on-demand. Suppose we have { t b } b L in an iteration of Algorithm 1 so that we need the values of β p ( t b , r b ) for b L . Let b be an element in arg max b L r b . The table has r b + 1 columns. By the criteria of selecting b in Algorithm 1 and by Lemmas 1(c) and (d),we have max b L t b = t b . From Figure 5, we know we have to calculate all rows of β ( t , r ) for t t b . Furthermore, the recursive relations on a row only depend on the previous row; thus, we need to prepare the values of B p in the next row so that we have the values to compute β p in the next row. As an example, Figure 6 illustrates the values we have prepared when t b = 2 .
Each entry in the table is modified at most twice during the two stages. Each assignment takes O ( 1 ) time. Therefore, the time and space complexities for building the table are both O ( M R ) , where R is the number of rows we want to construct. When restricted by the block size, we know that R t max L . The worst case is that we only receive one rank-M batch for the whole block, which is unlikely to occur. In this case, we have the worst case complexity O ( M t max L ) .
Note that we can use fixed-point numbers instead of floating point numbers for a more efficient table construction. Furthermore, the numerical values in the table are not important as long as the orders for any pair of values in the table are the same.

3.5. Throughput Evaluations

We now evaluate the performance of BAR in a feedbackless multi-hop network. Note that baseline recoding is a special case of BAR with block size 1. Our main goal here is to show the throughput gain of BAR among different block sizes. In the evaluation, all (recoded) packets of a batch are sent before sending those of another batch.
Let ( h 0 , h 1 , , h M ) be the incoming rank distribution of batches arriving at a network node. The normalized throughput at a network node is defined as the average rank of the received batches divided by the batch size, i.e., i = 0 M i h i / M . In our evaluations in this subsection, we set t max L = M | L | for every block L . That is, the source node transmits M packets per batch. We assume that every link in the line network has independent packet loss with the same packet loss rate p. In this topology, we set a sufficiently large t max b for every batch, say, t max b = t max L .
We first evaluate the normalized throughput with different batch sizes and packet loss rates. Figure 7 compares adaptive recoding (AR) and baseline recoding (BR) when we know the rank distribution of the batches arriving at each network node before the node applies BAR. In other words, Figure 7 shows the best possible throughput of AR. We compare the effect of block sizes later. We observe that
  • AR has a higher throughput than BR under the same setting;
  • the difference in throughput between AR and BR is larger when the batch size is smaller, the packet loss probability is larger, or the length of the line network is longer.
In terms of throughput, the percentage gains of AR over BR using M = 4 and p = 0.2 are 23.3 and 33.7% at the 20-th and 40-th hops, respectively. They become 43.8 and 70.3 % , respectively, when p = 0.3 .
Although the above figure shows that the throughput of BNC with AR maintains a good performance when the length of the line network is long, many applications use a much shorter line network. We zoom into the figure for the first 10 hops in Figure 8 for practical purposes.
Now, we consider the effect of different block sizes. Figure 9 shows the normalized throughput of different | L | and p with M = 8 . The first 10 hops in Figure 9 are zoomed in in Figure 10. We observe that
  • a larger | L | results a better throughput;
  • using | L | = 2 already gives a much larger throughput than using | L | = 1 ;
  • using | L | > 8 gives little extra gain in terms of throughput.
Next, we show the performance of the equal opportunity approximation scheme. Figure 11 compares the normalized throughput achieved by Algorithm 2 (AS) and the true optimal throughput (AR). We compare the best possible throughput of AR here, i.e., the same setting as in Figure 7. The first 10 hops in Figure 11 are zoomed in in Figure 12. We observe that
  • the approximation is close to the optimal solution;
  • the gap in the normalized throughput is smaller when the batch size is larger, the packet loss probability is smaller, or the length of the line network is shorter.

4. Impact of Inaccurate Channel Models

In this section, we first demonstrate that the throughput of BAR is insensitive to inaccurate channel models and packet loss rates. Then, we investigate the feedback design and show that although feedback can enhance the throughput, the benefit is insignificant. In other words, BAR works very well without the need of feedback.

4.1. Sensitivity of β p ( t , r )

We can see that our algorithms only depend on the order of the values of β p ( · , · ) ; therefore, it is possible that the optimal { t b } b L for an incorrect p is the same for a correct p. As shown in Figure 4, the boundaries 0 and 1 s are unaffected by p ( 0 , 1 ) . That is, we only need to investigate the stability of β p ( t , r ) for t r > 0 . We calculate values of β p ( t , r ) corrected to four digital places in Figure 13 for M = 4 , and p = 0.1 , 0.45 and their 1% relative changes. We can see that the order of the values are mostly the same when we slightly change p.
We can also check with the condition number [89] to verify the stability. Roughly speaking, the relative change in the function output is approximately equal to the condition number times the relative change in the function input. A small condition number of β p ( t , r ) means that the effect of the inaccurate p is small. As shown in Figure 13, the values of β p ( t , r ) drop quickly when t increases. In the view of the throughput, which is proportional to the sum of these values, we can tolerate a larger relative change, i.e., a larger condition number, when β p ( t , r ) is small. We calculate condition numbers of β p ( t , r ) in Figure 14 by the formula stated in Theorem 5.
Theorem 5.
Let p ( 0 , 1 ) and t r > 0 . The condition number of β p ( t , r ) with respect to p is p t r + 1 ( 1 p ) r 1 t ! I p ( t r + 1 , r ) ( t r ) ! ( r 1 ) ! , or equivalently, j = 0 r 1 ( 1 ) j r 1 j p t r + j + 1 j = 0 r 1 ( 1 ) j r 1 j p t r + j + 1 / ( t r + j + 1 ) .
Proof. 
See Appendix F. □

4.2. Impact of Inaccurate Channel Models

To demonstrate the impact of an inaccurate channel model, we consider three different channels to present our observations.
  • ch1: independent packet loss with constant loss rate p = 0.45 .
  • ch2: burst packet loss modelled by the GE model illustrated in Figure 2 with the parameters used in [70], namely p G B = p B G = p G = 0.1 , p B = 0.8 .
  • ch3: independent packet loss with varying loss rates p = 0.45 + 0.3 sin ( 2 π c / 1280 ) , where c is the number of transmitted batches.
All the three channels have the same average packet loss rate of 0.45 . The formula of ch3 is for demonstration purpose only.
We now demonstrate the impact of inaccurate p on the throughput. We consider a line network where all the links use the same channel (ch1, ch2, or ch3). In this topology, we set a sufficiently large t max b for every batch, say, t max b = t max L . Similar to the previous evaluation, all (recoded) packets of a batch are sent before sending those of another batch. Furthermore, we set t max L = M | L | for every block L .
In Figure 15 we plot the normalized throughput of the first 80 received blocks at the fourth hop where | L | = M = 4 or 8. We use BAR with (2) for each network although ch2 is a bursty channel. The black curves with BAR are the throughput of BAR where the loss rate is known. For ch1 and ch2, this loss rate p is a constant of 0.45 . The red and blue curves are the throughput of BAR when we guess p = 0.65 and 0.25 , respectively, which is ± 0.2 from the average loss rate of 0.45 . As there is no feedback, we do not change our guess on p for these curves. We can see that the throughput is actually very close to the corresponding black curves. This suggests that in the view of the throughput, BAR is not sensitive to p. Even with a wild guess on p, BAR still outperforms BR, as illustrated by the green curves. Regarding ch2, we also plot the orange curve with GE BAR, which is the throughput achieved by BAR with (3). We can see that the gap between the throughput achieved by BAR with (2) and (3) is very small. As a summary of our demonstration:
  • We can use BAR with (2) for bursty channels and the loss in throughput is insignificant.
  • BAR with an inaccurate constant p can achieve a throughput close to the one when we have the exact real-time loss rate.
  • We can see a significant throughput gain from BR by using BAR even with inaccurate channel models.

4.3. Feedback Design

Although an inaccurate p can give an acceptable throughput, we can further enhance the throughput by adapting the varying p values. To achieve this goal, we need to use feedback.
We adopt a simple feedback strategy which lets the next node return the number of received packets of the batches for the current node to estimate p. Although the next node does not know the number of lost packets per batch, it knows the number of received packets per batch. Therefore, we do not need to introduce more overhead to the transmitted packets by the current node.
When we estimate p, we have to know the number of packets lost during a certain time frame. If the time frame is too small, the estimation is too sensitive so the estimated p changes rapidly and unpredictably. If the time frame is too long, we captured too much out-dated information about the channel so the estimated p changes too slowly and may not be able to adapt to the real loss rate. Recall that the block size is not large as we want to keep the delay small. We use a block as an atomic unit of the time frame. The next node gives feedback on the number of received packets per block. The current node uses the feedback of the blocks in the time frame to estimate p. We perform an estimation of p per received feedback. In this way, the estimated p is the same for each block so that we can apply BAR with (2).
If the feedback is sent via a reliable side channel, then we can assume that the current node can always receive the feedback. However, if the feedback is sent via an unreliable channel, say, the reverse direction of the same channel the data packets were sent from, then we need to consider feedback loss. Let Λ be a set of blocks in a time frame with received feedback. We handle the case of feedback loss by considering the total number of packets transmitted for the blocks in Λ as the total number of packets transmitted during the time frame. In this way, we can also start the estimation before a node sends enough blocks to fill up a time frame. Suppose no feedback is received for every block in a time frame, then we reuse the previously estimated p for BAR.
At the beginning of the transmission, we have no feedback yet so we have no information to estimate p. To outperform BR without the knowledge of p, we can use the approximation of BAR given by Algorithm 2. Once we have received at least one feedback, we can then start estimating p.

4.4. Estimators

Let x and n be the total number of packets received by the current node and the total number of packets transmitted by the previous node, respectively, in a time frame for observation. That is, the number of packets lost in the time frame is n x . We introduce three types of estimators for our numerical evaluations.
(1) Maximum likelihood estimator (MLE): The MLE, denoted by p ^ MLE , estimates p by maximizing the likelihood function. p ^ MLE = ( n x ) / n is a well-known result which can be obtained via derivative tests. This form collides with the sample average, so by the law of large numbers, p ^ MLE p when n if p does not change over time.
(2) Minimax estimator: The minimax estimator achieves the smallest maximum risk among all estimators. With the popular mean squared error (MSE) as the risk function, it is a Bayes estimator with respect to the least favourable prior distribution. As studied in [90,91], such prior distribution is a beta distribution Beta ( n / 2 , n / 2 ) . The minimax estimator of p, denoted by p ^ MM , is the posterior mean, which is n 1 + n n x n + 1 1 + n 1 2 , or equivalently, n x + 0.5 n n + n .
(3) Weighted Bayesian update: Suppose the prior distribution is Beta ( a , b ) , where the hyperparameters can be interpreted as a pseudo-observation having a successes and b failures. Given a sample of s successes and f failures from a binomial distribution, the posterior distribution is Beta ( a + s , b + f ) . To fade out the old samples captured by the hyperparameters, we introduce a scaling factor 0 γ 1 and let the posterior distribution be Beta ( γ a + s , γ b + f ) . This factor can also prevent the hyperparameters from growing indefinitely. The estimation of p, denoted by p ^ Bayes , is the posterior mean with s = n x and f = x , which is γ a + n x γ ( a + b ) + n . To prevent a bias when there are insufficient samples, we select a non-informative prior as the initial hyperparameters. Specifically, we use the Jeffreys prior, which is Beta ( 1 / 2 , 1 / 2 ) .
We first show the estimation of p by different schemes in Figure 16. We use BAR with (2) and | L | = M = 4 . The size of the time frame is W blocks. For p ^ MLE and p ^ MM , the observations in the whole time frame have the same weight. For p ^ Bayes , the effect of each observation deceases exponentially faster. We consider an observation is out of the time frame when it is scaled into 10% of the original value. That is, we define the scaling factor by γ = 0.1 W . In each subplot, the black curve is the real-time p. The red and blue curves are for the estimation without and with feedback loss, respectively. In each case, the two curves are the 25 and 75% percentiles from 1000 runs, respectively.
We can see that a larger W has a slower response to the change in p in ch3. Among the estimators, p ^ Bayes has the fastest response speed as its observations in the time frame are not fairly weighted. Furthermore, although ch1 and ch2 have the same average loss rate, the estimation has a larger variance when the channel is bursty.

4.5. Throughput Evaluations

As discussed in Section 4.2, the guessed p values have an insignificant impact on the throughput. We now show the throughput achieved by the estimation schemes in Figure 17. The parameters for the networks and BAR are the same as in Section 4.2. We do not wildly guess p here so it is no surprise that we can achieve nearly the same throughput as when we know the real p for ch1 and ch2. If we look closely, we can see from Figure 15 that for ch3, there is a small gap between the throughput of BAR when we know the real-time p and the one of BAR when using a constant p. Although the estimation may not be accurate at all times, we can now adapt to the change in p to finally achieve a throughput nearly the same as when we know the real-time p. On the other hand, whether the feedback is lost or not, the plots shown in Figure 17 are basically the same.

5. Conclusions

We proposed BAR in this paper which can adapt to variations in the incoming channel condition. In a practical perspective, we discussed how to calculate the components of BAR and how to solve BAR efficiently. We also investigated the impact of an inaccurate channel model on the throughput achieved by BAR. Our evaluations showed that
  • BAR is insensitive to the channel model: guessing the loss rate still outperforms BR.
  • For bursty channels, the throughput achieved by BAR with an independent loss model is nearly identical to one with the real channel model. That is, we can use the independent loss model for BAR in practice and apply the techniques in this paper to reduce the computational costs of BAR.
  • Feedback can slightly enhance the throughput for channels with a dynamic loss rate. This suggests that BAR works very well without the need of feedback. On the other hand, feedback loss barely affects the throughput of BAR. Therefore, we can send the feedback through a lossy channel without the need of retransmission. Unless we need to use an accurate estimated loss rate in other applications, we can use MLE with a small time frame for BAR to reduce the computational time.
These encouraging results suggest that BAR is suitable to be deployed in real-world applications.
One drawback of our proposed scheme is that we need to change the default behaviour of some intermediate network nodes, which can be a practical problem in existing networks. In fact, this is a common issue for all network coding schemes. Some routers have hard-wired circuits to efficiently handle heavy traffic, so it is unfeasible to deploy other schemes on them without replacing the hardware. For these heavy-loaded nodes, one may consider producing a hardware to speed up the network coding operations, e.g., [92,93], inducing extra costs on the deployment. On the other hand, the protocol for BNC is not standardized yet, meaning two parties may adopt BAR with incompatible protocols, thus restricting the application of BNC in public networks. However, it is not easy to build a consensus on the protocol, because there are still many research directions to improve the performance of BNC so the protocol design is subject to change in the near future.

6. Patents

The algorithms in Section 3 are variants of those that can be found in the U.S. patent 10,425,192 granted on 24 September 2019 [94]. The linear programming-based algorithm for BAR in Appendix H can be found in the U.S. patent 11,452,003 granted on 20 September 2022 [95].

Author Contributions

Conceptualization, H.H.F.Y. and S.Y.; methodology, H.H.F.Y.; software, H.H.F.Y. and L.M.L.Y.; validation, H.H.F.Y., Q.Z. and K.H.N.; formal analysis, H.H.F.Y.; investigation, H.H.F.Y., S.Y. and Q.Z.; writing—original draft preparation, H.H.F.Y.; writing—review and editing, H.H.F.Y., S.Y. and Q.Z.; visualization, H.H.F.Y.; supervision, H.H.F.Y. and S.Y.; project administration, H.H.F.Y.; funding acquisition, S.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by NSFC under Grants 12141108 and 62171399.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Part of the work of Hoover H. F. Yin, Shanghao Yang and Qiaoqiao Zhou was conducted when they were with the Institute of Network Coding, The Chinese University of Hong Kong, Shatin, New Territories, Hong Kong. The work of Lily M. L. Yung was performed when she was with the Department of Computer Science and Engineering, The Chinese University of Hong Kong, Shatin, New Territories, Hong Kong. The work of Ka Hei Ng was performed when he was with the Department of Physics, The Chinese University of Hong Kong, Shatin, New Territories, Hong Kong. Some results in this paper were included in the thesis of Hoover H. F. Yin [96].

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BARBlockwise adaptive recoding
IoTInternet of Things
BNCBatched network coding
RLNCRandom linear network coding
LDPCLow-density parity-check code
BATS codeBatched sparse code
GE modelGilbert–Elliott model
MLEMaximum likelihood estimator
MSEMean squared error

Appendix A. Accuracy of the Approximation ζ j i , r δ j , min { i , r }

We demonstrate the accuracy of the approximation ζ j i , r δ j , min { i , r } by showing the percentage error of the expected rank function corrected to three decimal places when q = 2 8 , p = 0.2 and X t Binom ( t , 1 p ) in Table A1. That is, the table shows the values
100 i = 0 t t i ( 1 p ) i p t i ( min { r , i } j = 0 min { i , r } j ζ j i , r ) i = 0 t t i ( 1 p ) i p t i j = 0 min { i , r } j ζ j i , r
for different r and t.
From the table, we can see that only three pairs of ( r , t ) have percentage errors larger than 0.1%, where they occur when r , t 2 . For all the other cases, the percentage errors are less than 0.1%. Therefore, such an approximation is accurate enough for practical applications.
Table A1. Percentage error when approximating expected rank functions.
Table A1. Percentage error when approximating expected rank functions.
t r = 1 r = 2 r = 3 r = 4 r = 5 r = 6 r = 7 r = 8 r = 9 r = 10 r = 11 r = 12 r = 13 r = 14 r = 15 r = 16
10.392160.001530.000010.000000.000000.000000.000000.000000.000000.000000.000000.000000.000000.000000.000000.00000
20.131400.157410.000610.000000.000000.000000.000000.000000.000000.000000.000000.000000.000000.000000.000000.00000
30.038410.080320.083970.000330.000000.000000.000000.000000.000000.000000.000000.000000.000000.000000.000000.00000
40.010250.030910.057910.050420.000200.000000.000000.000000.000000.000000.000000.000000.000000.000000.000000.00000
50.002580.010240.027610.043980.032290.000130.000000.000000.000000.000000.000000.000000.000000.000000.000000.00000
60.000620.003080.010910.025020.034160.021550.000080.000000.000000.000000.000000.000000.000000.000000.000000.00000
70.000150.000870.003820.011470.022580.026850.014790.000060.000000.000000.000000.000000.000000.000000.000000.00000
80.000030.000230.001220.004570.011780.020230.021230.010360.000040.000000.000000.000000.000000.000000.000000.00000
90.000010.000060.000370.001650.005260.011820.017980.016860.007380.000030.000000.000000.000000.000000.000000.00000
100.000000.000020.000110.000550.002100.005850.011630.015850.013420.005320.000020.000000.000000.000000.000000.00000
110.000000.000000.000030.000170.000770.002570.006320.011250.013880.010700.003870.000020.000000.000000.000000.00000
120.000000.000000.000010.000050.000270.001030.003020.006650.010720.012080.008540.002840.000010.000000.000000.00000
130.000000.000000.000000.000020.000090.000380.001310.003440.006850.010090.010460.006820.002100.000010.000000.00000
140.000000.000000.000000.000000.000030.000130.000520.001600.003810.006930.009400.009010.005450.001560.000010.00000
150.000000.000000.000000.000000.000010.000040.000200.000690.001900.004120.006890.008660.007730.004360.001170.00000
160.000000.000000.000000.000000.000000.000010.000070.000280.000870.002190.004360.006770.007920.006600.003490.00088
170.000000.000000.000000.000000.000000.000000.000020.000100.000370.001060.002460.004540.006560.007180.005620.00279
180.000000.000000.000000.000000.000000.000000.000010.000040.000150.000480.001260.002710.004660.006290.006470.00476
190.000000.000000.000000.000000.000000.000000.000000.000010.000060.000200.000600.001470.002930.004710.005980.00579
200.000000.000000.000000.000000.000000.000000.000000.000000.000020.000080.000270.000730.001670.003110.004700.00563

Appendix B. Proof of Lemma 1

We have the following properties when a and b are positive integers and 0 x 1 ([79], Equations 8.17.20 and 8.17.21):
I x ( a , b ) I x ( a + 1 , b ) = a + b 1 a x a ( 1 x ) b ;
I x ( a , b + 1 ) I x ( a , b ) = a + b 1 b x a ( 1 x ) b .
Proof of Lemma 1(a) 
It is trivial for i = 0 . For i > 0 , recall the recursive formula of binomial coefficients ([79], Equation 1.2.7):
t + 1 i = t i 1 + t i , i = 1 , 2 , , t .
Applying the formula, we have
    B p ( t + 1 , i ) = t + 1 i ( 1 p ) i p t + 1 i = ( 1 p ) t i 1 ( 1 p ) i 1 p t ( i 1 ) + p t i ( 1 p ) i p t i = ( 1 p ) B p ( t , i 1 ) + p B p ( t , i ) .
Proof of Lemma 1(b). 
Case I: t < r . By (10) and (12), β p ( t + 1 , r ) 1 = β p ( t , r ) , and the equality holds if and only if t + 1 r 1 < r .
Case II: t r > 0 . By (6) and (A1),
    β p ( t , r ) β p ( t + 1 , r ) = I p ( t r + 1 , r ) I p ( t r + 2 , r ) = t t r + 1 p t r + 1 ( 1 p ) r > 0 .
Case III: t r = 0 . By (11), the equality always hold. □
Proof of Lemma 1(c). 
Case I: t < r . By (12), the equality always hold.
Case II: t r > 0 . By (6) and (A2),
    β p ( t + 1 , r + 1 ) β p ( t , r ) = I p ( t r + 1 , r + 1 ) I p ( t r + 1 , r ) = t r p t r + 1 ( 1 p ) r > 0 .
Case III: t r = 0 . By (10) and (11), β p ( t , r ) = 0 < β p ( t + 1 , r + 1 ) . □
Proof of Lemma 1(d). 
Case I: t = 1 . By definition, β p ( t , r + 1 ) = β p ( t , r ) = 1 .
Case II: t 0 . Recall that β p ( t , r ) is the partial sum of the probability mass of the binomial distribution Binom ( t , 1 p ) . By summing one more term, i.e., β p ( t , r + 1 ) , the partial sum must be larger than or equal to β p ( t , r ) . Note that B p ( t , i ) 0 when 0 i t , so the equality holds if and only if β p ( t , r ) = 1 and if t < r by (12). □
Proof of Lemma 1(e) and (f). 
Inductively by Lemma 1(b), we have
β p ( t a + u , r a ) β p ( t a , r a ) β p ( t a v , r a )
for all a Λ where u , v are non-negative integers such that t a v 1 . By (10),
0 min b Λ β p ( t b , r b ) β p ( t a , r a ) max b Λ β p ( t b , r b ) 1
for all a Λ . Combining (A3) and (A4), the proof is complete. □

Appendix C. Proof of Lemma 3

By Lemma 3(a), we have E ( r , t + 1 ) = E ( r , t ) + ( 1 p ) β p ( t , r ) . If t < r , we have
β p ( t , r ) = i = 0 r 1 B p ( t , i ) = 1 ,
which proves Lemma 3(a).
For Lemma 3(b), note that we have the initial condition
E ( r , 0 ) = B p ( 0 , 0 ) min { r , 0 } = 0 = ( 1 p ) j = 0 ( 0 ) 1 β p ( j , r ) .
We can evaluate Lemma 2 recursively and obtain the first equality in Lemma 3(b).
By Lemma 3(a), we can show that when t < r , we have
E ( r , t ) = t ( 1 p ) .
This implies that when t r , we have
E ( r , t ) = ( 1 p ) r + j = r t 1 β p ( j , r ) .
When t < r , the summation term j = r t 1 β p ( j , r ) in (A6) equals 0. So, we can combine (A5) and (A6) and give
E ( r , t ) = ( 1 p ) min { r , t } + j = r t 1 β p ( j , r ) .

Appendix D. Proof of Theorem 2

Suppose t m > t n for some r m < r n , i.e.,
t m > t n r n > r m .
We define
t b = t m if b = n , t n if b = m , t b otherwise
for all b L . We consider the difference of
    b L E ( r b , t b ) b L E ( r b , t b ) = [ E ( r m , t n ) + E ( r n , t m ) ] [ E ( r m , t m ) + E ( r n , t n ) ] = [ E ( r n , t m ) E ( r n , t n ) ] + [ E ( r m , t n ) E ( r m , t m ) ] = ( 1 p ) j = 0 t m 1 β p ( j , r n ) j = 0 t n 1 β p ( j , r n ) + j = 0 t n 1 β p ( j , r m ) j = 0 t m 1 β p ( j , r m ) ( A8 ) = ( 1 p ) j = t n t m 1 β p ( j , r n ) j = t n t m 1 β p ( j , r m ) = ( 1 p ) j = t n t m 1 β p ( j , r n ) β p ( j , r m ) > 0 , ( A9 )
where
  • (A8) follows Lemma 3(b);
  • (A9) follows Lemma 1(d) together with (A7).
The above result contradicts that { t b } b L solves (8), which gives us that t m t n for all r m < r n .
Next, we suppose t m = t n for some r m < r n , i.e.,
t m = t n r n > r m .
We define
t b = t n + 1 if b = n , t m 1 if b = m , t b otherwise
for all b L . Moreover, we compare the difference of
    b L E ( r i , t i ) b L E ( r i , t i ) = [ E ( r n , t n + 1 ) + E ( r m , t m 1 ) ] [ E ( r n , t n ) + E ( r m , t m ) ] = [ E ( r n , t n + 1 ) E ( r n , t n ) ] [ E ( r m , t m ) E ( r m , t m 1 ) ] = ( 1 p ) [ β p ( t n , r n ) β p ( t m 1 , r m ) ] ( A11 ) ( 1 p ) [ β p ( t m , r m + 1 ) β p ( t m 1 , r m ) ] ( A12 ) > 0 , ( A13 )
where
  • (A11) follows Lemma 3(a)
  • (A12) follows (A10) and Lemma 1(d)
  • (A13) follows Lemma 1(c) together with (A10).
This contradicts { t b } b L and solves (8). Therefore, we have t m t n for all r m < r n .
Combining the two cases, the proof is complete.

Appendix E. Performance Guarantee and Bounded Error of Algorithm 2

We start the discussion with the following theorem.
Theorem A1.
Let SOL and OPT be the solution given by Algorithm 2 and the optimal solution of (8), respectively, then
SOL ( 1 p ) OPT , OPT SOL ( 1 p ) b L j = r b + r b + | L | 1 β p ( j , r b ) ,
where = ( t max L b L r b ) / | L | and = .
Proof. 
We first show that the algorithm has a relative performance guarantee factor of 1 p . As stated in Theorem 3, when t max L b L r b , the algorithm guarantees an optimal solution. Therefore, we only consider t max L > b L r b . Let { t b } b L be the approximation given by the algorithm.
Note that any linear combinations of r-independent vectors cannot obtain more than r-independent vectors. Therefore, the expected rank of a batch at the next hop must be no larger than the rank of the batch at the current hop, and, be non-negative. That is,
0 E ( r b , t ) r b , t 0 , b L .
This gives a bound of the optimal solution by
0 OPT b L r b .
We consider the exact formula of the approximation:
    S O L = ( 1 p ) b L r b + ( 1 p ) b L j = r b t b 1 β p ( j , r b ) ( A16 ) ( 1 p ) b L r b ( A17 ) ( 1 p ) OPT , ( A18 )
where
  • (A16) is stated in Lemma 3(b)
  • (A17) holds as β p ( j , r b ) 0 for all j , r b , which is by (10);
  • (A18) follows the inequality (A15).
Lastly, we show the bounded error. Let { t b * } be a solution to (8). We write t b * = r b + b where b 0 for all b L . Note that the constraint of (8), i.e., b L t b * = t max L , suggests that
b t max L b L r b = | L | .
On the other hand, it is easy to see that the approximation must either give t b = r b + or t b = r b + + 1 . That is, we have t b r b + . By Lemma 3(b),we have
SOL ( 1 p ) b L r b + j = r b r b + 1 β p ( j , r b ) .
We consider the difference between OPT and SOL:
    OPT SOL ( 1 p ) b L j = r b r b + b 1 β p ( j , r b ) j = r b r b + 1 β p ( j , r b ) ( A21 ) = ( 1 p ) b L j = r b + , b > r b + b 1 β p ( j , r b ) j = r b + b , b < r b + 1 β p ( j , r b ) ( 1 p ) b L j = r b + , b > r b + b 1 β p ( j , r b ) ( A22 ) ( 1 p ) b L j = r b + r b + | L | 1 β p ( j , r b ) , ( A23 )
where
  • (A21) is the difference between the exact form of OPT by Lemma 3(b) after substituting the lower bound of SOL shown in (A20);
  • the condition b > in the summation of (A22) can be removed, as we have r b + b 1 < r b + if b ;
  • (A23) follows (A19) and the fact shown in (10) that the extra β p ( j , r b ) terms are non-negative.
The proof is done. □
If the relative performance guarantee factor of 1 p is tight, we need both equalities in (A17) and (A18) to hold. First, by (10), we know that β p ( j , r b ) is always non-negative. The equality in (A17) holds if and only if j = r b t b 1 β p ( j , r b ) = 0 for all b L . The sum equals 0 only when
  • r b = 0 and t b 0 according to (11); or
  • t b 1 < r b which forms an empty sum.
The equality in (A18) holds if and only if OPT = b L E ( r b , t b * ) = b L r b . Note that (A14) shows that E ( r b , t b * ) is upper bounded by r b . This implies that we need E ( r b , t b * ) = r b for all b L . When t b * r b , we can apply Lemma 3(a) to obtain E ( r b , t b * ) = ( 1 p ) t b * , which equals r b if and only if r b = 0 , as we assumed 0 < p < 1 in this paper. By Lemma 2, E ( r b , t ) is a monotonic increasing function in terms of t for all r b 0 . Therefore, when r b 0 we need t b * > r b , which implies that t max L > b L r b . Then, the approximation will also give t b > r b for some b L in this case, and the equality in (A17) does not hold.
That is, we have SOL = ( 1 p ) OPT only when r b = 0 for all b L . In this case, we have SOL = OPT = 0 . In practice, the probability of having r b = 0 for all b L is very small. Therefore, we can consider that the bound is not tight in most cases but it guarantees that the approximation is good when the packet loss probability is small.

Appendix F. Proof of Theorem 5

Let B ( a , b ; y ) : = 0 y x a 1 ( 1 x ) b 1 d x be the incomplete beta function. We have the beta function B ( a , b ) : = B ( a , b ; 1 ) .
From (6), we have β p ( t , r ) = I p ( t r + 1 , r ) = B ( t r + 1 , r ; p ) B ( t r + 1 , r ; 1 ) . By direct calculation, the condition number is
| p d β p ( t , r ) d p β p ( t , r ) = p d d p 0 p x t r ( 1 x ) r 1 d x B ( t r + 1 , r ; 1 ) I p ( t r + 1 , r )                                   = p t r + 1 ( 1 p ) r 1 B ( t r + 1 , r ; 1 ) I p ( t r + 1 , r )                                   = p t r + 1 ( 1 p ) r 1 0 p x t r ( 1 x ) r 1 d x                                   = p t r + 1 j = 0 r 1 ( 1 ) j r 1 j p j 0 p x t r j = 0 r 1 ( 1 ) j r 1 j x j d x                                   = j = 0 r 1 ( 1 ) j r 1 j p t r + j + 1 j = 0 r 1 ( 1 ) j r 1 j 0 p x t r + j d x                                   = j = 0 r 1 ( 1 ) j r 1 j p t r + j + 1 j = 0 r 1 ( 1 ) j r 1 j p t r + j + 1 / ( t r + j + 1 ) ,
where the absolute value disappears as both numerator and denominator are non-negative. The first form of the condition number can be obtained by substituting B ( t r + 1 , r ; 1 ) = ( t r ) ! ( r 1 ) ! t ! into (A24).

Appendix G. Corrupted Heaps

In this appendix, we explain why we can omit two of the heap updates in Algorithm 3.
Before we start, we need some mathematical descriptions of the optimal solutions of BAR. Then, we will introduce our lazy evaluation technique on a modified heap that we called the corrupted heap.
To simplify the notations, we redefine β p ( t b , r b ) as
β p ( t , r ) = 0 if t b t max b , 1 if t b min { r b , t max b } 1 , i = 0 r b 1 t b i ( 1 p ) i p t b i otherwise .
in this appendix. In other words, β p ( · , · ) is now a function of b. When t b t max b , β p ( t b , r b ) is the smallest possible value in the image of β p ( · , · ) . Therefore, the algorithms in this paper will not assign more recoded packets to the batch b.

Appendix G.1. Optimality Properties of BAR

First, we introduce the following theorem that states a condition for non-optimality (or optimality after taking contraposition).
Theorem A2.
Let { t b } b L be a feasible solution of (9). Then, { t b } b L is not an optimal solution of (9) if and only if there exists two distinct batches κ and ρ with t ρ 1 such that ( 1 p ) β p ( t κ , r κ ) > ( 1 p ) β p ( t ρ 1 , r ρ ) .
Proof. 
We first prove the sufficient condition. If { t b } b L does not solve (9), then it means that there exists another configuration { t b } b L which can give a higher objective value. Since b L t b = b L t b = t max L , there exists distinct κ , ρ L such that t κ > t κ and t ρ < t ρ . Note that t ρ 0 so we must have t ρ 1 . We define
Θ = { κ : t κ > t κ } and Φ = { ρ : t ρ < t ρ } ,
where
θ Θ ( t θ t θ ) = ϕ Φ ( t ϕ t ϕ ) > 0 .
Using the fact that { t b } b L gives a larger objective value and by Lemma 3(b),we have
θ Θ t = t θ t θ 1 ( 1 p ) β p ( t , r θ ) > ϕ Φ t = t ϕ t ϕ 1 ( 1 p ) β p ( t , r ϕ ) .
Now, we fix κ and ρ such that
κ arg max θ Θ β p ( t θ , r θ ) and ρ arg min ϕ Φ β p ( t ϕ 1 , r ϕ ) .
We have
    θ Θ ( t θ t θ ) ( 1 p ) β p ( t κ , r κ ) θ Θ ( t θ t θ ) ( 1 p ) β p ( t θ , r θ ) θ Θ t = t θ t θ 1 ( 1 p ) β p ( t , r θ ) ( A27 ) > ϕ Φ t = t ϕ t ϕ 1 ( 1 p ) β p ( t , r ϕ ) ( A28 ) ϕ Φ ( t ϕ t ϕ ) ( 1 p ) β p ( t ϕ 1 , r ϕ ) ( A29 ) ϕ Φ ( t ϕ t ϕ ) ( 1 p ) β p ( t ρ 1 , r ρ ) ,
where
  • (A27) and (A29) follows Lemma 1(b);
  • (A28) is the inequality shown in (A26).
Applying (A25), we have
( 1 p ) β p ( t κ , r κ ) > ( 1 p ) β p ( t ρ 1 , r ρ ) ,
which proves the sufficient condition.
Now we consider the necessary condition, where we have
( 1 p ) β p ( t κ , r κ ) > ( 1 p ) β p ( t ρ 1 , r ρ )
for some distinct κ and ρ . Let
t b = t κ + 1 if b = κ , t ρ 1 if b = ρ , t b otherwise
for all b L . Then, we consider the following:
    b L E ( r b , t b ) = b L { κ , ρ } E ( r b , t b ) + E ( r κ , t κ + 1 ) + E ( r ρ , t ρ 1 ) = b L { κ , ρ } E ( r b , t b ) + E ( r ρ , t ρ 1 ) + [ E ( r κ , t κ ) + ( 1 p ) β p ( t κ , r κ ) ] ( A31 ) > b L { κ , ρ } E ( r b , t b ) + E ( r κ , t κ ) + [ E ( r ρ , t ρ 1 ) + ( 1 p ) β p ( t ρ 1 , r ρ ) ] ( A32 ) = b L { ρ } E ( r b , t b ) + E ( r ρ , t ρ ) ( A33 ) = b L E ( r b , t b ) ,
where
  • (A31) and (A33) follow Lemma 3(a)
  • (A32) follows (A30).
meaning that { t b } b L is not an optimal solution of (9). □
Next, we define the sub-problems (A34) of (9) for k { 0 , 1 , , t max L } , which present the optimal substructure of (9).
max t b { 0 , 1 , 2 , } , b L b L E ( r b , t b ) s . t . b L t b = k t b t max b , b L .
Note that we assume b L t max b t max L k , or otherwise we should include more batches in the block.
We define a multiset Ω r that collects the value of ( 1 p ) β p ( t , r ) for all integers t 0 , i.e.,
Ω r = { ( 1 p ) β p ( t , r ) : t { 0 , 1 , 2 , } } .
By Lemma 3(b),we have E ( r , t ) = ( 1 p ) j = 0 t 1 β p ( j , r ) . As E ( r , t ) is concave with respect to t, β p ( t , r ) is a monotonic decreasing function on t (also stated in Lemma 1(b)).
This implies the following lemma.
Lemma A1.
E ( r , t ) equals the sum of the largest t elements in Ω r .
Proof. 
By Lemma 3(b), E ( r , t ) = ( 1 p ) j = 0 t 1 β p ( j , r ) . By Lemma 1(b), β p ( t , r ) is a monotonic decreasing function on t. By (10) and (12), we have β p ( 0 , r ) = 1 β p ( t , r ) for all positive integers t. Therefore, E ( r , t ) is the sum of the largest t elements in Ω r . □
We fix a block L and define a multiset
Ω : = b L Ω r b = { ( 1 p ) β p ( t b , r b ) : t b { 0 , 1 , 2 , , t max b 1 } , b L } .
If two batches a , b L have the same rank, i.e., r a = r b , then for all t, we have ( 1 p ) β p ( t , r a ) = ( 1 p ) β p ( t , r b ) . As Ω is a multiset, the duplicated values are not eliminated. Now, we have the following lemma to connect (A34) and Ω .
Lemma A2.
The optimal value of (A34) is the sum of the largest k elements in Ω.
Proof. 
Let { t b } b L solves (A34). We suppose the optimal value is not the sum of the largest k elements in Ω . However, Lemma A1 states that E ( r b , t b ) equals the sum of the largest t b elements in Ω r b for all b L . This means that there exists two distinct batches κ , ρ L with t κ t max b 1 and t ρ 1 such that ( 1 p ) β p ( t κ , r κ ) > ( 1 p ) β p ( t ρ 1 , r ρ ) .
By setting t max L = k , we can apply Theorem A2 which gives that { t b } b L is not an optimal solution of (A34). The proof is completed by contradiction. □
We define a multiset Ω t max L which is a collection of the largest t max L elements in Ω . By Lemma A2, ω Ω t max L ω is the optimal value of (9). For any non-optimal solution ( t b ( k ) ) b L , we define a multiset
k : = { ( 1 p ) β p ( t , r b ) : t { 0 , 1 , , t b ( k ) 1 } , b L } ,
where the value k { 0 , 1 , , t max L 1 } is the number of elements in Ω t max L which are also contained in k .

Appendix G.2. Lazy Evaluations

We consider an iteration in the last while loop in Algorithm 3. Suppose we choose to increase t b by 1 and decrease t a by 1.
Lemma A3.
If batch a is selected by the max-heap or batch b is selected by the min-heap in any future iteration, then the optimal solution is reached.
Proof. 
Suppose batch a with key A is selected by the max-heap in a future iteration. Note that A was once the smallest element in k for some k. Therefore, at the current state where k > k , every element in k must be no smaller than A . Equivalently, we have ( 1 p ) β p ( t κ , r κ ) ( 1 p ) β p ( t ρ 1 , r ρ ) for all κ , ρ L . By Theorem A2, the optimal solution is reached. The min-heap counterpart can be proved in a similar fashion. □
Suppose we omit the update for the batch ρ in the heap. We call the key of the batch ρ a corrupted key, or the key of the batch ρ is corrupted. A key which is not corrupted is called an uncorrupted key. A heap with corrupted keys is called a corrupted heap. In other words, the key of a batch is corrupted in a corrupted max-heap if and only if the same batch was once the minimum of the corresponding original min-heap, and vice versa. As a remark, we do not have a guaranteed maximum portion of corrupted keys as an input. Furthermore, we do not adopt the carpooling technique. This suggests that the heap here is not a soft heap [97].
Lemma A4.
If the root of a corrupted heap is a corrupted key, then the optimal solution is reached.
Proof. 
We only consider a corrupted max-heap in the proof. We can use similar arguments to show that a corrupted min-heap also works.
In a future iteration, suppose batch a is selected by the corrupted max-heap. We consider the real maximum in the original max-heap. There are three cases.
Case I: batch a is also the root of the original max-heap. As the key of a is corrupted, it means that the batch was once selected by the corresponding min-heap. By Lemma A3, the optimal solution is reached.
Case II: the root of the original max-heap is batch a where the key of a is also corrupted. Similar to Case I, batch a was once selected by the corresponding min-heap, and we can apply Lemma A3 to finish this case.
Case III: the root of the original max-heap is batch a where the key of a is not corrupted. In this case, the uncorrupted key of a is also in the corrupted max-heap. Note that the corrupted key of a is no larger than the actual key of a in the original max-heap. This means that the key of a, a and the corrupted key of a have the same value. It is equivalent to let the original max-heap select batch a, as every element in k must be no smaller than the key of a , where k represents the state of the current iteration. Then, the problem is reduced to Case I.
Combining the three cases, the proof is completed. □
Theorem A3.
The updates for batch a in the max-heap and batch b in the min-heap can be omitted.
Proof. 
When we omit the updates, the heap itself becomes a corrupted heap. We have to make sure that when a batch with corrupted key is selected, the termination condition of the algorithm is also met.
We can express the key of batch π in a corrupted max-heap and min-heap by β p ( t π + s π , r π ) and β p ( t π 1 u π , r π ) , respectively, where s π , u π are non-negative integers. When s π or u π is 0, the key is uncorrupted in the corresponding corrupted heap. By Lemma 1(b), we have
β p ( t π + s π , r π ) β p ( t π , r π ) , β p ( t π 1 , r π ) β p ( t π 1 u π , r π ) .
That is, the root of the corrupted max-heap is no larger than the root of the original max-heap. Similar for the min-heap. Mathematically, we have
max π L β p ( t π + s π , r π ) max π L β p ( t π , r π ) ,
min π L β p ( t π 1 , r π ) min π L β p ( t π 1 u π , r π ) .
Suppose a corrupted key is selected. By Lemma A4, we know that the optimal solution is reached. Therefore, we can apply the contrapositive of Theorem A2 and know that
( 1 p ) β p ( t κ , r κ ) ( 1 p ) β p ( t ρ 1 , r ρ )
for all κ , ρ L . We can omit the condition t ρ 1 because by (10) and (12), we have β p ( 1 , · ) = 1 β p ( · , · ) . The inequality (A37) is equivalent to
max π L β p ( t π , r π ) min π L β p ( t π 1 , r π ) .
We can mix this inequality with (A35) and (A36) to show that when a corrupted key is selected, we have
max π L β p ( t π + s π , r π ) min π L β p ( t π 1 u π , r π ) ,
which is the termination condition shown in Algorithm 3 after we replaced the heaps into corrupted heaps.
We just showed that once a corrupted key selected, the termination condition is reached. In the other words, before a corrupted key is selected, every previous selection must be an uncorrupted key. That is, the details inside the iterations are unaffected. If an uncorrupted key is selected where it also satisfies the termination condition, then no corrupted key is touched, and the corrupted heap still acts as a normal heap at this point.
The correctness of the algorithm when using a corrupted heap is proven. Moreover, we do not need to mark which key is corrupted. This is, we can omitted the mentioned heap updates for a normal heap. □
We do not need to mark down which key is corrupted while the algorithm still works, so we can simply omit the mentioned updates as lazy evaluations. As there are two heaps in algorithm, we can reduce from four to two heap updates.

Appendix H. Linear Programming Formulation of BAR

In [81], a distributionally robust optimization [98] for AR is formulated as a linear programming problem. It is based on an observation that when the expected rank function E ( r , t ) is concave with respect to t, we can reformulate it by
E ( r , t ) = min i { 0 , 1 , , ı ¯ } ( Δ r , i t + ξ r , i )
if we fix an artificial upper bound t ı ¯ , where Δ r , t : = E ( r , i + 1 ) E ( r , i ) and ξ r , i : = E ( r , i ) i Δ r , i . In (9), we implicitly have t t max L , so we can make use of this expression to write (9) as
max t b , e b 0 , b L b L e b s . t . b L t b = t max L t b t max b , b L e b E ( r b , i ) + ( E ( r b , i + 1 ) E ( r b , i ) ) ( t b i ) , b L , i { 0 , 1 , , t max L } ,
where t b is allowed to be a non-integer. A non-integer t b means that we first generate t b recoded packets, then we generate one more recoded packet with probability t b t b . Note that there are | L | t max L constraints for e b .
To turn such a non-deterministic solution into a deterministic one, we perform the following steps:
  • Collect the batches with non-integer recoded packets into a set S.
  • Calculate R = b S ( t b t b ) . Note that R is an integer for BAR.
  • For every b S , remove the fractional part of t b .
  • Randomly select R batches from S and add one recoded packet to each of these batches.
We have an integer R because b L t b = t max L . Furthermore, we have R < | S | . Referring to the idea of Algorithm 1, we have the same value of Δ r b , t b for all b S . After removing the fractional part of t b for all b S , it becomes the sub-problem (A34) (defined in Appendix G.1) with k = t max L R . The last step follows Algorithm 1 such that the output is a solution to (9) where t b for all b L are all integers.

References

  1. Luby, M. LT Codes. In Proceedings of the 43rd Annual IEEE Symposium on Foundations of Computer Science, Vancouver, BC, Canada, 19 November 2002; pp. 271–282. [Google Scholar]
  2. Shokrollahi, A. Raptor Codes. IEEE Trans. Inf. Theory 2006, 52, 2551–2567. [Google Scholar] [CrossRef]
  3. Maymounkov, P. Online Codes; Technical Report; New York University: New York, NY, USA, 2002. [Google Scholar]
  4. Ho, T.; Koetter, R.; Médard, M.; Karger, D.R.; Effros, M. The Benefits of Coding over Routing in a Randomized Setting. In Proceedings of the 2003 IEEE International Symposium on Information Theory (ISIT), Yokohama, Japan, 29 June–4 July 2003; p. 442. [Google Scholar]
  5. Ho, T.; Médard, M.; Koetter, R.; Karger, D.R.; Effros, M.; Shi, J.; Leong, B. A Random Linear Network Coding Approach to Multicast. IEEE Trans. Inf. Theory 2006, 52, 4413–4430. [Google Scholar] [CrossRef] [Green Version]
  6. Ahlswede, R.; Cai, N.; Li, S.Y.R.; Yeung, R.W. Network Information Flow. IEEE Trans. Inf. Theory 2000, 46, 1204–1216. [Google Scholar]
  7. Koetter, R.; Médard, M. An Algebraic Approach to Network Coding. IEEE/ACM Trans. Netw. 2003, 11, 782–795. [Google Scholar] [CrossRef] [Green Version]
  8. Li, S.Y.R.; Yeung, R.W.; Cai, N. Linear Network Coding. IEEE Trans. Inf. Theory 2003, 49, 371–381. [Google Scholar] [CrossRef] [Green Version]
  9. Wu, Y. A Trellis Connectivity Analysis of Random Linear Network Coding with Buffering. In Proceedings of the 2006 IEEE International Symposium on Information Theory (ISIT), Seattle, WA, USA, 9–14 July 2006; pp. 768–772. [Google Scholar]
  10. Lun, D.S.; Médard, M.; Koetter, R.; Effros, M. On coding for reliable communication over packet networks. Phys. Commun. 2008, 1, 3–20. [Google Scholar] [CrossRef] [Green Version]
  11. Chou, P.A.; Wu, Y.; Jain, K. Practical Network Coding. In Proceedings of the Annual Allerton Conference on Communication Control and Computing, Monticello, IL, USA, 1–3 October 2003; Volume 41, pp. 40–49. [Google Scholar]
  12. Pandi, S.; Gabriel, F.; Cabrera, J.A.; Wunderlich, S.; Reisslein, M.; Fitzek, F.H.P. PACE: Redundancy Engineering in RLNC for Low-Latency Communication. IEEE Access 2017, 5, 20477–20493. [Google Scholar] [CrossRef]
  13. Wunderlich, S.; Gabriel, F.; Pandi, S.; Fitzek, F.H.P.; Reisslein, M. Caterpillar RLNC (CRLNC): A Practical Finite Sliding Window RLNC Approach. IEEE Access 2017, 5, 20183–20197. [Google Scholar] [CrossRef]
  14. Lucani, D.E.; Pedersen, M.V.; Ruano, D.; Sørensen, C.W.; Fitzek, F.H.P.; Heide, J.; Geil, O.; Nguyen, V.; Reisslein, M. Fulcrum: Flexible Network Coding for Heterogeneous Devices. IEEE Access 2018, 6, 77890–77910. [Google Scholar] [CrossRef]
  15. Nguyen, V.; Tasdemir, E.; Nguyen, G.T.; Lucani, D.E.; Fitzek, F.H.P.; Reisslein, M. DSEP Fulcrum: Dynamic Sparsity and Expansion Packets for Fulcrum Network Coding. IEEE Access 2020, 8, 78293–78314. [Google Scholar] [CrossRef]
  16. Tasdemir, E.; Tömösközi, M.; Cabrera, J.A.; Gabriel, F.; You, D.; Fitzek, F.H.P.; Reisslein, M. SpaRec: Sparse Systematic RLNC Recoding in Multi-Hop Networks. IEEE Access 2021, 9, 168567–168586. [Google Scholar] [CrossRef]
  17. Tasdemir, E.; Nguyen, V.; Nguyen, G.T.; Fitzek, F.H.P.; Reisslein, M. FSW: Fulcrum sliding window coding for low-latency communication. IEEE Access 2022, 10, 54276–54290. [Google Scholar] [CrossRef]
  18. Fu, A.; Sadeghi, P.; Médard, M. Dynamic rate adaptation for improved throughput and delay in wireless network coded broadcast. IEEE/ACM Trans. Netw. 2014, 22, 1715–1728. [Google Scholar] [CrossRef]
  19. Chatzigeorgiou, I.; Tassi, A. Decoding delay performance of random linear network coding for broadcast. IEEE Trans. Veh. Technol. 2017, 66, 7050–7060. [Google Scholar] [CrossRef] [Green Version]
  20. Yazdani, N.; Lucani, D.E. Revolving codes: Overhead and computational complexity analysis. IEEE Commun. Lett. 2021, 25, 374–378. [Google Scholar] [CrossRef]
  21. Torres Compta, P.; Fitzek, F.H.P.; Lucani, D.E. Network Coding is the 5G Key Enabling Technology: Effects and Strategies to Manage Heterogeneous Packet Lengths. Trans. Emerg. Telecommun. Technol. 2015, 6, 46–55. [Google Scholar] [CrossRef]
  22. Torres Compta, P.; Fitzek, F.H.P.; Lucani, D.E. On the Effects of Heterogeneous Packet Lengths on Network Coding. In Proceedings of the European Wireless 2014, Barcelona, Spain, 14–16 May 2014; pp. 385–390. [Google Scholar]
  23. Taghouti, M.; Lucani, D.E.; Cabrera, J.A.; Reisslein, M.; Pedersen, M.V.; Fitzek, F.H.P. Reduction of Padding Overhead for RLNC Media Distribution with Variable Size Packets. IEEE Trans. Broadcast. 2019, 65, 558–576. [Google Scholar] [CrossRef]
  24. Taghouti, M.; Tömösközi, M.; Howeler, M.; Lucani, D.E.; Fitzek, F.H.; Bouallegue, A.; Ekler, P. Implementation of Network Coding with Recoding for Unequal-Sized and Header Compressed Traffic. In Proceedings of the 2019 IEEE Wireless Communications and Networking Conference (WCNC), Marrakech, Morocco, 15–19 April 2019. [Google Scholar]
  25. Schütz, B.; Aschenbruck, N. Packet-Preserving Network Coding Schemes for Padding Overhead Reduction. In Proceedings of the 2019 IEEE 44th Conference on Local Computer Networks (LCN), Osnabrueck, Germany, 14–17 October 2019; pp. 447–454. [Google Scholar]
  26. de Alwis, C.; Kodikara Arachchi, H.; Fernando, A.; Kondoz, A. Towards Minimising the Coefficient Vector Overhead in Random Linear Network Coding. In Proceedings of the 2013 IEEE International Conference on Acoustics, Speech and Signal Processing, Vancouver, BC, Canada, 26–31 May 2013; pp. 5127–5131. [Google Scholar]
  27. Silva, D. Minimum-Overhead Network Coding in the Short Packet Regime. In Proceedings of the 2012 International Symposium on Network Coding (NetCod), Boston, MA, USA, 29–30 June 2012; pp. 173–178. [Google Scholar]
  28. Gligoroski, D.; Kralevska, K.; Øverby, H. Minimal Header Overhead for Random Linear Network Coding. In Proceedings of the 2015 IEEE International Conference on Communication Workshop (ICCW), London, UK, 8–12 June 2015; pp. 680–685. [Google Scholar]
  29. Silva, D.; Zeng, W.; Kschischang, F.R. Sparse Network Coding with Overlapping Classes. In Proceedings of the 2009 Workshop on Network Coding, Theory, and Applications, Lausanne, Switzerland, 15–16 June 2009; pp. 74–79. [Google Scholar]
  30. Heidarzadeh, A.; Banihashemi, A.H. Overlapped Chunked Network Coding. In Proceedings of the 2010 IEEE Information Theory Workshop on Information Theory (ITW), Dublin, Ireland, 30 August–3 September 2010; pp. 1–5. [Google Scholar]
  31. Li, Y.; Soljanin, E.; Spasojevic, P. Effects of the Generation Size and Overlap on Throughput and Complexity in Randomized Linear Network Coding. IEEE Trans. Inf. Theory 2011, 57, 1111–1123. [Google Scholar] [CrossRef] [Green Version]
  32. Tang, B.; Yang, S.; Yin, Y.; Ye, B.; Lu, S. Expander graph based overlapped chunked codes. In Proceedings of the 2012 IEEE International Symposium on Information Theory (ISIT), Cambridge, MA, USA, 1–6 July 2012; pp. 2451–2455. [Google Scholar]
  33. Mahdaviani, K.; Ardakani, M.; Bagheri, H.; Tellambura, C. Gamma Codes: A Low-Overhead Linear-Complexity Network Coding Solution. In Proceedings of the 2012 International Symposium on Network Coding (NetCod), Cambridge, MA, USA, 29–30 June 2012; pp. 125–130. [Google Scholar]
  34. Mahdaviani, K.; Yazdani, R.; Ardakani, M. Linear-Complexity Overhead-Optimized Random Linear Network Codes. arXiv 2013, arXiv:1311.2123. [Google Scholar]
  35. Yang, S.; Tang, B. From LDPC to chunked network codes. In Proceedings of the 2014 IEEE Information Theory Workshop on Information Theory (ITW), Hobart, TAS, Australia, 2–5 November 2014; pp. 406–410. [Google Scholar]
  36. Tang, B.; Yang, S. An LDPC Approach for Chunked Network Codes. IEEE/ACM Trans. Netw. 2018, 26, 605–617. [Google Scholar] [CrossRef]
  37. Yang, S.; Yeung, R.W. Coding for a network coded fountain. In Proceedings of the 2011 IEEE International Symposium on Information Theory (ISIT), St. Petersburg, Russia, 31 July–5 August 2011; pp. 2647–2651. [Google Scholar]
  38. Yang, S.; Yeung, R.W. Batched Sparse Codes. IEEE Trans. Inf. Theory 2014, 60, 5322–5346. [Google Scholar] [CrossRef] [Green Version]
  39. Yang, S.; Yeung, R.W. BATS Codes: Theory and Practice; Synthesis Lectures on Communication Networks; Morgan & Claypool Publishers: San Rafael, CA, USA, 2017. [Google Scholar]
  40. Yang, S.; Ho, S.W.; Meng, J.; Yang, E.H. Capacity Analysis of Linear Operator Channels Over Finite Fields. IEEE Trans. Inf. Theory 2014, 60, 4880–4901. [Google Scholar] [CrossRef] [Green Version]
  41. Zhou, Q.; Yang, S.; Yin, H.H.F.; Tang, B. On BATS Codes with Variable Batch Sizes. IEEE Commun. Lett. 2017, 21, 1917–1920. [Google Scholar] [CrossRef]
  42. Huang, Q.; Sun, K.; Li, X.; Wu, D.O. Just FUN: A Joint Fountain Coding and Network Coding Approach to Loss-Tolerant Information Spreading. In Proceedings of the 15th ACM International Symposium on Mobile Ad Hoc Networking and Computing, Philadelphia, PA, USA, 11–14 August 2014; pp. 83–92. [Google Scholar]
  43. Yin, H.H.F.; Ng, K.H.; Wang, X.; Cao, Q. On the Minimum Delay of Block Interleaver for Batched Network Codes. In Proceedings of the 2019 IEEE International Symposium on Information Theory (ISIT), Paris, France, 7–12 July 2019; pp. 1957–1961. [Google Scholar]
  44. Yin, H.H.F.; Ng, K.H.; Wang, X.; Cao, Q.; Ng, L.K.L. On the Memory Requirements of Block Interleaver for Batched Network Codes. In Proceedings of the 2020 IEEE International Symposium on Information Theory (ISIT), Angeles, CA, USA, 21–26 June 2020; pp. 1658–1663. [Google Scholar]
  45. Zhou, Z.; Li, C.; Yang, S.; Guang, X. Practical Inner Codes for BATS Codes in Multi-Hop Wireless Networks. IEEE Trans. Veh. Technol. 2019, 68, 2751–2762. [Google Scholar] [CrossRef] [Green Version]
  46. Zhou, Z.; Kang, J.; Zhou, L. Joint BATS Code and Periodic Scheduling in Multihop Wireless Networks. IEEE Access 2020, 8, 29690–29701. [Google Scholar] [CrossRef]
  47. Yang, S.; Yeung, R.W.; Cheung, J.H.F.; Yin, H.H.F. BATS: Network Coding in Action. In Proceedings of the Annual Allerton Conference on Communication Control and Computing, Monticello, IL, USA, 30 September–3 October 2014; pp. 1204–1211. [Google Scholar]
  48. Tang, B.; Yang, S.; Ye, B.; Guo, S.; Lu, S. Near-Optimal One-Sided Scheduling for Coded Segmented Network Coding. IEEE Trans. Comput. 2016, 65, 929–939. [Google Scholar] [CrossRef]
  49. Yin, H.H.F.; Yang, S.; Zhou, Q.; Yung, L.M.L. Adaptive Recoding for BATS Codes. In Proceedings of the 2016 IEEE International Symposium on Information Theory (ISIT), Barcelona, Spain, 10–15 July 2016; pp. 2349–2353. [Google Scholar]
  50. Yin, H.H.F.; Tang, B.; Ng, K.H.; Yang, S.; Wang, X.; Zhou, Q. A Unified Adaptive Recoding Framework for Batched Network Coding. IEEE J. Sel. Areas Inf. Theory 2021, 2, 1150–1164. [Google Scholar] [CrossRef]
  51. Yin, H.H.F.; Ng, K.H. Impact of Packet Loss Rate Estimation on Blockwise Adaptive Recoding for Batched Network Coding. In Proceedings of the 2021 IEEE International Symposium on Information Theory (ISIT), Melbourne, VIC, Australia, 12–20 July 2021; pp. 1415–1420. [Google Scholar]
  52. Yin, H.H.F.; Ng, K.H.; Zhong, A.Z.; Yeung, R.W.; Yang, S.; Chan, I.Y.Y. Intrablock Interleaving for Batched Network Coding with Blockwise Adaptive Recoding. IEEE J. Sel. Areas Inf. Theory 2021, 2, 1135–1149. [Google Scholar] [CrossRef]
  53. Breidenthal, J.C. The Merits of Multi-Hop Communication in Deep Space. In Proceedings of the 2000 IEEE Aerospace Conference, Big Sky, MT, USA, 18–25 March 2000; Volume 1, pp. 211–222. [Google Scholar]
  54. Zhao, H.; Dong, G.; Li, H. Simplified BATS Codes for Deep Space Multihop Networks. In Proceedings of the 2016 IEEE Information Technology Networking, Electronic and Automation Control Conference, Chongqing, China, 20–22 May 2016; pp. 311–314. [Google Scholar]
  55. Yeung, R.W.; Dong, G.; Zhu, J.; Li, H.; Yang, S.; Chen, C. Space Communication and BATS Codes: A Marriage Made in Heaven. J. Deep. Space Explor. 2018, 5, 129–139. [Google Scholar]
  56. Sozer, E.M.; Stojanovic, M.; Proakis, J.G. Underwater Acoustic Networks. IEEE J. Ocean. Eng. 2000, 25, 72–83. [Google Scholar] [CrossRef]
  57. Yang, S.; Ma, J.; Huang, X. Multi-Hop Underwater Acoustic Networks Based on BATS Codes. In Proceedings of the 13th International Conference on Underwater Networks & Systems, Shenzhen, China, 3–5 December 2018; pp. 30:1–30:5. [Google Scholar]
  58. Sprea, N.; Bashir, M.; Truhachev, D.; Srinivas, K.V.; Schlegel, C.; Sacchi, C. BATS Coding for Underwater Acoustic Communication Networks. In Proceedings of the OCEANS 2019, Marseille, France, 17–20 June 2019; pp. 1–10. [Google Scholar]
  59. Yin, H.H.F.; Yeung, R.W.; Yang, S. A Protocol Design Paradigm for Batched Sparse Codes. Entropy 2020, 22, 790. [Google Scholar] [CrossRef] [PubMed]
  60. Yin, H.H.F.; Tahernia, M. On the Design of Timeout-Based Recoders for Batched Network Codes in the MU-MIMO Regime. In Proceedings of the 2022 IEEE Region 10 Conference (TENCON), Hong Kong, China, 1–4 November 2022. [Google Scholar]
  61. Qing, J.; Yin, H.H.F.; Yeung, R.W. Enhancing the Decoding Rates of BATS Codes by Learning with Guided Information. In Proceedings of the 2022 IEEE International Symposium on Information Theory (ISIT), Espoo, Finland, 26 June–1 July 2022; pp. 37–42. [Google Scholar]
  62. Yang, S.; Zhou, Q. Tree Analysis of BATS Codes. IEEE Commun. Lett. 2016, 20, 37–40. [Google Scholar] [CrossRef]
  63. Yang, S.; Ng, T.C.; Yeung, R.W. Finite-Length Analysis of BATS Codes. IEEE Trans. Inf. Theory 2018, 64, 322–348. [Google Scholar] [CrossRef] [Green Version]
  64. Yang, J.; Shi, Z.; Wang, C.; Ji, J. Design of Optimized Sliding-Window BATS Codes. IEEE Commun. Lett. 2019, 23, 410–413. [Google Scholar] [CrossRef]
  65. Xu, X.; Zeng, Y.; Guan, Y.L.; Yuan, L. Expanding-Window BATS Code for Scalable Video Multicasting Over Erasure Networks. IEEE Trans. Multimed. 2018, 20, 271–281. [Google Scholar] [CrossRef]
  66. Yin, H.H.F.; Wong, H.W.H.; Tahernia, M.; Qing, J. Packet Size Optimization for Batched Network Coding. In Proceedings of the 2022 IEEE International Symposium on Information Theory (ISIT), Espoo, Finland, 26 June–1 July 2022; pp. 1584–1589. [Google Scholar]
  67. Yang, S.; Yeung, R.W. Network Communication Protocol Design from the Perspective of Batched Network Coding. IEEE Commun. Mag. 2022, 60, 89–93. [Google Scholar] [CrossRef]
  68. Shokrollahi, A.; Luby, M. Raptor Codes. In Foundations and Trends in Communications and Information Theory; Now Publishers Inc.: Hanover, MA, USA, 2011; Volume 6. [Google Scholar]
  69. Shokrollahi, A.; Lassen, S.; Karp, R. Systems and Processes for Decoding Chain Reaction Codes through Inactivation. U.S. Patent 6,856,263, 15 February 2005. [Google Scholar]
  70. Xu, X.; Guan, Y.L.; Zeng, Y. Batched Network Coding with Adaptive Recoding for Multi-Hop Erasure Channels with Memory. IEEE Trans. Commun. 2018, 66, 1042–1052. [Google Scholar] [CrossRef]
  71. Yin, H.H.F.; Tahernia, M. Multi-Phase Recoding for Batched Network Coding. In Proceedings of the 2022 IEEE Information Theory Workshop on Information Theory (ITW), Mumbai, India, 6–9 November 2022; pp. 25–30. [Google Scholar]
  72. Ye, F.; Roy, S.; Wang, H. Efficient Data Dissemination in Vehicular Ad Hoc Networks. IEEE J. Sel. Areas Commun. (JSAC) 2012, 30, 769–779. [Google Scholar] [CrossRef]
  73. Yin, H.H.F.; Xu, X.; Ng, K.H.; Guan, Y.L.; Yeung, R.W. Packet Efficiency of BATS Coding on Wireless Relay Network with Overhearing. In Proceedings of the 2019 IEEE International Symposium on Information Theory (ISIT), Paris, France, 7–12 July 2019; pp. 1967–1971. [Google Scholar]
  74. Lucani, D.E.; Médard, M.; Stojanovic, M. Random Linear Network Coding for Time-Division Duplexing: Field Size Considerations. In Proceedings of the 2009 IEEE Global Telecommunications Conference, Honolulu, HI, USA, 30 November–4 December 2009; pp. 1–6. [Google Scholar]
  75. Yin, H.H.F.; Xu, X.; Ng, K.H.; Guan, Y.L.; Yeung, R.W. Analysis of Innovative Rank of Batched Network Codes for Wireless Relay Networks. In Proceedings of the 2021 IEEE Information Theory Workshop on Information Theory (ITW), Kanazawa, Japan, 17–21 October 2021. [Google Scholar]
  76. Zhang, C.; Tang, B.; Ye, B.; Lu, S. An efficient chunked network code based transmission scheme in wireless networks. In Proceedings of the 2017 IEEE International Conference on Communications (ICC), Paris, France, 21–25 May 2017; pp. 1–6. [Google Scholar]
  77. Gilbert, E.N. Capacity of a Burst-Noise Channel. Bell Syst. Tech. J. 1960, 39, 1253–1265. [Google Scholar] [CrossRef]
  78. Elliott, E.O. Estimates of Error Rates for Codes on Burst-Noise Channels. Bell Syst. Tech. J. 1963, 42, 1977–1997. [Google Scholar] [CrossRef]
  79. NIST Digital Library of Mathematical Functions. Release 1.1.0. Available online: http://dlmf.nist.gov/ (accessed on 15 December 2020).
  80. Galassi, M.; Davies, J.; Theiler, J.; Gough, B.; Jungman, G.; Booth, M.; Rossi, F. GNU Scientific Library Reference Manual, 3rd ed.; Network Theory Ltd.: London, UK, 2002. [Google Scholar]
  81. Wang, J.; Jia, Z.; Yin, H.H.F.; Yang, S. Small-Sample Inferred Adaptive Recoding for Batched Network Coding. In Proceedings of the 2021 IEEE International Symposium on Information Theory (ISIT), Melbourne, VIC, Australia, 12–20 July 2021; pp. 1427–1432. [Google Scholar]
  82. Dong, Y.; Jin, S.; Yang, S.; Yin, H.H.F. Network Utility Maximization for BATS Code Enabled Multihop Wireless Networks. In Proceedings of the 2020 IEEE International Conference on Communications (ICC), Dublin, Ireland, 7–11 June 2020. [Google Scholar]
  83. Dong, Y.; Jin, S.; Chen, Y.; Yang, S.; Yin, H.H.F. Utility Maximization for Multihop Networks Employing BATS Codes with Adaptive Recoding. IEEE J. Sel. Areas Inf. Theory 2021, 2, 1120–1134. [Google Scholar] [CrossRef]
  84. Fredman, M.L.; Tarjan, R.E. Fibonacci Heaps and Their Uses in Improved Network Optimization Algorithms. J. ACM (JACM) 1987, 34, 596–615. [Google Scholar] [CrossRef]
  85. Musser, D. Introspective Sorting and Selection Algorithms. Softw. Pract. Exp. 1997, 27, 983–993. [Google Scholar] [CrossRef]
  86. Hoare, C.A.R. Algorithm 65: Find. Commun. ACM 1961, 4, 321–322. [Google Scholar] [CrossRef]
  87. Blum, M.; Floyd, R.W.; Pratt, V.; Rivest, R.L.; Tarjan, R.E. Time Bounds for Selection. J. Comput. Syst. Sci. 1973, 7, 448–461. [Google Scholar] [CrossRef] [Green Version]
  88. Seward, H.H. Information Sorting in the Application of Electronic Digital Computers to Business Operations. Master’s Thesis, MIT Digital Computer Laboratory, Cambridge, UK, 1954. Report R-232. [Google Scholar]
  89. Higham, N.J. Accuracy and Stability of Numerical Algorithms: Second Edition; Other Titles in Applied Mathematics; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 2002. [Google Scholar]
  90. Hodges, J.L., Jr.; Lehmann, E.L. Some Problems in Minimax Point Estimation. Ann. Math. Stat. 1950, 21, 182–197. [Google Scholar] [CrossRef]
  91. Steinhaus, H. The Problem of Estimation. Ann. Math. Stat. 1957, 28, 633–648. [Google Scholar] [CrossRef]
  92. Qing, J.; Leong, P.H.W.; Yeung, R.W. Performance Analysis and Optimal Design of BATS Code: A Hardware Perspective. IEEE Trans. Veh. Technol. 2023, 1–14. [Google Scholar] [CrossRef]
  93. Yang, S.; Yeung, W.H.; Chao, T.I.; Lee, K.H.; Ho, C.I. Hardware Acceleration for Batched Sparse Codes. U.S. Patent 10,237,782, 19 March 2019. [Google Scholar]
  94. Yin, H.F.H.; Yang, S.; Yeung, W.H.R. Loss-Resilient Protocols for Communication Networks. U.S. Patent 10,425,192, 24 September 2019. [Google Scholar]
  95. Yin, H.F.H.; Ng, K.H.; Zhong, Z.; Yeung, R.W.H.; Yang, S. Compatible Packet Separation for Communication Networks. U.S. Patent 11,452,003, 20 September 2022. [Google Scholar]
  96. Yin, H.F.H. Recoding Optimizations in Batched Sparse Codes. Ph.D. Thesis, The Chinese University of Hong Kong, Hong Kong, China, 2019. [Google Scholar]
  97. Chazelle, B. The Soft Heap: An Approximate Priority Queue with Optimal Error Rate. J. ACM (JACM) 2000, 47, 16. [Google Scholar] [CrossRef]
  98. Gao, R.; Kleywegt, A.J. Distributionally Robust Stochastic Optimization with Wasserstein Distance. Math. Oper. Res. 2023, 48, 603–655. [Google Scholar] [CrossRef]
Figure 1. A three-hop line network. Network links only exist between two neighbouring nodes.
Figure 1. A three-hop line network. Network links only exist between two neighbouring nodes.
Entropy 25 01054 g001
Figure 2. A Gilbert–Elliott (GE) model. In each state, there is an independent event to decide whether a packet is lost or not.
Figure 2. A Gilbert–Elliott (GE) model. In each state, there is an independent event to decide whether a packet is lost or not.
Entropy 25 01054 g002
Figure 3. This figure illustrates the idea of modifying the output of an approximation scheme with two batches, where L = { 1 , 2 } , t max L r 1 + r 2 and t max 1 , t max 2 t max L . The red and blue dots represent the optimal and approximate solutions on the facet t 1 + t 2 = t max L , respectively. Algorithm 1 starts the search from an interior point ( r 1 , r 2 ) , while a modification approach starts the search from the blue dot.
Figure 3. This figure illustrates the idea of modifying the output of an approximation scheme with two batches, where L = { 1 , 2 } , t max L r 1 + r 2 and t max 1 , t max 2 t max L . The red and blue dots represent the optimal and approximate solutions on the facet t 1 + t 2 = t max L , respectively. Algorithm 1 starts the search from an interior point ( r 1 , r 2 ) , while a modification approach starts the search from the blue dot.
Entropy 25 01054 g003
Figure 4. The tabular appearance of the function β p ( t , r ) after introducing boundaries 0 and 1 s. The rows and columns correspond to t = 1 , 0 , 1 , and r = 0 , 1 , 2 , , respectively. The row above the line is β p ( 1 , · ) .
Figure 4. The tabular appearance of the function β p ( t , r ) after introducing boundaries 0 and 1 s. The rows and columns correspond to t = 1 , 0 , 1 , and r = 0 , 1 , 2 , , respectively. The row above the line is β p ( 1 , · ) .
Entropy 25 01054 g004
Figure 5. The figures illustrate the two stages of the table generation. The indices start from ( 1 , 0 ) . The first row has the index y = 1 , which is the row above the line. Compared to Figure 4, β p ( 0 , 1 ) = β p ( 1 , 2 ) = = 1 can be substituted in directly without using the relation (15). (a) The first stage of the table generation. The 1 and 0 s paddings are generated first. The solid and dashed arrows represent (13) and (14), respectively. (b) The second stage of the table generation. The 1 and 0 s paddings are kept. The arrows represent the recursive relation (15) with the B p function in-place.
Figure 5. The figures illustrate the two stages of the table generation. The indices start from ( 1 , 0 ) . The first row has the index y = 1 , which is the row above the line. Compared to Figure 4, β p ( 0 , 1 ) = β p ( 1 , 2 ) = = 1 can be substituted in directly without using the relation (15). (a) The first stage of the table generation. The 1 and 0 s paddings are generated first. The solid and dashed arrows represent (13) and (14), respectively. (b) The second stage of the table generation. The 1 and 0 s paddings are kept. The arrows represent the recursive relation (15) with the B p function in-place.
Entropy 25 01054 g005
Figure 6. The values prepared when t b = 2 . The asterisks represent the values that are not yet initialized.
Figure 6. The values prepared when t b = 2 . The asterisks represent the values that are not yet initialized.
Entropy 25 01054 g006
Figure 7. Adaptive recoding (AR) vs. baseline recoding (BR) in line networks of different lengths, batch sizes and packet loss rates.
Figure 7. Adaptive recoding (AR) vs. baseline recoding (BR) in line networks of different lengths, batch sizes and packet loss rates.
Entropy 25 01054 g007
Figure 8. The first 10 hops in Figure 7.
Figure 8. The first 10 hops in Figure 7.
Entropy 25 01054 g008
Figure 9. The effect of different block sizes with M = 8 .
Figure 9. The effect of different block sizes with M = 8 .
Entropy 25 01054 g009
Figure 10. The first 10 hops in Figure 9.
Figure 10. The first 10 hops in Figure 9.
Entropy 25 01054 g010
Figure 11. Approximation vs. optimal AR in line networks of different lengths, batch sizes and packet loss rates.
Figure 11. Approximation vs. optimal AR in line networks of different lengths, batch sizes and packet loss rates.
Entropy 25 01054 g011
Figure 12. The first 10 hops in Figure 11.
Figure 12. The first 10 hops in Figure 11.
Entropy 25 01054 g012
Figure 13. The values of β p ( t , r ) for r = 1 , 2 , 3 , 4 and t = 1 , 2 , with different p. The coloured numbers are the largest eight values smaller than 1.
Figure 13. The values of β p ( t , r ) for r = 1 , 2 , 3 , 4 and t = 1 , 2 , with different p. The coloured numbers are the largest eight values smaller than 1.
Entropy 25 01054 g013
Figure 14. The condition numbers of β p ( t , r ) for r = 1 , 2 , 3 , 4 and t = 1 , 2 , .
Figure 14. The condition numbers of β p ( t , r ) for r = 1 , 2 , 3 , 4 and t = 1 , 2 , .
Entropy 25 01054 g014
Figure 15. Throughput with inaccurate channel conditions.
Figure 15. Throughput with inaccurate channel conditions.
Entropy 25 01054 g015
Figure 16. The 25 and 75% percentiles of the estimation of p by different schemes where | L | = M = 4 in 1000 runs.
Figure 16. The 25 and 75% percentiles of the estimation of p by different schemes where | L | = M = 4 in 1000 runs.
Entropy 25 01054 g016
Figure 17. Throughput with estimated p via feedback where | L | = M = 4 .
Figure 17. Throughput with estimated p via feedback where | L | = M = 4 .
Entropy 25 01054 g017
Table 1. Terminology used for batched network coding (BNC).
Table 1. Terminology used for batched network coding (BNC).
TerminologyDescription
BatchA small set of coded packets.
Batch sizeThe number of coded packets in a batch.
Rank of a batchA measure of “useful” information (linearly independent packets) retained in the batch.
Expected rank of a batchThe expectation of the rank of the batch at the next network node.
Incoming rank distributionThe distribution of the ranks of the batches arriving at the current network node.
ThroughputThe expectation of the rank distribution of the batches arriving at the destination node.
RecodingThe network coding operations restricted to the packets belonging to the same batch.
Recoded packetsThe coded packets generated by recoding.
RecoderThe module that performs recoding.
Baseline recodingA strategy that generates the same number of recoded packets per batch.
Adaptive recodingA strategy that generates different number of recoded packets per batch.
BlockA set of batches.
Blockwise adaptive recodingApplying adaptive recoding block by block.
Table 2. Frequently used notations in this paper.
Table 2. Frequently used notations in this paper.
NotationDescription
MBatch size.
r b The rank of the batch b.
t b The number of recoded packets for batch b.
E ( r b , t b ) The expected rank of batch b when its rank is r b at the current node and t b recoded packets are sent.
( h 0 , , h M ) The incoming rank distribution.
pThe packet loss rate in the independent packet loss model.
L A block.
t max L The total number of recoded packets in block L .
t max b The maximum number of recoded packets allowed for batch b.
Binom ( n , p ) The binomial distribution.
B p ( t , i ) The probability mass function of the binomial distribution Binom ( t , 1 p ) .
β p ( t , r ) The sum of the first r probability masses of Binom ( t , 1 p ) .
Beta ( a , b ) The beta distribution.
I x ( a , b ) The regularized incomplete beta function.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Yin, H.H.F.; Yang, S.; Zhou, Q.; Yung, L.M.L.; Ng, K.H. BAR: Blockwise Adaptive Recoding for Batched Network Coding. Entropy 2023, 25, 1054. https://doi.org/10.3390/e25071054

AMA Style

Yin HHF, Yang S, Zhou Q, Yung LML, Ng KH. BAR: Blockwise Adaptive Recoding for Batched Network Coding. Entropy. 2023; 25(7):1054. https://doi.org/10.3390/e25071054

Chicago/Turabian Style

Yin, Hoover H. F., Shenghao Yang, Qiaoqiao Zhou, Lily M. L. Yung, and Ka Hei Ng. 2023. "BAR: Blockwise Adaptive Recoding for Batched Network Coding" Entropy 25, no. 7: 1054. https://doi.org/10.3390/e25071054

APA Style

Yin, H. H. F., Yang, S., Zhou, Q., Yung, L. M. L., & Ng, K. H. (2023). BAR: Blockwise Adaptive Recoding for Batched Network Coding. Entropy, 25(7), 1054. https://doi.org/10.3390/e25071054

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop