SphericalScattering.jl
SphericalScattering.jl copied to clipboard
Jump in electric field while testing against BEAST
While following the application example from your documentation I observed a jump in the electric field. If my points_cartNF is defined as
ϑ = range(0.0, stop=0.999999999999999π, length=18) # 10° steps ϕ = range(0.0, stop=0.9999999999999992π, length=36) # 10° steps P = [point(cos(φ) * sin(θ), sin(φ) * sin(θ), cos(θ)) for θ in ϑ, φ in ϕ] points_cartNF = P .* r
the test against BEAST results in 1.034% relative error.
With one additional 9 (moving towards π and 2π) I experienced a jump in the electric field at critical points.
The relative error increases over 100%.