From 659c287bfd1c78ed6b17d4d53d58119cd437bf97 Mon Sep 17 00:00:00 2001 From: Eugeniy Mikhailov Date: Fri, 8 Jul 2011 12:18:05 -0400 Subject: fix incendent angle calculaton --- face_beam_interaction.m | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/face_beam_interaction.m b/face_beam_interaction.m index 33345a2..3bc3e4a 100644 --- a/face_beam_interaction.m +++ b/face_beam_interaction.m @@ -6,6 +6,7 @@ function [is_face_hit, hit_position, hit_distance, new_beams] = face_beam_inter % beam.intensity - intensity of the beam % beam.face - if beam starts from face then its index is here % beam.polarization - can be either 1 for 's' or 2 for 'p' + % beam.status - could be 'reflected', 'refracted', 'incoming' % faces cell array of face structures % face - structure definiong the beam % face.vertex1 - [x,y] of the 1st point/vertex of the face @@ -97,7 +98,7 @@ function [is_face_hit, hit_position, hit_distance, new_beams] = face_beam_inter theta_i=pi-theta_i; end if( theta_i < -pi/2) - theta_i=-pi+theta_i; + theta_i=-pi-theta_i; end % angle of the normal with respect to horizon @@ -120,6 +121,7 @@ function [is_face_hit, hit_position, hit_distance, new_beams] = face_beam_inter beam_reflected.face=closest_face_index; beam_reflected.polarization = polarization; beam_reflected.intensity = beam.intensity * reflectivity; + beam_reflected.status = 'reflected'; new_beams{1} = beam_reflected; @@ -140,6 +142,7 @@ function [is_face_hit, hit_position, hit_distance, new_beams] = face_beam_inter beam_refracted.face=closest_face_index; beam_refracted.polarization = polarization; beam_refracted.intensity = beam.intensity * transmission; + beam_refracted.status = 'refracted'; new_beams{2} = beam_refracted; end end -- cgit v1.2.3