update for scimlops

src/composite_operators.jl

@@ -14,7 +14,7 @@ sparse1(A::DiffEqScaledOperator) = A.coeff * sparse1(A.op)
 # Linear Combination
 struct DiffEqOperatorCombination{T, O <: Tuple{Vararg{AbstractDiffEqLinearOperator{T}}},
                                  C <: AbstractVector{T}} <:
-       AbstractDiffEqCompositeOperator{T}
+       SciMLBase.AbstractDiffEqCompositeOperator{T}
     ops::O
     cache::C
     function DiffEqOperatorCombination(ops; cache = nothing)
@@ -77,7 +77,7 @@ end
 # Composition (A ∘ B)
 struct DiffEqOperatorComposition{T, O <: Tuple{Vararg{AbstractDiffEqLinearOperator{T}}},
                                  C <: Tuple{Vararg{AbstractVector{T}}}} <:
-       AbstractDiffEqCompositeOperator{T}
+       SciMLBase.AbstractDiffEqCompositeOperator{T}
     ops::O # stored in the order of application
     caches::C
     function DiffEqOperatorComposition(ops; caches = nothing)

src/derivative_operators/derivative_operator_functions.jl

@@ -184,7 +184,7 @@ end
 
 # A more efficient mul! implementation for compositions of operators which may include regular-grid, centered difference,
 # scalar coefficient, non-winding, DerivativeOperator, operating on a 2-D or 3-D AbstractArray
-function LinearAlgebra.mul!(x_temp::AbstractArray{T, 2}, A::AbstractDiffEqCompositeOperator,
+function LinearAlgebra.mul!(x_temp::AbstractArray{T, 2}, A::SciMLBase.AbstractDiffEqCompositeOperator,
                             M::AbstractArray{T, 2}) where {T}
 
     # opsA operators satisfy conditions for NNlib.conv! call, opsB operators do not
@@ -447,7 +447,7 @@ end
 
 # A more efficient mul! implementation for compositions of operators which may include regular-grid, centered difference,
 # scalar coefficient, non-winding, DerivativeOperator, operating on a 2-D or 3-D AbstractArray
-function LinearAlgebra.mul!(x_temp::AbstractArray{T, 3}, A::AbstractDiffEqCompositeOperator,
+function LinearAlgebra.mul!(x_temp::AbstractArray{T, 3}, A::SciMLBase.AbstractDiffEqCompositeOperator,
                             M::AbstractArray{T, 3}) where {T}
 
     # opsA operators satisfy conditions for NNlib.conv! call, opsB operators do not

src/derivative_operators/ghost_derivative_operator.jl

@@ -8,7 +8,7 @@ function *(L::AbstractDiffEqLinearOperator{T}, Q::AbstractBC{T}) where {T}
     return GhostDerivativeOperator{T, typeof(L), typeof(Q)}(L, Q)
 end
 
-function *(L::AbstractDiffEqCompositeOperator{T}, Q::AbstractBC{T}) where {T}
+function *(L::SciMLBase.AbstractDiffEqCompositeOperator{T}, Q::AbstractBC{T}) where {T}
     return sum(map(op -> op * Q, L.ops))
 end
 

src/docstrings.jl

@@ -24,7 +24,7 @@ calculate_weights
 
 Represent a finite-difference derivative operator.
 
-These operators implement the `DiffEqBase.AbstractDiffEqLinearOperator` interface. Therefore, `eltype` returns `T`.
+These operators implement the `SciMLBase.AbstractDiffEqLinearOperator` interface. Therefore, `eltype` returns `T`.
 
 These operators can be contracted over an arbitrary dimension, given by the type parameter `N`.
 
@@ -34,7 +34,7 @@ The finite-difference methods are defined for `DerivativeOperator`. In particula
 
 The key data are three stencils of coefficients, stored in `SVector`s. The interior stencil, stored in `stencil_coefs::S1`, is the normal one used in the interior of the grid. The others, `low_boundary_coefs::S2` and `high_boundary_coefs::S2`, are used where the normal stencil would jut out of the grid boundary. These can have a different length than the interior stencil, hence the two types.
 
-When the operator is applied by `mul!`, these stencils are multiplied by the vector in the `coefficients` field. It is set to `coeff_func.(1:len)` in `UpwindDifference`, which seems odd, because `len` is not the stencil length. Scalar multiplication is absorbed into this vector, as required by the `DiffEqBase.AbstractDiffEqLinearOperator` interface.
+When the operator is applied by `mul!`, these stencils are multiplied by the vector in the `coefficients` field. It is set to `coeff_func.(1:len)` in `UpwindDifference`, which seems odd, because `len` is not the stencil length. Scalar multiplication is absorbed into this vector, as required by the `SciMLBase.AbstractDiffEqLinearOperator` interface.
 
 The `coefficients` field appears to be a more general `DifferentialEquations` thing. There is an `update_coefficients!` method.
 

src/jacvec_operators.jl

@@ -1,6 +1,6 @@
 ### Operator Implementation
 mutable struct JacVecOperator{T, F, T1, T2, uType, P, tType, O} <:
-               DiffEqBase.AbstractDiffEqLinearOperator{T}
+               SciMLBase.AbstractDiffEqLinearOperator{T}
     f::F
     cache1::T1
     cache2::T2
@@ -151,7 +151,7 @@ end
 ### AnalyticalOperator Implementation
 
 mutable struct AnalyticalJacVecOperator{T, F, uType, P, tType, O} <:
-               DiffEqBase.AbstractDiffEqLinearOperator{T}
+               SciMLBase.AbstractDiffEqLinearOperator{T}
     f::F
     u::uType
     p::P

src/vecjac_operators.jl

@@ -1,5 +1,5 @@
 mutable struct VecJacOperator{T, F, T1, T2, uType, P, tType, O} <:
-               DiffEqBase.AbstractDiffEqLinearOperator{T}
+               SciMLBase.AbstractDiffEqLinearOperator{T}
     f::F
     cache1::T1
     cache2::T2