# Design: Sequence Decoder Generating LoDTensors

In tasks such as machine translation and visual captioning, a sequence decoder is necessary to generate sequences, one word at a time.

This documentation describes how to implement the sequence decoder as an operator.

## Beam Search based Decoder

The beam search algorithm is necessary when generating sequences. It is a heuristic search algorithm that explores the paths by expanding the most promising node in a limited set.

In the old version of PaddlePaddle, the C++ class RecurrentGradientMachine implements the general sequence decoder based on beam search, due to the complexity involved, the implementation relies on a lot of special data structures that are quite trivial and hard to be customized by users.

There are a lot of heuristic tricks in the sequence generation tasks, so the flexibility of sequence decoder is very important to users.

During the refactoring of PaddlePaddle, some new concepts are proposed such as: LoDTensor and TensorArray that can better support the sequence usage, and they can also help make the implementation of beam search based sequence decoder more transparent and modular .

For example, the RNN states, candidates IDs and probabilities of beam search can be represented all as LoDTensors; the selected candidate's IDs in each time step can be stored in a TensorArray, and Packed to the sentences translated.

## Changing LoD's absolute offset to relative offsets

The current LoDTensor is designed to store levels of variable-length sequences. It stores several arrays of integers where each represents a level.

The integers in each level represent the begin and end (not inclusive) offset of a sequence in the underlying tensor, let's call this format the absolute-offset LoD for clarity.

The absolute-offset LoD can retrieve any sequence very quickly but fails to represent empty sequences, for example, a two-level LoD is as follows

[[0, 3, 9]
[0, 2, 3, 3, 3, 9]]

The first level tells that there are two sequences: - the first's offset is [0, 3) - the second's offset is [3, 9)

while on the second level, there are several empty sequences that both begin and end at 3. It is impossible to tell how many empty second-level sequences exist in the first-level sequences.

There are many scenarios that rely on empty sequence representation, for example in machine translation or visual captioning, one instance has no translation or the empty candidate set for a prefix.

So let's introduce another format of LoD, it stores the offsets of the lower level sequences and is called relative-offset LoD.

For example, to represent the same sequences of the above data

[[0, 3, 6]
[0, 2, 3, 3, 3, 9]]


the first level represents that there are two sequences, their offsets in the second-level LoD is [0, 3) and [3, 5).

The second level is the same with the relative offset example because the lower level is a tensor. It is easy to find out the second sequence in the first-level LoD has two empty sequences.

The following examples are based on relative-offset LoD.

## Usage in a simple machine translation model

Let's start from a simple machine translation model that is simplified from the machine translation chapter to draw a blueprint of what a sequence decoder can do and how to use it.

The model has an encoder that learns the semantic vector from a sequence, and a decoder which uses the sequence encoder to generate new sentences.

Encoder

import paddle as pd

dict_size = 8000
source_dict_size = dict_size
target_dict_size = dict_size
word_vector_dim = 128
encoder_dim = 128
decoder_dim = 128
beam_size = 5
max_length = 120

# encoder
src_word_id = pd.data(
name='source_language_word',
type=pd.data.integer_value_sequence(source_dict_dim))
src_embedding = pd.embedding(size=source_dict_size, size=word_vector_dim)

src_word_vec = pd.lookup(src_embedding, src_word_id)

encoder_out_seq = pd.gru(input=src_word_vec, size=encoder_dim)

encoder_ctx = pd.last_seq(encoder_out_seq)
# encoder_ctx_proj is the learned semantic vector
encoder_ctx_proj = pd.fc(
encoder_ctx, size=decoder_dim, act=pd.activation.Tanh(), bias=None)


Decoder

def generate():
decoder = pd.while_loop()
with decoder.step():
decoder_mem = decoder.memory(init=encoder_ctx)  # mark the memory
generated_ids = decoder.memory() # TODO init to batch_size <s>s
generated_scores = decoder.memory() # TODO init to batch_size 1s or 0s

target_word = pd.lookup(trg_embedding, gendrated_ids)
# expand encoder_ctx's batch to fit target_word's lod
# for example
# decoder_mem.lod is
# [[0 1 3],
#  [0 1 3 6]]
# its tensor content is [a1 a2 a3 a4 a5]
# which means there are 2 sentences to translate
#   - the first sentence has 1 translation prefixes, the offsets are [0, 1)
#   - the second sentence has 2 translation prefixes, the offsets are [1, 3) and [3, 6)
# the target_word.lod is
# [[0, 1, 6]
#  [0, 2, 4, 7, 9 12]]
# which means 2 sentences to translate, each has 1 and 5 prefixes
# the first prefix has 2 candidates
# the following has 2, 3, 2, 3 candidates
# the encoder_ctx_expanded's content will be
# [a1 a1 a2 a2 a3 a3 a3 a4 a4 a5 a5 a5]
encoder_ctx_expanded = pd.lod_expand(encoder_ctx, target_word)
decoder_input = pd.fc(
act=pd.activation.Linear(),
input=[target_word, encoder_ctx_expanded],
size=3 * decoder_dim)
gru_out, cur_mem = pd.gru_step(
decoder_input, mem=decoder_mem, size=decoder_dim)
scores = pd.fc(
gru_out,
size=trg_dic_size,
bias=None,
act=pd.activation.Softmax())
# K is an config
topk_scores, topk_ids = pd.top_k(scores, K)

selected_ids, selected_generation_scores = decoder.beam_search(
topk_ids, topk_generated_scores)

# update the states
decoder_mem.update(cur_mem)  # tells how to update state
generated_ids.update(selected_ids)
generated_scores.update(selected_generation_scores)

decoder.output(selected_ids)
decoder.output(selected_generation_scores)

translation_ids, translation_scores = decoder()

The decoder.beam_search is an operator that, given the candidates and the scores of translations including the candidates, returns the result of the beam search algorithm.

In this way, users can customize anything on the input or output of beam search, for example:

1. Make the corresponding elements in topk_generated_scores zero or some small values, beam_search will discard this candidate.
2. Remove some specific candidate in selected_ids.
3. Get the final translation_ids, remove the translation sequence in it.

The implementation of sequence decoder can reuse the C++ class: RNNAlgorithm, so the python syntax is quite similar to that of an RNN.

Both of them are two-level LoDTensors:

• The first level represents batch_size of (source) sentences.
• The second level represents the candidate ID sets for translation prefix.

For example, 3 source sentences to translate, and has 2, 3, 1 candidates.

Unlike an RNN, in sequence decoder, the previous state and the current state have different LoD and shape, and an lod_expand operator is used to expand the LoD of the previous state to fit the current state.

For example, the previous state:

• LoD is [0, 1, 3][0, 2, 5, 6]
• content of tensor is a1 a2 b1 b2 b3 c1

the current state is stored in encoder_ctx_expanded:

• LoD is [0, 2, 7][0 3 5 8 9 11 11]
• the content is
• a1 a1 a1 (a1 has 3 candidates, so the state should be copied 3 times for each candidates)
• a2 a2
• b1 b1 b1
• b2
• b3 b3
• None (c1 has 0 candidates, so c1 is dropped)

The benefit from the relative offset LoD is that the empty candidate set can be represented naturally.

The status in each time step can be stored in TensorArray, and Packed to a final LoDTensor. The corresponding syntax is:

decoder.output(selected_ids)
decoder.output(selected_generation_scores)


The selected_ids are the candidate ids for the prefixes, and will be Packed by TensorArray to a two-level LoDTensor, where the first level represents the source sequences and the second level represents generated sequences.

Packing the selected_scores will get a LoDTensor that stores scores of each translation candidate.

Packing the selected_generation_scores will get a LoDTensor, and each tail is the probability of the translation.

## LoD and shape changes during decoding

According to the image above, the only phase that changes the LoD is beam search.

## Beam search design

The beam search algorithm will be implemented as one method of the sequence decoder and has 3 inputs:

1. topk_ids, the top K candidate ids for each prefix.
2. topk_scores, the corresponding scores for topk_ids
3. generated_scores, the score of the prefixes.

All of these are LoDTensors, so that the sequence affiliation is clear. Beam search will keep a beam for each prefix and select a smaller candidate set for each prefix.

It will return three variables:

1. selected_ids, the final candidate beam search function selected for the next step.
2. selected_scores, the scores for the candidates.
3. generated_scores, the updated scores for each prefix (with the new candidates appended).

## Introducing the LoD-based Pack and Unpack methods in TensorArray

The selected_ids, selected_scores and generated_scores are LoDTensors that exist at each time step, so it is natural to store them in arrays.

Currently, PaddlePaddle has a module called TensorArray which can store an array of tensors. It is better to store the results of beam search in a TensorArray.

The Pack and UnPack in TensorArray are used to pack tensors in the array to an LoDTensor or split the LoDTensor to an array of tensors. It needs some extensions to support the packing or unpacking an array of LoDTensors.