Comment extraire la matrice triangulaire supérieure w ou w/o diagonal dans la chapelle

Question

Les meilleures histoires commencent avec ma matrice A.

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

Je me souviens quelque part, vous pouvez obtenir le triangle supérieur, mais je ne peux pas le trouver dans le Chapel docs Pouvez-vous le mettre dans/exclure la diagonale?

Answer 1

Utilisation de la bibliothèque 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));