Cleanup

src/mueller.jl

+function müller_rec(::Linear; ϵ = 0.0, ϕ = 0.0, ΔG = 1.0, δα = 0.0, ψ = 0.0, s = +)
+    ± = s(1.0)
+    ∓ = -±
+    return [1             ∓*ΔG/2       2ϵ*cos(ϕ) 2ϵ*sin(ϕ)
+            ΔG/2          ±            0         ±*2δα
+            2ϵ*cos(ϕ + ψ) ±*2δα*sin(ψ) cos(ψ)    ∓*sin(ψ)
+            2ϵ*sin(ϕ + ψ) ∓*2δα*cos(ψ) sin(ψ)    ±*cos(ψ)]
+end
+
+function müller_rec(::Circular; ϵ = 0.0, ϕ = 0.0, ΔG = 1.0, δα = 0.0, ψ = 0.0, s = +)
+    ± = s(1.0)
+    ∓ = -±
+    return [1             ∓*2ϵ*sin(ϕ)  2ϵ*cos(ϕ)  ±*ΔG/2
+            ΔG/2          ∓*2δα            0         ±
+            2ϵ*cos(ϕ + ψ) ±*sin(ψ)      cos(ψ)   ±*2δα*sin(ψ)
+            2ϵ*sin(ϕ + ψ) ∓*cos(ψ)      sin(ψ)   ∓*2δα*cos(ψ)]
+end
+
+function müller_sky(PAaz)
+    cpa = cos(2*PAaz)
+    spa = sin(2*PAaz)
+    return [1    0    0   0
+            0   cpa  spa  0
+            0  -spa  cpa  0
+            0    0    0   1]
+end

src/test.jl

+using Plots
+
+pa = readdlm("pa")
+pa = reshape(pa', length(pa))
+i = readdlm("i")
+i = reshape(i', length(i))
+q = readdlm("q")
+q = reshape(q', length(q))
+u = readdlm("u")
+u = reshape(u', length(u))
+v = readdlm("v")
+v = reshape(v', length(v))
+
+plot(pa,i)