import yt
import time

yt.enable_parallelism()

# start timer
t1 = time.time()

# load data
ds = yt.load("Data_000001")

# heavy works
p1 = yt.ProjectionPlot(ds, "x", "density")
p2 = yt.ProjectionPlot(ds, "y", "density")
p3 = yt.ProjectionPlot(ds, "z", "density")
p1 = yt.ProjectionPlot(ds, "x", "density")
p2 = yt.ProjectionPlot(ds, "y", "density")
p3 = yt.ProjectionPlot(ds, "z", "density")

v, c = ds.find_max("density")
v, c = ds.find_max("density")
v, c = ds.find_max("density")

# end timer
t2 = time.time()

# show elapsed time
if yt.is_root():
    print("yt completed  %f sec" % (t2-t1))