diff options
| -rw-r--r-- | faraday_and_psr/pp_tasks.m | 6 | ||||
| -rw-r--r-- | 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); |