Missing umf_extract(UmfpackLU{Complex{Float64},Int64},) for lufact() applied to sparse complex matrices

[jasax]

Hello,

In the quest to get complex algebra working in julia in Win7-64 (following issue #7031), I turned to sparse matrices and lufact(). Apparently it was uncovered another glitch, which is the nonexistence of a method: "no method umf_extract(UmfpackLU{Complex{Float64},Int64},)"

So, it seems that no template for those types is present in umf_extract()...

Below is the test code and result.

Best Regards

Jose

N=4
J=0.0+1.0im

MR = rand(Float64,N,N)
MI = rand(Float64,N,N)

SpMR=sparse(MR)
SpMI=sparse(MI)
SpM=SpMR+J*SpMI
println("Sparse M \n",typeof(SpM))
println(SpM)

F=lufact(SpM)
println("F \n",typeof(F),"\n")

println(F[:L])
println(F[:U])
println(F[:p])
println(F[:q])
println(F[:Rs])

C:\Users...>julia bug4.jl
Sparse M
SparseMatrixCSC{Complex{Float64},Int64}
0.17324610206636093 + 0.556739526277584im 0.3636155270578807 + 0.779898130
0944341im 0.9462204048595892 + 0.0816916240409733im 0.32368934936414
506 + 0.7481843025943526im
0.9912417559106097 + 0.0051733113483452im 0.23271082818867894 + 0.50533292
47279121im 0.16763856188783266 + 0.07775290297860527im 0.73970414046621
24 + 0.243096819167115im
0.783090390333812 + 0.5863584102559694im 0.2521944895936421 + 0.201648277
29072243im 0.8912116317502279 + 0.13636989345048156im 0.27938051396603
88 + 0.4982856834741676im
0.4573921542170194 + 0.4777090518686091im 0.6704835515875494 + 0.388381970
5387357im 0.004182163343978029 + 0.7094696885852512im 0.76423637986060
87 + 0.12410491260558443im

F
UmfpackLU{Complex{Float64},Int64}

ERROR: no method umf_extract(UmfpackLU{Complex{Float64},Int64},)
in getindex at linalg/umfpack.jl:285
WARNING: backtraces on your platform are often misleading or partially incorrect

[ViralBShah]

Looking into this now.

[jasax]

Thanks, Viral.

Jose

On 6/3/14, Viral B. Shah [email protected] wrote:

Closed #7098 via 324acc0.

---

Reply to this email directly or view it on GitHub:
#7098 (comment)