% Restore a big raytrace (e.g. from raytrace_longrun) and plot the density. Resave the file with
% this new density information. This version uses density_ss to minimize artifacts in the density.
function density_ss_longrun(filename)

t0 = cputime;

filename_input = filename
load(filename);
filename = filename_input
t1 = cputime;
cpu_time_to_load_file = t1 - t0

ray_density_ss = density_ss(xs,ys,L,0.01);

t2 = cputime;
cpu_time_to_calculate_density = t2 - t1

clear t0 filename_input %Don't want these items saved.
save("-binary",filename);
print(strcat(filename,'.pdf'),'-dpdfwrite')

t3 = cputime;
cpu_time_to_save = t3 - t2