% Test of multiple density calculation methods
load("raytrace_longrun_data_N=1e5a.octavebinary")

% jit_enable(1,"local");

figure(1)
tic
foo=density(xs,ys,L,0.01);
toc



% Did density2 improve upon the efficiency of density.m?
figure(2)
tic
foo2=density2(xs,ys,L,0.01);
toc
% Conclusion: Miserable failure. density.m is actually faster!
