chapel - 如何在教堂中提取带或不带对角线的上三角矩阵

Question

最好的故事从我的矩阵开始A。

var A: [{1..4,1..4}] real = (
  (4, -30, 60, -35),
  (-30, 300, -675, 420),
  (60, -675, 1620, -1050),
  (-35, 420, -1050, 700)
  );

我记得在某个地方你可以得到上面的三角形，但我在Chapel 文档中找不到它你能把它设置为/排除对角线吗？

Answer 1

使用LinearAlgebra库：

use LinearAlgebra;

var A: [{1..4,1..4}] real = (
  (4, -30, 60, -35),
  (-30, 300, -675, 420),
  (60, -675, 1620, -1050),
  (-35, 420, -1050, 700)
  );

var upperTriangle = triu(A);

// Confirm it worked
writeln(isTriu(upperTriangle));