From e405091ddf823c49183f5a952f6a5b8e1e9ef679 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Tue, 22 Nov 2011 15:15:22 -0500 Subject: fix error about output fields with xi --- faraday_and_psr/pp_tasks.m | 6 +++--- squeezing_filter/pp_tasks.m | 5 ++--- 2 files changed, 5 insertions(+), 6 deletions(-) diff --git a/faraday_and_psr/pp_tasks.m b/faraday_and_psr/pp_tasks.m index a3919ba..0221897 100644 --- a/faraday_and_psr/pp_tasks.m +++ b/faraday_and_psr/pp_tasks.m @@ -20,9 +20,9 @@ for i=1:Nsteps Ep.right(i)=d.E_field_pos_freq.right; endfor -Ep_out.linear=(1-xi_linear).*Ep.linear; -Ep_out.left=(1-xi_left).*Ep.left; -Ep_out.right=(1-xi_right).*Ep.right; +Ep_out.left = exp(sqrt(1+1i*xi_left)) .*Ep.left; +Ep_out.right = exp(sqrt(1+1i*xi_right)) .*Ep.right; + Ep_out.x= (Ep_out.left + Ep_out.right)/sqrt(2); Ep_out.y= 1i*(Ep_out.left - Ep_out.right)/sqrt(2); diff --git a/squeezing_filter/pp_tasks.m b/squeezing_filter/pp_tasks.m index 6e5c692..c99b950 100644 --- a/squeezing_filter/pp_tasks.m +++ b/squeezing_filter/pp_tasks.m @@ -24,9 +24,8 @@ endfor xi_x=(xi_right+xi_left)/sqrt(2); xi_y=(1i*xi_right-1i*xi_left)/sqrt(2); -%Ep_out.linear=(1-xi_linear).*Ep.linear; -Ep_out.left=(1-xi_left).*Ep.left; -Ep_out.right=(1-xi_right).*Ep.right; +Ep_out.left = exp(sqrt(1+1i*xi_left)) .*Ep.left; +Ep_out.right = exp(sqrt(1+1i*xi_right)) .*Ep.right; Ep_out.x= (Ep_out.left + Ep_out.right)/sqrt(2); Ep_out.y= 1i*(Ep_out.left - Ep_out.right)/sqrt(2); -- cgit v1.2.3