Hello
I was going to investigate whether re-using the structure of a sparse
cholesky could save me some time.
Trying it, I got this:
julia> C = cholfact(K)
Base.SparseMatrix.CHOLMOD.Factor{Float64}
type: LLt
method: supernodal
maxnnz: 0
nnz: 4117190
julia> D = cholfact!(C, K)
ERROR: MethodError: `cholfact!` has no method matching cholfact!(::Base.
SparseMatrix.CHOLMOD.Factor{Float64}, ::SparseMatrixCSC{Float64,Int64})
help?> cholfact!
search: cholfact! cholfact
.. cholfact!(A [,LU=:U [,pivot=Val{false}]][;tol=-1.0]) -> Cholesky
``cholfact!`` is the same as :func:`cholfact`, but saves space by
overwriting the input ``A``, instead of creating a copy. ``cholfact!`` can
also reuse the symbolic factorization from a different matrix ``F`` with
the same structure when used as: ``cholfact!(F::CholmodFactor, A)``.
julia> VERSION
v"0.5.0-dev+749"
So I have a `Factor` object, but it looks like `cholfact!` wants a
`CholmodFactor`. Is this some kind of hangover from previous development,
or have I missed something?
Cheers
Matthew
