From ad7df20050f6e9d2e3f212d7043f32ea142c8606 Mon Sep 17 00:00:00 2001
From: Datseris <datseris.george@gmail.com>
Date: Sun, 19 Apr 2020 12:50:36 +0200
Subject: [PATCH] correparticlebeam

---
 src/billiards/particles.jl | 10 ++++++----
 1 file changed, 6 insertions(+), 4 deletions(-)

diff --git a/src/billiards/particles.jl b/src/billiards/particles.jl
index fb428622..832afb4f 100644
--- a/src/billiards/particles.jl
+++ b/src/billiards/particles.jl
@@ -141,15 +141,17 @@ end
 ## Aux
 ####################################################
 """
-    particlebeam(x0, y0, φ, N, dx, ω = nothing) → ps
+    particlebeam(x0, y0, φ, N, dx, ω = nothing, T = eltype(x0)) → ps
 Make `N` particles, all with direction `φ`, starting at `x0, y0`. The particles
 don't all have the same position, but are instead spread by up to `dx` in the
 direction normal to `φ`.
+
+The particle element type is `T`.
 """
-function particlebeam(x0, y0, φ, N, dx, ω = nothing, T = Float64)
-    n = sincos(φ)
+function particlebeam(x0, y0, φ, N, dx, ω = nothing, T = eltype(x0))
+    n = cossin(φ)
     xyφs = [
-    T.((x0 + i*dx*n[1]/N, y0 + i*dx*n[2]/N, φ)) for i in range(-N/2, N/2; length = N)
+    T.((x0 - i*dx*n[2]/N, y0 + i*dx*n[1]/N, φ)) for i in range(-N/2, N/2; length = N)
     ]
     if isnothing(ω)
         ps = [Particle(z...) for z in xyφs]