lambda = 0.124;

tilesize = 0.25;

beta = 2*pi/lambda;

c = 299792458;

t_closest = inf;

t_point = [];

t_found = 0;

t_cuts = 0;

if a.object_number > 0

for nobj = 1:a.object_number

obj1 = a.object_array(nobj);

if strcmp(obj1.rregion,current_ray.rregion)||strcmp(obj1.lregion,current_ray.rregion)||strcmp(current_ray.rregion,’Initial’Log in to see images!

f = 0;

[f,i_t_values,i_points,i_facets,i_angle] = find_intersection(current_ray,obj1);

if (f == 1)&&(current_ray.ignore ~= nobj)

if t_found == 0

t_closest = i_t_values(1);

t_point = i_points(1,Log in to see images!;

t_face = i_facets(1);

t_face2 = i_facets(2);

t_found = 1;

t_cuts = nobj;

elseif i_t_values(1) < t_closest

t_closest = i_t_values(1);

t_point = i_points(1,Log in to see images!;

t_face = i_facets(1);

t_face2 = i_facets(2);

t_found = t_found + 1;

t_cuts = nobj;

end

end

end

end

end

if t_found > 0

current_ray.decay = t_point;

current_ray.end_length = current_ray.starting_length + t_closest;

c_rel_power = current_ray.acgreat timesulated_starting_losses.*(lambda./(4.*pi.*current_ray.end_length.*tilesize)).^2;

if c_rel_power >= min_relative_power

chn = current_ray.children_number + 1;

current_ray.children_number = chn;

t_child = trans_child(current_ray,t_cuts,a.object_array(t_cuts),t_closest);

f_child(chn) = follow_ray(a, t_child,min_relative_power);

chn = chn + 1;

current_ray.children_number = chn;

r_child = refl_child(current_ray,t_cuts,a.object_array(t_cuts),t_face,t_closest);

f_child(chn) = follow_ray(a, r_child,min_relative_power);

current_ray.children_tree = f_child;

end

else

bounds = [];

[t_ray, err_code] = ray_out_of_bounds(current_ray, bounds);

current_ray = t_ray;

end

exit_ray = current_ray;