locoformer 0.0.7__py3-none-any.whl → 0.0.15__py3-none-any.whl

locoformer/locoformer.py

@@ -1,11 +1,23 @@
 from __future__ import annotations
 from functools import partial
 
+from pathlib import Path
+from contextlib import contextmanager
+from collections import namedtuple
+
+import numpy as np
+from numpy import ndarray
+from numpy.lib.format import open_memmap
+
+from beartype import beartype
+from beartype.door import is_bearable
+
 import torch
-from torch import nn, cat, stack, arange, Tensor, is_tensor
+from torch import nn, cat, stack, arange, Tensor, tensor, is_tensor, from_numpy
 import torch.nn.functional as F
 from torch.nn import Module, ModuleList, Linear, RMSNorm, Identity, Sequential
 from torch.utils._pytree import tree_map
+from torch.utils.data import Dataset, DataLoader
 
 import einx
 from einops import rearrange, einsum
@@ -31,20 +43,33 @@ def first(arr):
 def divisible_by(num, den):
     return (num % den) == 0
 
+# tensor helpers
+
+def log(t, eps = 1e-20):
+    return t.clamp_min(eps).log()
+
 def tree_map_tensor(x, fn):
     return tree_map(lambda t: t if not is_tensor(t) else fn(t), x)
 
 def detach_all(x):
     return tree_map_tensor(x, lambda t: t.detach())
 
-def combine_kv_cache(cache1, cache2):
-    combined_cache = []
+def pad_at_dim(
+    t,
+    pad: tuple[int, int],
+    dim = -1,
+    value = 0.
+):
+    if pad == (0, 0):
+        return t
 
-    for layer_cache1, layer_cache2 in zip(cache1, cache2):
-        next_cache = cat((layer_cache1, layer_cache2), dim = -2)
-        combined_cache.append(next_cache)
+    dims_from_right = (- dim - 1) if dim < 0 else (t.ndim - dim - 1)
+    zeros = ((0, 0) * dims_from_right)
+    return F.pad(t, (*zeros, *pad), value = value)
 
-    return combined_cache
+def calc_entropy(logits):
+    prob = logits.softmax(dim = -1)
+    return -(prob * log(prob)).sum(dim = -1)
 
 # generalized advantage estimate
 
@@ -52,7 +77,7 @@ def combine_kv_cache(cache1, cache2):
 def calc_gae(
     rewards,
     values,
-    masks,
+    masks = None,
     gamma = 0.99,
     lam = 0.95,
     use_accelerated = None
@@ -63,6 +88,9 @@ def calc_gae(
     values = F.pad(values, (0, 1), value = 0.)
     values, values_next = values[..., :-1], values[..., 1:]
 
+    if not exists(masks):
+        masks = torch.ones_like(values)
+
     delta = rewards + gamma * values_next * masks - values
     gates = gamma * lam * masks
 
@@ -114,8 +142,8 @@ def create_xl_mask(
         # handle intra-episodic attention if needed
 
         if exists(episode_ids):
-            q_episode = episodes[b, q + offset]
-            k_episode = episodes[b, k]
+            q_episode = episode_ids[b, q + offset]
+            k_episode = episode_ids[b, k]
 
             intra_episode_mask = q_episode == k_episode
             mask = mask & intra_episode_mask
@@ -146,6 +174,217 @@ def create_sliding_mask(
     create_kwargs = dict(device = device) if exists(device) else dict()
     return create_block_mask(sliding_mask, B = None, H = None, Q_LEN = seq_len, KV_LEN = kv_seq_len, _compile = True, **create_kwargs)
 
+# data
+
+def collate_var_time(data):
+
+    datum = first(data)
+    keys = datum.keys()
+
+    all_tensors = zip(*[datum.values() for datum in data])
+
+    collated_values = []
+
+    for key, tensors in zip(keys, all_tensors):
+
+        # the episode lens have zero dimension - think of a cleaner way to handle this later
+
+        if key != '_lens':
+
+            times = [t.shape[0] for t in tensors]
+            max_time = max(times)
+            tensors = [pad_at_dim(t, (0, max_time - t.shape[0]), dim = 0) for t in tensors]
+
+        collated_values.append(stack(tensors))
+
+    return dict(zip(keys, collated_values))
+
+class ReplayDataset(Dataset):
+    def __init__(
+        self,
+        folder: str | Path,
+        fields: tuple[str, ...] | None = None
+    ):
+        if isinstance(folder, str):
+            folder = Path(folder)
+
+        episode_lens = folder / 'episode_lens.npy'
+        self.episode_lens = open_memmap(str(episode_lens), mode = 'r')
+
+        # get indices of non-zero lengthed episodes
+
+        nonzero_episodes = self.episode_lens > 0
+        self.indices = np.arange(self.episode_lens.shape[-1])[nonzero_episodes]
+
+        # get all data files
+
+        filepaths = [*folder.glob('*.data.npy')]
+        assert len(filepaths) > 0
+
+        fieldname_to_filepath = {path.name.split('.')[0]: path for path in filepaths}
+
+        fieldnames_from_files = set(fieldname_to_filepath.keys())
+
+        fields = default(fields, fieldnames_from_files)
+
+        self.memmaps = dict()
+
+        for field in fields:
+            assert field in fieldnames_from_files, f'invalid field {field} - must be one of {fieldnames_from_files}'
+
+            path = fieldname_to_filepath[field]
+
+            self.memmaps[field] = open_memmap(str(path), mode = 'r')
+
+    def __len__(self):
+        return len(self.indices)
+
+    def __getitem__(self, idx):
+        episode_index = self.indices[idx]
+
+        episode_len = self.episode_lens[episode_index]
+
+        data = {field: from_numpy(memmap[episode_index, :episode_len].copy()) for field, memmap in self.memmaps.items()}
+
+        data['_lens'] = tensor(episode_len)
+
+        return data
+
+class ReplayBuffer:
+
+    @beartype
+    def __init__(
+        self,
+        folder: str | Path,
+        max_episodes: int,
+        max_timesteps: int,
+        fields: dict[
+            str,
+            str | tuple[str, int | tuple[int, ...]]
+        ]
+    ):
+
+        # folder for data
+
+        if not isinstance(folder, Path):
+            folder = Path(folder)
+            folder.mkdir(exist_ok = True)
+
+        self.folder = folder
+        assert folder.is_dir()
+
+        # keeping track of episode length
+
+        episode_lens = folder / 'episode_lens.npy'
+
+        self.episode_index = 0
+        self.timestep_index = 0
+
+        self.max_episodes = max_episodes
+        self.max_timesteps= max_timesteps
+
+        self.episode_lens = open_memmap(str(episode_lens), mode = 'w+', dtype = np.int32, shape = (max_episodes,))
+
+        # create the memmap for individual data tracks
+
+        self.shapes = dict()
+        self.dtypes = dict()
+        self.memmaps = dict()
+        self.fieldnames = set(fields.keys())
+
+        for field_name, field_info in fields.items():
+
+            # some flexibility
+
+            field_info = (field_info, ()) if isinstance(field_info, str) else field_info
+
+            dtype_str, shape = field_info
+            assert dtype_str in {'int', 'float', 'bool'}
+
+            dtype = dict(int = np.int32, float = np.float32, bool = np.bool_)[dtype_str]
+
+            # memmap file
+
+            filepath = folder / f'{field_name}.data.npy'
+            memmap = open_memmap(str(filepath), mode = 'w+', dtype = dtype, shape = (max_episodes, max_timesteps, *shape))
+
+            self.memmaps[field_name] = memmap
+            self.shapes[field_name] = shape
+            self.dtypes[field_name] = dtype
+
+        self.memory_namedtuple = namedtuple('Memory', list(fields.keys()))
+
+    def reset_(self):
+        self.episode_lens[:] = 0
+        self.episode_index = 0
+        self.timestep_index = 0
+
+    def advance_episode(self):
+        self.episode_index = (self.episode_index + 1) % self.max_episodes
+        self.timestep_index = 0
+
+    def flush(self):
+        self.episode_lens[self.episode_index] = self.timestep_index
+
+        for memmap in self.memmaps.values():
+            memmap.flush()
+
+        self.episode_lens.flush()
+
+    @contextmanager
+    def one_episode(self):
+
+        yield
+
+        self.flush()
+        self.advance_episode()
+
+    @beartype
+    def store_datapoint(
+        self,
+        episode_index: int,
+        timestep_index: int,
+        name: str,
+        datapoint: Tensor | ndarray
+    ):
+        assert 0 <= episode_index < self.max_episodes
+        assert 0 <= timestep_index < self.max_timesteps
+
+        if is_tensor(datapoint):
+            datapoint = datapoint.detach().cpu().numpy()
+
+        assert name in self.fieldnames, f'invalid field name {name} - must be one of {self.fieldnames}'
+
+        assert datapoint.shape == self.shapes[name], f'invalid shape {datapoint.shape} - shape must be {self.shapes[name]}'
+
+        self.memmaps[name][self.episode_index, self.timestep_index] = datapoint
+
+    def store(
+        self,
+        **data
+    ):
+        assert is_bearable(data, dict[str, Tensor | ndarray])
+
+        assert not self.timestep_index >= self.max_timesteps, 'you exceeded the `max_timesteps` set on the replay buffer'
+
+        for name, datapoint in data.items():
+
+            self.store_datapoint(self.episode_index, self.timestep_index, name, datapoint)
+
+        self.timestep_index += 1
+
+        return self.memory_namedtuple(**data)
+
+    def dataset(self) -> Dataset:
+        self.flush()
+
+        return ReplayDataset(self.folder)
+
+    def dataloader(self, batch_size, **kwargs) -> DataLoader:
+        self.flush()
+
+        return DataLoader(self.dataset(), batch_size = batch_size, collate_fn = collate_var_time, **kwargs)
+
 # transformer-xl with ppo
 
 class Attention(Module):
@@ -204,7 +443,6 @@ class Attention(Module):
         return_kv_cache = False,
     ):
         seq_len = tokens.shape[-2]
-        assert seq_len <= self.window_size
 
         device = tokens.device
 
@@ -365,7 +603,13 @@ class Locoformer(Module):
         embedder: Module,
         unembedder: Module,
         transformer: dict | TransformerXL,
-        value_network: Module | None = None
+        value_network: Module | None = None,
+        discount_factor = 0.999,
+        gae_lam = 0.95,
+        ppo_eps_clip = 0.2,
+        ppo_entropy_weight = 0.01,
+        ppo_value_clip = 0.4,
+        value_loss_weight = 0.5
     ):
         super().__init__()
 
@@ -382,6 +626,15 @@ class Locoformer(Module):
         self.fixed_window_size = transformer.fixed_window_size
         self.window_size = transformer.window_size
 
+        # ppo related
+
+        self.discount_factor = discount_factor
+        self.gae_lam = gae_lam
+        self.ppo_eps_clip = ppo_eps_clip
+        self.ppo_entropy_weight = ppo_entropy_weight
+        self.ppo_value_clip = ppo_value_clip
+        self.value_loss_weight = value_loss_weight
+
     @property
     def device(self):
         return next(self.parameters()).device
@@ -390,8 +643,95 @@ class Locoformer(Module):
         return self.unembedder.parameters()
 
     def critic_parameters(self):
+        if not exists(self.value_network):
+            return []
+
         return self.value_network.parameters()
 
+    def ppo(
+        self,
+        state,
+        action,
+        old_action_log_prob,
+        reward,
+        old_value,
+        mask,
+        actor_optim,
+        critic_optim
+    ):
+
+        (action_logits, value), _ = self.forward(state, return_values = True)
+        entropy = calc_entropy(action_logits)
+
+        action = rearrange(action, 'b t -> b t 1')
+        log_prob = action_logits.gather(-1, action)
+        log_prob = rearrange(log_prob, 'b t 1 -> b t')
+
+        # update actor, classic clipped surrogate loss
+
+        eps_clip = self.ppo_eps_clip
+        ratio = (log_prob - old_action_log_prob).exp()
+
+        returns = calc_gae(reward, old_value, lam = self.gae_lam, gamma = self.discount_factor)
+        advantage = returns - old_value
+
+        actor_loss = -torch.min(ratio * advantage, ratio.clamp(1. - eps_clip, 1. + eps_clip) * advantage)
+
+        actor_loss = actor_loss - self.ppo_entropy_weight * entropy
+
+        mean_actor_loss = actor_loss[mask].mean()
+        mean_actor_loss.backward(retain_graph = True)
+
+        # update critic
+
+        value_loss = F.mse_loss(returns, value, reduction = 'none')
+
+        value_clip = self.ppo_value_clip
+        clipped_value = old_value + (value - old_value).clamp(-value_clip, value_clip)
+        clipped_value_loss = F.mse_loss(returns, clipped_value, reduction = 'none')
+
+        critic_loss = torch.maximum(value_loss, clipped_value_loss) * self.value_loss_weight
+
+        mean_critic_loss = critic_loss[mask].mean()
+        mean_critic_loss.backward()
+
+        # optimizer update
+
+        actor_optim.step()
+        actor_optim.zero_grad()
+
+        critic_optim.step()
+        critic_optim.zero_grad()
+
+        # return losses for logging
+
+        return mean_actor_loss.detach(), mean_critic_loss.detach()
+
+    def wrap_env_functions(self, env):
+
+        def wrapped_reset(*args, **kwargs):
+            state, _ =  env.reset(*args, **kwargs)
+
+            if isinstance(state, ndarray):
+                state = from_numpy(state)
+
+            return state, _
+
+        def wrapped_step(action, *args, **kwargs):
+            out = env.step(action.item(), *args, **kwargs)
+
+            def transform_output(el):
+                if isinstance(el, ndarray):
+                    return from_numpy(el)
+                elif isinstance(el, (int, bool, float)):
+                    return tensor(el)
+                else:
+                    return el
+
+            return tree_map(transform_output, out)
+
+        return wrapped_reset, wrapped_step
+
     def get_stateful_forward(
         self,
         initial_states: Tensor | None = None,
@@ -463,6 +803,8 @@ class Locoformer(Module):
         return_values = False
     ):
 
+        state = state.to(self.device)
+
         tokens = self.embedder(state)
 
         embed, kv_cache = self.transformer(tokens, cache = cache, return_kv_cache = True)

{locoformer-0.0.7.dist-info → locoformer-0.0.15.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: locoformer
-Version: 0.0.7
+Version: 0.0.15
 Summary: LocoFormer
 Project-URL: Homepage, https://pypi.org/project/locoformer/
 Project-URL: Repository, https://github.com/lucidrains/locoformer
@@ -35,6 +35,7 @@ Classifier: Programming Language :: Python :: 3.9
 Classifier: Topic :: Scientific/Engineering :: Artificial Intelligence
 Requires-Python: >=3.9
 Requires-Dist: assoc-scan
+Requires-Dist: beartype
 Requires-Dist: einops>=0.8.0
 Requires-Dist: einx>=0.3.0
 Requires-Dist: rotary-embedding-torch

locoformer-0.0.15.dist-info/RECORD

@@ -0,0 +1,6 @@
+locoformer/__init__.py,sha256=XctsMGEZSR4mVl75fhds_1BtS5qGFiiItTDV7CmCt_I,45
+locoformer/locoformer.py,sha256=1jPK41G4HB1PEPtlusQxcrne489E-3QKXAULZ20FEZM,22740
+locoformer-0.0.15.dist-info/METADATA,sha256=IHtK7NvVQewYQ0GBB7v1KG90_H2Jakxir0MakUIA-jU,3218
+locoformer-0.0.15.dist-info/WHEEL,sha256=qtCwoSJWgHk21S1Kb4ihdzI2rlJ1ZKaIurTj_ngOhyQ,87
+locoformer-0.0.15.dist-info/licenses/LICENSE,sha256=1yCiA9b5nhslTavxPjsQAO-wpOnwJR9-l8LTVi7GJuk,1066
+locoformer-0.0.15.dist-info/RECORD,,

locoformer-0.0.7.dist-info/RECORD

@@ -1,6 +0,0 @@
-locoformer/__init__.py,sha256=XctsMGEZSR4mVl75fhds_1BtS5qGFiiItTDV7CmCt_I,45
-locoformer/locoformer.py,sha256=lJQs0CKr9iztF8tie1FRUVEItCt-IZbIILQqKcgK2sI,13142
-locoformer-0.0.7.dist-info/METADATA,sha256=PZ_phKV3t4Bha0GnUB5HPmE9w8A5fvNevsuN532Ls3s,3193
-locoformer-0.0.7.dist-info/WHEEL,sha256=qtCwoSJWgHk21S1Kb4ihdzI2rlJ1ZKaIurTj_ngOhyQ,87
-locoformer-0.0.7.dist-info/licenses/LICENSE,sha256=1yCiA9b5nhslTavxPjsQAO-wpOnwJR9-l8LTVi7GJuk,1066
-locoformer-0.0.7.dist-info/RECORD,,