#!/usr/bin/env python3
# -*- coding: utf-8 -*-
#
# Copyright (C) IBM Corporation 2018
#
# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at
#
# http://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and
# limitations under the License.
"""
repeat_reverse_recall_cl.py: Contains definition of repeat reverse recall problem with control markers and\
command lines
"""
__author__ = "Tomasz Kornuta & Vincent Marois"
import torch
import numpy as np
from miprometheus.problems.seq_to_seq.algorithmic.algorithmic_seq_to_seq_problem import AlgorithmicSeqToSeqProblem
[docs]class RepeatReverseRecallCommandLines(AlgorithmicSeqToSeqProblem):
"""
Class generating sequences of random bit-patterns and targets forcing the
system to learn repeated reverse recall problem.
1. There are 2 markers, indicating:
- beginning of storing/memorization,
- beginning of forward recalling from memory.
2. Additionally, there is a command line (3rd command bit) indicating whether given item is to be stored in\
memory (0) or recalled (1).
"""
[docs] def __init__(self, params):
"""
Constructor - stores parameters. Calls parent class ``AlgorithmicSeqToSeqProblem``\
initialization.
:param params: Dictionary of parameters (read from configuration ``.yaml`` file).
"""
# Set default number of bits for a given problem.
# This has to be done before calling base class constructor!
params.add_default_params({
'control_bits': 2,
'data_bits': 8 })
# Call parent constructor - sets e.g. the loss function, dtype.
# Additionally it extracts "standard" list of parameters for
# algorithmic tasks, like batch_size, numbers of bits, sequences etc.
super(RepeatReverseRecallCommandLines, self).__init__(params)
self.name = 'RepeatReverseRecallCommandLines'
assert self.control_bits >= 2, "Problem requires at least 2 control bits (currently %r)" % self.control_bits
assert self.data_bits >= 1, "Problem requires at least 1 data bit (currently %r)" % self.data_bits
# Min and max number of recalls.
self.params.add_default_params({
'min_recall_number': 1,
'max_recall_number': 5})
self.min_recall_number = params['min_recall_number']
self.max_recall_number = params['max_recall_number']
[docs] def generate_batch(self, batch_size):
"""
Generates a batch of samples of size ''batch_size'' on-the-fly.
.. note::
The sequence length is drawn randomly between ``self.min_sequence_length`` and \
``self.max_sequence_length``.
.. warning::
All the samples within the batch will have the same sequence lengt.
:param batch_size: Size of the batch to be returned.
:return: DataDict({'sequences', 'sequences_length', 'targets', 'masks', 'num_subsequences'}), with:
- sequences: [BATCH_SIZE, 2*SEQ_LENGTH+2, CONTROL_BITS+DATA_BITS]
- sequences_length: [BATCH_SIZE, 1] (the same random value between self.min_sequence_length and self.max_sequence_length)
- targets: [BATCH_SIZE, , 2*SEQ_LENGTH+2, DATA_BITS]
- masks: [BATCH_SIZE, 2*SEQ_LENGTH+2, 1]
- num_subsequences: [BATCH_SIZE, 1]
"""
# Store marker.
marker_start_main = np.zeros(self.control_bits)
marker_start_main[self.store_bit] = 1 # [1, 0, 0]
# Recall marker.
marker_start_aux = np.zeros(self.control_bits)
marker_start_aux[self.recall_bit] = 1 # [0, 1, 0]
# Define control lines.
ctrl_aux = np.zeros(self.control_bits)
if self.use_control_lines:
if self.control_bits >= 3:
if self.randomize_control_lines:
# Randomly pick one of the bits to be set.
ctrl_bit = np.random.randint(2, self.control_bits)
ctrl_aux[ctrl_bit] = 1
else:
# Set last.
ctrl_aux[self.control_bits - 1] = 1
# Else: no control lines!
# Set sequence length.
seq_length = np.random.randint(
self.min_sequence_length, self.max_sequence_length + 1)
# Number of recalls.
recall_number = np.random.randint(
self.min_recall_number, self.max_recall_number + 1)
# Generate batch of random bit sequences [BATCH_SIZE x SEQ_LENGTH X
# DATA_BITS]
bit_seq = np.random.binomial(
1, self.bias, (batch_size, seq_length, self.data_bits))
# 1. Generate inputs.
# Generate input: [BATCH_SIZE, 3*SEQ_LENGTH+3, CONTROL_BITS+DATA_BITS]
inputs = np.zeros([batch_size,
(recall_number + 1) * (seq_length + 1),
self.control_bits + self.data_bits],
dtype=np.float32)
# Set start main control marker.
inputs[:, 0, 0:self.control_bits] = np.tile(
marker_start_main, (batch_size, 1))
# Set bit sequence.
inputs[:, 1:seq_length + 1,
self.control_bits:self.control_bits + self.data_bits] = bit_seq
for r in range(recall_number):
# Set start aux serial recall control marker.
inputs[:,
(r + 1) * (seq_length + 1),
0:self.control_bits] = np.tile(marker_start_aux,
(batch_size,
1))
inputs[:,
(r + 1) * (seq_length + 1) + 1:(r + 2) * (seq_length + 1),
0:self.control_bits] = np.tile(ctrl_aux,
(batch_size,
seq_length,
1))
# 2. Generate targets.
# Generate target: [BATCH_SIZE, 3*SEQ_LENGTH+3, DATA_BITS] (only data
# bits!)
targets = np.zeros([batch_size,
(recall_number + 1) * (seq_length + 1),
self.data_bits],
dtype=np.float32)
# Set bit sequence for serial recall.
for r in range(recall_number):
targets[:, (r + 1) * (seq_length + 1) + 1:(r + 2) *
(seq_length + 1), :] = np.fliplr(bit_seq)
# 3. Generate mask.
# Generate target mask: [BATCH_SIZE, 3*SEQ_LENGTH+3, 1]
ptmasks = torch.zeros([batch_size, (recall_number + 1) *
(seq_length + 1), 1]).type(self.app_state.ByteTensor)
for r in range(recall_number):
ptmasks[:, (r + 1) * (seq_length + 1) +
1:(r + 2) * (seq_length + 1)] = 1
# Return data_dict.
data_dict = self.create_data_dict()
data_dict['sequences'] = torch.from_numpy(inputs).type(self.app_state.dtype)
data_dict['targets'] = torch.from_numpy(targets).type(self.app_state.dtype)
data_dict['masks'] = ptmasks
data_dict['sequences_length'] = torch.ones([batch_size,1]).type(torch.CharTensor) * seq_length
data_dict['num_subsequences'] = torch.ones([batch_size, 1]).type(torch.CharTensor)
return data_dict
if __name__ == "__main__":
""" Tests sequence generator - generates and displays a random sample"""
# "Loaded parameters".
from miprometheus.utils.param_interface import ParamInterface
params = ParamInterface()
params.add_config_params({#'control_bits': 4,
#'data_bits': 8,
# 'randomize_control_lines': False,
'min_sequence_length': 1,
'max_sequence_length': 10})
batch_size = 64
# Create problem object.
repeatreverserecallcl = RepeatReverseRecallCommandLines(params)
# get a sample
sample = repeatreverserecallcl[0]
print(repr(sample))
print('__getitem__ works.')
# wrap DataLoader on top
from torch.utils.data import DataLoader
def init_fn(worker_id):
np.random.seed(seed=worker_id)
problem = DataLoader(dataset=repeatreverserecallcl, batch_size=batch_size, collate_fn=repeatreverserecallcl.collate_fn,
shuffle=False, num_workers=0, worker_init_fn=init_fn)
# generate a batch
import time
s = time.time()
for i, batch in enumerate(problem):
#print('Batch # {} - {}'.format(i, type(batch)))
pass
print('Number of workers: {}'.format(problem.num_workers))
print('time taken to exhaust a dataset of size {}, with a batch size of {}: {}s'
.format(len(repeatreverserecallcl), batch_size, time.time() - s))
# Display single sample (0) from batch.
batch = next(iter(problem))
repeatreverserecallcl.show_sample(batch, 0)
print('Unit test completed.')