/usr/lib/python2.7/dist-packages/ufl/indexutils.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 | "Some utilities for working with index tuples."
# Copyright (C) 2008-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/>.
from collections import defaultdict
from ufl.indexing import Index
#--- Utility functions ---
def complete_shape(shape, default_dim):
"Complete shape tuple by replacing non-integers with a default dimension."
return tuple((s if isinstance(s, int) else default_dim) for s in shape)
def unique_indices(indices):
"""Return a tuple of all indices from input, with
no single index occuring more than once in output."""
s = set()
newindices = []
for i in indices:
if isinstance(i, Index):
if not i in s:
s.add(i)
newindices.append(i)
return tuple(newindices)
def repeated_indices(indices):
"Return tuple of indices occuring more than once in input."
ri = []
s = set()
for i in indices:
if isinstance(i, Index):
if i in s:
ri.append(i)
else:
s.add(i)
return tuple(ri)
def shared_indices(ai, bi):
"Return a tuple of indices occuring in both index tuples ai and bi."
bis = set(bi)
return tuple(i for i in unique_indices(ai) if i in bis)
def single_indices(indices):
"Return a tuple of all indices occuring exactly once in input."
count = defaultdict(int)
for i in indices:
count[i] += 1
return tuple(i for i in unique_indices(indices) if count[i] == 1)
|