JuliaLang/LinearAlgebra.jl

Linear Algebra operations with matrices of abstract element types are promoted incorrectly

Aberta

#287 aberto em 2 de dez. de 2015

 (2 comentários) (0 reação) (1 responsável)Julia (62 forks)github user discovery
bignumsbughelp wanted

Métricas do repositório

Stars
 (77 estrelas)
Métricas de merge de PR
 (Mesclagem média 8d 17h) (5 fundiu PRs em 30d)

Description

Follow up to JuliaLang/LinearAlgebra.jl#285.

Ex.

julia> a = Real[1.0 big(realmax(Float64))+1
       1.0 big(1.0)]
2x2 Array{Real,2}:
 1.0  1.797693134862315708145274237317043567980705675258449965989174768031572607800285e+308
 1.0  1.000000000000000000000000000000000000000000000000000000000000000000000000000000

julia> eigfact(a)
Base.LinAlg.Eigen{Float64,Float64,Array{Float64,2},Array{Float64,1}}([1.3407807929942597e154,-1.3407807929942597e154],2x2 Array{Float64,2}:
 1.0           -1.0
 7.45834e-155   7.45834e-155)

To be correct, this would have to widen the element type in a manner similar to map.

cc: @andreasnoack

Guia do colaborador