Parallel 3D Adaptive Mesh Refinement in Titanium

We describe a 3-dimensional adaptive mesh reenement Poisson solver. The complete program consists of about 3,500 lines of Titanium code and runs on both shared-memory and distributed-memory architectures. This paper focuses on the algorithm and on our experiences in writing AMR and tuning its performance.