提交 87850539 authored 作者: Iban Harlouchet's avatar Iban Harlouchet

flake8 for neighbourhoods.py

上级 c2cded9d
#!/usr/bin/python
"""WARNING: This code is not recommanded. It is not finished, it is """WARNING: This code is not recommanded. It is not finished, it is
slower then the version in sandbox/neighbours.py, and it do not work slower then the version in sandbox/neighbours.py, and it do not work
on the GPU. on the GPU.
...@@ -8,13 +7,13 @@ it cover more cases. But thoses cases aren't needed frequently, so you ...@@ -8,13 +7,13 @@ it cover more cases. But thoses cases aren't needed frequently, so you
probably don't want to use this version, go see neighbours.py!!!!!!! probably don't want to use this version, go see neighbours.py!!!!!!!
""" """
import theano
from theano import gof, Op, tensor, Variable, Apply
import numpy import numpy
from six.moves import xrange from six.moves import xrange
import six.moves.builtins as builtins import six.moves.builtins as builtins
import theano
from theano import gof, Op
class NeighbourhoodsFromImages(Op): class NeighbourhoodsFromImages(Op):
def __init__(self, n_dims_before, dims_neighbourhoods, def __init__(self, n_dims_before, dims_neighbourhoods,
...@@ -75,7 +74,7 @@ class NeighbourhoodsFromImages(Op): ...@@ -75,7 +74,7 @@ class NeighbourhoodsFromImages(Op):
""" """
self.n_dims_before = n_dims_before self.n_dims_before = n_dims_before
self.dims_neighbourhoods = dims_neighbourhoods self.dims_neighbourhoods = dims_neighbourhoods
if not strides is None: if strides is not None:
self.strides = strides self.strides = strides
else: else:
self.strides = dims_neighbourhoods self.strides = dims_neighbourhoods
...@@ -85,14 +84,6 @@ class NeighbourhoodsFromImages(Op): ...@@ -85,14 +84,6 @@ class NeighbourhoodsFromImages(Op):
self.code_string, self.code = self.make_py_code() self.code_string, self.code = self.make_py_code()
def _compute_neigh_strides(self):
neigh_strides = [1 for i in xrange(len(self.strides))]
cur_stride = 1
for i in xrange(len(self.strides)-1, -1, -1):
neigh_strides[i] = cur_stride
cur_stride *= self.dims_neighbourhoods[i]
return neigh_strides
def __eq__(self, other): def __eq__(self, other):
return type(self) == type(other) and \ return type(self) == type(other) and \
self.n_dims_before == other.n_dims_before and \ self.n_dims_before == other.n_dims_before and \
...@@ -108,8 +99,7 @@ class NeighbourhoodsFromImages(Op): ...@@ -108,8 +99,7 @@ class NeighbourhoodsFromImages(Op):
hash(self.ignore_border) hash(self.ignore_border)
def __str__(self): def __str__(self):
return '%s{%s,%s,%s,%s}' % \ return '%s{%s,%s,%s,%s}' % (self.__class__.__name__,
(self.__class__.__name__,
self.n_dims_before, self.n_dims_before,
self.dims_neighbourhoods, self.dims_neighbourhoods,
self.strides, self.strides,
...@@ -163,11 +153,11 @@ class NeighbourhoodsFromImages(Op): ...@@ -163,11 +153,11 @@ class NeighbourhoodsFromImages(Op):
x = theano.tensor.as_tensor_variable(x) x = theano.tensor.as_tensor_variable(x)
if self.inverse: if self.inverse:
# +1 in the inverse case # +1 in the inverse case
if x.type.ndim != (self.n_dims_before + \ if x.type.ndim != (self.n_dims_before +
len(self.dims_neighbourhoods) + 1): len(self.dims_neighbourhoods) + 1):
raise TypeError() raise TypeError()
else: else:
if x.type.ndim != (self.n_dims_before + \ if x.type.ndim != (self.n_dims_before +
len(self.dims_neighbourhoods)): len(self.dims_neighbourhoods)):
raise TypeError() raise TypeError()
return gof.Apply(self, [x], [x.type()]) return gof.Apply(self, [x], [x.type()])
...@@ -177,22 +167,24 @@ class NeighbourhoodsFromImages(Op): ...@@ -177,22 +167,24 @@ class NeighbourhoodsFromImages(Op):
z, = out z, = out
if self.inverse: if self.inverse:
# +1 in the inverse case # +1 in the inverse case
if len(x.shape) != (self.n_dims_before + \ if len(x.shape) != (self.n_dims_before +
len(self.dims_neighbourhoods) + 1): len(self.dims_neighbourhoods) + 1):
raise ValueError("Images passed as input don't match the " +\ raise ValueError("Images passed as input don't match the "
"dimensions passed when this (inversed) Apply node was created") "dimensions passed when this (inversed) "
"Apply node was created")
prod = 1 prod = 1
for dim in self.dims_neighbourhoods: for dim in self.dims_neighbourhoods:
prod *= dim prod *= dim
if x.shape[-1] != prod: if x.shape[-1] != prod:
raise ValueError(("Last dimension of neighbourhoods (%s) is not " +\ raise ValueError("Last dimension of neighbourhoods (%s) is not"
"the product of the neighbourhoods dimensions (%s)") % \ " the product of the neighbourhoods dimensions"
(str(x.shape[-1]), str(prod))) " (%s)" % (str(x.shape[-1]), str(prod)))
else: else:
if len(x.shape) != (self.n_dims_before + \ if len(x.shape) != (self.n_dims_before +
len(self.dims_neighbourhoods)): len(self.dims_neighbourhoods)):
raise ValueError("Images passed as input don't match the " +\ raise ValueError("Images passed as input don't match the "
"dimensions passed when this Apply node was created") "dimensions passed when this Apply node "
"was created")
if self.inverse: if self.inverse:
input_shape, num_strides = self.in_shape(x.shape) input_shape, num_strides = self.in_shape(x.shape)
...@@ -201,8 +193,6 @@ class NeighbourhoodsFromImages(Op): ...@@ -201,8 +193,6 @@ class NeighbourhoodsFromImages(Op):
input_shape = x.shape input_shape = x.shape
out_shape, num_strides = self.out_shape(input_shape) out_shape, num_strides = self.out_shape(input_shape)
neigh_strides = self._compute_neigh_strides()
if z[0] is None: if z[0] is None:
if self.inverse: if self.inverse:
z[0] = numpy.zeros(input_shape) z[0] = numpy.zeros(input_shape)
...@@ -228,7 +218,7 @@ class NeighbourhoodsFromImages(Op): ...@@ -228,7 +218,7 @@ class NeighbourhoodsFromImages(Op):
return code_before return code_before
def _py_innerloop(self, inner_dim_no): def _py_innerloop(self, inner_dim_no):
base_indent = ('\t' * (self.n_dims_before + inner_dim_no*2)) base_indent = ('\t' * (self.n_dims_before + inner_dim_no * 2))
code_before = base_indent + \ code_before = base_indent + \
"for stride_idx_%d in xrange(num_strides[%d]):\n" % \ "for stride_idx_%d in xrange(num_strides[%d]):\n" % \
(inner_dim_no, inner_dim_no) (inner_dim_no, inner_dim_no)
...@@ -237,36 +227,33 @@ class NeighbourhoodsFromImages(Op): ...@@ -237,36 +227,33 @@ class NeighbourhoodsFromImages(Op):
"dim_%d_offset = stride_idx_%d * self.strides[%d]\n" %\ "dim_%d_offset = stride_idx_%d * self.strides[%d]\n" %\
(inner_dim_no, inner_dim_no, inner_dim_no) (inner_dim_no, inner_dim_no, inner_dim_no)
code_before += base_indent + \ code_before += base_indent + \
"max_neigh_idx_%d = input_shape[%d] - dim_%d_offset\n" %\ "max_neigh_idx_%d = input_shape[%d] - dim_%d_offset\n" % \
(inner_dim_no, (inner_dim_no, self.n_dims_before + inner_dim_no, inner_dim_no)
self.n_dims_before+inner_dim_no, inner_dim_no)
code_before += base_indent + \ code_before += base_indent + \
("for neigh_idx_%d in xrange(min(max_neigh_idx_%d,"\ ("for neigh_idx_%d in xrange(min(max_neigh_idx_%d,"
+ " self.dims_neighbourhoods[%d])):\n") % \ " self.dims_neighbourhoods[%d])):\n") %\
(inner_dim_no, inner_dim_no, inner_dim_no) (inner_dim_no, inner_dim_no, inner_dim_no)
return code_before return code_before
def _py_flattened_idx(self): def _py_flattened_idx(self):
return "+".join(["neigh_strides[%d]*neigh_idx_%d" % (i, i) \ return "+".join(["neigh_strides[%d]*neigh_idx_%d" % (i, i)
for i in xrange(len(self.strides))]) for i in xrange(len(self.strides))])
def _py_assignment(self): def _py_assignment(self):
input_idx = "".join(["outer_idx_%d," % (i,) \ input_idx = "".join(["outer_idx_%d," % (i,)
for i in xrange(self.n_dims_before)]) for i in xrange(self.n_dims_before)])
input_idx += "".join(["dim_%d_offset+neigh_idx_%d," % \ input_idx += "".join(["dim_%d_offset+neigh_idx_%d," %
(i, i) for i in xrange(len(self.strides))]) (i, i) for i in xrange(len(self.strides))])
out_idx = "".join(\ out_idx = "".join(
["outer_idx_%d," % (i,) for i in \ ["outer_idx_%d," % (i,) for i in xrange(self.n_dims_before)] +
xrange(self.n_dims_before)] + \ ["stride_idx_%d," % (i,) for i in xrange(len(self.strides))])
["stride_idx_%d," % (i,) for i in \
xrange(len(self.strides))])
out_idx += self._py_flattened_idx() out_idx += self._py_flattened_idx()
#return_val = '\t' * (self.n_dims_before + len(self.strides)*2) # return_val = '\t' * (self.n_dims_before + len(self.strides)*2)
#return_val += "print "+input_idx+"'\\n',"+out_idx+"\n" # return_val += "print "+input_idx+"'\\n',"+out_idx+"\n"
return_val = '\t' * (self.n_dims_before + len(self.strides)*2) return_val = '\t' * (self.n_dims_before + len(self.strides) * 2)
if self.inverse: if self.inverse:
# remember z and x are inversed: # remember z and x are inversed:
...@@ -282,8 +269,9 @@ class NeighbourhoodsFromImages(Op): ...@@ -282,8 +269,9 @@ class NeighbourhoodsFromImages(Op):
class ImagesFromNeighbourhoods(NeighbourhoodsFromImages): class ImagesFromNeighbourhoods(NeighbourhoodsFromImages):
def __init__(self, n_dims_before, dims_neighbourhoods, def __init__(self, n_dims_before, dims_neighbourhoods,
strides=None, ignore_border=False): strides=None, ignore_border=False):
NeighbourhoodsFromImages.__init__(self, n_dims_before, dims_neighbourhoods, NeighbourhoodsFromImages.__init__(self, n_dims_before,
strides=strides, ignore_border=ignore_border, dims_neighbourhoods,
strides=strides,
ignore_border=ignore_border,
inverse=True) inverse=True)
# and that's all there is to it # and that's all there is to it
...@@ -137,7 +137,6 @@ whitelist_flake8 = [ ...@@ -137,7 +137,6 @@ whitelist_flake8 = [
"sandbox/rng_mrg.py", "sandbox/rng_mrg.py",
"sandbox/theano_object.py", "sandbox/theano_object.py",
"sandbox/scan.py", "sandbox/scan.py",
"sandbox/neighbourhoods.py",
"sandbox/fourier.py", "sandbox/fourier.py",
"sandbox/test_multinomial.py", "sandbox/test_multinomial.py",
"sandbox/minimal.py", "sandbox/minimal.py",
......
Markdown 格式
0%
您添加了 0 到此讨论。请谨慎行事。
请先完成此评论的编辑!
注册 或者 后发表评论