From 7f2f128a75c32982d4f7e0a9a8d1050ad9990e78 Mon Sep 17 00:00:00 2001
From: Dominic Kempf <dominic.kempf@iwr.uni-heidelberg.de>
Date: Fri, 7 Dec 2018 14:31:49 +0100
Subject: [PATCH] Remove unused quadrature code
---
python/dune/codegen/sumfact/quadrature.py | 57 -----------------------
1 file changed, 57 deletions(-)
diff --git a/python/dune/codegen/sumfact/quadrature.py b/python/dune/codegen/sumfact/quadrature.py
index b68da208..92a598d9 100644
--- a/python/dune/codegen/sumfact/quadrature.py
+++ b/python/dune/codegen/sumfact/quadrature.py
@@ -179,63 +179,6 @@ def constructor_quadrature_inames(dim):
return tuple(constructor_quad_iname(name, d, local_qps_per_dir[d]) for d in range(local_dimension()))
-def define_recursive_quadrature_weight(visitor, name, direction):
- info = visitor.current_info[1]
- if info is None:
- element = None
- else:
- element = info.element
- iname = quadrature_inames(element)[direction]
- temporary_variable(name, shape=(), shape_impl=())
-
- qps_per_dir = quadrature_points_per_direction()
- qp_bound = qps_per_dir[direction]
- instruction(expression=Product((recursive_quadrature_weight(direction + 1),
- Subscript(Variable(name_oned_quadrature_weights(qp_bound)),
- (Variable(iname),)
- ))
- ),
- assignee=Variable(name),
- forced_iname_deps=frozenset(quadrature_inames()[direction:]),
- forced_iname_deps_is_final=True,
- tags=frozenset({"quad"}),
- )
-
-
-def recursive_quadrature_weight(visitor, direction=0):
- if direction == local_dimension():
- return pymbolic_base_weight()
- else:
- name = 'weight_{}'.format(direction)
- define_recursive_quadrature_weight(visitor, name, direction)
- return Variable(name)
-
-
-def define_quadrature_position(name, local_index):
- # TODO: This whole function looks quite suspicious:
- # - qps_per_dir not used
- # - probably rename index to local_index
- #
- # This whole function is probably not used anymore since we use sum
- # factorization for evaluation of global coordinate, jacobian inverse
- # transposed and determinant of the jacobian of the geometry mapping.
- #
- # Needs to be fixed. Or removed. Added to my to do list ;)
- assert False
-
- qps_per_dir = quadrature_points_per_direction()
- local_qps_per_dir = local_quadrature_points_per_direction()
- qp_bound = local_qps_per_dir[index]
- expression = Subscript(Variable(name_oned_quadrature_points(qp_bound)),
- (Variable(quadrature_inames()[index]),))
- instruction(expression=expression,
- assignee=Subscript(Variable(name), (index,)),
- forced_iname_deps=frozenset(quadrature_inames(element)),
- forced_iname_deps_is_final=True,
- tags=frozenset({"quad"}),
- )
-
-
def additional_inames(visitor):
"""Return inames for iterating over ansatz space in the case of jacobians
"""
--
GitLab