Julia 中的稀疏复数矩阵

Question

在 Julia 中，我可以创建一个由零组成的稀疏矩阵：

julia> a = spzeros(2,2)
2x2 sparse matrix with 0 Float64 entries:

julia> a[1,1] = 1
1

julia> full(a)
2x2 Array{Float64,2}:
 1.0  0.0
 0.0  0.0

我可以创建一个复杂的矩阵：

julia> b = [ 1 ; im ]
2-element Array{Complex{Int64},1}:
 1+0im
 0+1im

如果我尝试将一个复数赋值给零的稀疏矩阵，我会得到一个错误：

julia> a[1,1] = im
ERROR: InexactError()
 in setindex! at sparse/sparsematrix.jl:1095

这与 Float64 参数化的 spzeros() 返回类型一致：

julia> typeof(a)
SparseMatrixCSC{Float64,Int64} (constructor with 1 method)

如何在 Julia 中创建复杂类型零的稀疏矩阵？

Answer 1

看看我们可以传递给什么 spzeros:

julia> methods(spzeros)
# 5 methods for generic function "spzeros":
spzeros(m::Integer,n::Integer) at sparse/sparsematrix.jl:406
spzeros(Tv::Type{T<:Top},m::Integer,n::Integer) at sparse/sparsematrix.jl:407
spzeros(Tv::Type{T<:Top},Ti::Type{T<:Top},m::Integer,n::Integer) at sparse/sparsematrix.jl:409
spzeros(m::Integer) at deprecated.jl:28
spzeros(Tv::Type{T<:Top},m::Integer) at deprecated.jl:28

我们看到我们应该能够将类型作为第一个参数传递：

julia> a = spzeros(Complex{Float64}, 2,2)
2x2 sparse matrix with 0 Complex{Float64} entries:

julia> full(a)
2x2 Array{Complex{Float64},2}:
 0.0+0.0im  0.0+0.0im
 0.0+0.0im  0.0+0.0im

julia> a[1,1] = 2+3.4im
2.0 + 3.4im

julia> a
2x2 sparse matrix with 1 Complex{Float64} entries:
    [1, 1]  =  2.0+3.4im

Julia 中的稀疏复数矩阵

matrix of sparse complex numbers in Julia

matrix

sparse-matrix

complex-numbers

julia