Define combining block sparse interfaces

Author: mtfishman

src/abstractblocksparsearray/wrappedabstractblocksparsearray.jl

@@ -29,8 +29,8 @@ const AnyAbstractBlockSparseVecOrMat{T,N} = Union{
   AnyAbstractBlockSparseVector{T},AnyAbstractBlockSparseMatrix{T}
 }
 
-function DerivableInterfaces.interface(::Type{<:AnyAbstractBlockSparseArray})
-  return BlockSparseArrayInterface()
+function DerivableInterfaces.interface(arrayt::Type{<:AnyAbstractBlockSparseArray})
+  return BlockSparseArrayInterface(interface(blocktype(arrayt)))
 end
 
 # a[1:2, 1:2]

src/blocksparsearrayinterface/blocksparsearrayinterface.jl

@@ -17,7 +17,12 @@ using BlockArrays:
   blocks,
   findblockindex
 using DerivableInterfaces:
-  DerivableInterfaces, @interface, AbstractArrayInterface, DefaultArrayInterface, zero!
+  DerivableInterfaces,
+  @interface,
+  AbstractArrayInterface,
+  DefaultArrayInterface,
+  interface,
+  zero!
 using LinearAlgebra: Adjoint, Transpose
 using SparseArraysBase:
   AbstractSparseArrayInterface,
@@ -125,7 +130,8 @@ function BlockSparseArrayInterface{N}(blockinterface::AbstractArrayInterface{N})
   return BlockSparseArrayInterface{N,typeof(blockinterface)}(blockinterface)
 end
 function BlockSparseArrayInterface{M,B}(::Val{N}) where {M,B<:AbstractArrayInterface{M},N}
-  return BlockSparseArrayInterface{N}(B(Val(N)))
+  B′ = B(Val(N))
+  return BlockSparseArrayInterface(B′)
 end
 function BlockSparseArrayInterface{N}() where {N}
   return BlockSparseArrayInterface{N}(DefaultArrayInterface{N}())
@@ -134,6 +140,14 @@ BlockSparseArrayInterface(::Val{N}) where {N} = BlockSparseArrayInterface{N}()
 BlockSparseArrayInterface{M}(::Val{N}) where {M,N} = BlockSparseArrayInterface{N}()
 BlockSparseArrayInterface() = BlockSparseArrayInterface{Any}()
 
+function DerivableInterfaces.combine_interface_rule(
+  interface1::AbstractBlockSparseArrayInterface,
+  interface2::AbstractBlockSparseArrayInterface,
+)
+  B = interface(blockinterface(interface1), blockinterface(interface2))
+  return BlockSparseArrayInterface(B)
+end
+
 @interface ::AbstractBlockSparseArrayInterface function BlockArrays.blocks(a::AbstractArray)
   return error("Not implemented")
 end