mirror of
https://github.com/ml-explore/mlx-examples.git
synced 2025-06-26 02:33:23 +08:00
340 lines
11 KiB
Python
340 lines
11 KiB
Python
import math
|
|
from dataclasses import dataclass, field
|
|
from typing import Optional, Tuple, Union
|
|
|
|
import mlx.core as mx
|
|
import mlx.nn as nn
|
|
|
|
from .base import BaseModelArgs
|
|
from .cache import Mamba2Cache
|
|
|
|
@dataclass
|
|
class ModelArgs(BaseModelArgs):
|
|
num_heads: int
|
|
head_dim: int
|
|
vocab_size: int
|
|
hidden_size: int
|
|
state_size: int
|
|
num_hidden_layers: int
|
|
layer_norm_epsilon: float
|
|
expand: int
|
|
conv_kernel: int
|
|
n_groups: int
|
|
use_bias: bool
|
|
use_conv_bias: bool
|
|
initializer_range: float
|
|
residual_in_fp32: bool
|
|
time_step_min: float
|
|
time_step_max: float
|
|
time_step_floor: float
|
|
rescale_prenorm_residual: bool
|
|
use_cache: bool
|
|
rms_norm: bool
|
|
chunk_size: int
|
|
tie_word_embeddings: bool
|
|
time_step_limit: Tuple[float, float] = field(default_factory=lambda: (0.0, float("inf")))
|
|
time_step_rank: Union[int, str] = "auto"
|
|
model_type: str = "mamba2"
|
|
|
|
def __post_init__(self):
|
|
if not hasattr(self, "intermediate_size"):
|
|
self.intermediate_size = int(self.expand * self.hidden_size)
|
|
if not hasattr(self, "head_dim"):
|
|
self.head_dim = self.hidden_size // self.num_heads
|
|
if self.time_step_rank == "auto":
|
|
self.time_step_rank = math.ceil(self.hidden_size / 16)
|
|
|
|
|
|
def selective_scan(x, A, B, C, chunk_size):
|
|
"""
|
|
Selective scan implementation for training.
|
|
|
|
Arguments
|
|
x: (batch, seqlen, n_heads, d_head)
|
|
A: (batch, seqlen, n_heads)
|
|
B: (batch, seqlen, n_heads, d_state)
|
|
C: (batch, seqlen, n_heads, d_state)
|
|
|
|
Return
|
|
y: (batch, seqlen, n_heads, d_head)
|
|
"""
|
|
assert x.shape[1] % chunk_size == 0
|
|
|
|
# Reshape into chunks
|
|
def chunk_reshape(m):
|
|
shape = list(m.shape)
|
|
shape[1:2] = [shape[1] // chunk_size, chunk_size]
|
|
return m.reshape(shape)
|
|
|
|
x, A, B, C = map(chunk_reshape, (x, A, B, C))
|
|
A = mx.transpose(A, [0, 3, 1, 2])
|
|
|
|
# Compute cumulative sums
|
|
A_cumsum = mx.cumsum(A, axis=-1)
|
|
|
|
# Process chunks
|
|
L = mx.exp(selective_cumsum(A))
|
|
Y_diag = mx.einsum('bclhn,bcshn,bhcls,bcshp->bclhp', C, B, L, x)
|
|
|
|
decay_states = mx.exp(A_cumsum[..., -1:] - A_cumsum)
|
|
states = mx.einsum('bclhn,bhcl,bclhp->bchpn', B, decay_states, x)
|
|
|
|
initial_states = mx.zeros_like(states[:, :1])
|
|
states = mx.concatenate([initial_states, states], axis=1)
|
|
decay_chunk = mx.exp(selective_cumsum(mx.pad(A_cumsum[..., -1], ((0,0), (0,0), (1,0)))))
|
|
new_states = mx.einsum('bhzc,bchpn->bzhpn', decay_chunk, states)
|
|
states = new_states[:, :-1]
|
|
|
|
state_decay_out = mx.exp(A_cumsum)
|
|
Y_off = mx.einsum('bclhn,bchpn,bhcl->bclhp', C, states, state_decay_out)
|
|
|
|
Y = (Y_diag + Y_off).reshape((-1, x.shape[1] * chunk_size, *Y_diag.shape[-2:]))
|
|
return Y
|
|
|
|
def selective_cumsum(x: mx.array) -> mx.array:
|
|
"""Stable selective cumulative sum calculation."""
|
|
T = x.shape[-1]
|
|
x = mx.repeat(x[..., None], T, axis=-1)
|
|
mask = mx.tril(mx.ones((T, T)), k=-1)
|
|
x = x * mask
|
|
x_cumsum = mx.cumsum(x, axis=-2)
|
|
mask = mx.tril(mx.ones((T, T)), k=0)
|
|
return mx.where(mask, x_cumsum, float('-inf'))
|
|
|
|
|
|
class Mamba2Block(nn.Module):
|
|
def __init__(self, args: ModelArgs):
|
|
super().__init__()
|
|
self.args = args
|
|
|
|
# Internal cache state
|
|
self.conv_state = None
|
|
self.ssm_state = None
|
|
|
|
# Project input to get various components
|
|
d_in_proj = (2 * args.intermediate_size + 2 * self.args.n_groups * args.state_size + args.num_heads)
|
|
self.in_proj = nn.Linear(
|
|
args.hidden_size,
|
|
d_in_proj,
|
|
bias=args.use_bias
|
|
)
|
|
|
|
# Convolution layer
|
|
conv_dim = args.intermediate_size + 2 * self.args.n_groups * args.state_size
|
|
self.conv1d = nn.Conv1d(
|
|
in_channels=conv_dim,
|
|
out_channels=conv_dim,
|
|
kernel_size=args.conv_kernel,
|
|
groups=conv_dim,
|
|
padding=args.conv_kernel - 1,
|
|
bias=args.use_conv_bias
|
|
)
|
|
|
|
# SSM parameters
|
|
dt_init_floor = math.log(args.time_step_floor)
|
|
self.dt_bias = mx.zeros((args.num_heads,)) * args.initializer_range
|
|
self.A_log = mx.zeros((args.num_heads,)) * args.initializer_range
|
|
self.D = mx.zeros((args.num_heads,)) * args.initializer_range
|
|
|
|
# Output projections
|
|
self.norm = nn.RMSNorm(args.intermediate_size, eps=args.layer_norm_epsilon)
|
|
self.out_proj = nn.Linear(args.intermediate_size, args.hidden_size, bias=args.use_bias)
|
|
|
|
def __call__(self, x: mx.array, cache=None) -> mx.array:
|
|
return self.forward_training(x) if x.shape[1] > 1 else self.forward_inference(x, cache)
|
|
|
|
def forward_training(self, u: mx.array) -> mx.array:
|
|
# Reset cache during training
|
|
self.cache = None
|
|
|
|
# Input projection and splitting
|
|
zxbcdt = self.in_proj(u)
|
|
z, xBC, dt = mx.split(
|
|
zxbcdt,
|
|
[
|
|
self.args.intermediate_size,
|
|
self.args.intermediate_size + 2 * self.args.state_size
|
|
],
|
|
axis=-1
|
|
)
|
|
|
|
# Time step processing
|
|
dt = mx.clip(
|
|
nn.softplus(dt + self.dt_bias),
|
|
self.args.time_step_min,
|
|
self.args.time_step_max
|
|
)
|
|
|
|
# Convolution processing
|
|
xBC_t = mx.transpose(xBC, [0, 2, 1])
|
|
conv_out = self.conv1d(xBC_t)
|
|
xBC = mx.transpose(conv_out, [0, 2, 1])[:, :u.shape[1]]
|
|
xBC = mx.sigmoid(xBC) * xBC # SiLU
|
|
|
|
# Split states
|
|
x, B, C = mx.split(
|
|
xBC,
|
|
[self.args.intermediate_size, self.args.state_size],
|
|
axis=-1
|
|
)
|
|
|
|
# Reshape for selective scan
|
|
x = x.reshape((-1, x.shape[1], self.args.num_heads, self.args.head_dim))
|
|
A = -mx.exp(self.A_log)
|
|
|
|
# Apply selective scan
|
|
y = selective_scan(
|
|
x * dt[..., None],
|
|
A * dt,
|
|
B[..., None, :],
|
|
C[..., None, :],
|
|
self.args.chunk_size
|
|
)
|
|
|
|
# Output processing
|
|
y = y + x * self.D[None, None, :, None]
|
|
y = y.reshape((-1, y.shape[1], self.args.intermediate_size))
|
|
y = self.norm(y, z)
|
|
y = self.out_proj(y)
|
|
|
|
return y
|
|
|
|
def forward_inference(self, u: mx.array, cache=None) -> mx.array:
|
|
"""Single token processing during inference."""
|
|
assert u.shape[1] == 1, "Inference mode expects single token"
|
|
|
|
batch_size = u.shape[0]
|
|
# Use provided cache or create new one
|
|
self.cache = cache if cache is not None else Mamba2Cache.get_cache(self.args, batch_size, None)
|
|
|
|
# Project input
|
|
zxbcdt = self.in_proj(mx.squeeze(u, 1))
|
|
parts = mx.split(
|
|
zxbcdt,
|
|
[
|
|
self.args.intermediate_size,
|
|
self.args.intermediate_size + 2 * self.args.state_size
|
|
],
|
|
axis=-1
|
|
)
|
|
z, xBC = parts[0], parts[1]
|
|
dt = zxbcdt[:, -self.args.num_heads:] # Extract dt separately
|
|
|
|
# Update convolution state and apply
|
|
conv_state = self.cache.update_conv_state(xBC)
|
|
xBC = mx.sum(
|
|
conv_state * mx.transpose(self.conv1d.weight, [1, 0, 2]),
|
|
axis=-1
|
|
)
|
|
if self.args.use_conv_bias:
|
|
xBC = xBC + self.conv1d.bias
|
|
xBC = mx.sigmoid(xBC) * xBC # SiLU
|
|
|
|
# Split states and ensure proper shapes
|
|
x_splits = mx.split(
|
|
xBC,
|
|
[self.args.intermediate_size, self.args.state_size],
|
|
axis=-1
|
|
)
|
|
x, B, C = x_splits[0], x_splits[1], x_splits[2]
|
|
|
|
# Process time steps - ensure proper broadcasting
|
|
dt = mx.reshape(dt, (batch_size, self.args.num_heads))
|
|
dt = mx.clip(
|
|
nn.softplus(dt + self.dt_bias[None, :]),
|
|
self.args.time_step_min,
|
|
self.args.time_step_max
|
|
)
|
|
|
|
# SSM step with explicit shapes
|
|
A = -mx.exp(self.A_log)
|
|
dA = mx.exp(dt * A[None, :]) # Shape: (batch_size, num_heads)
|
|
|
|
# Reshape x considering intermediate size
|
|
# x shape should be (batch_size * num_heads, head_dim)
|
|
x = mx.reshape(x, (batch_size, self.args.num_heads, -1))
|
|
assert x.shape[-1] == self.args.head_dim, f"Head dimension mismatch: {x.shape[-1]} vs {self.args.head_dim}"
|
|
|
|
# Reshape B and C for ssm computation
|
|
B = mx.reshape(B, (batch_size, -1)) # Should be (batch_size, state_size)
|
|
C = mx.reshape(C, (batch_size, -1)) # Should be (batch_size, state_size)
|
|
|
|
# Compute dBx with explicit shapes
|
|
dBx = mx.einsum('bh,bs,bhd->bhds', dt, B, x)
|
|
|
|
ssm_state = self.cache.update_ssm_state(dA, dBx)
|
|
|
|
y = mx.einsum('bhds,bs->bhd', ssm_state, C)
|
|
y = y + x * self.D[None, :, None]
|
|
y = mx.reshape(y, (batch_size, self.args.intermediate_size))
|
|
|
|
# Output processing
|
|
y = self.norm(y, z)
|
|
y = self.out_proj(y)
|
|
|
|
return mx.expand_dims(y, 1)
|
|
|
|
|
|
class ResidualBlock(nn.Module):
|
|
def __init__(self, args: ModelArgs):
|
|
super().__init__()
|
|
self.mixer = Mamba2Block(args)
|
|
self.norm = nn.RMSNorm(args.hidden_size)
|
|
|
|
def __call__(self, x: mx.array, cache=None) -> mx.array:
|
|
return self.mixer(self.norm(x), cache) + x
|
|
|
|
|
|
class Mamba2Model(nn.Module):
|
|
def __init__(self, args: ModelArgs):
|
|
super().__init__()
|
|
self.args = args
|
|
self.embeddings = nn.Embedding(args.vocab_size, args.hidden_size)
|
|
self.layers = [ResidualBlock(args) for _ in range(args.num_hidden_layers)]
|
|
self.norm_f = nn.RMSNorm(args.hidden_size, eps=args.layer_norm_epsilon)
|
|
|
|
def __call__(self, x: mx.array, cache=None) -> mx.array:
|
|
x = self.embeddings(x)
|
|
if cache is None:
|
|
cache = [None] * len(self.layers)
|
|
for layer, layer_cache in zip(self.layers, cache):
|
|
x = layer(x, layer_cache)
|
|
return self.norm_f(x)
|
|
|
|
|
|
class Model(nn.Module):
|
|
def __init__(self, args: ModelArgs):
|
|
super().__init__()
|
|
self.args = args
|
|
self.backbone = Mamba2Model(args)
|
|
|
|
if not args.tie_word_embeddings:
|
|
self.lm_head = nn.Linear(args.hidden_size, args.vocab_size, bias=False)
|
|
|
|
def __call__(self, inputs: mx.array, cache=None) -> mx.array:
|
|
B, T = inputs.shape
|
|
|
|
x = self.backbone(inputs, cache)
|
|
|
|
if self.args.tie_word_embeddings:
|
|
logits = self.backbone.embeddings.as_linear(x)
|
|
else:
|
|
logits = self.lm_head(x)
|
|
|
|
return logits
|
|
|
|
def make_cache(self, batch_size=1):
|
|
return [Mamba2Cache(
|
|
batch_size=batch_size,
|
|
intermediate_size=self.args.intermediate_size,
|
|
state_size=self.args.state_size,
|
|
conv_kernel=self.args.conv_kernel,
|
|
num_heads=self.args.num_heads,
|
|
head_dim=self.args.head_dim
|
|
) for _ in range(len(self.backbone.layers))]
|
|
|
|
def sanitize(self, weights):
|
|
for k, v in weights.items():
|
|
if "conv1d.weight" in k and v.ndim == 3:
|
|
weights[k] = v.moveaxis(2, 1)
|
|
return weights |