squash some more misc deprecations

Author: alexfikl

pytential/linalg/skeletonization.py

@@ -132,7 +132,9 @@ def prg():
             <> ioffset = starts[icluster]
             <> npoints = starts[icluster + 1] - ioffset
             result[icluster] = reduce(sum, i, waa[indices[i + ioffset]]) / npoints
-            """)
+            """,
+            lang_version=lp.MOST_RECENT_LANGUAGE_VERSION,
+            )
 
         return knl
 

pytential/symbolic/execution.py

@@ -193,11 +193,8 @@ def map_elementwise_max(self, expr):
     def map_ones(self, expr):
         discr = self.places.get_discretization(
                 expr.dofdesc.geometry, expr.dofdesc.discr_stage)
-        result = discr.empty(actx=self.array_context, dtype=discr.real_dtype)
-
-        for grp_ary in result:
-            grp_ary.fill(1)
-        return result
+        return self.array_context.np.ones_like(
+            self.array_context.thaw(discr.nodes()[0]))
 
     def map_node_coordinate_component(self, expr):
         discr = self.places.get_discretization(
@@ -468,7 +465,7 @@ def flatten(self, ary):
             ary = [ary]
 
         from arraycontext import flatten
-        result = self.array_context.empty(self.total_dofs, self.dtype)
+        result = self.array_context.zeros(self.total_dofs, self.dtype)
         for res_i, (start, end) in zip(ary, self.starts_and_ends):
             result[start:end] = flatten(res_i, self.array_context)
 

pytential/symbolic/pde/scalar.py

@@ -411,7 +411,8 @@ def Sp(density):
             else:
                 raise ValueError(f"no improved operator for '{self.kernel}' known")
         else:
-            DpS0u = sym.Dp(self.kernel, laplace_s_inv_sqrt_w_u, **kwargs)
+            DpS0u = sym.Dp(self.kernel, laplace_s_inv_sqrt_w_u,
+                           qbx_forced_limit=+1, **kwargs)
 
         if self.is_unique_only_up_to_constant():
             # The interior Neumann operator in this representation

test/test_matrix.py

@@ -52,9 +52,13 @@ def max_cluster_error(mat, clusters, mindex, p=None):
     error = -np.inf
     for i in range(mindex.nclusters):
         mat_i = mindex.cluster_take(mat, i, i)
+        norm_mat_i = la.norm(mat_i, ord=p)
+        if norm_mat_i < 1.0e-12:
+            norm_mat_i = 1.0
+
         error = max(
                 error,
-                la.norm(mat_i - clusters[i, i], ord=p) / la.norm(mat_i, ord=p)
+                la.norm(mat_i - clusters[i, i], ord=p) / norm_mat_i
                 )
 
     return error

test/test_symbolic.py

@@ -533,20 +533,23 @@ def test_derivative_with_spatial_constant():
     ambient_dim = 3
 
     from sumpy.kernel import LaplaceKernel
+    knl = LaplaceKernel(ambient_dim)
+    density = sym.var("sigma")
+
     sym.d_dx(ambient_dim,
             sym.SpatialConstant("kappa")
-            * sym.D(LaplaceKernel(ambient_dim), sym.Variable("sigma")))
+            * sym.D(knl, density, qbx_forced_limit="avg"))
 
     from sumpy.kernel import LaplaceKernel
     sym.d_dx(ambient_dim,
             (3+sym.SpatialConstant("kappa"))
-            * sym.D(LaplaceKernel(ambient_dim), sym.Variable("sigma")))
+            * sym.D(knl, density, qbx_forced_limit="avg"))
 
     from pytential.symbolic.mappers import _DerivativeTakerUnsupoortedProductError
     with pytest.raises(_DerivativeTakerUnsupoortedProductError):
         sym.d_dx(ambient_dim,
                 (3+sym.Variable("kappa"))
-                * sym.D(LaplaceKernel(ambient_dim), sym.Variable("sigma")))
+                * sym.D(knl, density, qbx_forced_limit="avg"))
 
 # }}}
 

test/test_tools.py

@@ -66,8 +66,6 @@ def test_gmres():
 # {{{ test_interpolatory_error_reporting
 
 def test_interpolatory_error_reporting(actx_factory):
-    logging.basicConfig(level=logging.INFO)
-
     actx = actx_factory()
 
     h = 0.2
@@ -97,7 +95,9 @@ def test_interpolatory_error_reporting(actx_factory):
     one = 1 + 0*vol_x[0]
     from meshmode.discretization import NoninterpolatoryElementGroupError
     with pytest.raises(NoninterpolatoryElementGroupError):
-        print("AREA", integral(vol_discr, one), 0.25**2*np.pi)
+        logger.info("AREA integral %g exact %g",
+                    actx.to_numpy(integral(vol_discr, one)).item(),
+                    0.25**2*np.pi)
 
 # }}}
 
@@ -154,7 +154,7 @@ def test_geometry_collection_caching(actx_factory):
     # construct a geometry collection
     from pytential import GeometryCollection
     places = GeometryCollection(dict(zip(sources, lpots)), auto_where=sources[0])
-    print(places.places)
+    logger.info("%s", places.places)
 
     # check on-demand refinement
     from pytential import bind, sym