Source code for pyro.advection_fv4.problems.smooth

"""Initialize a Gaussian profile (shifted so the minimum value is 1.0).  This
is smooth and the limiters should not kick in too much, so this can be used
for testing convergence."""

import numpy

from pyro.util import msg

DEFAULT_INPUTS = "inputs.smooth"

PROBLEM_PARAMS = {}


[docs] def init_data(my_data, rp): """ initialize the smooth advection problem """ if rp.get_param("driver.verbose"): msg.bold("initializing the smooth advection problem...") dens = my_data.get_var("density") xmin = my_data.grid.xmin xmax = my_data.grid.xmax ymin = my_data.grid.ymin ymax = my_data.grid.ymax xctr = 0.5*(xmin + xmax) yctr = 0.5*(ymin + ymax) dens[:, :] = 1.0 + numpy.exp(-60.0*((my_data.grid.x2d-xctr)**2 + (my_data.grid.y2d-yctr)**2))
[docs] def finalize(): """ print out any information to the user at the end of the run """