Skip to content

Issue with passed layout in default cholesky!_layout method for Symmetric matrices #244

@DanielVandH

Description

@DanielVandH

function cholesky!_layout(layout, axes, A::RealHermSymComplexHerm, ::CNoPivot; check::Bool = true)
C, info = _chol!(layout, A.data, A.uplo == 'U' ? UpperTriangular : LowerTriangular)
check && LinearAlgebra.checkpositivedefinite(info)
return Cholesky(C.data, A.uplo, info)
end

Should the call to _chol! be passing MemoryLayout(A.data) instead of simply layout? Something similar is done for reversecholesky here

https://github.com/JuliaLinearAlgebra/MatrixFactorizations.jl/blob/d3f30bb62fc282e9a3873a858f2823339888e143/src/reversecholesky.jl#L113-L117

which I think is more correct?

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions