Fix multivariate Dirichlet (#90)

Author: jishnub

src/FourierOperators.jl

@@ -365,7 +365,7 @@ Reverse(S::Fourier{D}) where {D} = ReverseWrapper(NegateEven(S,S))
 ## Multivariate
 
 
-for TYP in (:Fourier,:Laurent,:CosSpace,:SinSpace,:Taylor)
+for TYP in (:Fourier,)
     @eval begin
         function Dirichlet(S::TensorSpace{Tuple{$TYP{PeriodicSegment{T},R},PS}},k) where {PS,T,R}
             op = [Operator(I,S.spaces[1])⊗Evaluation(S.spaces[2],leftendpoint,k);
@@ -374,7 +374,7 @@ for TYP in (:Fourier,:Laurent,:CosSpace,:SinSpace,:Taylor)
         end
         function Dirichlet(S::TensorSpace{Tuple{PS,$TYP{PeriodicSegment{T},R}}},k) where {PS,T,R}
             op = [Evaluation(S.spaces[1],leftendpoint,k)⊗Operator(I,S.spaces[2]);
-                            Evaluation(S.spaces[1],rightendpoint,k)⊗ReverseOrientation(S[2]) ]
+                            Evaluation(S.spaces[1],rightendpoint,k)⊗ReverseOrientation(S.spaces[2]) ]
             DirichletWrapper(SpaceOperator(op,S,PiecewiseSpace(rangespace(op).spaces)),k)
         end
     end

test/runtests.jl

@@ -752,6 +752,18 @@ end
     @test f(1.1im) ≈ 2exp(-1.1^2)
 end
 
+@testset "multivariate Dirichlet" begin
+    f = Fun((x,y) -> cos(x)*sin(y), CosSpace()*Fourier())
+    g = Dirichlet() * f
+    @test f(0, pi/2) ≈ g(0, pi/2) atol=1e-10
+    @test f(2pi, pi/2) ≈ g(2pi, pi/2) atol=1e-10
+
+    f = Fun((x,y) -> cos(x)*cos(y), Fourier()*CosSpace())
+    g = Dirichlet() * f
+    @test f(pi/2, 0) ≈ g(pi/2, 0) atol=1e-10
+    @test f(pi/2, 2pi) ≈ g(pi/2, 2pi) atol=1e-10
+end
+
 @testset "show" begin
     for (spT, spstr) in Any[(Fourier, "Fourier"), (CosSpace, "CosSpace"),
             (SinSpace, "SinSpace"), (Laurent, "Laurent"), (Taylor, "Taylor")]