0

Is there a method to use the nalgebra glm methods on dynamically sized matrices The aim is to Find X such that X†SX = I: the current method is X=Us^{-1/2}U† where S=UsU†:

// in general these will be dynamically sized and read from a file
let N = 2;
let S = DMatrix::from_row_slice(N,N, &[1.0, 0.5017, 0.5017, 1.0]);
// S will always be symmetric
let decomp = SymmetricEigen::new(S);
let eigvals = decomp.eigenvalues;
// this is the code that does not work
let s = glm::inversesqrt(&eigvals);
JimGreen
  • 15
  • 4

0 Answers0