Results from the Brio-Wu shocktube in 3D, with an initial left state (ρ_L, P_L, v_x, v_y, v_z, B_x, B_y, B_z) = (1, 1, 0, 0, 0, 0.75, 1, 0) and right state (ρ_R, P_R, v_x, v_y, v_z, B_x, B_y, B_z) = (0.125, 0.1, 0, 0, 0, 0.75, −1, 0). The figure shows the active region of the shock after t = 0.2, which contains about n_x ≈ 300−400 particles across the x-direction. The blue line shows the reference solution and the black dots show the result from the GDSPH simulation, while red dots show the result from the TSPH simulation. There are minimal differences between the GDSPH and TSPH result.