diff --git a/grudge/symbolic/operators.py b/grudge/symbolic/operators.py
index b95575f6ee1dd21c91fb26d43c8f158f20f9c427..4ecd920c8d82ee429855183b319fd191ed524a9f 100644
--- a/grudge/symbolic/operators.py
+++ b/grudge/symbolic/operators.py
@@ -721,7 +721,7 @@ def norm(p, arg, dd=None):
 
     if p == 2:
         norm_squared = NodalSum(dd_in=dd)(
-                prim.fabs(arg * MassOperator()(arg)))
+                arg * MassOperator()(arg))
 
         if isinstance(norm_squared, np.ndarray):
             if len(norm_squared.shape) != 1:
diff --git a/test/test_grudge.py b/test/test_grudge.py
index 60f58ddf4f4accb6a11d536e14057a25c5da41a6..b3662b5fb5c85fdd25d1d2bb0248aef56d159301 100644
--- a/test/test_grudge.py
+++ b/test/test_grudge.py
@@ -669,7 +669,7 @@ def test_norm_obj_array(ctx_factory, p):
 
     norm_exact = w[0]
     logger.info("norm: %.5e %.5e", norm, norm_exact)
-    # assert abs(norm - norm_exact) < 1.0e-14
+    assert abs(norm - norm_exact) < 1.0e-14
 
     # }}}
 
@@ -680,7 +680,7 @@ def test_norm_obj_array(ctx_factory, p):
 
     norm_exact = np.sqrt(np.sum(w**2)) if p == 2 else np.max(w)
     logger.info("norm: %.5e %.5e", norm, norm_exact)
-    # assert abs(norm - norm_exact) < 1.0e-14
+    assert abs(norm - norm_exact) < 1.0e-14
 
     # }}}