cccode / models.py
WayneW's picture
Upload folder using huggingface_hub
705a8fd verified
raw
history blame
21.9 kB
# Copyright (c) Meta Platforms, Inc. and affiliates.
# All rights reserved.
# This source code is licensed under the license found in the
# LICENSE file in the root directory of this source tree.
# --------------------------------------------------------
# References:
# GLIDE: https://github.com/openai/glide-text2im
# MAE: https://github.com/facebookresearch/mae/blob/main/models_mae.py
# --------------------------------------------------------
import torch
import torch.nn as nn
import numpy as np
import math
from timm.models.vision_transformer import PatchEmbed, Attention, Mlp
def modulate(x, shift, scale):
return x * (1 + scale.unsqueeze(1)) + shift.unsqueeze(1)
#################################################################################
# Embedding Layers for Timesteps and Class Labels #
#################################################################################
class TimestepEmbedder(nn.Module):
"""
Embeds scalar timesteps into vector representations.
"""
def __init__(self, hidden_size, frequency_embedding_size=256):
super().__init__()
self.mlp = nn.Sequential(
nn.Linear(frequency_embedding_size, hidden_size, bias=True),
nn.SiLU(),
nn.Linear(hidden_size, hidden_size, bias=True),
)
self.frequency_embedding_size = frequency_embedding_size
@staticmethod
def timestep_embedding(t, dim, max_period=10000):
"""
Create sinusoidal timestep embeddings.
:param t: a 1-D Tensor of N indices, one per batch element.
These may be fractional.
:param dim: the dimension of the output.
:param max_period: controls the minimum frequency of the embeddings.
:return: an (N, D) Tensor of positional embeddings.
"""
# https://github.com/openai/glide-text2im/blob/main/glide_text2im/nn.py
half = dim // 2
freqs = torch.exp(
-math.log(max_period) * torch.arange(start=0, end=half, dtype=torch.float32) / half
).to(device=t.device)
args = t.float() * freqs[None]
embedding = torch.cat([torch.cos(args), torch.sin(args)], dim=-1)
if dim % 2:
embedding = torch.cat([embedding, torch.zeros_like(embedding[:, :1])], dim=-1)
return embedding
def forward(self, t):
t_freq = self.timestep_embedding(t, self.frequency_embedding_size)
t_emb = self.mlp(t_freq)
return t_emb
class ActionEmbedder(nn.Module):
"""
Embeds action xy into vector representations.
"""
def __init__(self, hidden_size, frequency_embedding_size=256):
super().__init__()
hsize = hidden_size//3
self.x_emb = TimestepEmbedder(hsize, frequency_embedding_size)
self.y_emb = TimestepEmbedder(hsize, frequency_embedding_size)
self.angle_emb = TimestepEmbedder(hidden_size -2*hsize, frequency_embedding_size)
def forward(self, xya):
return torch.cat([self.x_emb(xya[...,0:1]), self.y_emb(xya[...,1:2]), self.angle_emb(xya[...,2:3])], dim=-1)
#################################################################################
# Core AVCDiT Model #
#################################################################################
class AVCDiTBlock(nn.Module):
"""
A DiT block with adaptive layer norm zero (adaLN-Zero) conditioning and two modalities.
"""
def __init__(self, hidden_size, num_heads, mlp_ratio=4.0, mode="av", **block_kwargs):
super().__init__()
self.mode = mode
self.norm1 = nn.LayerNorm(hidden_size, elementwise_affine=False, eps=1e-6)
self.attn = Attention(hidden_size, num_heads=num_heads, qkv_bias=True, **block_kwargs)
self.norm2 = nn.LayerNorm(hidden_size, elementwise_affine=False, eps=1e-6)
self.norm_cond = nn.LayerNorm(hidden_size, elementwise_affine=False, eps=1e-6)
self.cttn = nn.MultiheadAttention(hidden_size, num_heads=num_heads, add_bias_kv=True, bias=True, batch_first=True, **block_kwargs)
self.adaLN_modulation = nn.Sequential(
nn.SiLU(),
nn.Linear(hidden_size, 11 * hidden_size, bias=True)
)
self.norm3 = nn.LayerNorm(hidden_size, elementwise_affine=False, eps=1e-6)
mlp_hidden_dim = int(hidden_size * mlp_ratio)
approx_gelu = lambda: nn.GELU(approximate="tanh")
if self.mode == "av" or self.mode == "v":
self.mlp_v = Mlp(in_features=hidden_size, hidden_features=mlp_hidden_dim, act_layer=approx_gelu, drop=0)
if self.mode == "av" or self.mode == "a":
self.mlp_a = Mlp(in_features=hidden_size, hidden_features=mlp_hidden_dim, act_layer=approx_gelu, drop=0)
# def forward(self, x_v, x_a, c, x_v_cond, x_a_cond, mode="av"):
def forward(self, *args):
if self.mode == "av":
x_v, x_a, c, x_v_cond, x_a_cond = args
shift_msa, scale_msa, gate_msa, shift_ca_xcond, scale_ca_xcond, shift_ca_x, scale_ca_x, gate_ca_x, shift_mlp, scale_mlp, gate_mlp = self.adaLN_modulation(c).chunk(11, dim=1)
_, v_token_num, _ = x_v.shape
x = torch.cat([x_v, x_a], dim=1)
x_cond = torch.cat([x_v_cond, x_a_cond], dim=1)
x = x + gate_msa.unsqueeze(1) * self.attn(modulate(self.norm1(x), shift_msa, scale_msa))
x_cond_norm = modulate(self.norm_cond(x_cond), shift_ca_xcond, scale_ca_xcond)
x = x + gate_ca_x.unsqueeze(1) * self.cttn(query=modulate(self.norm2(x), shift_ca_x, scale_ca_x), key=x_cond_norm, value=x_cond_norm, need_weights=False)[0]
x_v = x[:,:v_token_num,:]
x_a = x[:,v_token_num:,:]
x_v = x_v + gate_mlp.unsqueeze(1) * self.mlp_v(modulate(self.norm3(x_v), shift_mlp, scale_mlp))
x_a = x_a + gate_mlp.unsqueeze(1) * self.mlp_a(modulate(self.norm3(x_a), shift_mlp, scale_mlp))
return x_v, x_a
elif self.mode == "v":
x, c, x_cond = args
shift_msa, scale_msa, gate_msa, shift_ca_xcond, scale_ca_xcond, shift_ca_x, scale_ca_x, gate_ca_x, shift_mlp, scale_mlp, gate_mlp = self.adaLN_modulation(c).chunk(11, dim=1)
x = x + gate_msa.unsqueeze(1) * self.attn(modulate(self.norm1(x), shift_msa, scale_msa))
x_cond_norm = modulate(self.norm_cond(x_cond), shift_ca_xcond, scale_ca_xcond)
x = x + gate_ca_x.unsqueeze(1) * self.cttn(query=modulate(self.norm2(x), shift_ca_x, scale_ca_x), key=x_cond_norm, value=x_cond_norm, need_weights=False)[0]
x = x + gate_mlp.unsqueeze(1) * self.mlp_v(modulate(self.norm3(x), shift_mlp, scale_mlp))
return x
elif self.mode == "a":
x, c, x_cond = args
shift_msa, scale_msa, gate_msa, shift_ca_xcond, scale_ca_xcond, shift_ca_x, scale_ca_x, gate_ca_x, shift_mlp, scale_mlp, gate_mlp = self.adaLN_modulation(c).chunk(11, dim=1)
x = x + gate_msa.unsqueeze(1) * self.attn(modulate(self.norm1(x), shift_msa, scale_msa))
x_cond_norm = modulate(self.norm_cond(x_cond), shift_ca_xcond, scale_ca_xcond)
x = x + gate_ca_x.unsqueeze(1) * self.cttn(query=modulate(self.norm2(x), shift_ca_x, scale_ca_x), key=x_cond_norm, value=x_cond_norm, need_weights=False)[0]
x = x + gate_mlp.unsqueeze(1) * self.mlp_a(modulate(self.norm3(x), shift_mlp, scale_mlp))
return x
class FinalLayer(nn.Module):
"""
The final layer of DiT.
"""
def __init__(self, hidden_size, patch_size, out_channels):
super().__init__()
self.norm_final = nn.LayerNorm(hidden_size, elementwise_affine=False, eps=1e-6)
self.linear = nn.Linear(hidden_size, patch_size * patch_size * out_channels, bias=True)
self.adaLN_modulation = nn.Sequential(
nn.SiLU(),
nn.Linear(hidden_size, 2 * hidden_size, bias=True)
)
def forward(self, x, c):
shift, scale = self.adaLN_modulation(c).chunk(2, dim=1)
x = modulate(self.norm_final(x), shift, scale)
x = self.linear(x)
return x
class FinalLayer_audio(nn.Module):
def __init__(self, hidden_size, out_channels):
super().__init__()
self.norm_final = nn.LayerNorm(hidden_size, elementwise_affine=False, eps=1e-6)
self.linear = nn.Linear(hidden_size, out_channels, bias=True) # no patch²
self.adaLN_modulation = nn.Sequential(
nn.SiLU(),
nn.Linear(hidden_size, 2 * hidden_size, bias=True)
)
def forward(self, x, c):
# x: (B, N, hidden_size), c: (B, hidden_size)
shift, scale = self.adaLN_modulation(c).chunk(2, dim=1) # shape (B, hidden_size)
x = modulate(self.norm_final(x), shift, scale) # apply AdaLN
x = self.linear(x) # → (B, N, out_channels)
return x
class AVCDiT(nn.Module):
"""
Diffusion model with a Transformer backbone.
"""
def __init__(
self,
input_size=32,
context_size=2,
patch_size=2,
in_channels=4,
hidden_size=1152,
depth=28,
num_heads=16,
mlp_ratio=4.0,
learn_sigma=True,
num_patches_a=180,
mode="av",
):
super().__init__()
self.mode = mode
assert (self.mode=="av" or self.mode=="v" or self.mode=="a")
self.context_size = context_size
self.learn_sigma = learn_sigma
self.in_channels = in_channels
self.out_channels = in_channels * 2 if learn_sigma else in_channels
self.patch_size = patch_size
self.num_heads = num_heads
if self.mode == "av" or self.mode == "v":
self.x_embedder_v = PatchEmbed(input_size, patch_size, in_channels, hidden_size, bias=True)
num_patches_v = self.x_embedder_v.num_patches
self.pos_embed_v = nn.Parameter(torch.zeros(self.context_size + 1, num_patches_v, hidden_size), requires_grad=True) # for context and for predicted frame
self.final_layer = FinalLayer(hidden_size, patch_size, self.out_channels)
if self.mode == "av" or self.mode == "a":
self.x_embedder_a = nn.Conv1d(
in_channels=16,
out_channels=hidden_size, # [B]
kernel_size=1,
stride=1,
bias=True
) #TODO
self.pos_embed_a_cond = nn.Parameter(torch.zeros(self.context_size, num_patches_a, hidden_size), requires_grad=True)
self.pos_embed_a_pred = nn.Parameter(torch.zeros(1, num_patches_a+1, hidden_size), requires_grad=True)
self.final_layer_a = FinalLayer_audio(hidden_size=hidden_size, out_channels=32) # [B]
self.t_embedder = TimestepEmbedder(hidden_size)
self.y_embedder = ActionEmbedder(hidden_size)
# self.blocks = nn.ModuleList([AVCDiTBlock(hidden_size, num_heads, mlp_ratio=mlp_ratio) for _ in range(depth)])
self.blocks = nn.ModuleList([
AVCDiTBlock(hidden_size, num_heads, mlp_ratio=mlp_ratio, mode=self.mode)
for _ in range(depth)
])
self.time_embedder = TimestepEmbedder(hidden_size)
self.initialize_weights()
def initialize_weights(self):
# Initialize transformer layers:
def _basic_init(module):
if isinstance(module, nn.Linear):
torch.nn.init.xavier_uniform_(module.weight)
if module.bias is not None:
nn.init.constant_(module.bias, 0)
self.apply(_basic_init)
# Initialize (and freeze) pos_embed by sin-cos embedding:
if self.mode == "av" or self.mode == "v":
nn.init.normal_(self.pos_embed_v, std=0.02)
if self.mode == "av" or self.mode == "a":
nn.init.normal_(self.pos_embed_a_pred, std=0.02)
nn.init.normal_(self.pos_embed_a_cond, std=0.02)
# Initialize patch_embed like nn.Linear (instead of nn.Conv2d):
if self.mode == "av" or self.mode == "v":
w = self.x_embedder_v.proj.weight.data
nn.init.xavier_uniform_(w.view([w.shape[0], -1]))
nn.init.constant_(self.x_embedder_v.proj.bias, 0)
# Initialize x_embedder_a (Conv1d) like linear
if self.mode == "av" or self.mode == "a":
w = self.x_embedder_a.weight.data
nn.init.xavier_uniform_(w.view([w.shape[0], -1]))
nn.init.constant_(self.x_embedder_a.bias, 0)
# Initialize action embedding:
nn.init.normal_(self.y_embedder.x_emb.mlp[0].weight, std=0.02)
nn.init.normal_(self.y_embedder.x_emb.mlp[2].weight, std=0.02)
nn.init.normal_(self.y_embedder.y_emb.mlp[0].weight, std=0.02)
nn.init.normal_(self.y_embedder.y_emb.mlp[2].weight, std=0.02)
nn.init.normal_(self.y_embedder.angle_emb.mlp[0].weight, std=0.02)
nn.init.normal_(self.y_embedder.angle_emb.mlp[2].weight, std=0.02)
# Initialize timestep embedding MLP:
nn.init.normal_(self.t_embedder.mlp[0].weight, std=0.02)
nn.init.normal_(self.t_embedder.mlp[2].weight, std=0.02)
nn.init.normal_(self.time_embedder.mlp[0].weight, std=0.02)
nn.init.normal_(self.time_embedder.mlp[2].weight, std=0.02)
# Zero-out adaLN modulation layers in DiT blocks:
for block in self.blocks:
nn.init.constant_(block.adaLN_modulation[-1].weight, 0)
nn.init.constant_(block.adaLN_modulation[-1].bias, 0)
# Zero-out output layers:
if self.mode == "av" or self.mode == "v":
nn.init.constant_(self.final_layer.adaLN_modulation[-1].weight, 0)
nn.init.constant_(self.final_layer.adaLN_modulation[-1].bias, 0)
nn.init.constant_(self.final_layer.linear.weight, 0)
nn.init.constant_(self.final_layer.linear.bias, 0)
if self.mode == "av" or self.mode == "a":
nn.init.constant_(self.final_layer_a.adaLN_modulation[-1].weight, 0)
nn.init.constant_(self.final_layer_a.adaLN_modulation[-1].bias, 0)
nn.init.constant_(self.final_layer_a.linear.weight, 0)
nn.init.constant_(self.final_layer_a.linear.bias, 0)
def unpatchify(self, x):
"""
x: (N, T, patch_size**2 * C)
imgs: (N, H, W, C)
"""
c = self.out_channels
p = self.x_embedder_v.patch_size[0]
h = w = int(x.shape[1] ** 0.5)
assert h * w == x.shape[1]
x = x.reshape(shape=(x.shape[0], h, w, p, p, c))
x = torch.einsum('nhwpqc->nchpwq', x)
imgs = x.reshape(shape=(x.shape[0], c, h * p, h * p))
return imgs
# def forward(self, x_v, x_a, t, y, x_v_cond, x_a_cond, rel_t):
# def forward(self, *args):
def forward(self, *args, **kwargs):
"""
Forward pass of DiT.
x: (N, C, H, W) tensor of spatial inputs (images or latent representations of images)
t: (N,) tensor of diffusion timesteps
y: (N,) tensor of class labels
"""
if self.mode == "av":
if len(args) >= 7:
x_v, x_a, t, y, x_v_cond, x_a_cond, rel_t = args[:7]
else:
assert len(args) == 3, f"mode='v' expects 2 or 5 positional args, got {len(args)}"
x_v, x_a, t = args
y = kwargs["y"]
x_v_cond = kwargs["x_v_cond"]
x_a_cond = kwargs["x_a_cond"]
rel_t = kwargs["rel_t"]
x_v = self.x_embedder_v(x_v) + self.pos_embed_v[self.context_size:]
x_v_cond = self.x_embedder_v(x_v_cond.flatten(0, 1)).unflatten(0, (x_v_cond.shape[0], x_v_cond.shape[1])) + self.pos_embed_v[:self.context_size] # (N, T, D), where T = H * W / patch_size ** 2.flatten(1, 2)
x_v_cond = x_v_cond.flatten(1, 2)
x_a = self.x_embedder_a(x_a) # → (B, embed_dim, L')
x_a = x_a.transpose(1, 2) # → (B, L', embed_dim)
x_a = x_a + self.pos_embed_a_pred
x_a_cond = self.x_embedder_a(x_a_cond.flatten(0, 1)).transpose(1, 2).unflatten(0, (x_a_cond.shape[0], x_a_cond.shape[1])) + self.pos_embed_a_cond
x_a_cond = x_a_cond.flatten(1, 2)
t = self.t_embedder(t[..., None])
y = self.y_embedder(y)
time_emb = self.time_embedder(rel_t[..., None])
c = t + time_emb + y # if training on unlabeled data, dont add y.
for block in self.blocks:
x_v, x_a = block(x_v, x_a, c, x_v_cond, x_a_cond)
x_v = self.final_layer(x_v, c)
x_v = self.unpatchify(x_v)
x_a = self.final_layer_a(x_a, c)
x_a = x_a.transpose(1, 2)
return x_v, x_a
elif self.mode == "v":
if len(args) >= 5:
x, t, y, x_cond, rel_t = args[:5]
else:
assert len(args) == 2, f"mode='v' expects 2 or 5 positional args, got {len(args)}"
x, t = args
y = kwargs["y"]
x_cond = kwargs["x_cond"]
rel_t = kwargs["rel_t"]
x = self.x_embedder_v(x) + self.pos_embed_v[self.context_size:]
x_cond = self.x_embedder_v(x_cond.flatten(0, 1)).unflatten(0, (x_cond.shape[0], x_cond.shape[1])) + self.pos_embed_v[:self.context_size] # (N, T, D), where T = H * W / patch_size ** 2.flatten(1, 2)
x_cond = x_cond.flatten(1, 2)
t = self.t_embedder(t[..., None])
y = self.y_embedder(y)
time_emb = self.time_embedder(rel_t[..., None])
c = t + time_emb + y # if training on unlabeled data, dont add y.
for block in self.blocks:
x = block(x, c, x_cond)
x = self.final_layer(x, c)
x = self.unpatchify(x)
return x
elif self.mode == "a":
if len(args) >= 5:
x, t, y, x_cond, rel_t = args[:5]
else:
assert len(args) == 2, f"mode='v' expects 2 or 5 positional args, got {len(args)}"
x, t = args
y = kwargs["y"]
x_cond = kwargs["x_cond"]
rel_t = kwargs["rel_t"]
x = self.x_embedder_a(x) # → (B, embed_dim, L')
x = x.transpose(1, 2) # → (B, L', embed_dim)
x = x + self.pos_embed_a_pred # [REWARD]
x_cond = self.x_embedder_a(x_cond.flatten(0, 1)).transpose(1, 2).unflatten(0, (x_cond.shape[0], x_cond.shape[1])) + self.pos_embed_a_cond # [REWARD]
x_cond = x_cond.flatten(1, 2)
t = self.t_embedder(t[..., None])
y = self.y_embedder(y)
time_emb = self.time_embedder(rel_t[..., None])
c = t + time_emb + y # if training on unlabeled data, dont add y.
for block in self.blocks:
x = block(x, c, x_cond)
x = self.final_layer_a(x, c)
x = x.transpose(1, 2)
return x
#################################################################################
# Sine/Cosine Positional Embedding Functions #
#################################################################################
# https://github.com/facebookresearch/mae/blob/main/util/pos_embed.py
def get_2d_sincos_pos_embed(embed_dim, grid_size, cls_token=False, extra_tokens=0):
"""
grid_size: int of the grid height and width
return:
pos_embed: [grid_size*grid_size, embed_dim] or [1+grid_size*grid_size, embed_dim] (w/ or w/o cls_token)
"""
grid_h = np.arange(grid_size, dtype=np.float32)
grid_w = np.arange(grid_size, dtype=np.float32)
grid = np.meshgrid(grid_w, grid_h) # here w goes first
grid = np.stack(grid, axis=0)
grid = grid.reshape([2, 1, grid_size, grid_size])
pos_embed = get_2d_sincos_pos_embed_from_grid(embed_dim, grid)
if cls_token and extra_tokens > 0:
pos_embed = np.concatenate([np.zeros([extra_tokens, embed_dim]), pos_embed], axis=0)
return pos_embed
def get_2d_sincos_pos_embed_from_grid(embed_dim, grid):
assert embed_dim % 2 == 0
# use half of dimensions to encode grid_h
emb_h = get_1d_sincos_pos_embed_from_grid(embed_dim // 2, grid[0]) # (H*W, D/2)
emb_w = get_1d_sincos_pos_embed_from_grid(embed_dim // 2, grid[1]) # (H*W, D/2)
emb = np.concatenate([emb_h, emb_w], axis=1) # (H*W, D)
return emb
def get_1d_sincos_pos_embed_from_grid(embed_dim, pos):
"""
embed_dim: output dimension for each position
pos: a list of positions to be encoded: size (M,)
out: (M, D)
"""
assert embed_dim % 2 == 0
omega = np.arange(embed_dim // 2, dtype=np.float64)
omega /= embed_dim / 2.
omega = 1. / 10000**omega # (D/2,)
pos = pos.reshape(-1) # (M,)
out = np.einsum('m,d->md', pos, omega) # (M, D/2), outer product
emb_sin = np.sin(out) # (M, D/2)
emb_cos = np.cos(out) # (M, D/2)
emb = np.concatenate([emb_sin, emb_cos], axis=1) # (M, D)
return emb
#################################################################################
# AVCDiT Configs #
#################################################################################
def AVCDiT_XL_2(**kwargs):
return AVCDiT(depth=28, hidden_size=1152, patch_size=2, num_heads=16, **kwargs)
def AVCDiT_L_2(**kwargs):
return AVCDiT(depth=24, hidden_size=1024, patch_size=2, num_heads=16, **kwargs)
def AVCDiT_B_2(**kwargs):
return AVCDiT(depth=12, hidden_size=768, patch_size=2, num_heads=12, **kwargs)
def AVCDiT_S_2(**kwargs):
return AVCDiT(depth=12, hidden_size=384, patch_size=2, num_heads=6, **kwargs)
AVCDiT_models = {
'AVCDiT-XL/2': AVCDiT_XL_2,
'AVCDiT-L/2': AVCDiT_L_2,
'AVCDiT-B/2': AVCDiT_B_2,
'AVCDiT-S/2': AVCDiT_S_2
}