deepseek_v2.py 57.6 KB
Newer Older
1
# SPDX-License-Identifier: Apache-2.0
2
# SPDX-FileCopyrightText: Copyright contributors to the vLLM project
3

wangding zeng's avatar
wangding zeng committed
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
# Adapted from
# https://github.com/huggingface/transformers/blob/v4.28.0/src/transformers/models/llama/modeling_llama.py
# Copyright 2023 The vLLM team.
# Copyright 2023 DeepSeek-AI and the HuggingFace Inc. team. All rights reserved.
#
# This code is based on EleutherAI's GPT-NeoX library and the GPT-NeoX
# and OPT implementations in this library. It has been modified from its
# original forms to accommodate minor architectural differences compared
# to GPT-NeoX and OPT used by the Meta AI team that trained the model.
#
# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at
#
#     http://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and
# limitations under the License.
25
"""Inference-only DeepseekV2/DeepseekV3 model."""
26

27
28
import typing
from collections.abc import Callable, Iterable
29
from itertools import islice
wangding zeng's avatar
wangding zeng committed
30
31
32

import torch
from torch import nn
33
from transformers import DeepseekV2Config, DeepseekV3Config
wangding zeng's avatar
wangding zeng committed
34

35
from vllm._aiter_ops import rocm_aiter_ops
36
from vllm.compilation.decorators import support_torch_compile
37
38
39
40
41
42
43
44
from vllm.config import CacheConfig, ParallelConfig, VllmConfig, get_current_vllm_config
from vllm.distributed import (
    get_ep_group,
    get_pp_group,
    get_tensor_model_parallel_rank,
    get_tensor_model_parallel_world_size,
    tensor_model_parallel_all_gather,
)
45
from vllm.logger import init_logger
wangding zeng's avatar
wangding zeng committed
46
from vllm.model_executor.layers.activation import SiluAndMul
47
from vllm.model_executor.layers.attention import Attention
48
from vllm.model_executor.layers.attention_layer_base import AttentionLayerBase
49
from vllm.model_executor.layers.fused_moe import SharedFusedMoE
50
from vllm.model_executor.layers.layernorm import LayerNorm, RMSNorm
51
52
53
from vllm.model_executor.layers.linear import (
    ColumnParallelLinear,
    MergedColumnParallelLinear,
54
    QKVParallelLinear,
55
56
57
    ReplicatedLinear,
    RowParallelLinear,
)
wangding zeng's avatar
wangding zeng committed
58
from vllm.model_executor.layers.logits_processor import LogitsProcessor
59
from vllm.model_executor.layers.mla import MLAModules, MultiHeadLatentAttentionWrapper
60
from vllm.model_executor.layers.quantization import QuantizationConfig
61
from vllm.model_executor.layers.quantization.utils.fp8_utils import (
62
63
    per_token_group_quant_fp8,
)
wangding zeng's avatar
wangding zeng committed
64
from vllm.model_executor.layers.rotary_embedding import get_rope
65
from vllm.model_executor.layers.sparse_attn_indexer import SparseAttnIndexer
wangding zeng's avatar
wangding zeng committed
66
from vllm.model_executor.layers.vocab_parallel_embedding import (
67
68
69
    ParallelLMHead,
    VocabParallelEmbedding,
)
70
from vllm.model_executor.model_loader.weight_utils import (
71
72
73
    default_weight_loader,
    maybe_remap_kv_scale_name,
)
74
from vllm.model_executor.models.utils import sequence_parallel_chunk
75
from vllm.platforms import current_platform
76
from vllm.sequence import IntermediateTensors
77
from vllm.v1.attention.backend import AttentionBackend
78
79
80
from vllm.v1.attention.backends.mla.indexer import (
    DeepseekV32IndexerBackend,
)
81
from vllm.v1.kv_cache_interface import KVCacheSpec, MLAAttentionSpec
wangding zeng's avatar
wangding zeng committed
82

83
from .interfaces import MixtureOfExperts, SupportsEagle, SupportsLoRA, SupportsPP
84
85
86
87
88
89
90
from .utils import (
    PPMissingLayer,
    is_pp_missing_parameter,
    make_empty_intermediate_tensors_factory,
    make_layers,
    maybe_prefix,
)
91

92
93
logger = init_logger(__name__)

wangding zeng's avatar
wangding zeng committed
94

95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
class DeepseekAttention(nn.Module):
    """Normal MHA implementation used by Deepseek v1."""

    def __init__(
        self,
        vllm_config: VllmConfig,
        config: DeepseekV2Config | DeepseekV3Config,
        hidden_size: int,
        num_heads: int,
        max_position_embeddings: int = 8192,
        cache_config: CacheConfig | None = None,
        quant_config: QuantizationConfig | None = None,
        prefix: str = "",
        **kwargs,
    ) -> None:
        super().__init__()
        self.hidden_size = hidden_size
        tp_size = get_tensor_model_parallel_world_size()
        self.total_num_heads = num_heads
        assert self.total_num_heads % tp_size == 0
        self.num_heads = self.total_num_heads // tp_size
        self.total_num_kv_heads = config.num_key_value_heads
        if self.total_num_kv_heads >= tp_size:
            # Number of KV heads is greater than TP size, so we partition
            # the KV heads across multiple tensor parallel GPUs.
            assert self.total_num_kv_heads % tp_size == 0
        else:
            # Number of KV heads is less than TP size, so we replicate
            # the KV heads across multiple tensor parallel GPUs.
            assert tp_size % self.total_num_kv_heads == 0
        self.num_kv_heads = max(1, self.total_num_kv_heads // tp_size)
        self.head_dim = hidden_size // self.total_num_heads
        self.q_size = self.num_heads * self.head_dim
        self.kv_size = self.num_kv_heads * self.head_dim
        self.scaling = self.head_dim**-0.5
        self.max_position_embeddings = max_position_embeddings

        self.qkv_proj = QKVParallelLinear(
            hidden_size,
            self.head_dim,
            self.total_num_heads,
            self.total_num_kv_heads,
            bias=False,
            quant_config=quant_config,
        )

        self.o_proj = RowParallelLinear(
            self.total_num_heads * self.head_dim,
            hidden_size,
            bias=False,
            quant_config=quant_config,
        )

        self.rotary_emb = get_rope(
            self.head_dim,
            max_position=max_position_embeddings,
151
            rope_parameters=config.rope_parameters,
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
        )
        self.attn = Attention(
            self.num_heads,
            self.head_dim,
            self.scaling,
            num_kv_heads=self.num_kv_heads,
            cache_config=cache_config,
            quant_config=quant_config,
            prefix=f"{prefix}.attn",
        )

    def forward(
        self,
        positions: torch.Tensor,
        hidden_states: torch.Tensor,
    ) -> torch.Tensor:
        qkv, _ = self.qkv_proj(hidden_states)
        q, k, v = qkv.split([self.q_size, self.kv_size, self.kv_size], dim=-1)
        q, k = self.rotary_emb(positions, q, k)
        attn_output = self.attn(q, k, v)
        output, _ = self.o_proj(attn_output)
        return output


wangding zeng's avatar
wangding zeng committed
176
177
178
179
180
181
class DeepseekV2MLP(nn.Module):
    def __init__(
        self,
        hidden_size: int,
        intermediate_size: int,
        hidden_act: str,
182
        quant_config: QuantizationConfig | None = None,
wangding zeng's avatar
wangding zeng committed
183
        reduce_results: bool = True,
184
        is_sequence_parallel=False,
185
        prefix: str = "",
wangding zeng's avatar
wangding zeng committed
186
187
    ) -> None:
        super().__init__()
188
189
190
191
192

        # If is_sequence_parallel, the input and output tensors are sharded
        # across the ranks within the tp_group. In this case the weights are
        # replicated and no collective ops are needed.
        # Otherwise we use standard TP with an allreduce at the end.
wangding zeng's avatar
wangding zeng committed
193
        self.gate_up_proj = MergedColumnParallelLinear(
194
195
            hidden_size,
            [intermediate_size] * 2,
wangding zeng's avatar
wangding zeng committed
196
            bias=False,
197
            quant_config=quant_config,
198
            disable_tp=is_sequence_parallel,
199
200
201
202
203
204
205
206
207
208
209
            prefix=f"{prefix}.gate_up_proj",
        )
        self.down_proj = RowParallelLinear(
            intermediate_size,
            hidden_size,
            bias=False,
            quant_config=quant_config,
            reduce_results=reduce_results,
            disable_tp=is_sequence_parallel,
            prefix=f"{prefix}.down_proj",
        )
wangding zeng's avatar
wangding zeng committed
210
        if hidden_act != "silu":
211
212
213
            raise ValueError(
                f"Unsupported activation: {hidden_act}. Only silu is supported for now."
            )
wangding zeng's avatar
wangding zeng committed
214
215
216
217
218
219
220
221
222
223
224
225
        self.act_fn = SiluAndMul()

    def forward(self, x):
        gate_up, _ = self.gate_up_proj(x)
        x = self.act_fn(gate_up)
        x, _ = self.down_proj(x)
        return x


class DeepseekV2MoE(nn.Module):
    def __init__(
        self,
226
        config: DeepseekV2Config | DeepseekV3Config,
227
        parallel_config: ParallelConfig,
228
        quant_config: QuantizationConfig | None = None,
229
        prefix: str = "",
wangding zeng's avatar
wangding zeng committed
230
231
232
    ):
        super().__init__()
        self.tp_size = get_tensor_model_parallel_world_size()
233
234
        self.tp_rank = get_tensor_model_parallel_rank()

235
        self.routed_scaling_factor = getattr(config, "routed_scaling_factor", 1.0)
236
237

        self.ep_group = get_ep_group().device_group
238
        self.ep_rank = get_ep_group().rank_in_group
239
240
241
        self.ep_size = self.ep_group.size()
        self.n_routed_experts: int = config.n_routed_experts
        self.n_shared_experts: int = config.n_shared_experts
242

243
        self.is_sequence_parallel = parallel_config.use_sequence_parallel_moe
244

245
        if config.hidden_act != "silu":
246
247
248
249
250
251
252
253
254
255
256
257
            raise ValueError(
                f"Unsupported activation: {config.hidden_act}. "
                "Only silu is supported for now."
            )

        self.gate = ReplicatedLinear(
            config.hidden_size,
            config.n_routed_experts,
            bias=False,
            quant_config=None,
            prefix=f"{prefix}.gate",
        )
258
        if getattr(config, "topk_method", None) == "noaux_tc":
259
            self.gate.e_score_correction_bias = nn.Parameter(
260
261
                torch.empty(config.n_routed_experts, dtype=torch.float32)
            )
262
263
264
        else:
            self.gate.e_score_correction_bias = None

265
        # Load balancing settings.
266
267
        eplb_config = parallel_config.eplb_config
        self.enable_eplb = parallel_config.enable_eplb
268

269
        self.n_redundant_experts = eplb_config.num_redundant_experts
270
        self.n_logical_experts = self.n_routed_experts
271
        self.n_physical_experts = self.n_logical_experts + self.n_redundant_experts
272
273
        self.n_local_physical_experts = self.n_physical_experts // self.ep_size

274
275
276
277
        self.physical_expert_start = self.ep_rank * self.n_local_physical_experts
        self.physical_expert_end = (
            self.physical_expert_start + self.n_local_physical_experts
        )
278

279
        self.is_rocm_aiter_moe_enabled = rocm_aiter_ops.is_fused_moe_enabled()
280
281
282
283
        self.is_fusion_moe_shared_experts_enabled = (
            rocm_aiter_ops.is_fusion_moe_shared_experts_enabled()
        )
        if config.n_shared_experts is None or self.is_fusion_moe_shared_experts_enabled:
284
285
            self.shared_experts = None
        else:
286
            intermediate_size = config.moe_intermediate_size * config.n_shared_experts
287

wangding zeng's avatar
wangding zeng committed
288
289
290
291
292
            self.shared_experts = DeepseekV2MLP(
                hidden_size=config.hidden_size,
                intermediate_size=intermediate_size,
                hidden_act=config.hidden_act,
                quant_config=quant_config,
293
                is_sequence_parallel=self.is_sequence_parallel,
294
                reduce_results=False,
295
                prefix=f"{prefix}.shared_experts",
wangding zeng's avatar
wangding zeng committed
296
297
            )

298
299
300
301
302
303
304
305
        n_group = getattr(config, "n_group", 1)
        topk_group = getattr(config, "topk_group", 1)
        use_grouped_topk = True
        if (n_group, topk_group) == (1, 1):
            n_group = None
            topk_group = None
            use_grouped_topk = False

306
307
        self.experts = SharedFusedMoE(
            shared_experts=self.shared_experts,
308
            gate=self.gate,
309
310
311
312
313
314
315
            num_experts=config.n_routed_experts,
            top_k=config.num_experts_per_tok,
            hidden_size=config.hidden_size,
            intermediate_size=config.moe_intermediate_size,
            reduce_results=False,
            renormalize=config.norm_topk_prob,
            quant_config=quant_config,
316
317
318
            use_grouped_topk=use_grouped_topk,
            num_expert_group=n_group,
            topk_group=topk_group,
319
            prefix=f"{prefix}.experts",
320
            scoring_func=getattr(config, "scoring_func", "softmax"),
321
            # we do scaling outside, set factor to 1.0 to avoid double mul
322
323
            # aiter applies routed_scaling_factor internally
            routed_scaling_factor=1.0
324
            if not self.is_rocm_aiter_moe_enabled
325
            else self.routed_scaling_factor,
326
327
328
329
            e_score_correction_bias=self.gate.e_score_correction_bias,
            enable_eplb=self.enable_eplb,
            num_redundant_experts=self.n_redundant_experts,
            is_sequence_parallel=self.is_sequence_parallel,
330
            n_shared_experts=config.n_shared_experts
331
            if self.is_fusion_moe_shared_experts_enabled
332
            else None,
333
        )
334

wangding zeng's avatar
wangding zeng committed
335
336
337
    def forward(self, hidden_states: torch.Tensor) -> torch.Tensor:
        num_tokens, hidden_dim = hidden_states.shape
        hidden_states = hidden_states.view(-1, hidden_dim)
338
339
340
341
342
343

        # Chunk the hidden states so they aren't replicated across TP ranks.
        # This avoids duplicate computation in self.experts.
        # TODO: We can replace the all_reduce at the end of attn with a
        # reduce_scatter instead of chunking here.
        if self.is_sequence_parallel:
344
            hidden_states = sequence_parallel_chunk(hidden_states)
345

346
347
348
349
350
351
352
353
354
355
356
        if self.experts.is_internal_router:
            # In this case, the gate/router runs inside the FusedMoE class
            fused_moe_out = self.experts(
                hidden_states=hidden_states, router_logits=hidden_states
            )
        else:
            # router_logits: (num_tokens, n_experts)
            router_logits, _ = self.gate(hidden_states)
            fused_moe_out = self.experts(
                hidden_states=hidden_states, router_logits=router_logits
            )
357

358
359
360
        shared_output, final_hidden_states = fused_moe_out
        if self.shared_experts is None:
            assert shared_output is None
361
362
363
364

        # Fix FP16 overflow
        # See DeepseekV2DecoderLayer for more details.
        if hidden_states.dtype != torch.float16:
365
            if not self.is_rocm_aiter_moe_enabled:
366
                final_hidden_states *= self.routed_scaling_factor
367
368
        elif self.shared_experts is not None:
            assert shared_output is not None
369
            shared_output *= 1.0 / self.routed_scaling_factor
370
371
372
373

        if self.shared_experts is not None:
            assert shared_output is not None
            final_hidden_states += shared_output
374

375
376
        if self.is_sequence_parallel:
            final_hidden_states = tensor_model_parallel_all_gather(
377
378
                final_hidden_states, 0
            )
379
380
            final_hidden_states = final_hidden_states[:num_tokens]
        elif self.tp_size > 1:
381
382
383
            final_hidden_states = self.experts.maybe_all_reduce_tensor_model_parallel(
                final_hidden_states
            )
wangding zeng's avatar
wangding zeng committed
384
385
386
387
388
389

        return final_hidden_states.view(num_tokens, hidden_dim)


def yarn_get_mscale(scale: float = 1, mscale: float = 1) -> float:
    import math
390

wangding zeng's avatar
wangding zeng committed
391
392
393
394
395
    if scale <= 1:
        return 1.0
    return 0.1 * mscale * math.log(scale) + 1.0


396
397
398
399
400
401
402
403
404
405
def _get_llama_4_scaling(
    original_max_position_embeddings: int, scaling_beta: float, positions: torch.Tensor
) -> torch.Tensor:
    scaling = 1 + scaling_beta * torch.log(
        1 + torch.floor(positions / original_max_position_embeddings)
    )
    # Broadcast over num_heads and head_dim
    return scaling[..., None, None]


wangding zeng's avatar
wangding zeng committed
406
407
408
class DeepseekV2Attention(nn.Module):
    def __init__(
        self,
409
        vllm_config: VllmConfig,
410
        config: DeepseekV2Config | DeepseekV3Config,
wangding zeng's avatar
wangding zeng committed
411
412
413
414
415
416
417
418
        hidden_size: int,
        num_heads: int,
        qk_nope_head_dim: int,
        qk_rope_head_dim: int,
        v_head_dim: int,
        q_lora_rank: int,
        kv_lora_rank: int,
        max_position_embeddings: int = 8192,
419
420
421
        cache_config: CacheConfig | None = None,
        quant_config: QuantizationConfig | None = None,
        topk_indices_buffer: torch.Tensor | None = None,
422
        prefix: str = "",
wangding zeng's avatar
wangding zeng committed
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
    ) -> None:
        super().__init__()
        self.hidden_size = hidden_size
        self.qk_nope_head_dim = qk_nope_head_dim
        self.qk_rope_head_dim = qk_rope_head_dim
        self.qk_head_dim = qk_nope_head_dim + qk_rope_head_dim
        self.v_head_dim = v_head_dim
        self.q_lora_rank = q_lora_rank
        self.kv_lora_rank = kv_lora_rank
        self.num_heads = num_heads
        tp_size = get_tensor_model_parallel_world_size()
        assert num_heads % tp_size == 0
        self.num_local_heads = num_heads // tp_size
        self.scaling = self.qk_head_dim**-0.5
        self.max_position_embeddings = max_position_embeddings
438
439
        assert topk_indices_buffer is None, (
            "topk_indices_buffer is not \
440
        supported for DeepseekV2Attention"
441
        )
wangding zeng's avatar
wangding zeng committed
442
443

        if self.q_lora_rank is not None:
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
            self.q_a_proj = ReplicatedLinear(
                self.hidden_size,
                self.q_lora_rank,
                bias=False,
                quant_config=quant_config,
                prefix=f"{prefix}.q_a_proj",
            )
            self.q_a_layernorm = RMSNorm(self.q_lora_rank, eps=config.rms_norm_eps)
            self.q_b_proj = ColumnParallelLinear(
                q_lora_rank,
                self.num_heads * self.qk_head_dim,
                bias=False,
                quant_config=quant_config,
                prefix=f"{prefix}.q_b_proj",
            )
wangding zeng's avatar
wangding zeng committed
459
        else:
460
461
462
463
464
465
466
            self.q_proj = ColumnParallelLinear(
                self.hidden_size,
                self.num_heads * self.qk_head_dim,
                bias=False,
                quant_config=quant_config,
                prefix=f"{prefix}.q_proj",
            )
wangding zeng's avatar
wangding zeng committed
467

468
469
470
471
472
        self.kv_a_proj_with_mqa = ReplicatedLinear(
            self.hidden_size,
            self.kv_lora_rank + self.qk_rope_head_dim,
            bias=False,
            quant_config=quant_config,
473
474
475
            prefix=f"{prefix}.kv_a_proj_with_mqa",
        )
        self.kv_a_layernorm = RMSNorm(self.kv_lora_rank, eps=config.rms_norm_eps)
wangding zeng's avatar
wangding zeng committed
476
477
478
479
        self.kv_b_proj = ColumnParallelLinear(
            self.kv_lora_rank,
            self.num_heads * (self.qk_nope_head_dim + self.v_head_dim),
            bias=False,
480
            quant_config=quant_config,
481
482
            prefix=f"{prefix}.kv_b_proj",
        )
wangding zeng's avatar
wangding zeng committed
483
        # O projection.
484
485
486
487
488
489
490
        self.o_proj = RowParallelLinear(
            self.num_heads * self.v_head_dim,
            self.hidden_size,
            bias=False,
            quant_config=quant_config,
            prefix=f"{prefix}.o_proj",
        )
491
        if config.rope_parameters["rope_type"] != "default":
492
493
494
495
496
            config.rope_parameters["rope_type"] = (
                "deepseek_yarn"
                if config.rope_parameters.get("apply_yarn_scaling", True)
                else "deepseek_llama_scaling"
            )
497

498
499
500
        self.rotary_emb = get_rope(
            qk_rope_head_dim,
            max_position=max_position_embeddings,
501
            rope_parameters=config.rope_parameters,
502
503
            is_neox_style=False,
        )
wangding zeng's avatar
wangding zeng committed
504

505
506
507
508
        if (
            config.rope_parameters["rope_type"] != "default"
            and config.rope_parameters["rope_type"] == "deepseek_yarn"
        ):
509
510
            mscale_all_dim = config.rope_parameters.get("mscale_all_dim", False)
            scaling_factor = config.rope_parameters["factor"]
wangding zeng's avatar
wangding zeng committed
511
512
513
            mscale = yarn_get_mscale(scaling_factor, float(mscale_all_dim))
            self.scaling = self.scaling * mscale * mscale

514
515
516
517
518
519
520
521
522
        self.attn = Attention(
            self.num_local_heads,
            self.qk_head_dim,
            self.scaling,
            num_kv_heads=self.num_local_heads,
            cache_config=cache_config,
            quant_config=quant_config,
            prefix=f"{prefix}.attn",
        )
wangding zeng's avatar
wangding zeng committed
523
524
525
526
527

    def forward(
        self,
        positions: torch.Tensor,
        hidden_states: torch.Tensor,
528
        llama_4_scaling: torch.Tensor | None,
wangding zeng's avatar
wangding zeng committed
529
530
531
532
    ) -> torch.Tensor:
        if self.q_lora_rank is not None:
            q = self.q_a_proj(hidden_states)[0]
            q = self.q_a_layernorm(q)
533
            q = self.q_b_proj(q)[0].view(-1, self.num_local_heads, self.qk_head_dim)
wangding zeng's avatar
wangding zeng committed
534
        else:
535
536
537
538
            q = self.q_proj(hidden_states)[0].view(
                -1, self.num_local_heads, self.qk_head_dim
            )
        q_nope, q_pe = q.split([self.qk_nope_head_dim, self.qk_rope_head_dim], dim=-1)
wangding zeng's avatar
wangding zeng committed
539
        latent_cache = self.kv_a_proj_with_mqa(hidden_states)[0]
540
        kv_a, _ = latent_cache.split([self.kv_lora_rank, self.qk_rope_head_dim], dim=-1)
wangding zeng's avatar
wangding zeng committed
541
        latent_cache = latent_cache.unsqueeze(1)
542
        kv_a = self.kv_a_layernorm(kv_a)
wangding zeng's avatar
wangding zeng committed
543
        kv = self.kv_b_proj(kv_a)[0]
544
        kv = kv.view(-1, self.num_local_heads, self.qk_nope_head_dim + self.v_head_dim)
wangding zeng's avatar
wangding zeng committed
545
        k_nope, v = kv.split([self.qk_nope_head_dim, self.v_head_dim], dim=-1)
546
        k_pe = latent_cache[:, :, self.kv_lora_rank :]
547

wangding zeng's avatar
wangding zeng committed
548
        q_pe, k_pe = self.rotary_emb(positions, q_pe, k_pe)
549

550
        q[..., self.qk_nope_head_dim :] = q_pe
wangding zeng's avatar
wangding zeng committed
551
        k = torch.empty_like(q)
552
553
        k[..., : self.qk_nope_head_dim] = k_nope
        k[..., self.qk_nope_head_dim :] = k_pe
554
555
556
557
558

        # Apply llama 4 scaling if provided
        if llama_4_scaling is not None:
            q *= llama_4_scaling

559
560
        # padding value to qk_head_dim for alignment
        v = torch.nn.functional.pad(
561
562
            v, [0, self.qk_head_dim - self.v_head_dim], value=0
        ).view(-1, self.num_local_heads * self.qk_head_dim)
563
        attn_output = self.attn(q, k, v)
564
565
566
        attn_output = attn_output.view(-1, self.num_local_heads, self.qk_head_dim)[
            ..., : self.v_head_dim
        ].reshape(-1, self.num_local_heads * self.v_head_dim)
wangding zeng's avatar
wangding zeng committed
567
568
569
570
        output, _ = self.o_proj(attn_output)
        return output


571
class DeepseekV32IndexerCache(torch.nn.Module, AttentionLayerBase):
572
573
574
    def __init__(
        self, head_dim: int, dtype: torch.dtype, prefix: str, cache_config: CacheConfig
    ):
575
576
577
578
579
580
581
582
583
584
585
        super().__init__()
        self.kv_cache = [torch.tensor([])]
        self.head_dim = head_dim
        self.prefix = prefix
        self.cache_config = cache_config
        self.dtype = dtype
        compilation_config = get_current_vllm_config().compilation_config
        if prefix in compilation_config.static_forward_context:
            raise ValueError(f"Duplicate layer name: {prefix}")
        compilation_config.static_forward_context[prefix] = self

586
    def get_kv_cache_spec(self, vllm_config: VllmConfig) -> KVCacheSpec:
587
588
589
590
591
592
593
        return MLAAttentionSpec(  # Only has one vector instead of K + V
            block_size=self.cache_config.block_size,
            num_kv_heads=1,
            head_size=self.head_dim,
            dtype=self.dtype,
        )

594
    def forward(self): ...
595
596
597
598
599
600

    def get_attn_backend(self) -> AttentionBackend:
        return DeepseekV32IndexerBackend


class Indexer(nn.Module):
601
602
603
    def __init__(
        self,
        vllm_config: VllmConfig,
604
        config: DeepseekV2Config | DeepseekV3Config,
605
606
        hidden_size: int,
        q_lora_rank: int,
607
608
609
        quant_config: QuantizationConfig | None,
        cache_config: CacheConfig | None,
        topk_indices_buffer: torch.Tensor | None,
610
611
        prefix: str = "",
    ):
612
613
614
615
616
617
618
619
620
621
        super().__init__()
        self.vllm_config = vllm_config
        self.config = config
        # self.indexer_cfg = config.attn_module_list_cfg[0]["attn_index"]
        self.topk_tokens = config.index_topk
        self.n_head = config.index_n_heads  # 64
        self.head_dim = config.index_head_dim  # 128
        self.rope_dim = config.qk_rope_head_dim  # 64
        self.q_lora_rank = q_lora_rank  # 1536
        # no tensor parallel, just replicated
622
623
624
625
626
627
628
629
630
631
632
633
634
635
        self.wq_b = ReplicatedLinear(
            self.q_lora_rank,
            self.head_dim * self.n_head,
            bias=False,
            quant_config=quant_config,
            prefix=f"{prefix}.wq_b",
        )
        self.wk = ReplicatedLinear(
            hidden_size,
            self.head_dim,
            bias=False,
            quant_config=quant_config,
            prefix=f"{prefix}.wk",
        )
636
        self.k_norm = LayerNorm(self.head_dim, eps=1e-6)
637
        self.weights_proj = ReplicatedLinear(
638
639
640
641
642
            hidden_size,
            self.n_head,
            bias=False,
            quant_config=None,
            prefix=f"{prefix}.weights_proj",
643
        )
644
645
646
647
648
649
650
651
652
653
        self.softmax_scale = self.head_dim**-0.5

        self.scale_fmt = "ue8m0"
        self.quant_block_size = 128  # TODO: get from config
        self.topk_indices_buffer = topk_indices_buffer

        # NOTE: (zyongye) we use fp8 naive cache,
        #       where we store value in fp8 and scale in fp32
        #       per self.quant_block_size element
        self.k_cache = DeepseekV32IndexerCache(
654
            head_dim=self.head_dim + self.head_dim // self.quant_block_size * 4,
655
656
            dtype=torch.uint8,
            prefix=f"{prefix}.k_cache",
657
658
            cache_config=cache_config,
        )
659
660
        self.max_model_len = vllm_config.model_config.max_model_len
        self.prefix = prefix
661
662
        from vllm.v1.attention.backends.mla.indexer import get_max_prefill_buffer_size

663
        self.max_total_seq_len = get_max_prefill_buffer_size(vllm_config)
664
665
666
667
668
669
670
671
672
673
        self.indexer_op = SparseAttnIndexer(
            self.k_cache,
            self.quant_block_size,
            self.scale_fmt,
            self.topk_tokens,
            self.head_dim,
            self.max_model_len,
            self.max_total_seq_len,
            self.topk_indices_buffer,
        )
674

675
676
677
    def forward(
        self, hidden_states: torch.Tensor, qr: torch.Tensor, positions, rotary_emb
    ) -> torch.Tensor:
678
679
680
        q, _ = self.wq_b(qr)
        q = q.view(-1, self.n_head, self.head_dim)
        q_pe, q_nope = torch.split(
681
682
            q, [self.rope_dim, self.head_dim - self.rope_dim], dim=-1
        )
683
684
685
686

        k, _ = self.wk(hidden_states)
        k = self.k_norm(k)
        k_pe, k_nope = torch.split(
687
688
            k, [self.rope_dim, self.head_dim - self.rope_dim], dim=-1
        )
689
690

        q_pe, k_pe = rotary_emb(positions, q_pe, k_pe.unsqueeze(1))
691
692
693
694
695
        # Note: RoPE (NeoX) can introduce extra leading dimensions during compilation
        # so we need to reshape back to token-flattened shapes
        q_pe = q_pe.reshape(-1, self.n_head, self.rope_dim)
        k_pe = k_pe.reshape(-1, 1, self.rope_dim)

696
697
        # `rotary_emb` is shape-preserving; `q_pe` is already
        # [num_tokens, n_head, rope_dim].
698
699
        q = torch.cat([q_pe, q_nope], dim=-1)
        # `k_pe` is [num_tokens, 1, rope_dim] (MQA).
700
        k = torch.cat([k_pe.squeeze(-2), k_nope], dim=-1)
701
702
703

        # we only quant q here since k quant is fused with cache insertion
        q = q.view(-1, self.head_dim)
704
705
706
707
708
709
        q_fp8, q_scale = per_token_group_quant_fp8(
            q,
            self.quant_block_size,
            column_major_scales=False,
            use_ue8m0=self.scale_fmt is not None,
        )
710
711
712
713
        q_fp8 = q_fp8.view(-1, self.n_head, self.head_dim)
        q_scale = q_scale.view(-1, self.n_head, 1)

        weights, _ = self.weights_proj(hidden_states)
714
715
716
        weights = (
            weights.unsqueeze(-1) * q_scale * self.softmax_scale * self.n_head**-0.5
        )
717
718
        weights = weights.squeeze(-1)

719
        return self.indexer_op(hidden_states, q_fp8, k, weights)
720
721


722
723
724
725
class DeepseekV2MLAAttention(nn.Module):
    """
    Main reference: DeepseekV2 paper, and FlashInfer Implementation
    (https://arxiv.org/abs/2405.04434 and https://github.com/flashinfer-ai/flashinfer/pull/551).
726

727
728
        For more info see MLACommonImpl in:
        vllm/v1/attention/backends/mla/utils.py
729
730
731
732
    """

    def __init__(
        self,
733
        vllm_config: VllmConfig,
734
        config: DeepseekV2Config | DeepseekV3Config,
735
736
737
738
739
        hidden_size: int,
        num_heads: int,
        qk_nope_head_dim: int,
        qk_rope_head_dim: int,
        v_head_dim: int,
740
        q_lora_rank: int | None,
741
742
        kv_lora_rank: int,
        max_position_embeddings: int = 8192,
743
744
        cache_config: CacheConfig | None = None,
        quant_config: QuantizationConfig | None = None,
745
        prefix: str = "",
746
        topk_indices_buffer: torch.Tensor | None = None,
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
    ) -> None:
        super().__init__()
        self.hidden_size = hidden_size
        self.qk_nope_head_dim = qk_nope_head_dim
        self.qk_rope_head_dim = qk_rope_head_dim
        self.qk_head_dim = qk_nope_head_dim + qk_rope_head_dim
        self.v_head_dim = v_head_dim

        self.q_lora_rank = q_lora_rank
        self.kv_lora_rank = kv_lora_rank

        self.num_heads = num_heads
        tp_size = get_tensor_model_parallel_world_size()
        assert num_heads % tp_size == 0
        self.num_local_heads = num_heads // tp_size

        self.scaling = self.qk_head_dim**-0.5
        self.max_position_embeddings = max_position_embeddings

        if self.q_lora_rank is not None:
767
            self.fused_qkv_a_proj = MergedColumnParallelLinear(
768
769
770
771
                self.hidden_size,
                [self.q_lora_rank, self.kv_lora_rank + self.qk_rope_head_dim],
                bias=False,
                quant_config=quant_config,
772
                prefix=f"{prefix}.fused_qkv_a_proj",
773
774
                disable_tp=True,
            )
775
776
777
778
779
780
        else:
            self.kv_a_proj_with_mqa = ReplicatedLinear(
                self.hidden_size,
                self.kv_lora_rank + self.qk_rope_head_dim,
                bias=False,
                quant_config=quant_config,
781
782
                prefix=f"{prefix}.kv_a_proj_with_mqa",
            )
783
784

        if self.q_lora_rank is not None:
785
786
787
788
789
790
791
792
            self.q_a_layernorm = RMSNorm(self.q_lora_rank, eps=config.rms_norm_eps)
            self.q_b_proj = ColumnParallelLinear(
                self.q_lora_rank,
                self.num_heads * self.qk_head_dim,
                bias=False,
                quant_config=quant_config,
                prefix=f"{prefix}.q_b_proj",
            )
793
        else:
794
795
796
797
798
799
800
801
            self.q_proj = ColumnParallelLinear(
                self.hidden_size,
                self.num_heads * self.qk_head_dim,
                bias=False,
                quant_config=quant_config,
                prefix=f"{prefix}.q_proj",
            )
        self.kv_a_layernorm = RMSNorm(self.kv_lora_rank, eps=config.rms_norm_eps)
802
803
804
805
806
        self.kv_b_proj = ColumnParallelLinear(
            self.kv_lora_rank,
            self.num_heads * (self.qk_nope_head_dim + self.v_head_dim),
            bias=False,
            quant_config=quant_config,
807
808
809
810
811
812
813
814
815
            prefix=f"{prefix}.kv_b_proj",
        )
        self.o_proj = RowParallelLinear(
            self.num_heads * self.v_head_dim,
            self.hidden_size,
            bias=False,
            quant_config=quant_config,
            prefix=f"{prefix}.o_proj",
        )
816

817
        if config.rope_parameters["rope_type"] != "default":
818
819
820
821
822
823
            config.rope_parameters["rope_type"] = (
                "deepseek_yarn"
                if config.rope_parameters.get("apply_yarn_scaling", True)
                else "deepseek_llama_scaling"
            )

824
825
826
        self.rotary_emb = get_rope(
            qk_rope_head_dim,
            max_position=max_position_embeddings,
827
            rope_parameters=config.rope_parameters,
828
829
            is_neox_style=False,
        )
830
831
832
833
834

        if (
            config.rope_parameters["rope_type"] != "default"
            and config.rope_parameters["rope_type"] == "deepseek_yarn"
        ):
835
836
            mscale_all_dim = config.rope_parameters.get("mscale_all_dim", False)
            scaling_factor = config.rope_parameters["factor"]
837
838
839
            mscale = yarn_get_mscale(scaling_factor, float(mscale_all_dim))
            self.scaling = self.scaling * mscale * mscale

840
841
842
        self.is_v32 = hasattr(config, "index_topk")

        if self.is_v32:
843
844
845
            self.indexer_rope_emb = get_rope(
                qk_rope_head_dim,
                max_position=max_position_embeddings,
846
                rope_parameters=config.rope_parameters,
847
848
                is_neox_style=True,
            )
849
850
851
852
853
854
855
856
857
858
            self.indexer = Indexer(
                vllm_config,
                config,
                hidden_size,
                q_lora_rank,
                quant_config,
                cache_config,
                topk_indices_buffer,
                f"{prefix}.indexer",
            )
859
        else:
860
            self.indexer_rope_emb = None
861
862
            self.indexer = None

863
864
        mla_modules = MLAModules(
            kv_a_layernorm=self.kv_a_layernorm,
865
            kv_b_proj=self.kv_b_proj,
866
867
868
            rotary_emb=self.rotary_emb,
            o_proj=self.o_proj,
            fused_qkv_a_proj=self.fused_qkv_a_proj
869
870
            if self.q_lora_rank is not None
            else None,
871
            kv_a_proj_with_mqa=self.kv_a_proj_with_mqa
872
873
874
            if self.q_lora_rank is None
            else None,
            q_a_layernorm=self.q_a_layernorm if self.q_lora_rank is not None else None,
875
876
            q_b_proj=self.q_b_proj if self.q_lora_rank is not None else None,
            q_proj=self.q_proj if self.q_lora_rank is None else None,
877
            indexer=self.indexer,
878
            indexer_rotary_emb=self.indexer_rope_emb,
879
880
            is_sparse=self.is_v32,
            topk_indices_buffer=topk_indices_buffer,
881
        )
882

883
        self.mla_attn = MultiHeadLatentAttentionWrapper(
884
885
886
887
888
889
890
891
892
893
894
895
            self.hidden_size,
            self.num_local_heads,
            self.scaling,
            self.qk_nope_head_dim,
            self.qk_rope_head_dim,
            self.v_head_dim,
            self.q_lora_rank,
            self.kv_lora_rank,
            mla_modules,
            cache_config,
            quant_config,
            prefix,
896
897
898
899
900
901
        )

    def forward(
        self,
        positions: torch.Tensor,
        hidden_states: torch.Tensor,
902
        llama_4_scaling: torch.Tensor | None,
903
    ) -> torch.Tensor:
904
        return self.mla_attn(positions, hidden_states, llama_4_scaling)
905
906


wangding zeng's avatar
wangding zeng committed
907
class DeepseekV2DecoderLayer(nn.Module):
908
909
910
911
    def __init__(
        self,
        vllm_config: VllmConfig,
        prefix: str,
912
913
        config: DeepseekV2Config | None = None,
        topk_indices_buffer: torch.Tensor | None = None,
914
    ) -> None:
wangding zeng's avatar
wangding zeng committed
915
        super().__init__()
916

917
918
        if config is None:
            config = vllm_config.model_config.hf_config
919
920
921
922
923
        model_config = vllm_config.model_config
        cache_config = vllm_config.cache_config
        quant_config = vllm_config.quant_config
        parallel_config = vllm_config.parallel_config

wangding zeng's avatar
wangding zeng committed
924
        self.hidden_size = config.hidden_size
925
        max_position_embeddings = getattr(config, "max_position_embeddings", 8192)
926
        moe_layer_freq = getattr(config, "moe_layer_freq", 1)
927
928
        # DecoderLayers are created with `make_layers` which passes the prefix
        # with the layer's index.
929
        layer_idx = int(prefix.split(sep=".")[-1])
930
        self.layer_idx = layer_idx
931
932
933
934
935
936
937
938
939
940

        # verify MLA attention specific fields
        qk_nope_head_dim = getattr(config, "qk_nope_head_dim", 0)
        qk_rope_head_dim = getattr(config, "qk_rope_head_dim", 0)
        v_head_dim = getattr(config, "v_head_dim", 0)
        kv_lora_rank = getattr(config, "kv_lora_rank", 0)
        use_mha = config.model_type == "deepseek" or all(
            dim == 0 for dim in (qk_nope_head_dim, qk_rope_head_dim)
        )

941
942
        self.use_mha = use_mha

943
944
945
        if use_mha:
            attn_cls = DeepseekAttention
        elif model_config.use_mla:
946
947
948
949
            attn_cls = DeepseekV2MLAAttention
        else:
            attn_cls = DeepseekV2Attention
        self.self_attn = attn_cls(
950
            vllm_config=vllm_config,
wangding zeng's avatar
wangding zeng committed
951
952
953
            config=config,
            hidden_size=self.hidden_size,
            num_heads=config.num_attention_heads,
954
955
956
            qk_nope_head_dim=qk_nope_head_dim,
            qk_rope_head_dim=qk_rope_head_dim,
            v_head_dim=v_head_dim,
957
            q_lora_rank=config.q_lora_rank if hasattr(config, "q_lora_rank") else None,
958
            kv_lora_rank=kv_lora_rank,
wangding zeng's avatar
wangding zeng committed
959
960
961
            max_position_embeddings=max_position_embeddings,
            cache_config=cache_config,
            quant_config=quant_config,
962
            prefix=f"{prefix}.self_attn",
963
            topk_indices_buffer=topk_indices_buffer,
wangding zeng's avatar
wangding zeng committed
964
        )
965

966
967
968
        if (
            config.n_routed_experts is not None
            and layer_idx >= config.first_k_dense_replace
969
            and layer_idx % moe_layer_freq == 0
970
        ):
971
972
            self.mlp = DeepseekV2MoE(
                config=config,
973
                parallel_config=parallel_config,
974
975
976
                quant_config=quant_config,
                prefix=f"{prefix}.mlp",
            )
wangding zeng's avatar
wangding zeng committed
977
978
979
980
981
982
        else:
            self.mlp = DeepseekV2MLP(
                hidden_size=config.hidden_size,
                intermediate_size=config.intermediate_size,
                hidden_act=config.hidden_act,
                quant_config=quant_config,
983
                prefix=f"{prefix}.mlp",
wangding zeng's avatar
wangding zeng committed
984
            )
985
986
987
988
        self.input_layernorm = RMSNorm(config.hidden_size, eps=config.rms_norm_eps)
        self.post_attention_layernorm = RMSNorm(
            config.hidden_size, eps=config.rms_norm_eps
        )
989
        self.routed_scaling_factor = getattr(config, "routed_scaling_factor", 1.0)
wangding zeng's avatar
wangding zeng committed
990
991
992
993
994

    def forward(
        self,
        positions: torch.Tensor,
        hidden_states: torch.Tensor,
995
        residual: torch.Tensor | None,
996
        llama_4_scaling: torch.Tensor | None = None,
wangding zeng's avatar
wangding zeng committed
997
998
999
    ) -> torch.Tensor:
        # Self Attention
        if residual is None:
1000
            residual = hidden_states.clone()
wangding zeng's avatar
wangding zeng committed
1001
1002
            hidden_states = self.input_layernorm(hidden_states)
        else:
1003
            hidden_states, residual = self.input_layernorm(hidden_states, residual)
1004
1005
1006
1007
1008
1009
1010
1011

        attn_kwargs = {
            "positions": positions,
            "hidden_states": hidden_states,
        }
        if not self.use_mha:
            attn_kwargs["llama_4_scaling"] = llama_4_scaling
        hidden_states = self.self_attn(**attn_kwargs)
wangding zeng's avatar
wangding zeng committed
1012

1013
1014
1015
1016
        if (
            not isinstance(self.self_attn, DeepseekAttention)
            and hidden_states.dtype == torch.float16
        ):
1017
1018
1019
            # Fix FP16 overflow
            # We scale both hidden_states and residual before
            # rmsnorm, and rmsnorm result would not affect by scale.
1020
            hidden_states *= 1.0 / self.routed_scaling_factor
1021
1022
1023
            if self.layer_idx == 0:
                # The residual is shared by all layers, we only scale it on
                # first layer.
1024
                residual *= 1.0 / self.routed_scaling_factor
1025
1026

        # Fully Connected
1027
        hidden_states, residual = self.post_attention_layernorm(hidden_states, residual)
wangding zeng's avatar
wangding zeng committed
1028
        hidden_states = self.mlp(hidden_states)
1029

1030
        if isinstance(self.mlp, DeepseekV2MLP) and hidden_states.dtype == torch.float16:
1031
1032
1033
1034
1035
            # Fix FP16 overflow
            # Scaling the DeepseekV2MLP output, it is the input of
            # input_layernorm of next decoder layer.
            # The scaling of DeepseekV2MOE output would be done in the forward
            # of DeepseekV2MOE
1036
            hidden_states *= 1.0 / self.routed_scaling_factor
1037

wangding zeng's avatar
wangding zeng committed
1038
1039
1040
        return hidden_states, residual


1041
@support_torch_compile
wangding zeng's avatar
wangding zeng committed
1042
1043
1044
class DeepseekV2Model(nn.Module):
    fall_back_to_pt_during_load = False

1045
    def __init__(self, *, vllm_config: VllmConfig, prefix: str = ""):
wangding zeng's avatar
wangding zeng committed
1046
        super().__init__()
1047
1048
1049

        config = vllm_config.model_config.hf_config
        quant_config = vllm_config.quant_config
1050
        self.config = config
1051
        self.device = current_platform.device_type
1052

wangding zeng's avatar
wangding zeng committed
1053
        self.vocab_size = config.vocab_size
1054
1055
1056
1057
1058
1059
1060
        self.is_v32 = hasattr(config, "index_topk")
        if self.is_v32:
            topk_tokens = config.index_topk
            topk_indices_buffer = torch.empty(
                vllm_config.scheduler_config.max_num_batched_tokens,
                topk_tokens,
                dtype=torch.int32,
1061
                device=self.device,
1062
            )
1063
1064
        else:
            topk_indices_buffer = None
wangding zeng's avatar
wangding zeng committed
1065

1066
1067
1068
1069
        if get_pp_group().is_first_rank:
            self.embed_tokens = VocabParallelEmbedding(
                config.vocab_size,
                config.hidden_size,
1070
                quant_config=quant_config,
1071
1072
                prefix=f"{prefix}.embed_tokens",
            )
1073
1074
1075
1076
        else:
            self.embed_tokens = PPMissingLayer()
        self.start_layer, self.end_layer, self.layers = make_layers(
            config.num_hidden_layers,
1077
            lambda prefix: DeepseekV2DecoderLayer(
1078
                vllm_config, prefix, topk_indices_buffer=topk_indices_buffer
1079
1080
1081
            ),
            prefix=f"{prefix}.layers",
        )
1082
1083
1084
1085
1086

        if get_pp_group().is_last_rank:
            self.norm = RMSNorm(config.hidden_size, eps=config.rms_norm_eps)
        else:
            self.norm = PPMissingLayer()
1087
1088
1089
        self.make_empty_intermediate_tensors = make_empty_intermediate_tensors_factory(
            ["hidden_states", "residual"], config.hidden_size
        )
wangding zeng's avatar
wangding zeng committed
1090

1091
    def embed_input_ids(self, input_ids: torch.Tensor) -> torch.Tensor:
1092
1093
        return self.embed_tokens(input_ids)

wangding zeng's avatar
wangding zeng committed
1094
1095
    def forward(
        self,
1096
        input_ids: torch.Tensor | None,
wangding zeng's avatar
wangding zeng committed
1097
        positions: torch.Tensor,
1098
1099
1100
        intermediate_tensors: IntermediateTensors | None,
        inputs_embeds: torch.Tensor | None = None,
    ) -> torch.Tensor | IntermediateTensors:
1101
        if get_pp_group().is_first_rank:
1102
1103
1104
            if inputs_embeds is not None:
                hidden_states = inputs_embeds
            else:
1105
                hidden_states = self.embed_input_ids(input_ids)
1106
1107
1108
1109
1110
1111
            residual = None
        else:
            assert intermediate_tensors is not None
            hidden_states = intermediate_tensors["hidden_states"]
            residual = intermediate_tensors["residual"]

1112
1113
1114
1115
1116
1117
1118
1119
1120
1121
1122
1123
1124
1125
        # Compute llama 4 scaling once per forward pass if enabled
        llama_4_scaling_config = getattr(self.config, "llama_4_scaling", None)
        llama_4_scaling: torch.Tensor | None
        if llama_4_scaling_config is not None:
            llama_4_scaling = _get_llama_4_scaling(
                original_max_position_embeddings=llama_4_scaling_config[
                    "original_max_position_embeddings"
                ],
                scaling_beta=llama_4_scaling_config["beta"],
                positions=positions,
            )
        else:
            llama_4_scaling = None

1126
        for layer in islice(self.layers, self.start_layer, self.end_layer):
1127
1128
1129
            hidden_states, residual = layer(
                positions, hidden_states, residual, llama_4_scaling
            )
1130
1131

        if not get_pp_group().is_last_rank:
1132
1133
1134
            return IntermediateTensors(
                {"hidden_states": hidden_states, "residual": residual}
            )
1135

wangding zeng's avatar
wangding zeng committed
1136
1137
1138
        hidden_states, _ = self.norm(hidden_states, residual)
        return hidden_states

1139

1140
1141
1142
1143
1144
1145
1146
1147
1148
1149
1150
1151
1152
1153
1154
1155
1156
1157
1158
1159
1160
1161
1162
1163
1164
1165
1166
1167
1168
1169
1170
1171
1172
1173
1174
1175
1176
1177
1178
1179
1180
1181
class DeepseekV2MixtureOfExperts(MixtureOfExperts):
    moe_mlp_layers: list[DeepseekV2MoE]
    """
    List of MoE MLP layers in the model.
    """

    def extract_moe_parameters(self, example_moe: DeepseekV2MoE | None):
        if example_moe is None:
            self.num_moe_layers = 0
            self.num_expert_groups = 0
            self.num_logical_experts = 0
            self.num_physical_experts = 0
            self.num_local_physical_experts = 0
            self.num_routed_experts = 0
            self.num_shared_experts = 0
            self.num_redundant_experts = 0
            logger.warning("DeepSeekV2: No DeepseekV2MoE layer found in model.layers.")
        else:
            self.num_logical_experts = example_moe.n_logical_experts
            self.num_physical_experts = example_moe.n_physical_experts
            self.num_local_physical_experts = example_moe.n_local_physical_experts
            self.num_routed_experts = example_moe.n_routed_experts
            self.num_shared_experts = example_moe.n_shared_experts
            self.num_redundant_experts = example_moe.n_redundant_experts

    def update_physical_experts_metadata(
        self,
        num_physical_experts: int,
        num_local_physical_experts: int,
    ) -> None:
        assert self.num_local_physical_experts == num_local_physical_experts
        self.num_physical_experts = num_physical_experts
        self.num_local_physical_experts = num_local_physical_experts
        self.num_redundant_experts = num_physical_experts - self.num_logical_experts
        for moe in self.moe_mlp_layers:
            moe.n_local_physical_experts = num_local_physical_experts
            moe.n_physical_experts = num_physical_experts
            moe.n_redundant_experts = self.num_redundant_experts
            moe.experts.update_expert_map()


class DeepseekV2ForCausalLM(
1182
    nn.Module, SupportsPP, DeepseekV2MixtureOfExperts, SupportsLoRA, SupportsEagle
1183
):
1184
1185
1186
    packed_modules_mapping = {
        "gate_up_proj": ["gate_proj", "up_proj"],
    }
1187
    model_cls = DeepseekV2Model
1188
1189
1190
1191
1192
1193
1194

    def __init__(self, *, vllm_config: VllmConfig, prefix: str = ""):
        super().__init__()
        config = vllm_config.model_config.hf_config
        quant_config = vllm_config.quant_config
        self.config = config
        self.quant_config = quant_config
1195

1196
1197
1198
1199
1200
1201
1202
1203
1204
        qk_nope_head_dim = getattr(config, "qk_nope_head_dim", 0)
        qk_rope_head_dim = getattr(config, "qk_rope_head_dim", 0)
        self.use_mha = config.model_type == "deepseek" or all(
            dim == 0 for dim in (qk_nope_head_dim, qk_rope_head_dim)
        )

        if self.use_mha:
            self.packed_modules_mapping["qkv_proj"] = ["q_proj", "k_proj", "v_proj"]

1205
1206
1207
1208
        # `packed_modules_mapping` needs to be modified before
        # initializing DeepseekV2Model, as it is passed inplace to
        # quantization config init and may be used to select the
        # quant_method for relevant layers during initialization.
1209
1210
1211
        self.fuse_qkv_a_proj = (
            hasattr(config, "q_lora_rank") and config.q_lora_rank is not None
        )
1212
1213
1214
1215
1216
1217
        if self.fuse_qkv_a_proj:
            self.packed_modules_mapping["fused_qkv_a_proj"] = [
                "q_a_proj",
                "kv_a_proj_with_mqa",
            ]

1218
        self.model = self.model_cls(
1219
1220
            vllm_config=vllm_config, prefix=maybe_prefix(prefix, "model")
        )
1221
        if get_pp_group().is_last_rank:
1222
1223
1224
1225
1226
1227
            self.lm_head = ParallelLMHead(
                config.vocab_size,
                config.hidden_size,
                quant_config=quant_config,
                prefix=maybe_prefix(prefix, "lm_head"),
            )
1228
1229
1230
1231
        else:
            self.lm_head = PPMissingLayer()
        self.logits_processor = LogitsProcessor(config.vocab_size)
        self.make_empty_intermediate_tensors = (
1232
1233
            self.model.make_empty_intermediate_tensors
        )
1234
1235
1236
1237
1238
1239
1240
        # Set MoE hyperparameters
        self.num_moe_layers = (
            self.config.num_hidden_layers - self.config.first_k_dense_replace
        )
        self.set_moe_parameters()

    def set_moe_parameters(self):
1241
1242
        self.expert_weights = []

1243
        self.num_expert_groups = getattr(self.config, "n_group", 1)
1244

1245
1246
        self.moe_layers = []
        self.moe_mlp_layers = []
1247
        example_moe = None
1248
        for layer in self.model.layers:
1249
1250
1251
            if isinstance(layer, PPMissingLayer):
                continue

1252
1253
            assert isinstance(layer, DeepseekV2DecoderLayer)
            if isinstance(layer.mlp, DeepseekV2MoE):
1254
1255
                # Pick last one layer since the first ones may be dense layers.
                example_moe = layer.mlp
1256
                self.moe_mlp_layers.append(layer.mlp)
1257
1258
                self.moe_layers.append(layer.mlp.experts)

1259
        self.extract_moe_parameters(example_moe)
1260

1261
1262
    def embed_input_ids(self, input_ids: torch.Tensor) -> torch.Tensor:
        return self.model.embed_input_ids(input_ids)
1263
1264
1265

    def forward(
        self,
1266
        input_ids: torch.Tensor | None,
1267
        positions: torch.Tensor,
1268
1269
1270
        intermediate_tensors: IntermediateTensors | None = None,
        inputs_embeds: torch.Tensor | None = None,
    ) -> torch.Tensor | IntermediateTensors:
1271
1272
1273
        hidden_states = self.model(
            input_ids, positions, intermediate_tensors, inputs_embeds
        )
1274
1275
1276
1277
1278
        return hidden_states

    def compute_logits(
        self,
        hidden_states: torch.Tensor,
1279
    ) -> torch.Tensor | None:
1280
        logits = self.logits_processor(self.lm_head, hidden_states)
1281
1282
        return logits

1283
1284
1285
1286
    def get_expert_mapping(self) -> list[tuple[str, str, int, str]]:
        # Params for weights, fp8 weight scales, fp8 activation scales
        # (param_name, weight_name, expert_id, shard_id)
        return SharedFusedMoE.make_expert_params_mapping(
1287
            self,
1288
1289
1290
1291
1292
1293
1294
            ckpt_gate_proj_name="gate_proj",
            ckpt_down_proj_name="down_proj",
            ckpt_up_proj_name="up_proj",
            num_experts=self.config.n_routed_experts,
            num_redundant_experts=0,
        )

1295
    def load_weights(self, weights: Iterable[tuple[str, torch.Tensor]]) -> set[str]:
1296
1297
1298
        rocm_aiter_moe_shared_expert_enabled = (
            rocm_aiter_ops.is_fusion_moe_shared_experts_enabled()
        )
wangding zeng's avatar
wangding zeng committed
1299
1300
1301
1302
        stacked_params_mapping = [
            # (param_name, shard_name, shard_id)
            ("gate_up_proj", "gate_proj", 0),
            ("gate_up_proj", "up_proj", 1),
1303
1304
        ]
        mla_params_mapping = [
1305
1306
            ("fused_qkv_a_proj", "q_a_proj", 0),
            ("fused_qkv_a_proj", "kv_a_proj_with_mqa", 1),
wangding zeng's avatar
wangding zeng committed
1307
        ]
1308
1309
1310
1311
1312
1313
1314
1315
1316
        mha_params_mapping = [
            ("qkv_proj", "q_proj", "q"),
            ("qkv_proj", "k_proj", "k"),
            ("qkv_proj", "v_proj", "v"),
        ]
        if self.use_mha:
            stacked_params_mapping.extend(mha_params_mapping)
        else:
            stacked_params_mapping.extend(mla_params_mapping)
wangding zeng's avatar
wangding zeng committed
1317

1318
1319
        # Params for weights, fp8 weight scales, fp8 activation scales
        # (param_name, weight_name, expert_id, shard_id)
1320
        expert_params_mapping = SharedFusedMoE.make_expert_params_mapping(
1321
            self,
1322
1323
1324
            ckpt_gate_proj_name="gate_proj",
            ckpt_down_proj_name="down_proj",
            ckpt_up_proj_name="up_proj",
1325
1326
1327
            num_experts=self.config.n_routed_experts
            + (
                self.config.n_shared_experts
1328
                if rocm_aiter_moe_shared_expert_enabled
1329
1330
                else 0
            ),
1331
1332
            num_redundant_experts=self.num_redundant_experts,
        )
1333

wangding zeng's avatar
wangding zeng committed
1334
        params_dict = dict(self.named_parameters())
1335
        loaded_params: set[str] = set()
wangding zeng's avatar
wangding zeng committed
1336
        for name, loaded_weight in weights:
1337
1338
1339
            if "rotary_emb.inv_freq" in name:
                continue

1340
1341
1342
            spec_layer = get_spec_layer_idx_from_weight_name(self.config, name)
            if spec_layer is not None:
                continue  # skip spec decode layers for main model
1343

1344
1345
            is_fusion_moe_shared_experts_layer = (
                rocm_aiter_moe_shared_expert_enabled and ("mlp.shared_experts" in name)
1346
1347
            )

1348
            for param_name, weight_name, shard_id in stacked_params_mapping:
1349
                # Skip non-stacked layers and experts (experts handled below).
wangding zeng's avatar
wangding zeng committed
1350
1351
                if weight_name not in name:
                    continue
1352
1353
1354
1355
1356
1357
                # We have mlp.experts[0].gate_proj in the checkpoint.
                # Since we handle the experts below in expert_params_mapping,
                # we need to skip here BEFORE we update the name, otherwise
                # name will be updated to mlp.experts[0].gate_up_proj, which
                # will then be updated below in expert_params_mapping
                # for mlp.experts[0].gate_gate_up_proj, which breaks load.
1358
                if ("mlp.experts." in name) and name not in params_dict:
1359
                    continue
1360
                if is_fusion_moe_shared_experts_layer:
1361
                    continue
1362
                name_mapped = name.replace(weight_name, param_name)
1363
1364
1365

                # QKV fusion is optional, fall back to normal
                # weight loading if it's not enabled
1366
                # if go with fusion option, then update name
1367
1368
1369
                if (
                    param_name == "fused_qkv_a_proj"
                ) and name_mapped not in params_dict:
1370
                    continue
1371
1372
                else:
                    name = name_mapped
wangding zeng's avatar
wangding zeng committed
1373
1374
1375
                # Skip loading extra bias for GPTQ models.
                if name.endswith(".bias") and name not in params_dict:
                    continue
1376
1377
1378
1379

                if is_pp_missing_parameter(name, self):
                    continue

wangding zeng's avatar
wangding zeng committed
1380
1381
1382
1383
1384
                param = params_dict[name]
                weight_loader = param.weight_loader
                weight_loader(param, loaded_weight, shard_id)
                break
            else:
1385
                is_expert_weight = False
1386
1387
1388
1389
1390
1391
1392
1393
1394
1395

                # Special handling: when AITER fusion_shared_experts is enabled,
                # checkpoints may provide a single widened shared_experts tensor
                # without explicit expert indices
                # (e.g. ...mlp.shared_experts.gate_proj.weight).
                # For models with multiple shared experts, split that tensor
                # evenly into per-shared-expert slices and load them into
                # appended expert slots mlp.experts.{n_routed_experts + j}.*
                # accordingly.
                num_chunks = 1
1396
                if is_fusion_moe_shared_experts_layer:
1397
1398
1399
1400
                    num_chunks = getattr(self.config, "n_shared_experts", 1) or 1
                    # Determine split axis based on op type
                    # gate/up: ColumnParallel → split along dim 0
                    # down: RowParallel → split along dim 1
1401
1402
1403
1404
1405
                    split_dim = (
                        1
                        if ("down_proj.weight" in name and loaded_weight.ndim > 1)
                        else 0
                    )
1406
1407
1408
1409
                    total = loaded_weight.shape[split_dim]
                    assert total % num_chunks == 0, (
                        f"Shared expert weight dim {total} "
                        f"not divisible by num_chunks {num_chunks}"
1410
                    )
1411
1412
1413
1414
1415
1416
                    chunk_size = total // num_chunks

                for j in range(num_chunks):
                    chunk_name = name
                    weight_to_load = loaded_weight

1417
                    if is_fusion_moe_shared_experts_layer:
1418
1419
1420
1421
1422
                        chunk_slice = slice(j * chunk_size, (j + 1) * chunk_size)
                        if loaded_weight.ndim == 1:
                            weight_to_load = loaded_weight[chunk_slice]
                        elif split_dim == 0:
                            weight_to_load = loaded_weight[chunk_slice, :]
1423
                        else:
1424
                            weight_to_load = loaded_weight[:, chunk_slice]
1425
1426
1427
1428
1429
1430
1431
1432
1433
1434
1435
1436
1437
1438
1439
1440
1441
1442
1443
1444
1445
1446
1447
1448
1449
1450
1451
1452
1453
1454
1455
1456
1457
1458
1459
1460
1461
1462
1463
1464
1465
1466
                        # Synthesize an expert-style name so expert mapping
                        # can route it
                        chunk_name = name.replace(
                            "mlp.shared_experts",
                            f"mlp.experts.{self.config.n_routed_experts + j}",
                        )

                    # Use expert_params_mapping to locate the destination
                    # param and delegate to its expert-aware weight_loader
                    # with expert_id.
                    for mapping in expert_params_mapping:
                        param_name, weight_name, expert_id, shard_id = mapping
                        if weight_name not in chunk_name:
                            continue

                        # Anyway, this is an expert weight and should not be
                        # attempted to load as other weights later
                        is_expert_weight = True

                        # Do not modify `name` since the loop may continue here
                        # Instead, create a new variable
                        name_mapped = chunk_name.replace(weight_name, param_name)

                        if is_pp_missing_parameter(name_mapped, self):
                            continue

                        param = params_dict[name_mapped]
                        # We should ask the weight loader to return success or
                        # not here since otherwise we may skip experts with
                        # other available replicas.
                        weight_loader = typing.cast(
                            Callable[..., bool], param.weight_loader
                        )
                        success = weight_loader(
                            param,
                            weight_to_load,
                            name_mapped,
                            shard_id=shard_id,
                            expert_id=expert_id,
                            return_success=True,
                        )
                        if success:
1467
                            if not is_fusion_moe_shared_experts_layer:
1468
1469
1470
1471
1472
1473
1474
1475
1476
1477
1478
1479
1480
1481
1482
1483
1484
1485
1486
1487
1488
1489
1490
1491
1492
1493
1494
1495
                                name = name_mapped
                            else:
                                loaded_params.add(name_mapped)
                            break
                    else:
                        if is_expert_weight:
                            # We've checked that this is an expert weight
                            # However it's not mapped locally to this rank
                            # So we simply skip it
                            continue

                        # Skip loading extra bias for GPTQ models.
                        if name.endswith(".bias") and name not in params_dict:
                            continue

                        # Remapping the name of FP8 kv-scale.
                        name = maybe_remap_kv_scale_name(name, params_dict)
                        if name is None:
                            continue

                        if is_pp_missing_parameter(name, self):
                            continue

                        param = params_dict[name]
                        weight_loader = getattr(
                            param, "weight_loader", default_weight_loader
                        )
                        weight_loader(param, loaded_weight)
1496
            if not is_fusion_moe_shared_experts_layer:
1497
                loaded_params.add(name)
1498

1499
        return loaded_params
1500
1501


1502
1503
1504
1505
class DeepseekForCausalLM(DeepseekV2ForCausalLM):
    pass


1506
1507
class DeepseekV3ForCausalLM(DeepseekV2ForCausalLM):
    pass
1508
1509


1510
1511
# Compatibility with
# https://huggingface.co/deepseek-ai/DeepSeek-V3-Base/blob/main/configuration_deepseek.py
1512
def get_spec_layer_idx_from_weight_name(
1513
1514
    config: DeepseekV2Config | DeepseekV3Config, weight_name: str
) -> int | None:
1515
1516
1517
1518
    if (
        hasattr(config, "num_nextn_predict_layers")
        and config.num_nextn_predict_layers > 0
    ):
1519
1520
        layer_idx = config.num_hidden_layers
        for i in range(config.num_nextn_predict_layers):
1521
            if weight_name.startswith(f"model.layers.{layer_idx + i}."):
1522
1523
                return layer_idx + i
    return None