diff --git a/firedrake/pointquery_utils.py b/firedrake/pointquery_utils.py index 739703506c..2653db12c6 100644 --- a/firedrake/pointquery_utils.py +++ b/firedrake/pointquery_utils.py @@ -146,15 +146,12 @@ def to_reference_coords_newton_step(ufl_coordinate_element, parameters, x0_dtype x0_expr = builder._coefficient(x0, "x0") loopy_args = [ lp.GlobalArg( - "C", dtype=ScalarType, shape=(numpy.prod(Cexpr.shape, dtype=int),)), + "C", dtype=ScalarType, shape=(numpy.prod(Cexpr.shape, dtype=int),) + ), lp.GlobalArg( - "x0", dtype=x0_dtype, shape=(numpy.prod(x0_expr.shape, dtype=int),))] - - loopy_args = [ - lp.GlobalArg( - "C", dtype=ScalarType_c, shape=(numpy.prod(Cexpr.shape, dtype=int),)), - lp.GlobalArg( - "x0", dtype=ScalarType, shape=(numpy.prod(x0_expr.shape, dtype=int),))] + "x0", dtype=x0_dtype, shape=(numpy.prod(x0_expr.shape, dtype=int),) + ), + ] dim = cell.topological_dimension() point = gem.Variable('X', (dim,))