Updated polynomial solution for Yukawa

Author: SonglanNie

docs/src/examples/cavity2d_scattering.jl

@@ -56,7 +56,7 @@ println("Number of quadrature points: ", length(Q))
 ## Setup the integral operators
 op = Inti.Helmholtz(; dim = 2, k)
 S, D = Inti.single_double_layer(;
-    op
+    op,
     target = Q,
     source = Q,
     correction = (method = :none,),

src/vdim.jl

@@ -357,8 +357,14 @@ function polynomial_solution(op::Helmholtz, p::ElementaryPDESolutions.Polynomial
     return ElementaryPDESolutions.convert_coefs(P, Float64)
 end
 
+function polynomial_solution(op::Yukawa, p::ElementaryPDESolutions.Polynomial)
+    k = im*op.λ
+    P = ElementaryPDESolutions.solve_helmholtz(p; k)
+    return ElementaryPDESolutions.convert_coefs(P, Float64)
+end
+
 function neumann_trace(
-    ::Union{Laplace,Helmholtz},
+    ::Union{Laplace,Helmholtz, Yukawa},
     P::ElementaryPDESolutions.Polynomial{N,T},
 ) where {N,T}
     return _normal_derivative(P)