diff --git a/pytential/qbx/__init__.py b/pytential/qbx/__init__.py
index f3d4007e2268b7ecd352ad5d4d8d4bae5231932d..2ab1488a02198022c39bec95566cd9d77f6b90f2 100644
--- a/pytential/qbx/__init__.py
+++ b/pytential/qbx/__init__.py
@@ -379,25 +379,6 @@ class QBXLayerPotentialSource(LayerPotentialSourceBase):
 
     # }}}
 
-    @memoize_method
-    def weights_and_area_elements(self):
-        import pytential.symbolic.primitives as p
-        from pytential.symbolic.execution import bind
-        with cl.CommandQueue(self.cl_context) as queue:
-            # fine_density_discr is not guaranteed to be usable for
-            # interpolation/differentiation. Use density_discr to find
-            # area element instead, then upsample that.
-
-            area_element = self.resampler(queue,
-                    bind(
-                        self.density_discr,
-                        p.area_element(self.ambient_dim, self.dim)
-                        )(queue))
-
-            qweight = bind(self.fine_density_discr, p.QWeight())(queue)
-
-            return (area_element.with_queue(queue)*qweight).with_queue(None)
-
     # {{{ helpers for symbolic operator processing
 
     def preprocess_optemplate(self, name, discretizations, expr):