import sys
from pyro.mesh import patch
[docs]
def init_data(my_data, rp):
""" an init routine for unit testing """
del rp # this problem doesn't use runtime params
# make sure that we are passed a valid patch object
if not isinstance(my_data, patch.CellCenterData2d):
print("ERROR: patch invalid in sedov.py")
print(my_data.__class__)
sys.exit()
# get the density, momenta, and energy as separate variables
dens = my_data.get_var("density")
# initialize the components, remember, that ener here is rho*eint
# + 0.5*rho*v**2, where eint is the specific internal energy
# (erg/g)
dens[:, :] = 1.0
[docs]
def finalize():
""" print out any information to the user at the end of the run """