Open jondeuce opened 5 years ago
Also, this does seem to work for other factorizations, e.g. lu
:
julia> B = lu(sparse(Matrix(I,10,10)));
julia> isapprox(B\@views(z[:]), z)
true
My julia version info is following, in case it's relevant:
julia> versioninfo()
Julia Version 1.1.0
Commit 80516ca202 (2019-01-21 21:24 UTC)
Platform Info:
OS: Linux (x86_64-pc-linux-gnu)
CPU: AMD Phenom(tm) II X6 1100T Processor
WORD_SIZE: 64
LIBM: libopenlibm
LLVM: libLLVM-6.0.1 (ORCJIT, amdfam10)
Environment:
JULIA_NUM_THREADS = 24
JULIA_EDITOR = code
Nice to have a PR here.
I believe CHOLMOD can handle strides, I'll look
The following errors for
cholesky
factorizations:while the same operation works without the view:
I've chased the issue down to being a simple missing definition in
SuiteSparse
, where we have the following (incholmod.jl
):So, the analogous definitions for complex
StridedVector
s andStridedMatrix
s are missing. I can submit the PR to base myself, but I wanted to open an issue first in case there is a reason they aren't defined.