Skip to content

get determinant of operator array #62

Description

@vpuri3
using SciMLOperators, LinearAlgebra
D = DiagonalOperator(rand(4))
A = SciMLOperators.AbstractSciMLOperator[D D; D D]
det(A)
ERROR: MethodError: no method matching SciMLOperators.AbstractSciMLOperator(::ScalarOperator{Bool, typeof(SciMLOperators.DEFAULT_UPDATE_FUNC)})
Stacktrace:
 [1] oneunit(#unused#::Type{SciMLOperators.AbstractSciMLOperator})
   @ Base ./number.jl:358
 [2] lutype(T::Type)
   @ LinearAlgebra /Applications/Julia-1.8.app/Contents/Resources/julia/share/julia/stdlib/v1.8/LinearAlgebra/src/lu.jl:204
 [3] lu(A::Matrix{SciMLOperators.AbstractSciMLOperator}, pivot::RowMaximum; check::Bool)
   @ LinearAlgebra /Applications/Julia-1.8.app/Contents/Resources/julia/share/julia/stdlib/v1.8/LinearAlgebra/src/lu.jl:279
 [4] det(A::Matrix{SciMLOperators.AbstractSciMLOperator})
   @ LinearAlgebra /Applications/Julia-1.8.app/Contents/Resources/julia/share/julia/stdlib/v1.8/LinearAlgebra/src/generic.jl:1549
 [5] top-level scope
   @ REPL[1]:1

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Fields

    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions