search

CalculustJL / CalculustCore.jl

Ecosystem for writing partial differential equation solvers
MIT License
14 stars 0 forks source link

need a notion of `size(::AbstractOperator)` to use `IterativeSolvers.jl` #5

Closed vpuri3 closed 2 years ago

vpuri3 commented 2 years ago

do we need to do size checking for mul!, *?

vpuri3 commented 2 years ago

^ yes

vpuri3 commented 2 years ago

instead of defining a separate IdentityOp, NullOp, just define a ScaleOp and set

const IdentityOp{D,N} = ScaleOp{D,N}(true)
const NullOp{D,N} = ScaleOp{D,N}(false)