def get_particles(with_boundary=True, **kwargs): adke, left, right = get_data.get_shock_tube_data( nl=nl, nr=nr, xl=xl, xr=xr, pl=pl, pr=pr, rhol=rhol, rhor=rhor, ul=ul, ur=ur, g1=g1, g2=g2, h0=h0, gamma=1.4 ) adke.m[:] = m0 return [adke]
def get_particles(with_boundary=False, **kwargs): adke, left, right = data.get_shock_tube_data( nl=nl, nr=nr, xl=xl, xr=xr, pl=pl, pr=pr, rhol=rhol, rhor=rhor, ul=ul, ur=ur, g1=g1, g2=g2, h0=h0, gamma=gamma ) if with_boundary: return [adke, left, right] else: return [adke]