如何使用犰狳得到归一化拉普拉斯矩阵

Question

任何人都可以告诉我Armadillo线性代数包中以下操作的等价物  L = D ^ -0.5 * A * D ^ -0.5

一般如何在犰狳中计算A ^ n或A ^ -0.5，其中A是方阵

Answer 1

我可以想到一种方法来做到这一点

 mat K1,K2;
 K1.load(argv[1],auto_detect);
 colvec c = sum(K1,1);
 mat D = diagmat(c);
 mat D1 = pow(inv(D),0.5);
 mat I(10,10);I.eye();
 mat L = I - D1*K1*D1;

还有其他更简单的方法吗？