function [] = test_poisson3d()
% Quick and dirty test of the poisson_3d function

m = 58;
n = 58;
q = 58;

% Random right hand side
rho = rand(m,n,q);

delta = 1;
disp(['poisson_3d: Domain Size (m,n,q): ' num2str(m) ' x ' num2str(n) ' x ' num2str(q)])
tic
poisson_3d(rho,delta);
disp(['Elapsed time using poisson_3d: ' num2str(toc) 's'])