我必须计算大小为n
的{{1}}个对称矩阵的秩。但是,这样做时,我意识到对称矩阵的秩的计算时间比普通矩阵的秩的计算时间慢得多(请参见下面的代码)。
m
由于这个原因,我只计算了大小为# Part I
m = 77;
n = 20_000;
A = rand(m,m);
SA = Symmetric(A);
@btime rank($A); #result: 357.101 μs (9 allocations: 92.77 KiB)
@btime rank($SA); #result: 2.915 ms (13 allocations: 75.69 KiB)
的{{1}}个正方形矩阵的秩。完成该工作所需的计算时间约为8-9(s)。但是我对此不满意。
n
那么,有人知道如何在第二部分中加快我的代码的速度吗?预先谢谢你!