Computing the Pseudoinverse of Specific Toeplitz Matrices Using Rank-One Updates

Table 1

Comparison of the rank-one updates method and the BP method.

Method	Sizes , ℓ,	CPU time
BP	50, 15, 10	0.0025	2.0011e − 15	2.0385e − 14	3.773e − 15	1.6463e − 14
SR1	50, 15, 10	0.0490	2.2762e − 14	5.5984e − 13	1.2594e − 14	4.0623e − 13
BP	50, 20, 10	0.0028	6.9626e − 16	2.8774e − 14	3.4226e − 15	2.7706e − 14
SR1	50, 20, 10	0.0364	6.16e − 14	2.0705e − 12	1.3341e − 14	1.3718e − 12
BP	50, 15, 500	0.0029	8.0717e − 16	1.6354e − 14	3.7149e − 15	7.0839e − 15
SR1	50, 15, 500	0.0501	1.9666e − 14	4.8611e − 13	8.979e − 15	3.5898e − 13
BP	50, 20, 500	0.0016	6.8218e − 16	2.1407e − 14	3.2663e − 15	6.0107e − 15
SR1	50, 20, 500	0.0362	2.0844e − 14	5.3265e − 13	7.9594e − 15	4.1980e − 13
BP	250, 15, 500	0.0059	3.4724e − 15	2.3693e − 13	1.5467e − 14	7.4894e − 14
SR1	250, 15, 500	23.9762	7.81e − 12	3.6574e − 11	1.2707e − 13	2.2523e − 11
BP	400, 20, 500	0.0135	4.731e − 15	3.7769e − 13	2.6494e − 14	9.909e − 14
SR1	400, 20, 500	320.0578	3.8409e − 11	1.4246e − 10	2.2311e − 13	8.0923e − 11