import yt
import time

yt.enable_parallelism()

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

# start timer
t1 = time.time()

# heavy works
p1 = yt.ProjectionPlot(ds, "x", "density")

# end timer
t2 = time.time()

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