Starting to set up 3d poisson kernel.

I have not tested it yet.
Added 3d indexing to jacobi_3d kernel
Corrected major stupid error in 2D jacobi that slowed convergence down to a crawl.

Increased the problem size in run_jacobi_2d. Before the fix it did not finish in over a minute. Now it finishes in just a second.