multiplication

Author: jishnub

src/Spaces/ProductSpaceOperators.jl

@@ -220,18 +220,24 @@ end
 ## Derivative
 
 #TODO: do in @calculus_operator?
+_spacename(::SumSpace) = SumSpace
+_spacename(::PiecewiseSpace) = PiecewiseSpace
+
+function InterlaceOperator_Diagonal(t, S)
+    allbanded = all(isbanded, t)
+    ds, rs = S, _spacename(S)(map(rangespace, t))
+    D = Diagonal(convert_vector_or_svector(t))
+    bw = interlace_bandwidths(D, ds, rs, allbanded)
+    InterlaceOperator(D, ds, rs, bw)
+end
 
 for (Op,OpWrap) in ((:Derivative,:DerivativeWrapper),(:Integral,:IntegralWrapper))
     _Op = Symbol(:_, Op)
     @eval begin
         @inline function $_Op(S::PiecewiseSpace, k::Number)
             assert_integer(k)
             t = map(s->$Op(s,k),components(S))
-            allbanded = all(isbanded, t)
-            ds, rs = S, PiecewiseSpace(map(rangespace, t))
-            D = Diagonal(convert_vector_or_svector(t))
-            bw = interlace_bandwidths(D, ds, rs, allbanded)
-            O = InterlaceOperator(D, ds, rs, bw)
+            O = InterlaceOperator_Diagonal(t, S)
             $OpWrap(O,k)
         end
         @inline function $_Op(S::ArraySpace, k::Number)
@@ -257,11 +263,7 @@ end
     # mixed bases.
     if typeof(canonicaldomain(S))==typeof(domain(S))
         t = map(s->Derivative(s,k),components(S))
-        allbanded = all(isbanded, t)
-        D = Diagonal(convert_vector_or_svector(t))
-        ds, rs = S, SumSpace(map(rangespace, t))
-        bw = interlace_bandwidths(D, ds, rs, allbanded)
-        O = InterlaceOperator(D, ds, rs, bw)
+        O = InterlaceOperator_Diagonal(t, S)
         DerivativeWrapper(O,k)
     else
         DefaultDerivative(S,k)
@@ -297,8 +299,7 @@ function Multiplication(f::Fun{<:PiecewiseSpace}, sp::PiecewiseSpace)
     p=perm(domain(f).domains,domain(sp).domains)  # sort f
     vf=components(f)[p]
     t = map(Multiplication,vf,sp.spaces)
-    D = Diagonal(convert_vector_or_svector(t))
-    O = InterlaceOperator(D, PiecewiseSpace)
+    O = InterlaceOperator_Diagonal(t, sp)
     MultiplicationWrapper(f, O)
 end
 
@@ -307,8 +308,7 @@ Multiplication(f::Fun{SumSpace{SV1,D,R1}},sp::SumSpace{SV2,D,R2}) where {SV1,SV2
 
 function Multiplication(f::Fun, sp::SumSpace)
     t = map(s->Multiplication(f,s),components(sp))
-    D = Diagonal(convert_vector_or_svector(t))
-    O = InterlaceOperator(D, SumSpace)
+    O = InterlaceOperator_Diagonal(t, sp)
     MultiplicationWrapper(f, O)
 end