Sharpclaw3d

This pull request is linked to current pull request on PyClaw. We have implemented a 3D SharpClaw solver. However in the process we discovered inconsistencies in the sequence calls for rpnx, x=1,2,3. We have standardized them to follow the following order:

ixyz,maxnx,num_eqn,num_waves,num_aux,num_ghost,mx,ql,qr,auxl,auxr,wave,s,amdq,apdq

The suggested changes plus the on-going pull request in PyClaw yield the following results:

test_1d_acoustics ... ok
test_1d_acoustics ... ok
test_1d_acoustics ... ok
test_1d_acoustics ... ok
test_1d_acoustics ... ok
test_1d_acoustics ... ok
test_1d_acoustics ... ok
test_1d_acoustics ... ok
test_1d_acoustics ... ok
test_1d_acoustics ... ok
test_2d_acoustics ... ok
test_2d_acoustics ... ok
test_2d_acoustics ... ok
test_2d_acoustics ... ok
Test variable-coefficient 2D acoustics ... ok
Test variable-coefficient 2D acoustics ... ok
Test variable-coefficient 2D acoustics ... ok
Test variable-coefficient 2D acoustics ... ok
Test I/O on variable-coefficient 2D acoustics application ... ok
Test I/O on variable-coefficient 2D acoustics application ... ok
test_3d_acoustics ... ok
test_3d_acoustics ... ok
test_3d_acoustics ... ok
test_3d_acoustics ... ok