flashinfer / wrapperfe59e7

import torch
import flashinfer

# GQA group_size=3 (24 qo_heads / 8 kv_heads) is not a power-of-2 and is
# unsupported by FlashInfer kernels.  Work-around: expand KV heads from 8
# to 24 (repeat_interleave x3) so group_size=1 (MHA), which is
# mathematically equivalent.

_WORKSPACE_SIZE_BYTES = 128 * 1024 * 1024
_workspace_cache = {}
_wrapper_cache = {}
_plan_state = {}


def _get_workspace(device):
    key = str(device)
    buf = _workspace_cache.get(key)
    if buf is None:
        buf = torch.empty(_WORKSPACE_SIZE_BYTES, dtype=torch.uint8, device=device)
        _workspace_cache[key] = buf
    return buf


def _get_wrapper(key, device):
    w = _wrapper_cache.get(key)
    if w is None:
        w = flashinfer.BatchPrefillWithRaggedKVCacheWrapper(_get_workspace(device), kv_layout="NHD")
        _wrapper_cache[key] = w
    return w


def run(q, k, v, qo_indptr, kv_indptr, sm_scale):
    total_q, num_qo_heads, head_dim = q.shape
    total_kv, num_kv_heads, _ = k.shape
    batch_size = qo_indptr.shape[0] - 1
    device = q.device
    group_size = num_qo_heads // num_kv_heads  # 3
    # Expand KV heads: [total_kv, 8, 128] -> [total_kv, 24, 128]
    k_exp = k.repeat_interleave(group_size, dim=1)
    v_exp = v.repeat_interleave(group_size, dim=1)
    expanded_heads = num_qo_heads  # 24
    wkey = (str(device), num_qo_heads, expanded_heads, head_dim, q.dtype, k.dtype)
    wrapper = _get_wrapper(wkey, device)
    state = _plan_state.get(wkey)
    needs_plan = (
        state is None
        or state["total_q"] != total_q
        or state["total_kv"] != total_kv
        or state["qo_ptr"] != qo_indptr.data_ptr()
        or state["kv_ptr"] != kv_indptr.data_ptr()
    )
    if needs_plan:
        wrapper.plan(
            qo_indptr=qo_indptr,
            kv_indptr=kv_indptr,
            num_qo_heads=num_qo_heads,
            num_kv_heads=expanded_heads,
            head_dim_qk=head_dim,
            causal=True,
            sm_scale=float(sm_scale),
            q_data_type=q.dtype,
            kv_data_type=k.dtype,
        )
        _plan_state[wkey] = {"total_q": total_q, "total_kv": total_kv, "qo_ptr": qo_indptr.data_ptr(), "kv_ptr": kv_indptr.data_ptr()}
    output, lse = wrapper.run(q, k_exp, v_exp, return_lse=True)
    return output, lse