Skip to content
KernelIndex
Search⌘K

submission 753984

guojun21 · python · License unknown

Use it

Vendorable · source mirrored · license unknownView source →

No package. Vendor the mirrored source: 233 lines, June 9 Researcher Reciprocity License v1.0.

submission_v177_single_4wg.py
curl "https://kernelindex.com/api/v1/implementations/kernelbot-amd-moe-mxfp4-753984?include=source"
interfacepython
Compatibility
measured onAMD Instinct MI355X
declared hardwareAMD Instinct MI355X
architecturesgfx950
dtypesbf16, fp32, fp8_e8m0, int32, mxfp4

Benchmark evidence

1 measurement across 1 GPU, fastest first.

Operation / workload
Hardware
Latency
Rank
Observed
AMD MXFP4 MoEsuite of 7 cases
AMD Instinct MI355X
120.3µs
#37 of 782
2026-04-07

Reported · How evidence levels are derived →

Source and license

sourceavailable
revision digestsha256:ddb45f973ec257ffc2a8f3f5a558d0831c31408557538b07f48c6126df1baff0
license declaredunknown
license concludedunknown
authorsguojun21
imported2026-08-15

Techniques

Extracted from the mirrored source by pattern, never inferred. Each row cites its line.

fp4"stage": 2, "a_dtype": "fp4", "b_dtype": "fp4", "out_dtype": "bf16",

Kernel source

submission_v177_single_4wg.py233 lines
#!POPCORN leaderboard amd-moe-mxfp4
#!POPCORN gpu MI355X

"""
v177: Single 4WG module strategy - ALL ksplit=0 shapes use 4WG_M128.
E=257 bs=512 uses 4WG_M128 with block_m=64 (same module as E=33 shapes).
Avoids 4WG_M32 to stay within 12min compilation budget.
"""

import os
import functools
import torch
import triton
from task import input_t, output_t

import aiter
from aiter import ActivationType, QuantType, dtypes
from aiter.fused_moe import get_2stage_cfgs, get_padded_M, get_inter_dim
import aiter.fused_moe as _fused_moe_module
import aiter.ops.flydsl.moe_kernels as _flydsl_moe_kernels
from aiter.ops.triton._triton_kernels.quant.fused_mxfp4_quant import (
    _fused_dynamic_mxfp4_quant_moe_sort_kernel,
)

_flydsl_moe_kernels._KERNEL_PARAMS["flydsl_moe2_afp4_wfp4_bf16_t16x128x128_atomic"] = {
    "stage": 2, "a_dtype": "fp4", "b_dtype": "fp4", "out_dtype": "bf16",
    "tile_m": 16, "tile_n": 128, "tile_k": 128, "mode": "atomic", "MPerBlock": 16,
}

_CUSTOM_CONFIGS = {}

def _make_key(token, inter_dim, expert):
    return (
        256, token, 7168, inter_dim, expert, 9,
        "ActivationType.Silu", "torch.bfloat16",
        "torch.float4_e2m1fn_x2", "torch.float4_e2m1fn_x2",
        "QuantType.per_1x32", True, False,
    )

_4WG = "moe_ck2stages_gemm1_256x128x128x128_1x4_MulABScaleShuffled_v3_Nswizzle0_Quant3_MulRoutedWeight0_silu_FP4X2_FP4X2_B16"
_FLY = "flydsl_moe2_afp4_wfp4_bf16_t16x128x128_atomic"

# E=257: bs=16/128 ksplit=2 (sparse, bf16), bs=512 ksplit=0 (fp4 with shared 4WG_M128)
_CUSTOM_CONFIGS[_make_key(16, 256, 257)] = {
    "block_m": 16, "ksplit": 2, "kernelName1": "", "kernelName2": "", "run_1stage": False,
}
_CUSTOM_CONFIGS[_make_key(128, 256, 257)] = {
    "block_m": 16, "ksplit": 2, "kernelName1": "", "kernelName2": "", "run_1stage": False,
}
_CUSTOM_CONFIGS[_make_key(512, 256, 257)] = {
    "block_m": 64, "ksplit": 0,
    "kernelName1": _4WG, "kernelName2": _FLY, "run_1stage": False,
}

# E=33: bs=16 ksplit=2, bs=128+ ksplit=0 (all use same 4WG_M128)
_CUSTOM_CONFIGS[_make_key(16, 512, 33)] = {
    "block_m": 32, "ksplit": 2, "kernelName1": "", "kernelName2": "", "run_1stage": False,
}
_CUSTOM_CONFIGS[_make_key(128, 512, 33)] = {
    "block_m": 64, "ksplit": 0,
    "kernelName1": _4WG, "kernelName2": _FLY, "run_1stage": False,
}
_CUSTOM_CONFIGS[_make_key(512, 512, 33)] = {
    "block_m": 64, "ksplit": 0,
    "kernelName1": _4WG, "kernelName2": _FLY, "run_1stage": False,
}
_CUSTOM_CONFIGS[_make_key(512, 2048, 33)] = {
    "block_m": 64, "ksplit": 0,
    "kernelName1": _4WG, "kernelName2": _FLY, "run_1stage": False,
}

# --- Buffer caches ---
_buf = {}
_qbuf = {}

def _get_sorting_bufs(M, E, topk, model_dim, block_m, device):
    key = (M, E, topk, model_dim, block_m)
    if key not in _buf:
        max_pad = M * topk + E * block_m - topk
        max_blk = (max_pad + block_m - 1) // block_m
        _buf[key] = {
            "sid": torch.empty(max_pad, dtype=dtypes.i32, device=device),
            "sw": torch.empty(max_pad, dtype=dtypes.fp32, device=device),
            "se": torch.empty(max_blk, dtype=dtypes.i32, device=device),
            "nv": torch.empty(2, dtype=dtypes.i32, device=device),
            "out": torch.empty((M, model_dim), dtype=torch.bfloat16, device=device),
            "a2": torch.empty((M, topk, 0), dtype=torch.bfloat16, device=device),
        }
    return _buf[key]

def _get_a2(M, topk, inter_dim, device):
    key = ("a2", M, topk, inter_dim)
    if key not in _buf:
        _buf[key] = torch.empty((M, topk, inter_dim), dtype=torch.bfloat16, device=device)
    return _buf[key]

def _quant_prealloc(x, sorted_ids, num_valid_ids, token_num, topk, block_m, device):
    M, N = x.shape
    QBS = 32
    BLK_Mx = 128
    BLK_M, BLK_N = 32, 8
    BLK_M_u32, BLK_N_u32 = 16, 4

    scaleN = triton.cdiv(N, QBS)
    M_o = sorted_ids.shape[0]

    qk = (M, N, M_o, topk)
    if qk not in _qbuf:
        _qbuf[qk] = {
            "fp4": torch.empty((M, N // 2), dtype=torch.uint8, device=device),
            "bs": torch.empty(
                (triton.cdiv(M_o, BLK_M), triton.cdiv(scaleN, BLK_N),
                 BLK_N_u32, BLK_M_u32, 4),
                dtype=torch.uint8, device=device,
            ),
        }
    qb = _qbuf[qk]

    num_pid = triton.cdiv(M, BLK_Mx) * scaleN + triton.cdiv(
        M_o, BLK_M
    ) * triton.cdiv(scaleN, BLK_N)

    _fused_dynamic_mxfp4_quant_moe_sort_kernel[(num_pid,)](
        x, qb["fp4"], sorted_ids, num_valid_ids, qb["bs"],
        M, N, scaleN,
        *x.stride(), *qb["fp4"].stride(), *qb["bs"].stride(),
        token_num=token_num, M_i=M, N_i=scaleN,
        MXFP4_QUANT_BLOCK_SIZE=QBS, BLOCK_SIZE_Mx=BLK_Mx,
        BLOCK_SIZE_M=BLK_M // 2, BLOCK_SIZE_N=BLK_N // 2,
        TOPK=topk,
    )

    return (
        qb["fp4"].view(dtypes.fp4x2),
        qb["bs"].view(dtypes.fp8_e8m0).view(-1, scaleN),
    )


_injected = False

def _inject():
    global _injected
    if _injected:
        return
    _injected = True
    if _fused_moe_module.cfg_2stages is None:
        import pandas as pd
        from aiter.jit.core import AITER_CONFIGS
        tf = AITER_CONFIGS.AITER_CONFIG_FMOE_FILE
        if os.path.exists(tf):
            _IDX = [
                "cu_num", "token", "model_dim", "inter_dim", "expert", "topk",
                "act_type", "dtype", "q_dtype_a", "q_dtype_w", "q_type",
                "use_g1u1", "doweight_stage1",
            ]
            df = pd.read_csv(tf)
            if "_tag" in df.columns:
                df = df[df["_tag"].fillna("") == ""]
            _fused_moe_module.cfg_2stages = df.set_index(_IDX).to_dict("index")
        else:
            _fused_moe_module.cfg_2stages = {}
    _fused_moe_module.cfg_2stages.update(_CUSTOM_CONFIGS)


def custom_kernel(data: input_t) -> output_t:
    (
        hidden_states, _w1r, _w2r, _w1sr, _w2sr,
        w1, w2, w1s, w2s,
        topk_weights, topk_ids, config,
    ) = data

    _inject()

    M = hidden_states.shape[0]
    topk = topk_ids.shape[1]
    device = hidden_states.device
    dhp = config["d_hidden_pad"]
    dep = config["d_expert_pad"]
    hidden_pad = dhp - config["d_hidden"]
    intermediate_pad = dep - config["d_expert"]

    E, model_dim, inter_dim = get_inter_dim(w1.shape, w2.shape)
    padded_M = get_padded_M(M)

    metadata = get_2stage_cfgs(
        padded_M, model_dim, inter_dim, E, topk,
        torch.bfloat16, dtypes.fp4x2, dtypes.fp4x2,
        QuantType.per_1x32, True, ActivationType.Silu,
        False, hidden_pad, intermediate_pad, True,
    )
    block_m = int(metadata.block_m)

    b = _get_sorting_bufs(M, E, topk, model_dim, block_m, device)
    aiter.moe_sorting_fwd(
        topk_ids, topk_weights,
        b["sid"], b["sw"], b["se"], b["nv"], b["out"],
        E, block_m, None, None, 0,
    )

    w1sv = w1s.view(dtypes.fp8_e8m0)
    w2sv = w2s.view(dtypes.fp8_e8m0)
    a2_buf = _get_a2(M, topk, inter_dim, device)

    if metadata.ksplit > 1:
        a1 = hidden_states.to(torch.bfloat16)
        a2 = metadata.stage1(
            a1, w1, w2, b["sid"], b["se"], b["nv"], a2_buf, topk,
            block_m=block_m, a1_scale=None, w1_scale=w1sv, sorted_weights=None,
        )
        metadata.stage2(
            a2, w1, w2, b["sid"], b["se"], b["nv"], b["out"], topk,
            w2_scale=w2sv, a2_scale=None, block_m=block_m, sorted_weights=b["sw"],
        )
    else:
        a1, a1s = _quant_prealloc(
            hidden_states, b["sid"], b["nv"], M, 1, block_m, device,
        )
        a2 = metadata.stage1(
            a1, w1, w2, b["sid"], b["se"], b["nv"], a2_buf, topk,
            block_m=block_m, a1_scale=a1s, w1_scale=w1sv, sorted_weights=None,
        )
        a2_flat = a2.view(-1, inter_dim)
        a2q, a2s = _quant_prealloc(
            a2_flat, b["sid"], b["nv"], M, topk, block_m, device,
        )
        a2q = a2q.view(M, topk, -1)
        metadata.stage2(
            a2q, w1, w2, b["sid"], b["se"], b["nv"], b["out"], topk,
            w2_scale=w2sv, a2_scale=a2s, block_m=block_m, sorted_weights=b["sw"],
        )

    return b["out"]
scrolls · 233 lines total

Source code from GPU Mode and the KernelBot dataset · June 9 Researcher Reciprocity License v1.0

Changes from previous submission

Against this author's previous submission submission 732358.

⋯ 1 unchanged lines
#!POPCORN gpu MI355X
"""
- v168: Pre-allocate quantization output buffers (x_fp4, blockscale_e8m0_sorted)
- for both stage1 and stage2 quant calls. Inline the fused_dynamic_mxfp4_quant_moe_sort
- Triton kernel launch with cached output tensors to eliminate 4 torch.empty allocations
- per forward pass on CK 2-stage shapes.
+ v177: Single 4WG module strategy - ALL ksplit=0 shapes use 4WG_M128.
+ E=257 bs=512 uses 4WG_M128 with block_m=64 (same module as E=33 shapes).
+ Avoids 4WG_M32 to stay within 12min compilation budget.
"""
+
import os
import functools
import torch
import triton
- from typing import Dict, Tuple, Optional
from task import input_t, output_t
import aiter
from aiter import ActivationType, QuantType, dtypes
- from aiter.fused_moe import (
- get_2stage_cfgs, get_padded_M, get_inter_dim,
- ck_moe_stage1, cktile_moe_stage1, cktile_moe_stage2,
- _flydsl_stage2_wrapper,
- )
+ from aiter.fused_moe import get_2stage_cfgs, get_padded_M, get_inter_dim
import aiter.fused_moe as _fused_moe_module
import aiter.ops.flydsl.moe_kernels as _flydsl_moe_kernels
from aiter.ops.triton._triton_kernels.quant.fused_mxfp4_quant import (
_fused_dynamic_mxfp4_quant_moe_sort_kernel,
)
- from aiter.utility import fp4_utils
- # Register FlyDSL tile_k=128 kernels
- _flydsl_moe_kernels._KERNEL_PARAMS["flydsl_moe2_afp4_wfp4_bf16_t32x128x128_atomic"] = {
- "stage": 2, "a_dtype": "fp4", "b_dtype": "fp4", "out_dtype": "bf16",
- "tile_m": 32, "tile_n": 128, "tile_k": 128, "mode": "atomic", "MPerBlock": 32,
- }
- _flydsl_moe_kernels._KERNEL_PARAMS["flydsl_moe2_afp4_wfp4_bf16_t32x256x128_atomic"] = {
- "stage": 2, "a_dtype": "fp4", "b_dtype": "fp4", "out_dtype": "bf16",
- "tile_m": 32, "tile_n": 256, "tile_k": 128, "mode": "atomic", "MPerBlock": 32,
- }
- _flydsl_moe_kernels._KERNEL_PARAMS["flydsl_moe2_afp4_wfp4_bf16_t16x256x128_atomic"] = {
- "stage": 2, "a_dtype": "fp4", "b_dtype": "fp4", "out_dtype": "bf16",
- "tile_m": 16, "tile_n": 256, "tile_k": 128, "mode": "atomic", "MPerBlock": 16,
- }
_flydsl_moe_kernels._KERNEL_PARAMS["flydsl_moe2_afp4_wfp4_bf16_t16x128x128_atomic"] = {
"stage": 2, "a_dtype": "fp4", "b_dtype": "fp4", "out_dtype": "bf16",
"tile_m": 16, "tile_n": 128, "tile_k": 128, "mode": "atomic", "MPerBlock": 16,
}
- # Shape configs
_CUSTOM_CONFIGS = {}
def _make_key(token, inter_dim, expert):
⋯ 4 unchanged lines
"QuantType.per_1x32", True, False,
)
- _4WG_STAGE1_M128 = "moe_ck2stages_gemm1_256x128x128x128_1x4_MulABScaleShuffled_v3_Nswizzle0_Quant3_MulRoutedWeight0_silu_FP4X2_FP4X2_B16"
- _4WG_STAGE1_M32 = "moe_ck2stages_gemm1_256x32x128x128_1x4_MulABScaleShuffled_v3_Nswizzle0_Quant3_MulRoutedWeight0_silu_FP4X2_FP4X2_B16"
- _FLYDSL_STAGE2_M16_K128 = "flydsl_moe2_afp4_wfp4_bf16_t16x128x128_atomic"
+ _4WG = "moe_ck2stages_gemm1_256x128x128x128_1x4_MulABScaleShuffled_v3_Nswizzle0_Quant3_MulRoutedWeight0_silu_FP4X2_FP4X2_B16"
+ _FLY = "flydsl_moe2_afp4_wfp4_bf16_t16x128x128_atomic"
- # E=33 shapes
+ # E=257: bs=16/128 ksplit=2 (sparse, bf16), bs=512 ksplit=0 (fp4 with shared 4WG_M128)
+ _CUSTOM_CONFIGS[_make_key(16, 256, 257)] = {
+ "block_m": 16, "ksplit": 2, "kernelName1": "", "kernelName2": "", "run_1stage": False,
+ }
+ _CUSTOM_CONFIGS[_make_key(128, 256, 257)] = {
+ "block_m": 16, "ksplit": 2, "kernelName1": "", "kernelName2": "", "run_1stage": False,
+ }
+ _CUSTOM_CONFIGS[_make_key(512, 256, 257)] = {
+ "block_m": 64, "ksplit": 0,
+ "kernelName1": _4WG, "kernelName2": _FLY, "run_1stage": False,
+ }
+
+ # E=33: bs=16 ksplit=2, bs=128+ ksplit=0 (all use same 4WG_M128)
_CUSTOM_CONFIGS[_make_key(16, 512, 33)] = {
- "block_m": 32, "ksplit": 2, "kernelName1": "", "kernelName2": "",
- "run_1stage": False,
+ "block_m": 32, "ksplit": 2, "kernelName1": "", "kernelName2": "", "run_1stage": False,
}
_CUSTOM_CONFIGS[_make_key(128, 512, 33)] = {
"block_m": 64, "ksplit": 0,
- "kernelName1": _4WG_STAGE1_M128, "kernelName2": _FLYDSL_STAGE2_M16_K128,
- "run_1stage": False,
+ "kernelName1": _4WG, "kernelName2": _FLY, "run_1stage": False,
}
_CUSTOM_CONFIGS[_make_key(512, 512, 33)] = {
"block_m": 64, "ksplit": 0,
- "kernelName1": _4WG_STAGE1_M128, "kernelName2": _FLYDSL_STAGE2_M16_K128,
- "run_1stage": False,
+ "kernelName1": _4WG, "kernelName2": _FLY, "run_1stage": False,
}
_CUSTOM_CONFIGS[_make_key(512, 2048, 33)] = {
"block_m": 64, "ksplit": 0,
- "kernelName1": _4WG_STAGE1_M128, "kernelName2": _FLYDSL_STAGE2_M16_K128,
- "run_1stage": False,
+ "kernelName1": _4WG, "kernelName2": _FLY, "run_1stage": False,
}
- # E=257 shapes
- _CUSTOM_CONFIGS[_make_key(16, 256, 257)] = {
- "block_m": 16, "ksplit": 2, "kernelName1": "", "kernelName2": "",
- "run_1stage": False,
- }
- _CUSTOM_CONFIGS[_make_key(128, 256, 257)] = {
- "block_m": 16, "ksplit": 2, "kernelName1": "", "kernelName2": "",
- "run_1stage": False,
- }
- _CUSTOM_CONFIGS[_make_key(512, 256, 257)] = {
- "block_m": 32, "ksplit": 0,
- "kernelName1": _4WG_STAGE1_M32, "kernelName2": _FLYDSL_STAGE2_M16_K128,
- "run_1stage": False,
- "use_non_temporal_load": True,
- }
+ # --- Buffer caches ---
+ _buf = {}
+ _qbuf = {}
- # Pre-allocated buffer cache
- _buffer_cache = {}
+ def _get_sorting_bufs(M, E, topk, model_dim, block_m, device):
+ key = (M, E, topk, model_dim, block_m)
+ if key not in _buf:
+ max_pad = M * topk + E * block_m - topk
+ max_blk = (max_pad + block_m - 1) // block_m
+ _buf[key] = {
+ "sid": torch.empty(max_pad, dtype=dtypes.i32, device=device),
+ "sw": torch.empty(max_pad, dtype=dtypes.fp32, device=device),
+ "se": torch.empty(max_blk, dtype=dtypes.i32, device=device),
+ "nv": torch.empty(2, dtype=dtypes.i32, device=device),
+ "out": torch.empty((M, model_dim), dtype=torch.bfloat16, device=device),
+ "a2": torch.empty((M, topk, 0), dtype=torch.bfloat16, device=device),
+ }
+ return _buf[key]
- def _get_or_alloc_sorting_buffers(M, E, topk, model_dim, block_size_M, device):
- """Pre-allocate moe_sorting output buffers."""
- key = ("sort", M, E, topk, model_dim, block_size_M)
- if key in _buffer_cache:
- return _buffer_cache[key]
-
- max_num_tokens_padded = int(M * topk + E * block_size_M - topk)
- max_num_m_blocks = int((max_num_tokens_padded + block_size_M - 1) // block_size_M)
-
- bufs = {
- "sorted_ids": torch.empty(max_num_tokens_padded, dtype=dtypes.i32, device=device),
- "sorted_weights": torch.empty(max_num_tokens_padded, dtype=dtypes.fp32, device=device),
- "sorted_expert_ids": torch.empty(max_num_m_blocks, dtype=dtypes.i32, device=device),
- "num_valid_ids": torch.empty(2, dtype=dtypes.i32, device=device),
- "moe_buf": torch.empty((M, model_dim), dtype=torch.bfloat16, device=device),
- }
- _buffer_cache[key] = bufs
- return bufs
-
- def _get_or_alloc_a2(M, topk, inter_dim, device):
- """Pre-allocate a2 intermediate buffer."""
+ def _get_a2(M, topk, inter_dim, device):
key = ("a2", M, topk, inter_dim)
- if key in _buffer_cache:
- return _buffer_cache[key]
- buf = torch.empty((M, topk, inter_dim), dtype=torch.bfloat16, device=device)
- _buffer_cache[key] = buf
- return buf
+ if key not in _buf:
+ _buf[key] = torch.empty((M, topk, inter_dim), dtype=torch.bfloat16, device=device)
+ return _buf[key]
- def _get_or_alloc_quant_buffers(M, N, sorted_ids_len, topk, device):
- """Pre-allocate quantization output buffers for fused_dynamic_mxfp4_quant_moe_sort."""
- MXFP4_QUANT_BLOCK_SIZE = 32
- BLOCK_SIZE_M, BLOCK_SIZE_N = 32, 8
- BLOCK_SIZE_M_u32, BLOCK_SIZE_N_u32 = 16, 4
-
- key = ("quant", M, N, sorted_ids_len, topk)
- if key in _buffer_cache:
- return _buffer_cache[key]
-
- x_fp4 = torch.empty((M, N // 2), dtype=torch.uint8, device=device)
- scaleN = triton.cdiv(N, MXFP4_QUANT_BLOCK_SIZE)
- M_o = sorted_ids_len
- N_o = scaleN
-
- blockscale_e8m0_sorted = torch.empty(
- (
- triton.cdiv(M_o, BLOCK_SIZE_M),
- triton.cdiv(N_o, BLOCK_SIZE_N),
- BLOCK_SIZE_N_u32,
- BLOCK_SIZE_M_u32,
- 4,
- ),
- dtype=torch.uint8,
- device=device,
- )
-
- bufs = {"x_fp4": x_fp4, "blockscale": blockscale_e8m0_sorted}
- _buffer_cache[key] = bufs
- return bufs
-
- def _quant_prealloc(x, sorted_ids, num_valid_ids, token_num, topk, block_size, device):
- """Inline fused_dynamic_mxfp4_quant_moe_sort with pre-allocated output buffers."""
+ def _quant_prealloc(x, sorted_ids, num_valid_ids, token_num, topk, block_m, device):
M, N = x.shape
- MXFP4_QUANT_BLOCK_SIZE = 32
- BLOCK_SIZE_Mx = 128
- BLOCK_SIZE_M, BLOCK_SIZE_N = 32, 8
+ QBS = 32
+ BLK_Mx = 128
+ BLK_M, BLK_N = 32, 8
+ BLK_M_u32, BLK_N_u32 = 16, 4
- scaleN = triton.cdiv(N, MXFP4_QUANT_BLOCK_SIZE)
- M_i, N_i = M, scaleN
+ scaleN = triton.cdiv(N, QBS)
M_o = sorted_ids.shape[0]
- # Get pre-allocated buffers
- qbufs = _get_or_alloc_quant_buffers(M, N, M_o, topk, device)
- x_fp4 = qbufs["x_fp4"]
- blockscale_e8m0_sorted = qbufs["blockscale"]
+ qk = (M, N, M_o, topk)
+ if qk not in _qbuf:
+ _qbuf[qk] = {
+ "fp4": torch.empty((M, N // 2), dtype=torch.uint8, device=device),
+ "bs": torch.empty(
+ (triton.cdiv(M_o, BLK_M), triton.cdiv(scaleN, BLK_N),
+ BLK_N_u32, BLK_M_u32, 4),
+ dtype=torch.uint8, device=device,
+ ),
+ }
+ qb = _qbuf[qk]
- num_pid = triton.cdiv(M, BLOCK_SIZE_Mx) * scaleN + triton.cdiv(
- M_o, BLOCK_SIZE_M
- ) * triton.cdiv(N_i, BLOCK_SIZE_N)
+ num_pid = triton.cdiv(M, BLK_Mx) * scaleN + triton.cdiv(
+ M_o, BLK_M
+ ) * triton.cdiv(scaleN, BLK_N)
_fused_dynamic_mxfp4_quant_moe_sort_kernel[(num_pid,)](
- x,
- x_fp4,
- sorted_ids,
- num_valid_ids,
- blockscale_e8m0_sorted,
- M,
- N,
- scaleN,
- *x.stride(),
- *x_fp4.stride(),
- *blockscale_e8m0_sorted.stride(),
- token_num=token_num,
- M_i=M_i,
- N_i=N_i,
- MXFP4_QUANT_BLOCK_SIZE=MXFP4_QUANT_BLOCK_SIZE,
- BLOCK_SIZE_Mx=BLOCK_SIZE_Mx,
- BLOCK_SIZE_M=BLOCK_SIZE_M // 2,
- BLOCK_SIZE_N=BLOCK_SIZE_N // 2,
+ x, qb["fp4"], sorted_ids, num_valid_ids, qb["bs"],
+ M, N, scaleN,
+ *x.stride(), *qb["fp4"].stride(), *qb["bs"].stride(),
+ token_num=token_num, M_i=M, N_i=scaleN,
+ MXFP4_QUANT_BLOCK_SIZE=QBS, BLOCK_SIZE_Mx=BLK_Mx,
+ BLOCK_SIZE_M=BLK_M // 2, BLOCK_SIZE_N=BLK_N // 2,
TOPK=topk,
)
return (
- x_fp4.view(dtypes.fp4x2),
- blockscale_e8m0_sorted.view(dtypes.fp8_e8m0).view(-1, scaleN),
+ qb["fp4"].view(dtypes.fp4x2),
+ qb["bs"].view(dtypes.fp8_e8m0).view(-1, scaleN),
)
_injected = False
- def _inject_configs():
+ def _inject():
global _injected
if _injected:
return
_injected = True
-
if _fused_moe_module.cfg_2stages is None:
import pandas as pd
from aiter.jit.core import AITER_CONFIGS
- tune_file = AITER_CONFIGS.AITER_CONFIG_FMOE_FILE
- if os.path.exists(tune_file):
- _INDEX_COLS = [
+ tf = AITER_CONFIGS.AITER_CONFIG_FMOE_FILE
+ if os.path.exists(tf):
+ _IDX = [
"cu_num", "token", "model_dim", "inter_dim", "expert", "topk",
"act_type", "dtype", "q_dtype_a", "q_dtype_w", "q_type",
"use_g1u1", "doweight_stage1",
]
- df = pd.read_csv(tune_file)
+ df = pd.read_csv(tf)
if "_tag" in df.columns:
df = df[df["_tag"].fillna("") == ""]
- _fused_moe_module.cfg_2stages = df.set_index(_INDEX_COLS).to_dict("index")
+ _fused_moe_module.cfg_2stages = df.set_index(_IDX).to_dict("index")
else:
_fused_moe_module.cfg_2stages = {}
-
_fused_moe_module.cfg_2stages.update(_CUSTOM_CONFIGS)
- # Monkeypatch get_2stage_cfgs to support use_non_temporal_load from config
- _original_get_2stage_cfgs = _fused_moe_module.get_2stage_cfgs
- @functools.lru_cache(maxsize=2048)
- def _patched_get_2stage_cfgs(
- token, model_dim, inter_dim, expert, topk,
- dtype, q_dtype_a, q_dtype_w, q_type, use_g1u1,
- activation, doweight_stage1, hidden_pad, intermediate_pad, is_shuffled=True,
- ):
- metadata = _original_get_2stage_cfgs(
- token, model_dim, inter_dim, expert, topk,
- dtype, q_dtype_a, q_dtype_w, q_type, use_g1u1,
- activation, doweight_stage1, hidden_pad, intermediate_pad, is_shuffled,
- )
- from aiter.jit.utils.chip_info import get_cu_num
- cu_num = get_cu_num()
- keys = (
- cu_num, token, model_dim, inter_dim, expert, topk,
- str(activation), str(dtype), str(q_dtype_a), str(q_dtype_w),
- str(q_type), use_g1u1, doweight_stage1,
- )
- cfg = _fused_moe_module.cfg_2stages.get(keys)
- if cfg and cfg.get("use_non_temporal_load") is not None:
- nt = cfg["use_non_temporal_load"]
- old_s1 = metadata.stage1
- if hasattr(old_s1, 'func') and old_s1.func is not None:
- if 'use_non_temporal_load' in (old_s1.keywords or {}):
- new_kw = dict(old_s1.keywords)
- new_kw['use_non_temporal_load'] = nt
- metadata = _fused_moe_module.MOEMetadata(
- functools.partial(old_s1.func, **{k: v for k, v in new_kw.items()}),
- metadata.stage2,
- metadata.block_m,
- metadata.ksplit,
- metadata.run_1stage,
- metadata.has_bias,
- nt,
- )
- old_s2 = metadata.stage2
- if old_s2 and hasattr(old_s2, 'keywords') and 'use_non_temporal_load' in (old_s2.keywords or {}):
- new_kw2 = dict(old_s2.keywords)
- new_kw2['use_non_temporal_load'] = nt
- metadata = _fused_moe_module.MOEMetadata(
- metadata.stage1,
- functools.partial(old_s2.func, **{k: v for k, v in new_kw2.items()}),
- metadata.block_m,
- metadata.ksplit,
- metadata.run_1stage,
- metadata.has_bias,
- nt,
- )
- return metadata
-
- _fused_moe_module.get_2stage_cfgs = _patched_get_2stage_cfgs
-
-
def custom_kernel(data: input_t) -> output_t:
(
- hidden_states, gate_up_weight, down_weight,
- gate_up_weight_scale, down_weight_scale,
- gate_up_weight_shuffled, down_weight_shuffled,
- gate_up_weight_scale_shuffled, down_weight_scale_shuffled,
+ hidden_states, _w1r, _w2r, _w1sr, _w2sr,
+ w1, w2, w1s, w2s,
topk_weights, topk_ids, config,
) = data
- _inject_configs()
+ _inject()
- hidden_pad = config["d_hidden_pad"] - config["d_hidden"]
- intermediate_pad = config["d_expert_pad"] - config["d_expert"]
-
M = hidden_states.shape[0]
topk = topk_ids.shape[1]
- device = topk_ids.device
- w1 = gate_up_weight_shuffled
- w2 = down_weight_shuffled
- E, model_dim, inter_dim = get_inter_dim(w1.shape, w2.shape)
+ device = hidden_states.device
+ dhp = config["d_hidden_pad"]
+ dep = config["d_expert_pad"]
+ hidden_pad = dhp - config["d_hidden"]
+ intermediate_pad = dep - config["d_expert"]
+ E, model_dim, inter_dim = get_inter_dim(w1.shape, w2.shape)
padded_M = get_padded_M(M)
+
metadata = get_2stage_cfgs(
padded_M, model_dim, inter_dim, E, topk,
torch.bfloat16, dtypes.fp4x2, dtypes.fp4x2,
QuantType.per_1x32, True, ActivationType.Silu,
False, hidden_pad, intermediate_pad, True,
)
+ block_m = int(metadata.block_m)
- block_size_M = int(metadata.block_m)
-
- # === Pre-allocated moe_sorting ===
- bufs = _get_or_alloc_sorting_buffers(M, E, topk, model_dim, block_size_M, device)
- sorted_ids = bufs["sorted_ids"]
- sorted_weights = bufs["sorted_weights"]
- sorted_expert_ids = bufs["sorted_expert_ids"]
- num_valid_ids = bufs["num_valid_ids"]
- moe_out = bufs["moe_buf"]
-
+ b = _get_sorting_bufs(M, E, topk, model_dim, block_m, device)
aiter.moe_sorting_fwd(
topk_ids, topk_weights,
- sorted_ids, sorted_weights, sorted_expert_ids, num_valid_ids, moe_out,
- E, int(block_size_M), None, None, 0,
+ b["sid"], b["sw"], b["se"], b["nv"], b["out"],
+ E, block_m, None, None, 0,
)
- # === Inline 2-stage pipeline ===
- token_num = M
+ w1sv = w1s.view(dtypes.fp8_e8m0)
+ w2sv = w2s.view(dtypes.fp8_e8m0)
+ a2_buf = _get_a2(M, topk, inter_dim, device)
if metadata.ksplit > 1:
- # cktile_moe path: bf16 activations, no fp4 quant
a1 = hidden_states.to(torch.bfloat16)
- a1_scale = None
- w1_scale_view = gate_up_weight_scale_shuffled.view(dtypes.fp8_e8m0)
- w2_scale_view = down_weight_scale_shuffled.view(dtypes.fp8_e8m0)
-
a2 = metadata.stage1(
- a1, w1, w2,
- sorted_ids, sorted_expert_ids, num_valid_ids,
- _get_or_alloc_a2(M, topk, inter_dim, device), # pre-allocated
- topk,
- block_m=block_size_M,
- a1_scale=a1_scale,
- w1_scale=w1_scale_view,
- sorted_weights=None,
+ a1, w1, w2, b["sid"], b["se"], b["nv"], a2_buf, topk,
+ block_m=block_m, a1_scale=None, w1_scale=w1sv, sorted_weights=None,
)
-
- # cktile_moe stage2: a2 is bf16, no inter-stage requant
- a2_scale = None
metadata.stage2(
- a2, w1, w2,
- sorted_ids, sorted_expert_ids, num_valid_ids,
- moe_out, topk,
- w2_scale=w2_scale_view,
- a2_scale=a2_scale,
- block_m=block_size_M,
- sorted_weights=sorted_weights,
+ a2, w1, w2, b["sid"], b["se"], b["nv"], b["out"], topk,
+ w2_scale=w2sv, a2_scale=None, block_m=block_m, sorted_weights=b["sw"],
)
else:
- # CK 2-stage path: fp4 activation quant with pre-allocated buffers
- w1_scale_view = gate_up_weight_scale_shuffled.view(dtypes.fp8_e8m0)
- w2_scale_view = down_weight_scale_shuffled.view(dtypes.fp8_e8m0)
-
- # Stage 1: quant activations + gate_up GEMM + SwiGLU
- a1, a1_scale = _quant_prealloc(
- hidden_states, sorted_ids, num_valid_ids,
- token_num, 1, block_size_M, device,
+ a1, a1s = _quant_prealloc(
+ hidden_states, b["sid"], b["nv"], M, 1, block_m, device,
)
-
- a2 = _get_or_alloc_a2(M, topk, inter_dim, device)
a2 = metadata.stage1(
- a1, w1, w2,
- sorted_ids, sorted_expert_ids, num_valid_ids,
- a2, topk,
- block_m=block_size_M,
- a1_scale=a1_scale,
- w1_scale=w1_scale_view,
- sorted_weights=None,
+ a1, w1, w2, b["sid"], b["se"], b["nv"], a2_buf, topk,
+ block_m=block_m, a1_scale=a1s, w1_scale=w1sv, sorted_weights=None,
)
-
- # Inter-stage requant: bf16 -> fp4 with pre-allocated buffers
a2_flat = a2.view(-1, inter_dim)
- a2_quant, a2_scale = _quant_prealloc(
- a2_flat, sorted_ids, num_valid_ids,
- token_num, topk, block_size_M, device,
+ a2q, a2s = _quant_prealloc(
+ a2_flat, b["sid"], b["nv"], M, topk, block_m, device,
)
- a2_quant = a2_quant.view(token_num, topk, -1)
-
- # Stage 2: down GEMM + weighted reduction
+ a2q = a2q.view(M, topk, -1)
metadata.stage2(
- a2_quant, w1, w2,
- sorted_ids, sorted_expert_ids, num_valid_ids,
- moe_out, topk,
- w2_scale=w2_scale_view,
- a2_scale=a2_scale,
- block_m=block_size_M,
- sorted_weights=sorted_weights,
+ a2q, w1, w2, b["sid"], b["se"], b["nv"], b["out"], topk,
+ w2_scale=w2sv, a2_scale=a2s, block_m=block_m, sorted_weights=b["sw"],
)
- return moe_out
+ return b["out"]
scrolls · 519 diff lines total

Best evidence level for this revision: reported

JSON