import math


def two_strong_shocks(left, right, g):

	"""
	Pressure at the interface of the Riemann problem in the case of two strong shocks
	Input:
	left - Primitive variables on the left side
	right - Primitive variables on the right side
	g - Adiabatic index
	"""
	
	
	dl = left.Density
	vl = left.Velocity
	dr = right.Density
	vr = right.Velocity
	return (dl*dr*(1 + g)*(vl - vr)**2)/(2.*(math.sqrt(dl) + math.sqrt(dr))**2)