deepseek_v2.py 60.2 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."""
王敏's avatar
王敏 committed
26
27
import os
import re
28
from vllm import envs
zhuwenwen's avatar
zhuwenwen committed
29

30
31
import typing
from collections.abc import Callable, Iterable
32
from itertools import islice
wangding zeng's avatar
wangding zeng committed
33
34
35

import torch
from torch import nn
36
from transformers import DeepseekV2Config, DeepseekV3Config
wangding zeng's avatar
wangding zeng committed
37

38
from vllm._aiter_ops import rocm_aiter_ops
39
from vllm.attention.layer import Attention
40
from vllm.compilation.decorators import support_torch_compile
41
42
43
44
45
46
47
48
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,
)
49
from vllm.logger import init_logger
wangding zeng's avatar
wangding zeng committed
50
from vllm.model_executor.layers.activation import SiluAndMul
51
from vllm.model_executor.layers.attention_layer_base import AttentionLayerBase
52
from vllm.model_executor.layers.fused_moe import SharedFusedMoE
53
from vllm.model_executor.layers.layernorm import LayerNorm, RMSNorm
54
55
56
from vllm.model_executor.layers.linear import (
    ColumnParallelLinear,
    MergedColumnParallelLinear,
57
    QKVParallelLinear,
58
59
60
    ReplicatedLinear,
    RowParallelLinear,
)
wangding zeng's avatar
wangding zeng committed
61
from vllm.model_executor.layers.logits_processor import LogitsProcessor
62
from vllm.model_executor.layers.mla import MLAModules, MultiHeadLatentAttentionWrapper
63
from vllm.model_executor.layers.quantization import QuantizationConfig
64
from vllm.model_executor.layers.quantization.utils.fp8_utils import (
65
66
    per_token_group_quant_fp8,
)
wangding zeng's avatar
wangding zeng committed
67
from vllm.model_executor.layers.rotary_embedding import get_rope
68
from vllm.model_executor.layers.sparse_attn_indexer import SparseAttnIndexer
wangding zeng's avatar
wangding zeng committed
69
from vllm.model_executor.layers.vocab_parallel_embedding import (
70
71
72
    ParallelLMHead,
    VocabParallelEmbedding,
)
73
from vllm.model_executor.model_loader.weight_utils import (
74
75
76
    default_weight_loader,
    maybe_remap_kv_scale_name,
)
77
from vllm.model_executor.models.utils import sequence_parallel_chunk
78
from vllm.platforms import current_platform
79
from vllm.sequence import IntermediateTensors
80

81
from vllm.v1.attention.backend import AttentionBackend
82
83
84
from vllm.v1.attention.backends.mla.indexer import (
    DeepseekV32IndexerBackend,
)
85
from vllm.v1.kv_cache_interface import KVCacheSpec, MLAAttentionSpec
wangding zeng's avatar
wangding zeng committed
86

87
from .interfaces import MixtureOfExperts, SupportsEagle, SupportsLoRA, SupportsPP
88
89
90
91
92
93
94
from .utils import (
    PPMissingLayer,
    is_pp_missing_parameter,
    make_empty_intermediate_tensors_factory,
    make_layers,
    maybe_prefix,
)
95

96
from vllm import _custom_ops as ops
zhuwenwen's avatar
zhuwenwen committed
97
from vllm.utils import W8a8GetCacheJSON
98
99
100

logger = init_logger(__name__)

wangding zeng's avatar
wangding zeng committed
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
151
152
153
154
155
156
157
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,
158
            rope_parameters=config.rope_parameters,
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
        )
        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
182
183
184
185
186
187
188

class DeepseekV2MLP(nn.Module):
    def __init__(
        self,
        hidden_size: int,
        intermediate_size: int,
        hidden_act: str,
189
        quant_config: QuantizationConfig | None = None,
wangding zeng's avatar
wangding zeng committed
190
        reduce_results: bool = True,
191
        is_sequence_parallel=False,
192
        prefix: str = "",
wangding zeng's avatar
wangding zeng committed
193
194
    ) -> None:
        super().__init__()
195
196
197
198
199

        # 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
200
        self.gate_up_proj = MergedColumnParallelLinear(
201
202
            hidden_size,
            [intermediate_size] * 2,
wangding zeng's avatar
wangding zeng committed
203
            bias=False,
204
            quant_config=quant_config,
205
            disable_tp=is_sequence_parallel,
206
207
208
209
210
            prefix=f"{prefix}.gate_up_proj",
        )
        self.down_proj = RowParallelLinear(
            intermediate_size,
            hidden_size,
wangding zeng's avatar
wangding zeng committed
211
            bias=False,
212
            quant_config=quant_config,
213
            reduce_results=reduce_results,
214
            disable_tp=is_sequence_parallel,
215
216
            prefix=f"{prefix}.down_proj",
        )
wangding zeng's avatar
wangding zeng committed
217
        if hidden_act != "silu":
218
219
220
            raise ValueError(
                f"Unsupported activation: {hidden_act}. Only silu is supported for now."
            )
wangding zeng's avatar
wangding zeng committed
221
222
        self.act_fn = SiluAndMul()

223
224
225
226
227
    def forward(self, x):
        gate_up, _ = self.gate_up_proj(x)
        x = self.act_fn(gate_up)
        x, _ = self.down_proj(x)
        return x
wangding zeng's avatar
wangding zeng committed
228
229
230
231
232


class DeepseekV2MoE(nn.Module):
    def __init__(
        self,
233
        config: DeepseekV2Config | DeepseekV3Config,
234
        parallel_config: ParallelConfig,
235
        quant_config: QuantizationConfig | None = None,
236
        prefix: str = "",
wangding zeng's avatar
wangding zeng committed
237
238
239
    ):
        super().__init__()
        self.tp_size = get_tensor_model_parallel_world_size()
240
241
        self.tp_rank = get_tensor_model_parallel_rank()

242
        self.routed_scaling_factor = getattr(config, "routed_scaling_factor", 1.0)
243
244

        self.ep_group = get_ep_group().device_group
245
        self.ep_rank = get_ep_group().rank_in_group
246
247
248
        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
249

250
        self.is_sequence_parallel = parallel_config.use_sequence_parallel_moe
251

252
        if config.hidden_act != "silu":
253
254
255
256
257
258
259
260
261
262
263
264
            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",
        )
265
        if getattr(config, "topk_method", None) == "noaux_tc":
zhuwenwen's avatar
zhuwenwen committed
266
267
268
269
270
271
272
273
            if envs.VLLM_ENABLE_MOE_FUSED_GATE:
                # avoid moe_fused_gate precision error
                self.gate.e_score_correction_bias = nn.Parameter(
                torch.empty(config.n_routed_experts))
            else:
                self.gate.e_score_correction_bias = nn.Parameter(
                    torch.empty(config.n_routed_experts, dtype=torch.float32)
                )
274
275
276
        else:
            self.gate.e_score_correction_bias = None

277
        # Load balancing settings.
278
279
        eplb_config = parallel_config.eplb_config
        self.enable_eplb = parallel_config.enable_eplb
280

281
        self.n_redundant_experts = eplb_config.num_redundant_experts
282
        self.n_logical_experts = self.n_routed_experts
283
        self.n_physical_experts = self.n_logical_experts + self.n_redundant_experts
284
285
        self.n_local_physical_experts = self.n_physical_experts // self.ep_size

286
287
288
289
        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
        )
290

291
        self.is_rocm_aiter_moe_enabled = rocm_aiter_ops.is_fused_moe_enabled()
292
293
294
295
        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:
296
297
            self.shared_experts = None
        else:
298
            intermediate_size = config.moe_intermediate_size * config.n_shared_experts
299

wangding zeng's avatar
wangding zeng committed
300
301
302
303
304
            self.shared_experts = DeepseekV2MLP(
                hidden_size=config.hidden_size,
                intermediate_size=intermediate_size,
                hidden_act=config.hidden_act,
                quant_config=quant_config,
305
                is_sequence_parallel=self.is_sequence_parallel,
306
                reduce_results=False,
307
                prefix=f"{prefix}.shared_experts",
wangding zeng's avatar
wangding zeng committed
308
309
            )

310
311
        self.experts = SharedFusedMoE(
            shared_experts=self.shared_experts,
312
            gate=self.gate,
313
314
315
316
317
318
319
320
            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,
            use_grouped_topk=True,
321
322
            num_expert_group=getattr(config, "n_group", 1),
            topk_group=getattr(config, "topk_group", 1),
323
            prefix=f"{prefix}.experts",
324
            scoring_func=getattr(config, "scoring_func", "softmax"),
325
            # we do scaling outside, set factor to 1.0 to avoid double mul
326
327
            # aiter applies routed_scaling_factor internally
            routed_scaling_factor=1.0
328
            if not self.is_rocm_aiter_moe_enabled
329
            else self.routed_scaling_factor,
330
331
332
333
            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,
334
            n_shared_experts=config.n_shared_experts
335
            if self.is_fusion_moe_shared_experts_enabled
336
            else None,
337
        )
338

339
    def forward(self, hidden_states: torch.Tensor,
zhuwenwen's avatar
zhuwenwen committed
340
341
                rms_weight: torch.Tensor | None = None,
                residual: torch.Tensor | None = None
342
                ) -> torch.Tensor:
wangding zeng's avatar
wangding zeng committed
343
344
        num_tokens, hidden_dim = hidden_states.shape
        hidden_states = hidden_states.view(-1, hidden_dim)
345
346
347
348
349
350

        # 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:
351
            hidden_states = sequence_parallel_chunk(hidden_states)
352

353
354
355
356
357
        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
            )
358
        else:
359
360
361
362
363
            # 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
            )
364

365
366
367
        shared_output, final_hidden_states = fused_moe_out
        if self.shared_experts is None:
            assert shared_output is None
368
369
370

        # Fix FP16 overflow
        # See DeepseekV2DecoderLayer for more details.
zhuwenwen's avatar
zhuwenwen committed
371
        if hidden_states.dtype != torch.float16:
372
            if not self.is_rocm_aiter_moe_enabled:
373
                final_hidden_states *= self.routed_scaling_factor
374
375
        elif self.shared_experts is not None:
            assert shared_output is not None
376
            shared_output *= 1.0 / self.routed_scaling_factor
377
378
379
380

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

382
383
        if self.is_sequence_parallel:
            final_hidden_states = tensor_model_parallel_all_gather(
384
385
                final_hidden_states, 0
            )
386
387
            final_hidden_states = final_hidden_states[:num_tokens]
        elif self.tp_size > 1:
388
389
390
            final_hidden_states = self.experts.maybe_all_reduce_tensor_model_parallel(
                final_hidden_states
            )
wangding zeng's avatar
wangding zeng committed
391
392
393
394
395
396

        return final_hidden_states.view(num_tokens, hidden_dim)


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

wangding zeng's avatar
wangding zeng committed
398
399
400
401
402
    if scale <= 1:
        return 1.0
    return 0.1 * mscale * math.log(scale) + 1.0


403
404
405
406
407
408
409
410
411
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
412
413
414
415

class DeepseekV2Attention(nn.Module):
    def __init__(
        self,
416
        vllm_config: VllmConfig,
417
        config: DeepseekV2Config | DeepseekV3Config,
wangding zeng's avatar
wangding zeng committed
418
419
420
421
422
423
424
425
        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,
426
427
428
        cache_config: CacheConfig | None = None,
        quant_config: QuantizationConfig | None = None,
        topk_indices_buffer: torch.Tensor | None = None,
429
        prefix: str = "",
wangding zeng's avatar
wangding zeng committed
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
    ) -> 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
445
446
        assert topk_indices_buffer is None, (
            "topk_indices_buffer is not \
447
        supported for DeepseekV2Attention"
448
        )
wangding zeng's avatar
wangding zeng committed
449
450

        if self.q_lora_rank is not None:
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
            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
466
        else:
467
468
469
470
471
472
473
            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
474

475
476
477
478
479
        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,
480
481
482
            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
483
484
485
486
        self.kv_b_proj = ColumnParallelLinear(
            self.kv_lora_rank,
            self.num_heads * (self.qk_nope_head_dim + self.v_head_dim),
            bias=False,
487
            quant_config=quant_config,
488
489
            prefix=f"{prefix}.kv_b_proj",
        )
wangding zeng's avatar
wangding zeng committed
490
        # O projection.
491
492
493
494
495
496
497
        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",
        )
498
        if config.rope_parameters["rope_type"] != "default":
499
500
501
502
503
            config.rope_parameters["rope_type"] = (
                "deepseek_yarn"
                if config.rope_parameters.get("apply_yarn_scaling", True)
                else "deepseek_llama_scaling"
            )
504

505
506
507
        self.rotary_emb = get_rope(
            qk_rope_head_dim,
            max_position=max_position_embeddings,
508
            rope_parameters=config.rope_parameters,
509
510
            is_neox_style=False,
        )
wangding zeng's avatar
wangding zeng committed
511

512
513
514
515
        if (
            config.rope_parameters["rope_type"] != "default"
            and config.rope_parameters["rope_type"] == "deepseek_yarn"
        ):
516
517
            mscale_all_dim = config.rope_parameters.get("mscale_all_dim", False)
            scaling_factor = config.rope_parameters["factor"]
wangding zeng's avatar
wangding zeng committed
518
519
520
            mscale = yarn_get_mscale(scaling_factor, float(mscale_all_dim))
            self.scaling = self.scaling * mscale * mscale

521
522
523
524
525
526
527
528
529
        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
530
531
532
533
534

    def forward(
        self,
        positions: torch.Tensor,
        hidden_states: torch.Tensor,
535
        llama_4_scaling: torch.Tensor | None,
wangding zeng's avatar
wangding zeng committed
536
537
538
539
    ) -> torch.Tensor:
        if self.q_lora_rank is not None:
            q = self.q_a_proj(hidden_states)[0]
            q = self.q_a_layernorm(q)
540
            q = self.q_b_proj(q)[0].view(-1, self.num_local_heads, self.qk_head_dim)
wangding zeng's avatar
wangding zeng committed
541
        else:
542
543
544
545
            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
546
        latent_cache = self.kv_a_proj_with_mqa(hidden_states)[0]
547
        kv_a, _ = latent_cache.split([self.kv_lora_rank, self.qk_rope_head_dim], dim=-1)
wangding zeng's avatar
wangding zeng committed
548
        latent_cache = latent_cache.unsqueeze(1)
549
        kv_a = self.kv_a_layernorm(kv_a)
wangding zeng's avatar
wangding zeng committed
550
        kv = self.kv_b_proj(kv_a)[0]
551
        kv = kv.view(-1, self.num_local_heads, self.qk_nope_head_dim + self.v_head_dim)
wangding zeng's avatar
wangding zeng committed
552
        k_nope, v = kv.split([self.qk_nope_head_dim, self.v_head_dim], dim=-1)
553
        k_pe = latent_cache[:, :, self.kv_lora_rank :]
554

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

557
        q[..., self.qk_nope_head_dim :] = q_pe
wangding zeng's avatar
wangding zeng committed
558
        k = torch.empty_like(q)
559
560
        k[..., : self.qk_nope_head_dim] = k_nope
        k[..., self.qk_nope_head_dim :] = k_pe
561
562
563
564
565

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

566
567
        # padding value to qk_head_dim for alignment
        v = torch.nn.functional.pad(
568
569
            v, [0, self.qk_head_dim - self.v_head_dim], value=0
        ).view(-1, self.num_local_heads * self.qk_head_dim)
570
        attn_output = self.attn(q, k, v)
571
572
573
        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
574
575
576
577
        output, _ = self.o_proj(attn_output)
        return output


578
class DeepseekV32IndexerCache(torch.nn.Module, AttentionLayerBase):
579
580
581
    def __init__(
        self, head_dim: int, dtype: torch.dtype, prefix: str, cache_config: CacheConfig
    ):
582
583
584
585
586
587
588
589
590
591
592
        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

593
    def get_kv_cache_spec(self, vllm_config: VllmConfig) -> KVCacheSpec:
594
595
596
597
598
599
600
        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,
        )

601
    def forward(self): ...
602
603
604
605
606
607

    def get_attn_backend(self) -> AttentionBackend:
        return DeepseekV32IndexerBackend


class Indexer(nn.Module):
608
609
610
    def __init__(
        self,
        vllm_config: VllmConfig,
611
        config: DeepseekV2Config | DeepseekV3Config,
612
613
        hidden_size: int,
        q_lora_rank: int,
614
615
616
        quant_config: QuantizationConfig | None,
        cache_config: CacheConfig | None,
        topk_indices_buffer: torch.Tensor | None,
617
618
        prefix: str = "",
    ):
619
620
621
622
623
624
625
626
627
628
        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
629
630
631
632
633
634
635
636
637
638
639
640
641
642
        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",
        )
643
        self.k_norm = LayerNorm(self.head_dim, eps=1e-6)
644
        self.weights_proj = ReplicatedLinear(
645
646
647
648
649
            hidden_size,
            self.n_head,
            bias=False,
            quant_config=None,
            prefix=f"{prefix}.weights_proj",
650
        )
651
652
653
654
655
656
657
658
659
660
        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(
661
662
            head_dim=self.head_dim + self.head_dim // self.quant_block_size * 4 if torch.cuda.get_device_properties("cuda").gcnArchName.split(':')[0] == "gfx938" else self.head_dim,
            dtype=torch.uint8  if torch.cuda.get_device_properties("cuda").gcnArchName.split(':')[0] == "gfx938" else torch.bfloat16,
663
            prefix=f"{prefix}.k_cache",
664
665
            cache_config=cache_config,
        )
666
667
        self.max_model_len = vllm_config.model_config.max_model_len
        self.prefix = prefix
668
669
        from vllm.v1.attention.backends.mla.indexer import get_max_prefill_buffer_size

670
        self.max_total_seq_len = get_max_prefill_buffer_size(vllm_config)
671
672
673
674
675
676
677
678
679
680
        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,
        )
681

682
683
684
    def forward(
        self, hidden_states: torch.Tensor, qr: torch.Tensor, positions, rotary_emb
    ) -> torch.Tensor:
685
686
687
        q, _ = self.wq_b(qr)
        q = q.view(-1, self.n_head, self.head_dim)
        q_pe, q_nope = torch.split(
688
689
            q, [self.rope_dim, self.head_dim - self.rope_dim], dim=-1
        )
690
691
692
693

        k, _ = self.wk(hidden_states)
        k = self.k_norm(k)
        k_pe, k_nope = torch.split(
694
695
            k, [self.rope_dim, self.head_dim - self.rope_dim], dim=-1
        )
696
697

        q_pe, k_pe = rotary_emb(positions, q_pe, k_pe.unsqueeze(1))
698
699
700
701
702
        # 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)

703
704
        # `rotary_emb` is shape-preserving; `q_pe` is already
        # [num_tokens, n_head, rope_dim].
705
706
        q = torch.cat([q_pe, q_nope], dim=-1)
        # `k_pe` is [num_tokens, 1, rope_dim] (MQA).
707
        k = torch.cat([k_pe.squeeze(-2), k_nope], dim=-1)
708
709

        # we only quant q here since k quant is fused with cache insertion
710
711
712
713
714
715
716
717
718
719
        if not current_platform.is_rocm() or torch.cuda.get_device_properties("cuda").gcnArchName.split(':')[0] == "gfx938":
            q = q.view(-1, self.head_dim)
            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,
            )
            q_fp8 = q_fp8.view(-1, self.n_head, self.head_dim)
            q_scale = q_scale.view(-1, self.n_head, 1)
720
721

        weights, _ = self.weights_proj(hidden_states)
722
723
724
725
726
        if not current_platform.is_rocm() or torch.cuda.get_device_properties("cuda").gcnArchName.split(':')[0] == "gfx938":
            weights = (
                weights.unsqueeze(-1) * q_scale * self.softmax_scale * self.n_head**-0.5
            )
            weights = weights.squeeze(-1)
727

728
        return self.indexer_op(hidden_states, q_fp8, k, weights)
729
730


731
732
733
734
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).
735

736
737
        For more info see MLACommonImpl in:
        vllm/v1/attention/backends/mla/utils.py
738
739
740
741
    """

    def __init__(
        self,
742
        vllm_config: VllmConfig,
743
        config: DeepseekV2Config | DeepseekV3Config,
744
745
746
747
748
        hidden_size: int,
        num_heads: int,
        qk_nope_head_dim: int,
        qk_rope_head_dim: int,
        v_head_dim: int,
749
        q_lora_rank: int | None,
750
751
        kv_lora_rank: int,
        max_position_embeddings: int = 8192,
752
753
        cache_config: CacheConfig | None = None,
        quant_config: QuantizationConfig | None = None,
754
        prefix: str = "",
755
        topk_indices_buffer: torch.Tensor | None = None,
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
    ) -> 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:
776
            self.fused_qkv_a_proj = MergedColumnParallelLinear(
777
778
779
780
                self.hidden_size,
                [self.q_lora_rank, self.kv_lora_rank + self.qk_rope_head_dim],
                bias=False,
                quant_config=quant_config,
781
                prefix=f"{prefix}.fused_qkv_a_proj",
782
783
                disable_tp=True,
            )
784
785
786
787
788
789
        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,
790
791
                prefix=f"{prefix}.kv_a_proj_with_mqa",
            )
792
793

        if self.q_lora_rank is not None:
794
795
796
797
798
799
800
801
            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",
            )
802
        else:
803
804
805
806
807
808
809
810
            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)
811
812
813
814
815
        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,
816
817
818
819
820
821
822
823
824
            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",
        )
825

826
        if config.rope_parameters["rope_type"] != "default":
827
828
829
830
831
832
            config.rope_parameters["rope_type"] = (
                "deepseek_yarn"
                if config.rope_parameters.get("apply_yarn_scaling", True)
                else "deepseek_llama_scaling"
            )

833
834
835
        self.rotary_emb = get_rope(
            qk_rope_head_dim,
            max_position=max_position_embeddings,
836
            rope_parameters=config.rope_parameters,
837
838
            is_neox_style=False,
        )
839
840
841
842
843

        if (
            config.rope_parameters["rope_type"] != "default"
            and config.rope_parameters["rope_type"] == "deepseek_yarn"
        ):
844
845
            mscale_all_dim = config.rope_parameters.get("mscale_all_dim", False)
            scaling_factor = config.rope_parameters["factor"]
846
847
848
            mscale = yarn_get_mscale(scaling_factor, float(mscale_all_dim))
            self.scaling = self.scaling * mscale * mscale

zhuwenwen's avatar
zhuwenwen committed
849
        self.is_v32 = hasattr(config, "index_topk")
850
851

        if self.is_v32:
852
853
854
            self.indexer_rope_emb = get_rope(
                qk_rope_head_dim,
                max_position=max_position_embeddings,
855
                rope_parameters=config.rope_parameters,
856
857
                is_neox_style=True,
            )
858
859
860
861
862
863
864
865
866
867
            self.indexer = Indexer(
                vllm_config,
                config,
                hidden_size,
                q_lora_rank,
                quant_config,
                cache_config,
                topk_indices_buffer,
                f"{prefix}.indexer",
            )
868
        else:
869
            self.indexer_rope_emb = None
870
871
            self.indexer = None

872
873
        mla_modules = MLAModules(
            kv_a_layernorm=self.kv_a_layernorm,
874
            kv_b_proj=self.kv_b_proj,
875
876
877
            rotary_emb=self.rotary_emb,
            o_proj=self.o_proj,
            fused_qkv_a_proj=self.fused_qkv_a_proj
878
879
            if self.q_lora_rank is not None
            else None,
880
            kv_a_proj_with_mqa=self.kv_a_proj_with_mqa
881
882
883
            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,
884
885
            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,
886
            indexer=self.indexer,
887
            indexer_rotary_emb=self.indexer_rope_emb,
888
889
            is_sparse=self.is_v32,
            topk_indices_buffer=topk_indices_buffer,
890
        )
891

892
        self.mla_attn = MultiHeadLatentAttentionWrapper(
893
894
895
896
897
898
899
900
901
902
903
904
            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,
905
906
907
908
909
910
        )

    def forward(
        self,
        positions: torch.Tensor,
        hidden_states: torch.Tensor,
911
        llama_4_scaling: torch.Tensor | None,
912
    ) -> torch.Tensor:
913
        return self.mla_attn(positions, hidden_states, llama_4_scaling)
914
915


wangding zeng's avatar
wangding zeng committed
916
class DeepseekV2DecoderLayer(nn.Module):
917
918
919
920
    def __init__(
        self,
        vllm_config: VllmConfig,
        prefix: str,
921
922
        config: DeepseekV2Config | None = None,
        topk_indices_buffer: torch.Tensor | None = None,
923
    ) -> None:
wangding zeng's avatar
wangding zeng committed
924
        super().__init__()
925

926
927
        if config is None:
            config = vllm_config.model_config.hf_config
928
929
930
931
932
        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
933
        self.hidden_size = config.hidden_size
934
        max_position_embeddings = getattr(config, "max_position_embeddings", 8192)
935
        moe_layer_freq = getattr(config, "moe_layer_freq", 1)
936
937
        # DecoderLayers are created with `make_layers` which passes the prefix
        # with the layer's index.
938
        layer_idx = int(prefix.split(sep=".")[-1])
939
        self.layer_idx = layer_idx
940
941
942
943
944
945
946
947
948
949

        # 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)
        )

950
951
        self.use_mha = use_mha

952
953
954
        if use_mha:
            attn_cls = DeepseekAttention
        elif model_config.use_mla:
955
956
957
958
            attn_cls = DeepseekV2MLAAttention
        else:
            attn_cls = DeepseekV2Attention
        self.self_attn = attn_cls(
959
            vllm_config=vllm_config,
wangding zeng's avatar
wangding zeng committed
960
961
962
            config=config,
            hidden_size=self.hidden_size,
            num_heads=config.num_attention_heads,
963
964
965
            qk_nope_head_dim=qk_nope_head_dim,
            qk_rope_head_dim=qk_rope_head_dim,
            v_head_dim=v_head_dim,
966
            q_lora_rank=config.q_lora_rank if hasattr(config, "q_lora_rank") else None,
967
            kv_lora_rank=kv_lora_rank,
wangding zeng's avatar
wangding zeng committed
968
969
970
            max_position_embeddings=max_position_embeddings,
            cache_config=cache_config,
            quant_config=quant_config,
971
            prefix=f"{prefix}.self_attn",
972
            topk_indices_buffer=topk_indices_buffer,
wangding zeng's avatar
wangding zeng committed
973
        )
974

975
976
977
        if (
            config.n_routed_experts is not None
            and layer_idx >= config.first_k_dense_replace
978
            and layer_idx % moe_layer_freq == 0
979
        ):
980
981
            self.mlp = DeepseekV2MoE(
                config=config,
982
                parallel_config=parallel_config,
983
984
985
                quant_config=quant_config,
                prefix=f"{prefix}.mlp",
            )
wangding zeng's avatar
wangding zeng committed
986
987
988
989
990
991
        else:
            self.mlp = DeepseekV2MLP(
                hidden_size=config.hidden_size,
                intermediate_size=config.intermediate_size,
                hidden_act=config.hidden_act,
                quant_config=quant_config,
992
                prefix=f"{prefix}.mlp",
wangding zeng's avatar
wangding zeng committed
993
            )
994
995
996
997
        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
        )
998
        self.routed_scaling_factor = getattr(config, "routed_scaling_factor", 1.0)
wangding zeng's avatar
wangding zeng committed
999
1000
1001
1002
1003

    def forward(
        self,
        positions: torch.Tensor,
        hidden_states: torch.Tensor,
1004
        residual: torch.Tensor | None,
1005
        llama_4_scaling: torch.Tensor | None = None,
wangding zeng's avatar
wangding zeng committed
1006
    ) -> torch.Tensor:
zhuwenwen's avatar
zhuwenwen committed
1007
1008
1009
1010
1011
1012
1013
        # Self Attention
        # Fix residual FP16 overflow
        residual_fix_overflow = False
        if residual is None:
            residual = hidden_states.clone()
            hidden_states = self.input_layernorm(hidden_states)
            residual_fix_overflow = True
wangding zeng's avatar
wangding zeng committed
1014
        else:
zhuwenwen's avatar
zhuwenwen committed
1015
            hidden_states, residual = self.input_layernorm(hidden_states, residual)
1016
            
1017

1018
1019
1020
1021
1022
1023
1024
        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
1025

1026
1027
1028
1029
        if (
            not isinstance(self.self_attn, DeepseekAttention)
            and hidden_states.dtype == torch.float16
        ):
1030
1031
1032
            # Fix FP16 overflow
            # We scale both hidden_states and residual before
            # rmsnorm, and rmsnorm result would not affect by scale.
1033
            hidden_states *= 1.0 / self.routed_scaling_factor
1034
1035
1036
            if self.layer_idx == 0:
                # The residual is shared by all layers, we only scale it on
                # first layer.
1037
                residual *= 1.0 / self.routed_scaling_factor
1038
1039

        # Fully Connected
1040
        hidden_states, residual = self.post_attention_layernorm(hidden_states, residual)
wangding zeng's avatar
wangding zeng committed
1041
        hidden_states = self.mlp(hidden_states)
1042

1043
        if isinstance(self.mlp, DeepseekV2MLP) and hidden_states.dtype == torch.float16:
1044
1045
1046
1047
1048
            # 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
1049
            hidden_states *= 1.0 / self.routed_scaling_factor
1050

wangding zeng's avatar
wangding zeng committed
1051
1052
1053
        return hidden_states, residual


1054
@support_torch_compile
wangding zeng's avatar
wangding zeng committed
1055
1056
1057
class DeepseekV2Model(nn.Module):
    fall_back_to_pt_during_load = False

1058
    def __init__(self, *, vllm_config: VllmConfig, prefix: str = ""):
wangding zeng's avatar
wangding zeng committed
1059
        super().__init__()
1060
1061
1062

        config = vllm_config.model_config.hf_config
        quant_config = vllm_config.quant_config
1063
        self.config = config
1064
        self.device = current_platform.device_type
1065

wangding zeng's avatar
wangding zeng committed
1066
        self.vocab_size = config.vocab_size
zhuwenwen's avatar
zhuwenwen committed
1067
        self.is_v32 = hasattr(config, "index_topk")
1068
1069
1070
1071
1072
1073
        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,
1074
                device=self.device,
1075
            )
1076
1077
        else:
            topk_indices_buffer = None
wangding zeng's avatar
wangding zeng committed
1078

1079
1080
1081
1082
        if get_pp_group().is_first_rank:
            self.embed_tokens = VocabParallelEmbedding(
                config.vocab_size,
                config.hidden_size,
1083
                quant_config=quant_config,
1084
1085
                prefix=f"{prefix}.embed_tokens",
            )
1086
1087
1088
1089
        else:
            self.embed_tokens = PPMissingLayer()
        self.start_layer, self.end_layer, self.layers = make_layers(
            config.num_hidden_layers,
1090
            lambda prefix: DeepseekV2DecoderLayer(
1091
                vllm_config, prefix, topk_indices_buffer=topk_indices_buffer
1092
1093
1094
            ),
            prefix=f"{prefix}.layers",
        )
1095
1096
1097
1098
1099

        if get_pp_group().is_last_rank:
            self.norm = RMSNorm(config.hidden_size, eps=config.rms_norm_eps)
        else:
            self.norm = PPMissingLayer()
1100
1101
1102
        self.make_empty_intermediate_tensors = make_empty_intermediate_tensors_factory(
            ["hidden_states", "residual"], config.hidden_size
        )
wangding zeng's avatar
wangding zeng committed
1103

1104
    def embed_input_ids(self, input_ids: torch.Tensor) -> torch.Tensor:
1105
1106
        return self.embed_tokens(input_ids)

wangding zeng's avatar
wangding zeng committed
1107
1108
    def forward(
        self,
1109
        input_ids: torch.Tensor | None,
wangding zeng's avatar
wangding zeng committed
1110
        positions: torch.Tensor,
1111
1112
1113
        intermediate_tensors: IntermediateTensors | None,
        inputs_embeds: torch.Tensor | None = None,
    ) -> torch.Tensor | IntermediateTensors:
1114
        if get_pp_group().is_first_rank:
1115
1116
1117
            if inputs_embeds is not None:
                hidden_states = inputs_embeds
            else:
1118
                hidden_states = self.embed_input_ids(input_ids)
1119
1120
1121
1122
1123
1124
            residual = None
        else:
            assert intermediate_tensors is not None
            hidden_states = intermediate_tensors["hidden_states"]
            residual = intermediate_tensors["residual"]

1125
1126
1127
1128
1129
1130
1131
1132
1133
1134
1135
1136
1137
1138
        # 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

1139
        for layer in islice(self.layers, self.start_layer, self.end_layer):
1140
1141
1142
            hidden_states, residual = layer(
                positions, hidden_states, residual, llama_4_scaling
            )
1143
1144

        if not get_pp_group().is_last_rank:
1145
1146
1147
            return IntermediateTensors(
                {"hidden_states": hidden_states, "residual": residual}
            )
1148

wangding zeng's avatar
wangding zeng committed
1149
1150
1151
1152
        hidden_states, _ = self.norm(hidden_states, residual)
        return hidden_states


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
1182
1183
1184
1185
1186
1187
1188
1189
1190
1191
1192
1193
1194
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(
1195
    nn.Module, SupportsPP, DeepseekV2MixtureOfExperts, SupportsLoRA, SupportsEagle
1196
):
1197
1198
1199
    packed_modules_mapping = {
        "gate_up_proj": ["gate_proj", "up_proj"],
    }
1200
    model_cls = DeepseekV2Model
wangding zeng's avatar
wangding zeng committed
1201

1202
    def __init__(self, *, vllm_config: VllmConfig, prefix: str = ""):
wangding zeng's avatar
wangding zeng committed
1203
        super().__init__()
1204
1205
        config = vllm_config.model_config.hf_config
        quant_config = vllm_config.quant_config
1206
1207
1208
1209
1210
1211
1212

        self.quant_method = None
        if quant_config is not None:
            self.quant_method = quant_config.get_name()
            os.environ['LLAMA_NN'] = '0'
            os.environ['LM_NN'] = '0'

wangding zeng's avatar
wangding zeng committed
1213
1214
        self.config = config
        self.quant_config = quant_config
王敏's avatar
王敏 committed
1215

1216
1217
1218
1219
1220
1221
1222
1223
1224
        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"]

1225
1226
1227
1228
        # `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.
1229
1230
1231
        self.fuse_qkv_a_proj = (
            hasattr(config, "q_lora_rank") and config.q_lora_rank is not None
        )
1232
1233
1234
1235
1236
1237
        if self.fuse_qkv_a_proj:
            self.packed_modules_mapping["fused_qkv_a_proj"] = [
                "q_a_proj",
                "kv_a_proj_with_mqa",
            ]

1238
        self.model = self.model_cls(
1239
1240
            vllm_config=vllm_config, prefix=maybe_prefix(prefix, "model")
        )
1241
        if get_pp_group().is_last_rank:
1242
1243
1244
1245
1246
1247
            self.lm_head = ParallelLMHead(
                config.vocab_size,
                config.hidden_size,
                quant_config=quant_config,
                prefix=maybe_prefix(prefix, "lm_head"),
            )
1248
1249
        else:
            self.lm_head = PPMissingLayer()
wangding zeng's avatar
wangding zeng committed
1250
        self.logits_processor = LogitsProcessor(config.vocab_size)
1251
        self.make_empty_intermediate_tensors = (
1252
1253
            self.model.make_empty_intermediate_tensors
        )
1254
1255
1256
1257
1258
1259
1260
        # 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):
1261
1262
        self.expert_weights = []

1263
        self.num_expert_groups = getattr(self.config, "n_group", 1)
1264

1265
1266
        self.moe_layers = []
        self.moe_mlp_layers = []
1267
        example_moe = None
1268
        for layer in self.model.layers:
1269
1270
1271
            if isinstance(layer, PPMissingLayer):
                continue

1272
1273
            assert isinstance(layer, DeepseekV2DecoderLayer)
            if isinstance(layer.mlp, DeepseekV2MoE):
1274
1275
                # Pick last one layer since the first ones may be dense layers.
                example_moe = layer.mlp
1276
                self.moe_mlp_layers.append(layer.mlp)
1277
1278
                self.moe_layers.append(layer.mlp.experts)

1279
        self.extract_moe_parameters(example_moe)
王敏's avatar
王敏 committed
1280
        self.use_llama_nn = os.environ.get('LLAMA_NN') == '1'
1281
        self.use_awq_pad = os.environ.get('AWQ_PAD') == '1'
1282
        self.tritonsingleton= W8a8GetCacheJSON() 
zhuwenwen's avatar
zhuwenwen committed
1283
        self.tritonsingleton.topk = self.config.num_experts_per_tok
1284
        self.tritonsingleton.quant_method=self.quant_method 
1285

1286
1287
    def embed_input_ids(self, input_ids: torch.Tensor) -> torch.Tensor:
        return self.model.embed_input_ids(input_ids)
1288

wangding zeng's avatar
wangding zeng committed
1289
1290
    def forward(
        self,
1291
        input_ids: torch.Tensor | None,
wangding zeng's avatar
wangding zeng committed
1292
        positions: torch.Tensor,
1293
1294
1295
        intermediate_tensors: IntermediateTensors | None = None,
        inputs_embeds: torch.Tensor | None = None,
    ) -> torch.Tensor | IntermediateTensors:
1296
1297
1298
        hidden_states = self.model(
            input_ids, positions, intermediate_tensors, inputs_embeds
        )
wangding zeng's avatar
wangding zeng committed
1299
1300
        return hidden_states

1301
1302
1303
    def compute_logits(
        self,
        hidden_states: torch.Tensor,
1304
    ) -> torch.Tensor | None:
1305
        logits = self.logits_processor(self.lm_head, hidden_states)
wangding zeng's avatar
wangding zeng committed
1306
1307
        return logits

1308
1309
1310
1311
    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(
1312
            self,
1313
1314
1315
1316
1317
1318
            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,
        )
1319

1320
    def load_weights(self, weights: Iterable[tuple[str, torch.Tensor]]) -> set[str]:
1321
1322
1323
        rocm_aiter_moe_shared_expert_enabled = (
            rocm_aiter_ops.is_fusion_moe_shared_experts_enabled()
        )
wangding zeng's avatar
wangding zeng committed
1324
1325
1326
1327
        stacked_params_mapping = [
            # (param_name, shard_name, shard_id)
            ("gate_up_proj", "gate_proj", 0),
            ("gate_up_proj", "up_proj", 1),
1328
1329
        ]
        mla_params_mapping = [
1330
1331
            ("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
1332
        ]
1333
1334
1335
1336
1337
1338
1339
1340
1341
        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
1342

1343
1344
        # Params for weights, fp8 weight scales, fp8 activation scales
        # (param_name, weight_name, expert_id, shard_id)
1345
        expert_params_mapping = SharedFusedMoE.make_expert_params_mapping(
1346
            self,
王敏's avatar
王敏 committed
1347
1348
1349
            ckpt_gate_proj_name="gate_proj",
            ckpt_down_proj_name="down_proj",
            ckpt_up_proj_name="up_proj",
1350
1351
1352
            num_experts=self.config.n_routed_experts
            + (
                self.config.n_shared_experts
1353
                if rocm_aiter_moe_shared_expert_enabled
1354
1355
                else 0
            ),
1356
1357
            num_redundant_experts=self.num_redundant_experts,
        )
1358

wangding zeng's avatar
wangding zeng committed
1359
        params_dict = dict(self.named_parameters())
1360
        loaded_params: set[str] = set()
wangding zeng's avatar
wangding zeng committed
1361
1362
1363
        for name, loaded_weight in weights:
            if "rotary_emb.inv_freq" in name:
                continue
1364

1365
1366
1367
            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
1368

1369
1370
            is_fusion_moe_shared_experts_layer = (
                rocm_aiter_moe_shared_expert_enabled and ("mlp.shared_experts" in name)
1371
1372
            )

1373
            for param_name, weight_name, shard_id in stacked_params_mapping:
1374
                # Skip non-stacked layers and experts (experts handled below).
wangding zeng's avatar
wangding zeng committed
1375
1376
                if weight_name not in name:
                    continue
1377
1378
1379
1380
1381
1382
                # 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.
1383
                if ("mlp.experts." in name) and name not in params_dict:
1384
                    continue
1385
                if is_fusion_moe_shared_experts_layer:
1386
                    continue
1387
                name_mapped = name.replace(weight_name, param_name)
1388
1389
1390

                # QKV fusion is optional, fall back to normal
                # weight loading if it's not enabled
1391
                # if go with fusion option, then update name
1392
1393
1394
                if (
                    param_name == "fused_qkv_a_proj"
                ) and name_mapped not in params_dict:
1395
                    continue
1396
1397
                else:
                    name = name_mapped
wangding zeng's avatar
wangding zeng committed
1398
1399
1400
                # Skip loading extra bias for GPTQ models.
                if name.endswith(".bias") and name not in params_dict:
                    continue
1401
1402
1403
1404

                if is_pp_missing_parameter(name, self):
                    continue

wangding zeng's avatar
wangding zeng committed
1405
1406
1407
1408
1409
                param = params_dict[name]
                weight_loader = param.weight_loader
                weight_loader(param, loaded_weight, shard_id)
                break
            else:
1410
1411
                is_expert_weight = False

1412
1413
1414
1415
1416
1417
1418
1419
1420
                # 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
1421
                if is_fusion_moe_shared_experts_layer:
1422
1423
1424
1425
                    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
1426
1427
1428
1429
1430
                    split_dim = (
                        1
                        if ("down_proj.weight" in name and loaded_weight.ndim > 1)
                        else 0
                    )
1431
1432
1433
1434
                    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}"
1435
                    )
1436
1437
1438
1439
1440
1441
                    chunk_size = total // num_chunks

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

1442
                    if is_fusion_moe_shared_experts_layer:
1443
1444
1445
1446
1447
                        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, :]
1448
                        else:
1449
                            weight_to_load = loaded_weight[:, chunk_slice]
1450
1451
1452
1453
1454
1455
1456
1457
1458
1459
1460
1461
1462
1463
1464
1465
1466
1467
1468
1469
1470
1471
1472
1473
1474
1475
1476
1477
1478
1479
1480
1481
1482
1483
1484
1485
1486
1487
1488
1489
1490
1491
                        # 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:
1492
                            if not is_fusion_moe_shared_experts_layer:
1493
1494
1495
1496
1497
1498
1499
1500
1501
1502
1503
1504
1505
1506
1507
1508
1509
1510
1511
1512
1513
1514
                                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
1515

1516
                        param = params_dict[name]
1517
1518
1519
1520
                        weight_loader = getattr(
                            param, "weight_loader", default_weight_loader
                        )
                        weight_loader(param, loaded_weight)
1521
            if not is_fusion_moe_shared_experts_layer:
1522
                loaded_params.add(name)
1523
                
1524
1525
1526
1527
1528
1529
1530
1531
1532
1533
1534
1535
1536
1537
1538
1539
1540
1541
1542
1543
1544
1545
1546
1547
1548
1549
1550
1551
1552
        if self.use_llama_nn and self.quant_method is None:
            lay_key_words = [
                "self_attn.q_proj.weight",
                "self_attn.q_a_proj.weight",
                "self_attn.q_b_proj.weight",
                "self_attn.kv_a_proj_with_mqa.weight",
                "self_attn.kv_b_proj.weight",
                "self_attn.o_proj.weight",
                "mlp.gate_up_proj.weight",
                "mlp.down_proj.weight",
                "mlp.gate.weight",
                "shared_experts.gate_up_proj.weight",
                "shared_experts.down_proj.weight",
                "lm_head.weight",
            ]

            combined_words = "|".join(lay_key_words)
            
            for layername in loaded_params:
                weight = params_dict[layername]
                matches = re.findall(combined_words, layername)
                if matches:
                    _weight = torch.zeros_like(weight.data)
                    ori_shape =_weight.shape
                    
                    ops.trans_w16_gemm(_weight, weight.data, _weight.shape[0], _weight.shape[1])
                    weight.data.copy_(_weight)
                    
                    weight.data=weight.data.reshape(ori_shape[1],-1)
zhuwenwen's avatar
zhuwenwen committed
1553
            
1554
        return loaded_params
1555
1556


1557
1558
1559
1560
class DeepseekForCausalLM(DeepseekV2ForCausalLM):
    pass


1561
1562
class DeepseekV3ForCausalLM(DeepseekV2ForCausalLM):
    pass
1563
1564


1565
1566
# Compatibility with
# https://huggingface.co/deepseek-ai/DeepSeek-V3-Base/blob/main/configuration_deepseek.py
1567
def get_spec_layer_idx_from_weight_name(
1568
1569
    config: DeepseekV2Config | DeepseekV3Config, weight_name: str
) -> int | None:
1570
1571
1572
1573
    if (
        hasattr(config, "num_nextn_predict_layers")
        and config.num_nextn_predict_layers > 0
    ):
1574
1575
        layer_idx = config.num_hidden_layers
        for i in range(config.num_nextn_predict_layers):
1576
            if weight_name.startswith(f"model.layers.{layer_idx + i}."):
1577
1578
                return layer_idx + i
    return None