From: Maximum-parsimony haplotype frequencies inference based on a joint constrained sparse representation of pooled DNA
1
Set k = 0
2
Set ρ > 0
3
Set X0 = 0, Z0 = 0, U 1 0 = 0 , U 2 0 = 0
4
U4 = (I−HT(I+H HT)−1H)
5
Repeat
6
k = k + 1
7
u 3 T = 1 T U 4 H H U 2 k − A + U 1 k − Z k − 1 T 1 T U 4 1
8
X k + 1 = U 4 U 1 k − Z k + H T ( U 2 k − A ) − 1 u 3 T
9
Z k + 1 =max Shrink X k + 1 + 1 ρ U 1 k , 1 ρ , 0
10
U 1 k + 1 = U 1 k + X k + 1 − Z k + 1
11
U 2 k + 1 = U 2 k +H X k + 1 −A
12
until convergence