Bug fix in the construction of the compressible divergence matrix: the multiplication by density_bdy in the surface integral was accidentally deleted. We now have detwei_bdy*density_bdy*face_val_at_quad(nvfrac, sele). This doesn't affect any existing compressible multiphase test cases as we are not integrating by parts, but this might explain why the shock tube test was not working properly with P0-P1.
Bug fix in the construction of the compressible divergence matrix: the multiplication by density_bdy in the surface integral was accidentally deleted. We now have detwei_bdy*density_bdy*face_val_at_quad(nvfrac, sele). This doesn't affect any existing compressible multiphase test cases as we are not integrating by parts, but this might explain why the shock tube test was not working properly with P0-P1.

jrper/fluidity,jrper/fluidity,iakovos-panourgias/fluidity,iakovos-panourgias/fluidity,iakovos-panourgias/fluidity,rjferrier/fluidity,jrper/fluidity,rjferrier/fluidity,jjo31/ATHAM-Fluidity,rjferrier/fluidity,jrper/fluidity,jjo31/ATHAM-Fluidity,jjo31/ATHAM-Fluidity,rjferrier/fluidity,jjo31/ATHAM-Fluidity,jjo31/ATHAM-Fluidity,jrper/fluidity,rjferrier/fluidity,iakovos-panourgias/fluidity,iakovos-panourgias/fluidity