Completed
Push — master ( dccd0d...37cade )
by Raphael
01:33
created

RecurrentLayer.get_step_inputs()   C

Complexity

Conditions 7

Size

Total Lines 22

Duplication

Lines 0
Ratio 0 %

Importance

Changes 1
Bugs 1 Features 0
Metric Value
cc 7
dl 0
loc 22
rs 5.7894
c 1
b 1
f 0
1
#!/usr/bin/env python
2
# -*- coding: utf-8 -*-
3
4
from . import NeuralLayer
5
from var import NeuralVariable
6
from deepy.utils import build_activation, FLOATX, XavierGlorotInitializer, OrthogonalInitializer, Scanner, neural_computation
7
import numpy as np
8
import theano.tensor as T
9
from abc import ABCMeta, abstractmethod
10
11
OUTPUT_TYPES = ["sequence", "one"]
12
INPUT_TYPES = ["sequence", "one"]
13
14
15
16
class RecurrentLayer(NeuralLayer):
17
    __metaclass__ = ABCMeta
18
19
    def __init__(self, name, state_names, hidden_size=100, input_type="sequence", output_type="sequence",
20
                 inner_init=None, outer_init=None,
21
                 gate_activation='sigmoid', activation='tanh',
22
                 steps=None, backward=False, mask=None,
23
                 additional_input_dims=None):
24
        super(RecurrentLayer, self).__init__(name)
25
        self.state_names = state_names
26
        self.main_state = state_names[0]
27
        self.hidden_size = hidden_size
28
        self._gate_activation = gate_activation
29
        self._activation = activation
30
        self.gate_activate = build_activation(self._gate_activation)
31
        self.activate = build_activation(self._activation)
32
        self._input_type = input_type
33
        self._output_type = output_type
34
        self.inner_init = inner_init if inner_init else OrthogonalInitializer()
35
        self.outer_init = outer_init if outer_init else XavierGlorotInitializer()
36
        self._steps = steps
37
        self._mask = mask.tensor if type(mask) == NeuralVariable else mask
38
        self._go_backwards = backward
39
        self.additional_input_dims = additional_input_dims if additional_input_dims else []
40
41
        if input_type not in INPUT_TYPES:
42
            raise Exception("Input type of {} is wrong: {}".format(name, input_type))
43
        if output_type not in OUTPUT_TYPES:
44
            raise Exception("Output type of {} is wrong: {}".format(name, output_type))
45
46
    @neural_computation
47
    def step(self, step_inputs):
48
        new_states = self.compute_new_state(step_inputs)
49
50
        # apply mask for each step for any `output_type`
51
        if step_inputs.get("mask"):
52
            mask = step_inputs["mask"].dimshuffle(0, 'x')
53
            for state_name in new_states:
54
                new_states[state_name] = new_states[state_name] * mask + step_inputs[state_name] * (1 - mask)
55
56
        return new_states
57
58
    @abstractmethod
59
    def compute_new_state(self, step_inputs):
60
        """
61
        :type step_inputs: dict
62
        :rtype: dict
63
        """
64
65
    @abstractmethod
66
    def merge_inputs(self, input_var, additional_inputs=None):
67
        """
68
        Merge inputs and return a map, which will be passed to core_step.
69
        :type input_var: T.var
70
        :param additional_inputs: list
71
        :rtype: dict
72
        """
73
74
    @abstractmethod
75
    def prepare(self):
76
        pass
77
78
    @neural_computation
79
    def get_initial_states(self, input_var):
80
        """
81
        :type input_var: T.var
82
        :rtype: dict
83
        """
84
        initial_states = {}
85
        for state in self.state_names:
86
            initial_states[state] = T.alloc(np.cast[FLOATX](0.), input_var.shape[0], self.hidden_size)
87
        return initial_states
88
89
    @neural_computation
90
    def get_step_inputs(self, input_var, states=None, mask=None, additional_inputs=None):
91
        """
92
        :type input_var: T.var
93
        :rtype: dict
94
        """
95
        step_inputs = {}
96
        if self._input_type == "sequence":
97
            if not additional_inputs:
98
                additional_inputs = []
99
            if mask:
100
                step_inputs['mask'] = mask.dimshuffle(1, 0)
101
            step_inputs.update(self.merge_inputs(input_var, additional_inputs=additional_inputs))
102
        else:
103
            # step_inputs["mask"] = mask.dimshuffle((1,0)) if mask else None
104
            if additional_inputs:
105
                step_inputs.update(self.merge_inputs(None, additional_inputs=additional_inputs))
106
        if states:
107
            for name in self.state_names:
108
                step_inputs[name] = states[name]
109
110
        return step_inputs
111
112
    def compute(self, input_var, mask=None, additional_inputs=None, steps=None, backward=False):
113
        if additional_inputs and not self.additional_input_dims:
114
            self.additional_input_dims = map(lambda var: var.dim(), additional_inputs)
115
        return super(RecurrentLayer, self).compute(input_var, mask=mask, additional_inputs=additional_inputs, steps=steps, backward=backward)
116
117
    def compute_tensor(self, input_var, mask=None, additional_inputs=None, steps=None, backward=False):
118
        # prepare parameters
119
        backward = backward if backward else self._go_backwards
120
        steps = steps if steps else self._steps
121
        mask = mask if mask else self._mask
122
        if mask and self._input_type == "one":
123
            raise Exception("Mask only works with sequence input")
124
        # get initial states
125
        init_state_map = self.get_initial_states(input_var)
126
        # get input sequence map
127
        if self._input_type == "sequence":
128
            # Move middle dimension to left-most position
129
            # (sequence, batch, value)
130
            input_var = input_var.dimshuffle((1,0,2))
131
            seq_map = self.get_step_inputs(input_var, mask=mask, additional_inputs=additional_inputs)
132
        else:
133
            init_state_map[self.main_state] = input_var
134
            seq_map = self.get_step_inputs(None, mask=mask, additional_inputs=additional_inputs)
135
        # scan
136
        retval_map, _ = Scanner(
137
            self.step,
138
            sequences=seq_map,
139
            outputs_info=init_state_map,
140
            n_steps=steps,
141
            go_backwards=backward
142
        ).compute()
143
        # return main states
144
        main_states = retval_map[self.main_state]
145
        if self._output_type == "one":
146
            return main_states[-1]
147
        elif self._output_type == "sequence":
148
            main_states = main_states.dimshuffle((1,0,2)) # ~ batch, time, size
149
            # if mask: # ~ batch, time
150
            #     main_states *= mask.dimshuffle((0, 1, 'x'))
151
            return main_states
152
153
154
class RNN(RecurrentLayer):
155
156
    def  __init__(self, hidden_size, **kwargs):
157
        kwargs["hidden_size"] = hidden_size
158
        super(RNN, self).__init__("RNN", ["state"], **kwargs)
159
160
    @neural_computation
161
    def compute_new_state(self, step_inputs):
162
        xh_t, h_tm1 = map(step_inputs.get, ["xh_t", "state"])
163
        if not xh_t:
164
            xh_t = 0
165
166
        h_t = self.activate(xh_t + T.dot(h_tm1, self.W_h) + self.b_h)
167
168
        return {"state": h_t}
169
170
    @neural_computation
171
    def merge_inputs(self, input_var, additional_inputs=None):
172
        if not additional_inputs:
173
            additional_inputs = []
174
        all_inputs = ([input_var] if input_var else []) + additional_inputs
175
        h_inputs = []
176
        for x, weights in zip(all_inputs, self.input_weights):
177
            wi, = weights
178
            h_inputs.append(T.dot(x, wi))
179
        merged_inputs = {
180
            "xh_t": sum(h_inputs)
181
        }
182
        return merged_inputs
183
184
    def prepare(self):
185
        self.output_dim = self.hidden_size
186
187
        self.W_h = self.create_weight(self.hidden_size, self.hidden_size, "h", initializer=self.outer_init)
188
        self.b_h = self.create_bias(self.hidden_size, "h")
189
190
        self.register_parameters(self.W_h, self.b_h)
191
192
        self.input_weights = []
193
        if self._input_type == "sequence":
194
            normal_input_dims = [self.input_dim]
195
        else:
196
            normal_input_dims = []
197
198
        all_input_dims = normal_input_dims + self.additional_input_dims
199
        for i, input_dim in enumerate(all_input_dims):
200
            wi = self.create_weight(input_dim, self.hidden_size, "wi_{}".format(i+1), initializer=self.outer_init)
201
            weights = [wi]
202
            self.input_weights.append(weights)
203
            self.register_parameters(*weights)