A Simple Proof of Threshold Saturation for Coupled Scalar Recursions
Henry D. Pfister
joint work with Arvind Yedla, Yung-Yih Jian, and Phong S. Nguyen
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions - - PowerPoint PPT Presentation
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions Henry D. Pfister joint work with Arvind Yedla, Yung-Yih Jian, and Phong S. Nguyen Stanford University Electrical Engineering August 17th, 2012 LDPC Codes Spatial Coupling
joint work with Arvind Yedla, Yung-Yih Jian, and Phong S. Nguyen
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 2 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 3 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 4 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Regular (l❀ r): H has l ones per column and r ones per row ✎ Irregular: number of ones given by degree distribution (✕❀ ✚) ✎ Introduced by Gallager in 1960, but largely forgotten until 1995
✎ An edge connects check node i to bit node j if Hij = 1 ✎ Naturally leads to message-passing iterative (MPI) decoding
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 5 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Low-complexity message-passing decoder introduced by Gallager ✎ Local inference assuming all input messages are independent
✎ Tracks distribution of messages during iterative decoding ✎ BP noise threshold can be computed via DE ✎ Long codes decode almost surely if DE predicts success
✎ Optimum decoder that chooses the most likely codeword ✎ Infeasible in practice due to enormous number of codewords ✎ MAP noise threshold can be bounded using EXIT curves
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 6 / 33 LDPC Codes Spatial Coupling Simple Proof
✎
✎
✎ Msgs passed along edges in phases: bit-to-check and check-to-bit ✎ Each output message depends only on the input messages ✎ All messages are either correct value or erasure
✎ Bits pass the correct value unless all other inputs are erased ✎ Checks pass the correct value only if all other inputs are correct
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 7 / 33 LDPC Codes Spatial Coupling Simple Proof
1
2
✎ Illustrates decoding from the perspective of a single bit-node ✎ For long random LDPC codes, the graph is typically a tree ✎ Allows density evolution to track message erasure probability ✎ If x❂y are erasure prob. of bit/check output messages, then y y y y y4 x x x x 1 − (1 − x)4
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 8 / 33 LDPC Codes Spatial Coupling Simple Proof
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.1 0.2 0.3 0.4 0.5 0.6 0.7 Erasure probability xi Erasure probability xi+1 Density Evolution for a (3,4) LDPC Code Iter = 16 = 0.600 ε
✎ Easily computed numerically for each code ensemble
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 9 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Ex. hBP(✧) = (x∞(✧))4 ✎ Equals 0 below BP thresh ✎ Upper bounds MAP EXIT
✎ Equals 0 below MAP thresh ✎ Area underneath equals rate
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 10 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Magnetism caused by alignment of electron spins ✛i ✷ {+1❀ −1} ✎ The system energy in an external field is modeled by
✎ In equilibrium, the configuration probability is approximated by
✎ Spin systems are mathematically similar to binary inference ✎ Pairwise correlations in a binary vector controlled by Jij ✎ Observations encoded into the local magnetic fields hi ✎ The minimium-energy configuration is maximum a posteriori
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 11 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Inverse temperature ☞ = 1❂T scales coupling and field strength ✎ At high temperature (☞ → 0), spin system resembles a liquid ✎ At low temperature (☞ → ∞), it can freeze into a ground state ✎ The transition can be very complicated
✎ Code defined using generalized coupling coefficients J☛ ✎ Codewords are ordered crystalline structures ✎ Field hi is a function of Yi and channel parameter ✎ System is a supercooled liquid between BP and MAP threshold ✎ Correct answer (crystalline state) has minimum energy w.h.p. ✎ Spontaneous crystallization (i.e., decoding) does not occur w.h.p.
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 12 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 13 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 14 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 14 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 14 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 14 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 14 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 15 / 33 LDPC Codes Spatial Coupling Simple Proof
5 10 15 20 25 30 35 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45
Protograph Bit Node Message Erasure Probability
i
w−1
w−1
i+j −k
l−1
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 15 / 33 LDPC Codes Spatial Coupling Simple Proof
5 10 15 20 25 30 35 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45
Protograph Bit Node Message Erasure Probability
i
w−1
w−1
i+j −k
l−1
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 15 / 33 LDPC Codes Spatial Coupling Simple Proof
5 10 15 20 25 30 35 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45
Protograph Bit Node Message Erasure Probability
i
w−1
w−1
i+j −k
l−1
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 15 / 33 LDPC Codes Spatial Coupling Simple Proof
5 10 15 20 25 30 35 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45
Protograph Bit Node Message Erasure Probability
i
w−1
w−1
i+j −k
l−1
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 15 / 33 LDPC Codes Spatial Coupling Simple Proof
5 10 15 20 25 30 35 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45
Protograph Bit Node Message Erasure Probability
i
w−1
w−1
i+j −k
l−1
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 15 / 33 LDPC Codes Spatial Coupling Simple Proof
5 10 15 20 25 30 35 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45
Protograph Bit Node Message Erasure Probability
i
w−1
w−1
i+j −k
l−1
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 15 / 33 LDPC Codes Spatial Coupling Simple Proof
5 10 15 20 25 30 35 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45
Protograph Bit Node Message Erasure Probability
i
w−1
w−1
i+j −k
l−1
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 16 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ As L → ∞, code rate increases to original rate ✎ BP Threshold decreases to MAP threshold of original
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 16 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ As L → ∞, code rate increases to original rate ✎ BP Threshold decreases to MAP threshold of original
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 16 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ As L → ∞, code rate increases to original rate ✎ BP Threshold decreases to MAP threshold of original
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 16 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ As L → ∞, code rate increases to original rate ✎ BP Threshold decreases to MAP threshold of original
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 16 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ As L → ∞, code rate increases to original rate ✎ BP Threshold decreases to MAP threshold of original
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 16 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ As L → ∞, code rate increases to original rate ✎ BP Threshold decreases to MAP threshold of original
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 17 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 18 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Show BP and MAP thresholds can be computed from potential
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 19 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Established many properties and tools used by later approaches
✎ Our use of potential functions was motivated by this paper
✎ Using a vector potential function in the continuous limit
✎ Informally reported a proof based on the continuous limit
✎ Can our simplified approach work here?
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 19 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Established many properties and tools used by later approaches
✎ Our use of potential functions was motivated by this paper
✎ Using a vector potential function in the continuous limit
✎ Informally reported a proof based on the continuous limit
✎ Can our simplified approach work here? yes (see Allerton 2012)
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 20 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Single recursion converges to 0 for ✧ below BP threshold
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 21 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 22 / 33 LDPC Codes Spatial Coupling Simple Proof
0 f (z; ✧)dz and G(x) =
0 g(z)dz.
s = sup {✧ ✷ [0❀ 1] | U ✵(x; ✧) ✕ 0 ✽ x ✷ [0❀ 1]}
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 23 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 23 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 23 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 23 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 23 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 23 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 23 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 23 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 24 / 33 LDPC Codes Spatial Coupling Simple Proof
✿✿✿ ✿✿✿
✙0 ✙ ✵ −L −2 −1 1 2 L
✿✿✿ ✿✿✿ ✿✿✿ ✿✿✿
✙ ✙ ✵
✿✿✿ ✿✿✿
✙ ✙ ✵
✿✿✿ ✿✿✿
✙ ✙ ✵
✿✿✿ ✿✿✿
✙ ✙ ✵
✿✿✿ ✿✿✿
✙ ✙ ✵
✿✿✿ ✿✿✿
✙ ✙ ✵ ✿✿✿ ✿✿✿ ✿✿✿ ✿✿✿ ✙
✿✿✿
✙
✿✿✿
✙
✿✿✿
✙
✿✿✿
✙ ✵
✿✿✿
✙ ✵
✿✿✿
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 24 / 33 LDPC Codes Spatial Coupling Simple Proof
✿✿✿ ✿✿✿
✙0 ✙ ✵ −L −2 −1 1 2 L
✿✿✿ ✿✿✿ ✿✿✿ ✿✿✿
✙−L ✙ ✵
−L
✿✿✿ ✿✿✿
✙−2 ✙ ✵
−2
✿✿✿ ✿✿✿
✙−1 ✙ ✵
−1
✿✿✿ ✿✿✿
✙1 ✙ ✵
1
✿✿✿ ✿✿✿
✙2 ✙ ✵
2
✿✿✿ ✿✿✿
✙L ✙ ✵
L
✿✿✿ ✿✿✿ ✿✿✿ ✿✿✿ ✙
✿✿✿
✙
✿✿✿
✙
✿✿✿
✙
✿✿✿
✙ ✵
✿✿✿
✙ ✵
✿✿✿
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 24 / 33 LDPC Codes Spatial Coupling Simple Proof
✿✿✿ ✿✿✿
✙0 ✙ ✵ −L −2 −1 1 2 L
✿✿✿ ✿✿✿ ✿✿✿ ✿✿✿
✙−L ✙ ✵
−L
✿✿✿ ✿✿✿
✙−2 ✙ ✵
−2
✿✿✿ ✿✿✿
✙−1 ✙ ✵
−1
✿✿✿ ✿✿✿
✙1 ✙ ✵
1
✿✿✿ ✿✿✿
✙2 ✙ ✵
2
✿✿✿ ✿✿✿
✙L ✙ ✵
L
✿✿✿ ✿✿✿ ✿✿✿ ✿✿✿ − L−2 − L−1 L+1 L+2 ✙−
L − 2
✿✿✿
✙−
L − 1
✿✿✿
✙L
+ 1
✿✿✿
✙L
+ 2
✿✿✿
✙ ✵
L + 1
✿✿✿
✙ ✵
L + 2
✿✿✿
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 25 / 33 LDPC Codes Spatial Coupling Simple Proof
i
w−1
w−1
i+j −k
i
i
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 25 / 33 LDPC Codes Spatial Coupling Simple Proof
i
w−1
w−1
i+j −k
i
2 ❝}
i
i
0 ❀ i ❃ ❜w−1 2 ❝
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 26 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 26 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 26 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 26 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 27 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 28 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 28 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 28 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 28 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 28 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 28 / 33 LDPC Codes Spatial Coupling Simple Proof
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 29 / 33 LDPC Codes Spatial Coupling Simple Proof
✵
✵✵
✵✵
✷ ❀
✵✵
✷ ❀
✵✵
❀
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 29 / 33 LDPC Codes Spatial Coupling Simple Proof
✵
✵✵
✵✵
✷ ❀
✵✵
✷ ❀
✵✵
❀
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 29 / 33 LDPC Codes Spatial Coupling Simple Proof
✵
✵✵
✵✵
✷ ❀
✵✵
✷ ❀
✵✵
❀
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 29 / 33 LDPC Codes Spatial Coupling Simple Proof
t✷[0❀1] ❦U ✵✵(x(t); ✧)❦∞ ❦Sx − x❦∞
t✷[0❀1] ❦U ✵✵(x(t); ✧)❦∞
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 30 / 33 LDPC Codes Spatial Coupling Simple Proof
✵
✎
✵
✎
✵
✎
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 30 / 33 LDPC Codes Spatial Coupling Simple Proof
✎
✵
✎
✵
✎
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 30 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Implies [x − A⊺f (Ag (x) ; ✧)]i g ✵(xi) ❃ 0 for some i ✷ ▲
✎
✵
✎
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 30 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Implies [x − A⊺f (Ag (x) ; ✧)]i g ✵(xi) ❃ 0 for some i ✷ ▲
✎
✵
✎
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 30 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Implies [x − A⊺f (Ag (x) ; ✧)]i g ✵(xi) ❃ 0 for some i ✷ ▲
✎ But, g ✵(x) ❃ 0 implies [A⊺f (Ag (x) ; ✧)]i ❁ xi ✎
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 30 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Implies [x − A⊺f (Ag (x) ; ✧)]i g ✵(xi) ❃ 0 for some i ✷ ▲
✎ But, g ✵(x) ❃ 0 implies [A⊺f (Ag (x) ; ✧)]i ❁ xi ✎ Implies one step of DE decreases xi and contradicts “x is a f.p.”
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 30 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Implies [x − A⊺f (Ag (x) ; ✧)]i g ✵(xi) ❃ 0 for some i ✷ ▲
✎ But, g ✵(x) ❃ 0 implies [A⊺f (Ag (x) ; ✧)]i ❁ xi ✎ Implies one step of DE decreases xi and contradicts “x is a f.p.”
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 30 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Implies [x − A⊺f (Ag (x) ; ✧)]i g ✵(xi) ❃ 0 for some i ✷ ▲
✎ But, g ✵(x) ❃ 0 implies [A⊺f (Ag (x) ; ✧)]i ❁ xi ✎ Implies one step of DE decreases xi and contradicts “x is a f.p.”
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 31 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ ✎
✎ ✎
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 31 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Irregular LDPC on the BEC and ISI channels with erasures ✎
✎ ✎
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 31 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Irregular LDPC on the BEC and ISI channels with erasures ✎ A GLDPC product code on the BEC and BSC w/binary msgs
✎ ✎
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 31 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Irregular LDPC on the BEC and ISI channels with erasures ✎ A GLDPC product code on the BEC and BSC w/binary msgs
✎ ✎
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 31 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Irregular LDPC on the BEC and ISI channels with erasures ✎ A GLDPC product code on the BEC and BSC w/binary msgs
✎ A noisy Slepian-Wolf problem with erasure noise ✎
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 31 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Irregular LDPC on the BEC and ISI channels with erasures ✎ A GLDPC product code on the BEC and BSC w/binary msgs
✎ A noisy Slepian-Wolf problem with erasure noise ✎ Multiple-access channel with erasures
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 32 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Powerful technique for designing codes and decoders ✎ Related to the statistical physics of supercooled liquids and
✎ Simple proof of threshold saturation for scalar recursions
✎ ✎ ✎ ✎
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 32 / 33 LDPC Codes Spatial Coupling Simple Proof
✎ Powerful technique for designing codes and decoders ✎ Related to the statistical physics of supercooled liquids and
✎ Simple proof of threshold saturation for scalar recursions
✎ Optimized LDPC codes are not universal ✎ But suboptimal decoding is the main problem ✎ Spatially-coupled joint decoders appear to be universal ✎ Observed in general and now proven for some erasure models
A Simple Proof of Threshold Saturation for Coupled Scalar Recursions 33 / 33 LDPC Codes Spatial Coupling Simple Proof