/usr/lib/python2.7/dist-packages/ufl/checks.py is in python-ufl 1.4.0-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 | """Utility functions for checking properties of expressions."""
# Copyright (C) 2013-2014 Martin Sandve Alnes
#
# This file is part of UFL.
#
# UFL is free software: you can redistribute it and/or modify
# it under the terms of the GNU Lesser General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# UFL is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU Lesser General Public License for more details.
#
# You should have received a copy of the GNU Lesser General Public License
# along with UFL. If not, see <http://www.gnu.org/licenses/>.
#
# Modified by Anders Logg, 2008-2009
from ufl.expr import Expr
def is_python_scalar(expression):
"Return True iff expression is of a Python scalar type."
return isinstance(expression, (int, float))
def is_ufl_scalar(expression):
"""Return True iff expression is scalar-valued,
but possibly containing free indices."""
return isinstance(expression, Expr) and not expression.shape()
def is_true_ufl_scalar(expression):
"""Return True iff expression is scalar-valued,
with no free indices."""
return isinstance(expression, Expr) and \
not (expression.shape() or expression.free_indices())
def is_globally_constant(expr):
"""Check if an expression is globally constant, which
includes spatially independent constant coefficients that
are not known before assembly time."""
# TODO: This does not consider gradients of coefficients, so false negatives are possible.
from ufl.algorithms.traversal import traverse_unique_terminals
from ufl.argument import Argument
from ufl.coefficient import Coefficient
for e in traverse_unique_terminals(expr):
if isinstance(e, Argument):
return False
if isinstance(e, Coefficient) and e.element().family() != "Real":
return False
if not e.is_cellwise_constant():
return False
# All terminals passed constant check
return True
def is_scalar_constant_expression(expr):
"""Check if an expression is a globally constant scalar expression."""
if is_python_scalar(expr):
return True
if expr.shape() != ():
return False
return is_globally_constant(expr)
|