The following code doesn't make a lot of sense, is there a reason why the
backslash doesnt pivot for sparse matrices?
*julia> *A =[1.0 0.0 1.0
0.0 0.0 1.0
1.0 1.0 0.0]
3x3 Array{Float64,2}:
1.0 0.0 1.0
0.0 0.0 1.0
1.0 1.0 0.0
*julia>* inv(A)
3x3 Array{Float64,2}:
1.0 -1.0 -0.0
-1.0 1.0 1.0
0.0 1.0 0.0
*julia>* A\ones(3)
3-element Array{Float64,1}:
0.0
1.0
1.0
*julia>* sparse(A)\ones(3)
ERROR: ArgumentError: matrix has one or more zero pivots
in ldltfact at sparse/cholmod.jl:1246
in ldltfact at sparse/cholmod.jl:1253
in factorize at sparse/linalg.jl:849
in \ at linalg/generic.jl:326
