From 1d7ce21a0b91b433de4d2a37ea25926c824d784e Mon Sep 17 00:00:00 2001 From: Miklos Homolya Date: Fri, 21 Oct 2016 10:34:41 +0100 Subject: [PATCH] cause not to trigger coneoproject/COFFEE#97 --- tsfc/coffee.py | 7 +------ 1 file changed, 1 insertion(+), 6 deletions(-) diff --git a/tsfc/coffee.py b/tsfc/coffee.py index 3cae8ed9..399dd013 100644 --- a/tsfc/coffee.py +++ b/tsfc/coffee.py @@ -182,11 +182,6 @@ def statement_evaluate(leaf, parameters): return coffee.Block(ops, open_scope=False) elif isinstance(expr, gem.Constant): assert parameters.declare[leaf] - # Take all axes except the last one - axes = tuple(range(len(expr.array.shape) - 1)) - nz_indices, = expr.array.any(axis=axes).nonzero() - nz_bounds = tuple([(i, 0)] for i in expr.array.shape[:-1]) - nz_bounds += ([(max(nz_indices) - min(nz_indices) + 1, min(nz_indices))],) table = numpy.array(expr.array) # FFC uses one less digits for rounding than for printing epsilon = eval("1e-%d" % (parameters.precision - 1)) @@ -195,7 +190,7 @@ def statement_evaluate(leaf, parameters): table[abs(table + 1.0) < epsilon] = -1.0 table[abs(table - 0.5) < epsilon] = 0.5 table[abs(table + 0.5) < epsilon] = -0.5 - init = coffee.SparseArrayInit(table, parameters.precision, nz_bounds) + init = coffee.ArrayInit(table, parameters.precision) return coffee.Decl(SCALAR_TYPE, _decl_symbol(expr, parameters), init,