如何使用犰狳

Question

誰能告訴我什麼是以下操作的犰狳線性代數包相當於 L到得到正規拉普拉斯矩陣= d^-0.5 * A * d^-0.5

一般如何來計算^ 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; 

有沒有其他更簡單的方法？