Skip to content

About issymmetric(REAL_SYM) #112

@briochemc

Description

@briochemc

Could marking real structurally symmetric matrices (REAL_SYM) as symmetric cause any issues? Specifically this line:

LinearAlgebra.issymmetric(v::MatrixType) = v in (REAL_SYM, REAL_SYM_POSDEF, REAL_SYM_INDEF, COMPLEX_STRUCT_SYM,

where I understand REAL_SYM (type 1 in Intel docs) as "the sparsity pattern is symmetric", and thus a REAL_SYM matrix may not be symmetric.

I'm worried Pardiso.jl might ccall with the wrong parameters and thus run an incorrect branch. It's very likely I'm just confused and this is irrelevant, so please don't hesitate to close this issue if that's the case! 😅

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