p = h3d_params;

s = h3d_create_struct(p.angular_cell_size_deg, p.rho_cell_size, -p.max_norm, p.max_norm);


for face=1:6
level = 0.1 + face/8;

for i=1:s.cube_ncells
for j=1:s.cube_ncells
	
	if mod(i+j,2) == 0
		color = 1;
	else
		color = level;
	end
	
	s.cube_hs(face,i,j) = color;
end
end


end


figure
h3d_draw_sphere_hs(s, 300);
axis('equal');

print('-depsc', 'sphere.eps');




 p.num_alpha_steps = 200;
 p.num_beta_steps = 1000;
 weight = 0;
 amplitude = pi;

axis = [0;0;1];
support = [1;0;0];

cyl1 = h3d_project_hs_on_cylinder(s, axis,support, p.num_alpha_steps, p.num_beta_steps, amplitude,weight, 0);


axis = normalize([0;0.7;0.3]);
support = [1;0;0];

cyl2 = h3d_project_hs_on_cylinder(s, axis,support, p.num_alpha_steps, p.num_beta_steps, amplitude,weight, 0);

cyls = {cyl1, cyl2};

for i=1:numel(cyls)
	figure
	imagesc(cyls{i});
	set(gca,'XTickLabel',[])
	set(gca,'XTick',[])
	set(gca,'YTickLabel',[])
	set(gca,'YTick',[]);
	print('-depsc', sprintf('cyl%d.eps', i));
end