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

Table 2

Comparison of the combination of AlgSR1 and AlgBP methods and the block partitioning method.

Method	Sizes , ℓ,	CPU time
BP	50, 15, 10	0.0037	2.0011e − 15	2.0385e − 14	3.773e − 15	1.6463e − 14
HSR1 + BP	50, 15, 10	0.0592	2.6951e − 14	2.9305e − 13	1.3171e − 13	4.0769e − 14
BP	50, 20, 10	3.4243e − 04	6.9626e − 16	2.8774e − 14	3.4226e − 15	2.7706e − 14
HSR1 + BP	50, 20, 10	0.0305	5.3084e − 14	4.1354e − 13	3.1072e − 13	9.9898e − 14
BP	50, 15, 500	0.0019	8.0717e − 16	1.6354e − 14	3.7149e − 15	7.0839e − 15
HSR1 + BP	50, 15, 500	0.0411	2.1453e − 14	2.2518e − 13	7.0401e − 14	3.2624e − 14
BP	50, 20, 500	3.2256e − 04	6.8218e − 16	2.1407e − 14	3.2663e − 15	6.0107e − 15
HSR1 + BP	50, 20, 500	0.0288	6.029e − 15	1.3385e − 13	2.6565e − 14	1.3927e − 14
BP	250, 15, 500	0.0048	3.4724e − 15	2.3693e − 13	1.5467e − 14	7.4894e − 14
HSR1 + BP	250, 15, 500	21.9978	1.8868e − 12	3.4807e − 10	1.154e − 11	4.0825e − 12
BP	400, 20, 500	0.0132	4.731e − 15	3.7769e − 13	2.6494e − 14	9.909e − 14
HSR1 + BP	400, 20, 500	306.1949	7.3251e − 12	2.7048e − 09	6.4634e − 11	2.0001e − 11