A New Parallel Method for Binary Black Hole Simulations

Algorithm 1

Submesh based Parallel Method.

(1) procedure SPM_Evolution

(2)   Input: meshes and boundary conditions

(3)   evolve one step on current level

(4)   if (total_level > 1) call recursive_step(1)

(5)   analyze the results

(6)   regrid if needed

(7)   Output: values on all the mesh grid points

(8) end procedure

(9) procedure  recursive_step()

(10)   for  ; ;   do

(11)     evolve one step on current level

(12)     if (level < total_level − 1) then

(13)      call recursive_step(level + 1)

(14)     end if

(15)     exchange data with related processes

(16)   end for

(17) end procedure