diff --git a/examples/dagrt-fusion.py b/examples/dagrt-fusion.py
index d8b6dcc7d8fe118f5a017b7512b6c70eeccdbb92..9ab70e2c535c64e62bdc17add11353b872ef0b2d 100755
--- a/examples/dagrt-fusion.py
+++ b/examples/dagrt-fusion.py
@@ -1116,21 +1116,23 @@ else:
 
 def problem_stats(order=3):
     cl_ctx = cl.create_some_context()
+    queue = cl.CommandQueue(cl_ctx)
+    actx = PyOpenCLArrayContext(queue)
 
     with open_output_file("grudge-problem-stats.txt") as outf:
         _, dg_discr_2d = get_strong_wave_op_with_discr_direct(
-            cl_ctx, dims=2, order=order)
+            actx, dims=2, order=order)
         print("Number of 2D elements:", dg_discr_2d.mesh.nelements, file=outf)
         vol_discr_2d = dg_discr_2d.discr_from_dd("vol")
-        dofs_2d = {group.nunit_nodes for group in vol_discr_2d.groups}
+        dofs_2d = {group.nunit_dofs for group in vol_discr_2d.groups}
         from pytools import one
         print("Number of DOFs per 2D element:", one(dofs_2d), file=outf)
 
         _, dg_discr_3d = get_strong_wave_op_with_discr_direct(
-            cl_ctx, dims=3, order=order)
+            actx, dims=3, order=order)
         print("Number of 3D elements:", dg_discr_3d.mesh.nelements, file=outf)
         vol_discr_3d = dg_discr_3d.discr_from_dd("vol")
-        dofs_3d = {group.nunit_nodes for group in vol_discr_3d.groups}
+        dofs_3d = {group.nunit_dofs for group in vol_discr_3d.groups}
         from pytools import one
         print("Number of DOFs per 3D element:", one(dofs_3d), file=outf)