使用分块方法对矩阵进行对角化。
示例代码:
// 初始化矩阵 int n = 1000; mat A = randu(n, n);
// 分块参数 int b = 128;
// 分块对角化 for (int i = 0; i < n; i += b) { int j = i + b < n ? i + b : n; mat diagA = A.submat(i, i, j-1, j-1); eig_sym(diagA); }
// 对角化后的矩阵 eig_sym(A);
上一篇:Armadillo并发链式操作挂起
下一篇:ArmadilloC++ifft性能较差