Comment obtenir la matrice laplacian normalisée en utilisant Armadillo

Question

Quelqu'un peut-il s'il vous plaît me dire ce qui est l'équivalent de l'opération suivante dans le paquet d'algèbre linéaire Armadillo L = D^-0,5 * A * D^-0,5

En général comment pour calculer a^n ou a^-0,5 dans Armadillo où a est une matrice carrée

Answer 1

je peux penser à une façon de le faire

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; 

est-il une autre façon plus simple?