Research Article

Smoothing and Regularization with Modified Sparse Approximate Inverses

Table 4

Optimal reconstruction error 𝑥 2 ̃ 𝑥 2 2 for the 1D blur operator 𝐻 4 1 and the original data 𝑥 2 . The problem has size 𝑛 = 1 0 3 and is affected by random noise of order 0.1%. 𝜌 = 1 0 2 .

Regularization method Optimal value Reached at iter.

Tikhonov ( 𝛾 = 0 . 3 ) 0.1286
CG ( 𝐻 4 1 ) 0.1240 8

PCG ( 𝐻 4 1 ) 𝜌 𝑅 3 𝐺 1 no boundary correction 0.1210 130
𝜏 = 0 . 2 5 0.0433 58
𝜏 = 0 . 2 5 , middle correction 0.0263 52