Currently the defaults in pyclaw are solver.dimensional_split = False and solver.transverse_waves=1/2 in 2D/3D. These don't match the defaults in Clawpack, which are to use the unsplit algorithm with transverse_waves=2. It might be worthwhile to make these consistent.