julia - Julia 主对角线下的总和

Question

如何在 julia 中不使用矩阵中的主对角线的情况下在主对角线下求和？

我正在使用 sum=tril(a)-1 但它在 julia 中不起作用。

我知道我需要一个面具，但我不知道如何使用它。

Answer 1

您正在寻找 LinearAlgebra module ，它是标准库的一部分，包含一个 tril 函数:

julia> using LinearAlgebra

julia> A = ones(5, 5)
5×5 Array{Float64,2}:
 1.0  1.0  1.0  1.0  1.0
 1.0  1.0  1.0  1.0  1.0
 1.0  1.0  1.0  1.0  1.0
 1.0  1.0  1.0  1.0  1.0
 1.0  1.0  1.0  1.0  1.0

julia> tril(A, -1)
5×5 Array{Float64,2}:
 0.0  0.0  0.0  0.0  0.0
 1.0  0.0  0.0  0.0  0.0
 1.0  1.0  0.0  0.0  0.0
 1.0  1.0  1.0  0.0  0.0
 1.0  1.0  1.0  1.0  0.0

julia> sum(tril(A, -1))
10.0