A high-performance attention mechanism that computes softmax normalization in a single streaming pass using running accumulators (online softmax). The design achieves O(N) memory, strong numerical stability via log-sum-exp, and competitive throughput with modern flash attention baselines.

This repository provides:

• A production-oriented StreamAttention module for drop-in use in transformer models
• A fused online softmax kernel in Triton with safe PyTorch SDPA fallbacks
• A FlashAttention-3 baseline wrapper
• Benchmark and accuracy CLIs for reproducible comparisons
• Utilities for memory profiling and KV-cache compression
• Optional integration helpers for Hugging Face models

• Added a single-sweep Triton backward path (streaming dQ/dK/dV using saved lse) that mirrors the forward pass, including masks, dropout, and ALiBi.
• Triton forward now supports boolean/additive masks, dropout, deterministic Philox seeding, and ALiBi bias without SDPA fallback.
• Expanded tests/docs covering mask/dropout/ALiBi parity plus deterministic mode usage.

The project depends on PyTorch (CUDA optional) and Triton (optional for the custom fused kernel).

# Editable install (preferred for development)
pip install -e .

# Optional extras
pip install -e .[hf]        # Hugging Face integration helpers

Notes:

• If the environment is system-managed, you may need to pass --break-system-packages to pip or use a virtual environment.
• Triton is optional. If Triton is not available, the implementation falls back to PyTorch SDPA (flash backend on CUDA where available).

import torch
from stream_attention import StreamAttention, StreamAttentionConfig

# Configure the attention
config = StreamAttentionConfig(
    num_heads=32,
    head_dim=128,
    tile_size_q=128,
    tile_size_k=64,
    use_qkv_projections=True
)

# Create the module
attention = (StreamAttention(config).cuda() if torch.cuda.is_available() else StreamAttention(config))

# Hidden-states path (with internal Q,K,V projections)
batch_size, seq_len = 2, 1024
hidden_dim = config.num_heads * config.head_dim
x = torch.randn(batch_size, seq_len, hidden_dim, device=attention.attention.device, dtype=(torch.float16 if torch.cuda.is_available() else torch.float32))

with torch.no_grad():
    y = attention(hidden_states=x)
print(y.shape)

# Explicit Q,K,V path (supports attn_mask/dropout/ALiBi in Triton, falling back to SDPA when needed)
q = torch.randn(batch_size, seq_len, config.num_heads, config.head_dim, device=x.device, dtype=x.dtype)
k = torch.randn_like(q)
v = torch.randn_like(q)
with torch.no_grad():
    # Example boolean mask [B, S_k]
    key_padding_mask = torch.ones(batch_size, seq_len, dtype=torch.bool, device=x.device)
    key_padding_mask[:, -16:] = False  # mask out last 16 positions
    y_qkv = attention(q, k, v, causal=True, attention_mask=key_padding_mask)
print(y_qkv.shape)

• Purpose: High-level module. Accepts either hidden_states ([B, T, H*D]) or explicit (query, key, value) ([B, T, H, D]).
• Signature (selected):
  • forward(query, key, value, causal: bool=True, return_lse: bool=False, attention_mask: Optional[Tensor]=None, dropout_p: float=0.0, alibi_slopes: Optional[Tensor]=None, deterministic: Optional[bool]=None) -> Tensor (and lse if requested)
  • benchmark(seq_len: int, batch_size: int=1, warmup: int=10, iterations: int=100) -> metrics dict
  • set_deterministic(enabled: bool, seed: Optional[int]=None) -> control deterministic dropout/mask behavior
• Autograd: The Triton kernel now runs a single-sweep backward pass (streaming dQ/dK/dV using the saved lse) covering masks, dropout, and ALiBi. PyTorch SDPA is used only when Triton is unavailable.

Use create_stream_attention to obtain an attention layer with a familiar nn.MultiheadAttention interface. Triton kernels are used automatically when available, otherwise PyTorch's SDPA backend is selected:

import torch
from stream_attention import create_stream_attention

mha = create_stream_attention(embed_dim=512, num_heads=8, batch_first=True)
x = torch.randn(2, 16, 512)
out, _ = mha(x, x, x)

• Purpose: Baseline using PyTorch SDPA with the flash backend on CUDA, falling back gracefully on CPU.
• Signature (selected):
  • forward(query, key, value, causal: bool=True) → Tensor
  • benchmark(...) → metrics dict

Selected fields (see source for full set):

• num_heads, head_dim
• tile_size_q, tile_size_k
• use_fp16, use_qkv_projections, qkv_bias, use_layer_norm, dropout
• enable_distributed
• max_sequence_length
• Ring/Star attention parameters for long-context variants
• .from_env() reads STREAM_ATTENTION_* variables; see .env.example

Two CLIs are provided:

# Performance comparison across sequence lengths
stream-attention-benchmark --seq 512 1024 2048 4096 --batch 1 --heads 8 --dim 64 --warmup 10 --iters 50

# Accuracy sanity check on random inputs
stream-attention-test --seq 1024 --batch 2 --heads 8 --dim 64 --dtype fp16

Behavior and methodology:

• On CUDA, the baseline uses PyTorch SDPA with the flash backend (FlashAttention-3 path). On CPU, both implementations use SDPA in fp32.
• Metrics reported: per-iteration latency, estimated TFLOPS, and approximate bandwidth based on tensor traffic. Measurements are averaged after warmup.
• The fused kernel uses Triton on CUDA for the forward pass; when gradients are required and dropout_p == 0, the streaming backward (single sweep over K/V) is invoked. If dropout is active during training, the module falls back to SDPA for gradient computation.
• For reproducibility, fix random seeds, pin CUDA clocks if applicable, and isolate runs. Actual performance depends on GPU architecture, drivers, and PyTorch/Triton versions.

Example output (format):

SeqLen  Fused(ms)  Fused(TF)  FA3(ms)  FA3(TF)  FA3/Fused(ms)
  1024      0.650     90.12     0.700     83.70          1.08

The fused kernel streams over K/V tiles while maintaining per-query running statistics for online softmax, avoiding materialization of the full attention matrix.

• Tiling:

  • Queries are processed in blocks of TILE_M (configurable via autotune); keys/values are streamed in tiles of TILE_N.
  • Head dimension D is processed vectorized. The kernel uses tl.dot(q, tl.trans(k)) to compute QK^T for each tile.

• Online softmax with log-sum-exp:

  • Maintain running_max[m], acc_num[m, :], acc_den[m] for each query row m.
  • For each K/V tile:
    • tile_max = max_j qk[m, j], new_max = max(running_max[m], tile_max)
    • Rescale previous accumulators by exp(running_max - new_max)
    • exp_qk = exp(qk - new_max)
    • acc_num += exp_qk @ V_tile, acc_den += sum(exp_qk, axis=tile)
    • running_max = new_max
  • Final output: acc_num / acc_den[:, None]. Log-sum-exp lse = running_max + log(acc_den) may be returned for analysis.

• Autotuning:

  • Multiple Triton configurations for TILE_M/TILE_N, warps, and stages are provided.
  • Kernel grid: (ceil_div(M, TILE_M), batch, heads)

• Numerical stability:

  • The log-sum-exp re-centering preserves stability across tiles and long sequences.
  • On CPU, fp16/bf16 inputs are upcast to fp32 where necessary.

• Backward:

  • The Triton path is presently forward-oriented. For training, the module selects PyTorch SDPA, which supports autograd.

• Ring Attention (stream_attention/core/ring_attention.py): partitions sequences across devices with overlapped communication and computation. Suitable for near-infinite contexts with linear memory scaling.
• Star Attention (stream_attention/core/star_attention.py): two-phase approach (context encoding, then query processing with aggregated attention). Reduces end-to-end latency for long sequences while preserving accuracy.

Both modules follow numerically stable aggregation strategies (log-sum-exp weighted combining) and can be paired with KV compression.

• KV-cache compression strategies: importance-based, chunk-based, quantized, and hybrid.
• Gradient checkpointing helpers for sequential modules.
• MemoryProfiler to snapshot and report peak allocations.

Example:

from stream_attention.utils.memory import create_kv_compressor
comp = create_kv_compressor('hybrid')
ck, cv, stats = comp.compress(keys, values, compression_ratio=8.0)
print(stats)

Helpers are provided to replace attention modules in HF models:

from transformers import AutoModel
from stream_attention.integration.hf import replace_llama_attention
from stream_attention import StreamAttentionConfig

model = AutoModel.from_pretrained("meta-llama/Llama-2-7b-hf")
cfg = StreamAttentionConfig(num_heads=model.config.num_attention_heads, head_dim=model.config.hidden_size // model.config.num_attention_heads)
num_replaced = replace_llama_attention(model, cfg)
print(f"Replaced {num_replaced} attention modules")

• PyTorch 2.0+
• CUDA: fp16 and bf16 supported via SDPA (flash backend where available). Triton kernel requires CUDA.
• CPU: falls back to SDPA with fp32 compute for correctness and stability.
• Distributed: query sharding is supported in the fused module for multi-GPU; ring/star provides long-context strategies.

• Native RoPE / relative positional bias fusion in the Triton kernel (forward + backward)
• Advanced pipelining (warp specialization, asynchronous staging) and Hopper-specific paths (WGMMA/TMA)
• Extended autotune coverage across architectures and sequence regimes
• Optional FP8 path with block-wise scaling

• Benchmarks: stream-attention-benchmark CLI
• Accuracy checks: stream-attention-test CLI
• Examples: examples/ directory includes basic usage, integrations, and long-context runs
• Environment variables: see .env.example; StreamAttentionConfig.from_env() can bootstrap configuration

Apache License. See LICENSE for details.