base.py 43.8 KB
Newer Older
1
# Copyright (c) 2022-2024, NVIDIA CORPORATION & AFFILIATES. All rights reserved.
2
3
4
5
#
# See LICENSE for license information.

"""Base modules and utilities for TransformerEngine PyTorch API"""
6
import io
7
8
9
import os
import pickle
import warnings
10
11
12
import socket
import fcntl
import struct
13
from abc import ABC, abstractmethod
14
from typing import Dict, Generator, List, Optional, Tuple, Union
15
16
17
18
19
from contextlib import contextmanager

import torch
import torch.nn.functional as F

20
import transformer_engine_torch as tex
21
from ._common import _ParameterInitMeta
22
from ..export import is_in_onnx_export_mode
23
24
25
from ..fp8 import (
    get_default_fp8_recipe,
    get_fp8_te_dtype,
26
    FP8GlobalStateManager,
27
28
29
30
31
)
from ..distributed import (
    gather_along_first_dim,
    is_fp8_activation_recompute_enabled,
    in_fp8_activation_recompute_phase,
32
    _fsdp_gather_tensors,
33
34
35
36
37
38
39
)
from ..cpp_extensions import (
    fp8_cast_transpose_fused,
    fp8_cast_transpose_bgrad_fused,
    cast_to_fp8,
)
from ..constants import dist_group_type
40
from ..float8_tensor import Float8Tensor
41

42
43
__all__ = ["initialize_ub", "destroy_ub"]

44
45
46
_2X_ACC_FPROP = False
_2X_ACC_DGRAD = True
_2X_ACC_WGRAD = True
47
_multi_stream_cublas_workspace = []
48
49
50
_cublas_workspace = None
_ub_communicators = None
_NUM_MAX_UB_STREAMS = 3
51
layers_atomic_ring_exchange = []
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70


def get_cublas_workspace_size_bytes() -> None:
    """Return 32 MiB if using hopper, 4 MiB for all other architectures."""
    if torch.cuda.get_device_properties(torch.cuda.current_device()).major >= 9:
        return 33_554_432
    return 4_194_304


def get_workspace() -> torch.Tensor:
    """Returns workspace for cublas."""
    global _cublas_workspace
    if _cublas_workspace is None:
        _cublas_workspace = torch.empty(
            get_cublas_workspace_size_bytes(), dtype=torch.uint8, device="cuda"
        )
    return _cublas_workspace


71
72
73
74
def get_multi_stream_cublas_workspace() -> List[torch.Tensor]:
    """Returns workspace for multi-stream cublas."""
    global _multi_stream_cublas_workspace
    if not _multi_stream_cublas_workspace:
75
        for _ in range(tex._num_cublas_streams):
76
77
78
79
80
81
            _multi_stream_cublas_workspace.append(
                torch.empty(get_cublas_workspace_size_bytes(), dtype=torch.uint8, device="cuda")
            )
    return _multi_stream_cublas_workspace


82
83
def initialize_ub(
    shape: list,
84
    tp_size: int,
85
    use_fp8: bool = False,
86
    dtype: torch.dtype = torch.bfloat16,
87
    ub_cfgs: Optional[dict] = None,
88
    bootstrap_backend: Union[str, torch.distributed.Backend] = None,
89
90
) -> None:
    """Initialize communicators for TP comm overlap using userbuffers."""
91
92
93
94
95
96
    if not tex.device_supports_multicast():
        assert bool(os.getenv("UB_SKIPMC", "0")), (
            "CUDA device, driver and/or toolkit version does not support comm+GEMM overlap with "
            + "CUDA Multicast. Launch app with UB_SKIPMC=1 to try CUDA IPC instead."
        )

97
98
99
    global _ub_communicators
    assert _ub_communicators is None, "UB communicators are already initialized."
    _ub_communicators = {}
100
101
102
103
104
105
106
107
108
109

    if tex.ubuf_built_with_mpi():
        # Userbuffers will ignore all these values when it is built with MPI, so these are just
        # placeholders based on an assumption that tp_size covers all devices in a physical node.
        assert torch.distributed.is_mpi_available()
        mpi_group = torch.distributed.new_group(backend="mpi")
        world_rank = torch.distributed.get_rank(mpi_group)
        world_size = torch.distributed.get_world_size(mpi_group)
        local_rank = world_rank % tp_size
        local_size = tp_size
110
        self_node_idx = world_rank // tp_size
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
        num_nodes = world_size // tp_size
        ub_callbacks = tex.UbufBootstrapCallbacks()
    else:
        assert (
            torch.distributed.is_initialized()
        ), "torch.distributed must be initialized before Userbuffers"
        if bootstrap_backend is None:
            bootstrap_backend = "nccl"
            if torch.distributed.is_gloo_available():
                bootstrap_backend = "gloo"
            elif torch.distributed.is_mpi_available():
                bootstrap_backend = "mpi"
        else:
            assert bootstrap_backend in ["gloo", "mpi", "nccl"]

        world_group = torch.distributed.new_group(backend=bootstrap_backend)
        world_rank = torch.distributed.get_rank(world_group)
        world_size = torch.distributed.get_world_size(world_group)

        # Construct an intra-node communicator based on global ranks that share the same hostname
        # NOTE: If the user specified a valid network interface for NCCL or GLOO, use the host
        #       address on that interface instead of the hostname. This can help avoid issues when
        #       different hosts have the same hostname on Kubernetes clusters.
        hostname = socket.gethostname()
        ifname = os.getenv(
            "NVTE_UB_SOCKET_IFNAME",
            os.getenv("NCCL_SOCKET_IFNAME", os.getenv("GLOO_SOCKET_IFNAME")),
        )

        if ifname is not None:
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
            # Make sure the ifname found in the environment is a valid network interface
            if ifname in [name for _, name in socket.if_nameindex()]:
                s = socket.socket(socket.AF_INET, socket.SOCK_DGRAM)
                try:
                    hostname = socket.inet_ntoa(
                        fcntl.ioctl(
                            s.fileno(), 0x8915, struct.pack("256s", ifname[:15].encode("UTF-8"))
                        )[20:24]
                    )
                except OSError as err:
                    raise OSError(f"Invalid network interface: {ifname}") from err
                finally:
                    s.close()
            else:
                ifname_warning = (
                    f"'{ifname}' is not a valid network interface! `te.initialize_ub()` will"
                    " attempt to "
                    + "detect ranks on the same node by matching 'socket.gethostname()', which is "
                    + "known to fail on virtual clusters like Kubernetes. If Userbuffers "
                    + "initialization fails, please set the 'NVTE_UB_SOCKET_IFNAME' variable in "
                    + "your environment to the correct network interface."
162
                )
163
                warnings.warn(ifname_warning, UserWarning)
164
165
166

        hostnames = [None for _ in range(world_size)]
        torch.distributed.all_gather_object(hostnames, hostname, world_group)
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
        unique_hosts = []
        for host in hostnames:
            if host not in unique_hosts:
                unique_hosts.append(host)
        num_nodes = len(unique_hosts)

        if num_nodes > 1:
            ranks_per_node_list = [[] for _ in range(num_nodes)]
            self_node_idx = -1
            for i, host in enumerate(hostnames):
                node_idx = unique_hosts.index(host)
                ranks_per_node_list[node_idx].append(i)
                if host == hostname:
                    self_node_idx = node_idx
            assert self_node_idx >= 0, "Internal TE error!"

            intra_node_group, _ = torch.distributed.new_subgroups_by_enumeration(
                ranks_per_node_list, backend=bootstrap_backend
            )
            local_rank = torch.distributed.get_rank(intra_node_group)
            local_size = torch.distributed.get_world_size(intra_node_group)
            intra_node_ranks = torch.distributed.get_process_group_ranks(intra_node_group)

        else:
            self_node_idx = 0
192
193
194
195
196
            intra_node_group = world_group
            local_rank = world_rank
            local_size = world_size
            intra_node_ranks = list(range(world_size))

197
198
        if world_rank == 0:
            print(f"!!! [UB] Number of physical nodes: {num_nodes}\n", end="", flush=True)
199
200
        if local_rank == 0:
            print(
201
                f"!!! [UB] Global ranks on node {self_node_idx}: {intra_node_ranks}\n",
202
203
204
205
206
                end="",
                flush=True,
            )

        ub_callbacks = tex.UbufBootstrapCallbacks(world_group, intra_node_group)
207
208
209
210
211
212

    # Increase the workspace by the number of maximum concurrent streams
    global _cublas_workspace
    _cublas_workspace = get_workspace().repeat(_NUM_MAX_UB_STREAMS)

    # Default buffer precision: AllGather buffers use fp8 when using fp8 recipe
213
    layers_all_gather_overlap = [
214
215
216
217
218
219
        "qkv_fprop",
        "qkv_dgrad",
        "proj_dgrad",
        "fc1_fprop",
        "fc1_dgrad",
        "fc2_dgrad",
220
    ]
221
    layers_reduce_scatter_overlap = ["proj_fprop", "fc2_fprop", "qkv_wgrad", "fc1_wgrad"]
Jaemin Choi's avatar
Jaemin Choi committed
222
    dgrad_reduce_scatter_overlap = ["qkv_dgrad", "fc1_dgrad"]
223
224
    # Default overlap methods for layers
    methods = {
225
226
227
        "ring_exchange": ["qkv_fprop", "fc1_fprop", "proj_dgrad", "fc2_dgrad"],
        "pipeline": ["proj_fprop", "fc2_fprop"],
        "bulk": ["qkv_dgrad", "qkv_wgrad", "fc1_dgrad", "fc1_wgrad"],
228
229
    }

230
    # AG-RS overlap pairs of layers forming a tensor-parallel block
231
232
    ag_rs_pairs = {"qkv_fprop": "proj_fprop", "fc1_fprop": "fc2_fprop"}
    rs_ag_pairs = {v: k for k, v in ag_rs_pairs.items()}
233
234
235
    global layers_atomic_ring_exchange
    layers_atomic_ring_exchange = []

236
237
238
239
240
241
    def get_method(name):
        for method, names in methods.items():
            if name in names:
                return method
        raise KeyError(f"Given layer name {name} does not exist.")

242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
    def get_default_config(name):
        method = get_method(name)
        is_reduce_scatter = name in layers_reduce_scatter_overlap
        default_cfg = {
            "method": method,
            "is_reduce_scatter": is_reduce_scatter,
            "num_sm": 1 if method == "ring_exchange" else 16,
            "cga_size": 1 if method == "ring_exchange" else 2,
            "set_sm_margin": False,
            "num_splits": 4 if method == "pipeline" else tp_size,
            "aggregate": False,
            "atomic_gemm": False,
            "use_ce": True,
            "fp8_buf": name in layers_all_gather_overlap,
        }
        return default_cfg

259
260
261
    def add_ub(
        name: str,
        method: str,
262
        is_reduce_scatter: int,
263
264
265
        num_sm: int = 16,
        cga_size: int = 2,
        set_sm_margin: int = 0,
266
        num_splits: int = 0,
267
        aggregate: int = 0,
268
        atomic_gemm: int = 0,
269
        use_ce: bool = True,
270
        fp8_buf: bool = False,
271
    ) -> None:
272
273
274
275
276
        if atomic_gemm:
            warnings.warn(
                "Atomic GEMM uses a beta API from cublas and is not tested for all use cases."
            )
            assert use_fp8, "Atomic GEMM overlap supported only for FP8 GEMM."
277
            if method == "bulk":
278
                warnings.warn(
279
                    f"At {name}, atoimic GEMM not is supported for a bulk overlap."
280
281
282
                    "Defaulting to `atomic_gemm=False`."
                )
                atomic_gemm = 0
283
        if not is_reduce_scatter and method == "pipeline":
284
            raise ValueError(
285
                f"At {name}, `pipeline` overlap method is not supported for AllGather."
286
            )
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
        # Check if both AG and RS overlaps use `atomic GEMM`` + `p2p ring-exchange`.
        # Using atomic GEMM + p2p ring-exchange in only one of the pair breaks functionality.
        global layers_atomic_ring_exchange
        if atomic_gemm and method == "ring_exchange" and name in ag_rs_pairs:
            layers_atomic_ring_exchange += [name, ag_rs_pairs[name]]
        if name in rs_ag_pairs:
            assert_message = (
                f"At {name}, atomic AG-GEMM overlap with `ring_exchange` shuffles GEMM chunk "
                "outputs, and  RS-GEMM overlap un-suffle them. When one of the GEMM-AG and "
                "GEMM-RS overlaps forming a TP block (e.g., qkv_fprop and proj_fprop) uses "
                "`atomic gemm` and `ring_exhcnage`, its pair must use the same overlap config "
                "for functionality."
            )
            if name in layers_atomic_ring_exchange:
                assert atomic_gemm and method == "ring_exchange", assert_message
            else:
                if atomic_gemm and method == "ring_exchange":
                    assert rs_ag_pairs[name] in layers_atomic_ring_exchange, assert_message

306
        sample_buffer = torch.empty(
307
308
309
            shape, dtype=torch.uint8 if (use_fp8 and fp8_buf) else dtype, device="cuda"
        )
        if method == "ring_exchange":
310
            ub_obj = tex.UbufP2PCommOverlap(
311
                sample_buffer,  # Sample userbuffer
312
                world_rank,  # World rank
313
                world_size,  # World size
314
315
                local_rank,  # Rank within the node
                local_size,  # Number of ranks/GPUs per node
316
                self_node_idx,  # Node ID
317
318
                num_nodes,  # Number of nodes
                tp_size,  # Tensor-parallel group size (may be different than local_size)
319
320
321
322
323
                num_sm,  # Number of communication SMs
                cga_size,  # CGA cluster size
                set_sm_margin,  # Set SM margin
                aggregate,  # Aggregate 2X GEMM chunks
                _NUM_MAX_UB_STREAMS,  # Max concurrent GEMM streams
324
325
326
327
                is_reduce_scatter,  # Overlap with reduce scatter
                atomic_gemm,  # Use a single GEMM with atomic-counters
                use_ce,  # Use copy engine for P2P communications
                ub_callbacks,
328
            )
329
330
        else:
            ub_obj = tex.UbufCommOverlap(
331
                sample_buffer,  # Sample userbuffer
332
                world_rank,  # World rank
333
                world_size,  # World size
334
335
                local_rank,  # Rank within the node
                local_size,  # Number of ranks/GPUs per node
336
                self_node_idx,  # Node ID
337
338
                num_nodes,  # Number of nodes
                tp_size,  # Tensor-parallel group size (may be different than local_size)
339
340
341
342
343
                num_sm,  # Number of communication SMs
                cga_size,  # CGA cluster size
                num_splits,  # Number of communication splits
                set_sm_margin,  # Set SM margin
                _NUM_MAX_UB_STREAMS,  # Max concurrent GEMM streams
344
345
                atomic_gemm,  # Use a single GEMM with atomic-counters
                ub_callbacks,
346
            )
347
348
        _ub_communicators[name] = ub_obj

Jaemin Choi's avatar
Jaemin Choi committed
349
350
    if ub_cfgs is not None:
        for name in dgrad_reduce_scatter_overlap:
351
352
            if name in ub_cfgs and "method" in ub_cfgs[name] and ub_cfgs[name]["method"] != "bulk":
                wgrad_name = name.replace("dgrad", "wgrad")
Jaemin Choi's avatar
Jaemin Choi committed
353
354
                assert wgrad_name not in ub_cfgs
                layers_reduce_scatter_overlap.remove(wgrad_name)
355
                layers_all_gather_overlap.remove(name)
Jaemin Choi's avatar
Jaemin Choi committed
356
                layers_reduce_scatter_overlap.append(name)
357
358
359
                methods["bulk"].remove(name)
                new_method = ub_cfgs[name]["method"]
                methods[new_method].append(name)
Jaemin Choi's avatar
Jaemin Choi committed
360

361
    for name in methods["ring_exchange"] + methods["pipeline"] + methods["bulk"]:
362
        ub_cfg = get_default_config(name)
363
        if ub_cfgs is not None and name in ub_cfgs:
364
            fp8_buf = (name in layers_all_gather_overlap) or (
365
                ub_cfgs[name].get("fp8_buf", False) and name in methods["pipeline"]
366
            )
367
368
369
            ub_cfg.update(ub_cfgs[name])
            ub_cfg["fp8_buf"] = fp8_buf
        add_ub(name, **ub_cfg)
370
371
372
373
374
375
376
377


def get_ub(name: str):
    """Get userbuffer communicator corresponding to give key."""
    assert _ub_communicators is not None, "UB manager is not initialized."
    assert name in _ub_communicators, f"UB for {name} is not registered."
    return _ub_communicators[name]

378

379
380
381
382
383
384
385
def destroy_ub():
    """Destroy all allocated userbuffer communicators."""
    global _ub_communicators
    _ub_communicators = None
    global layers_atomic_ring_exchange
    layers_atomic_ring_exchange = []

386
387
388
389
390
391
392
393
394
395
396

class TransformerEngineBaseModule(torch.nn.Module, ABC):
    """Base TE module."""

    def __init__(self) -> None:
        super().__init__()
        assert torch.cuda.is_available(), "TransformerEngine needs CUDA."
        self.fp8_initialized = False
        self.fp8 = False
        self.fp8_calibration = False
        self.fp8_meta = {}
397
        self.fp8_meta["fp8_checkpoint"] = False
398
399
400
401
402
403
        self.fp8_meta["fp8_group"] = None
        self.fp8_meta["recipe"] = get_default_fp8_recipe()
        self.fp8_meta_tensors_initialized = False
        self.tp_group = None
        self.tp_size = 1
        self.sequence_parallel = False
404
405
        self.param_init_meta = {}
        self.primary_weights_in_fp8 = FP8GlobalStateManager.with_fp8_parameters()
406
407
        self.fsdp_wrapped = False
        self.fsdp_group = None
408
        self._fp8_workspaces: Dict[str, Float8Tensor] = {}
409

410
411
412
413
414
415
416
417
418
419
420
421
    def adjust_amax_history_length(self, length: int, fwd: Optional[bool] = None) -> None:
        """Increase or decrease size of amax history based on given `length`.

        .. warning::
            This changes the underlying amax memory location.
        """
        if fwd is None:
            fp8_meta_tensor_keys = ("scaling_fwd", "scaling_bwd")
        else:
            fp8_meta_tensor_keys = ("scaling_fwd" if fwd else "scaling_bwd",)

        for meta_key in fp8_meta_tensor_keys:
422
423
424
            if meta_key not in self.fp8_meta:
                # Handles non-parameter FP8 modules, e.g. DPA.
                continue
425
426
427
428
429
            curr_len = self.fp8_meta[meta_key].amax_history.shape[0]
            if length == curr_len:
                continue
            if length < curr_len:
                self.fp8_meta[meta_key].amax_history = (
430
431
                    self.fp8_meta[meta_key].amax_history[:length].clone()
                )
432
433
434
435
436
437
438
439
            elif length > curr_len:
                extra_rows = length - curr_len
                self.fp8_meta[meta_key].amax_history = F.pad(
                    self.fp8_meta[meta_key].amax_history, pad=(0, 0, 0, extra_rows)
                )

            # Update the global buffers with new amax and history pointers.
            if FP8GlobalStateManager.get_buffer_info() in self.fp8_meta:
440
441
442
                fwd_pos, fwd_key, bwd_pos, bwd_key = self.fp8_meta[
                    FP8GlobalStateManager.get_buffer_info()
                ]
443
444
445
446
447
                for pos, buffer_key in zip((fwd_pos, bwd_pos), (fwd_key, bwd_key)):
                    if buffer_key in FP8GlobalStateManager.global_amax_buffer:
                        assert (
                            buffer_key in FP8GlobalStateManager.global_amax_history_buffer
                        ), "TE internal error during amax history change."
448
449
450
                        FP8GlobalStateManager.global_amax_buffer[buffer_key][pos] = self.fp8_meta[
                            meta_key
                        ].amax_history[0]
451
                        FP8GlobalStateManager.global_amax_history_buffer[buffer_key][pos] = (
452
453
                            self.fp8_meta[meta_key].amax_history
                        )
454

455
456
457
458
459
460
    def set_meta_tensor(self, fwd: bool) -> None:
        """Init scales and amaxes for fwd | bwd."""
        fp8_meta_tensor_key = "scaling_fwd" if fwd else "scaling_bwd"

        if self.fp8_meta_tensors_initialized:
            # Handle changed amax history size.
461
            self.adjust_amax_history_length(self.fp8_meta["recipe"].amax_history_len, fwd=fwd)
462
463
464
465
            return

        # Max. number of fp8 tensors per GEMM = 3 (input, weight, output) for fwd and
        # 2 (grad_output and grad_input) for bwd
466
        num_fp8_tensors = self.fp8_meta["num_gemms"] * 3 if fwd else self.fp8_meta["num_gemms"] * 2
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487

        self.fp8_meta[fp8_meta_tensor_key] = tex.FP8TensorMeta()
        self.fp8_meta[fp8_meta_tensor_key].scale = torch.ones(
            num_fp8_tensors, dtype=torch.float32, device="cuda"
        )
        self.fp8_meta[fp8_meta_tensor_key].scale_inv = torch.ones(
            num_fp8_tensors, dtype=torch.float32, device="cuda"
        )
        self.fp8_meta[fp8_meta_tensor_key].amax_history = torch.zeros(
            self.fp8_meta["recipe"].amax_history_len,
            num_fp8_tensors,
            dtype=torch.float32,
            device="cuda",
        )

    def init_fp8_meta_tensors(self) -> None:
        """Init scales and amaxes."""
        self.set_meta_tensor(True)
        self.set_meta_tensor(False)
        self.fp8_meta_tensors_initialized = True

488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
    def get_fp8_meta_tensors(self) -> None:
        """Get scales and amaxes."""
        fwd_key, bwd_key = "scaling_fwd", "scaling_bwd"
        if fwd_key not in self.fp8_meta or bwd_key not in self.fp8_meta:
            return None

        fp8_meta_tensors = {fwd_key: [], bwd_key: []}
        with torch.no_grad():
            for key in (fwd_key, bwd_key):
                fp8_meta_tensors[key].append(self.fp8_meta[key].scale.clone())
                fp8_meta_tensors[key].append(self.fp8_meta[key].scale_inv.clone())
                fp8_meta_tensors[key].append(self.fp8_meta[key].amax_history.clone())
        return fp8_meta_tensors

    def reset_fp8_meta_tensors(self, fp8_meta_tensors=None) -> None:
        """Reset scales and amaxes."""
504

505
506
507
508
509
        def reset(key):
            if key in self.fp8_meta:
                if fp8_meta_tensors is None:
                    self.fp8_meta[key].scale.copy_(torch.ones_like(self.fp8_meta[key].scale))
                    self.fp8_meta[key].scale_inv.copy_(
510
511
                        torch.ones_like(self.fp8_meta[key].scale_inv)
                    )
512
                    self.fp8_meta[key].amax_history.copy_(
513
514
                        torch.zeros_like(self.fp8_meta[key].amax_history)
                    )
515
516
517
518
519
                else:
                    assert key in fp8_meta_tensors, "Cannot reset fp8 tensors."
                    self.fp8_meta[key].scale.copy_(fp8_meta_tensors[key][0])
                    self.fp8_meta[key].scale_inv.copy_(fp8_meta_tensors[key][1])
                    self.fp8_meta[key].amax_history.copy_(fp8_meta_tensors[key][2])
520

521
522
523
524
        with torch.no_grad():
            reset("scaling_fwd")
            reset("scaling_bwd")

525
526
527
    def get_extra_state(self) -> torch.Tensor:
        """Save before checkpointing."""
        state = None
528

529
        fp8_checkpoint = self.fp8_meta["fp8_checkpoint"] or self.fp8 or self.fp8_calibration
530
531

        if fp8_checkpoint:
532
            state = {}
533
534
535
536
537
538
            state["scale_fwd"] = self.fp8_meta["scaling_fwd"].scale
            state["scale_inv_fwd"] = self.fp8_meta["scaling_fwd"].scale_inv
            state["amax_history_fwd"] = self.fp8_meta["scaling_fwd"].amax_history
            state["scale_bwd"] = self.fp8_meta["scaling_bwd"].scale
            state["scale_inv_bwd"] = self.fp8_meta["scaling_bwd"].scale_inv
            state["amax_history_bwd"] = self.fp8_meta["scaling_bwd"].amax_history
539
540
541
542

            # Store other pickelable values.
            extra = {}
            for k, v in self.fp8_meta.items():
543
544
545
                if k != "buffer_index_and_autocast_key" and isinstance(
                    v, (bool, int, float, str, tuple, list)
                ):
546
547
548
                    extra[k] = v
            state["extra_fp8_variables"] = extra

549
550
551
552
553
        if is_in_onnx_export_mode():
            state_serialized = torch.frombuffer(pickle.dumps(state), dtype=torch.uint8)
        else:
            state_serialized = io.BytesIO()
            torch.save(state, state_serialized)
554

555
        return state_serialized
556
557
558
559
560
561
562
563

    def set_extra_state(self, state: torch.Tensor) -> None:
        """Load previous state."""
        if state is None:
            return

        if isinstance(state, torch.Tensor):
            state = pickle.loads(state.detach().cpu().numpy().tobytes())
564
565
        elif isinstance(state, io.BytesIO):
            state.seek(0)
566
            state = torch.load(state, map_location="cuda")
567
568
        else:
            raise RuntimeError("Unsupported checkpoint format.")
569
570
571

        if state is None:
            return
572
573
574
575
576
577
578
579
580
581
582
583
584

        # Load extra items.
        self.fp8_meta.update(state["extra_fp8_variables"])
        self.fp8_meta["recipe"].amax_history_len = state["amax_history_fwd"].shape[0]
        if "global_fp8_buffer_pos_fwd_recompute" in self.fp8_meta:
            del self.fp8_meta["global_fp8_buffer_pos_fwd_recompute"]

        # Initialize before loading.
        self.init_fp8_meta_tensors()
        self.fp8_meta["scaling_fwd"].scale.copy_(state["scale_fwd"])
        self.fp8_meta["scaling_fwd"].amax_history.copy_(state["amax_history_fwd"])
        self.fp8_meta["scaling_bwd"].scale.copy_(state["scale_bwd"])
        self.fp8_meta["scaling_bwd"].amax_history.copy_(state["amax_history_bwd"])
585
586
        self.fp8_meta["scaling_fwd"].scale_inv.copy_(state["scale_inv_fwd"])
        self.fp8_meta["scaling_bwd"].scale_inv.copy_(state["scale_inv_bwd"])
587
588
589
590
591
592
593
594
595

    def set_activation_dtype(self, inp: torch.Tensor) -> None:
        """Get activation data type for AMP."""
        # Native AMP (`torch.autocast`) gets highest priority
        if torch.is_autocast_enabled():
            self.activation_dtype = torch.get_autocast_gpu_dtype()
            return

        # All checks after this have already been performed once, thus skip
596
        if hasattr(self, "activation_dtype") and self.activation_dtype == inp.dtype:
597
598
            return

599
600
601
602
603
604
605
606
        dtype = inp.dtype
        for name, param in self.named_parameters():
            if param is not None:
                assert dtype == param.dtype, (
                    "Data types for parameters must match when outside of autocasted region. "
                    f" Found input dtype: {dtype} and {name!r} dtype: {param.dtype}"
                )
        self.activation_dtype = dtype
607
608

    def set_tensor_parallel_group(self, tp_group: Union[dist_group_type, None]) -> None:
609
610
611
612
613
614
615
616
617
        """
        Set the tensor parallel group for the given
        module before executing the forward pass.

        Parameters
        ----------
        tp_group : ProcessGroup, default = `None`
                  tensor parallel process group.
        """
618
619
620
        self.tp_group = tp_group
        self.tp_group_initialized = True

621
622
623
    def _get_fp8_params(self) -> Union[List[torch.Tensor], None]:
        """returns the FP8 weights."""
        fp8_params = []
624
        for param in self.parameters(recurse=False):
625
626
627
628
629
630
            if isinstance(param, Float8Tensor) and param.requires_grad:
                fp8_params.append(param)
        if len(fp8_params) == 0:
            return None
        return fp8_params

631
632
    # This routine is shared across FP8 and FP8_calibration paths so should not actually
    # assume FP8 execution.
633
    def init_fp8_metadata(self, num_gemms: int = 1) -> None:
634
        """Initialize fp8 related metadata and tensors during fprop."""
635
        self.fp8_parameters = FP8GlobalStateManager.with_fp8_parameters()
636
637
        self.fp8 = FP8GlobalStateManager.is_fp8_enabled()
        self.fp8_calibration = FP8GlobalStateManager.is_fp8_calibration()
638
        self.fp8_meta["fp8_checkpoint"] = self.fp8 or self.fp8_calibration
639

640
641
642
643
        if self.fp8_parameters and not self.fp8_initialized:
            self.fp8_meta["num_gemms"] = num_gemms
            self.init_fp8_meta_tensors()

644
645
        if self.fp8 or self.fp8_calibration:
            # FP8 init has already been run and recipe is the same, don't do anything.
646
647
648
649
            if (
                self.fp8_initialized
                and FP8GlobalStateManager.get_fp8_recipe() == self.fp8_meta["recipe"]
            ):
650
651
652
                return

            # Set FP8, recipe, and other FP8 metadata
653
            self.fp8_meta["recipe"] = FP8GlobalStateManager.get_fp8_recipe()
654
            self.fp8_meta["num_gemms"] = num_gemms
655
            self.fp8_meta["fp8_group"] = FP8GlobalStateManager.get_fp8_group()
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672

            # Set FP8_MAX per tensor according to recipe
            self.fp8_meta["fp8_max_fwd"] = self.fp8_meta["recipe"].fp8_format.value.max_fwd
            self.fp8_meta["fp8_max_bwd"] = self.fp8_meta["recipe"].fp8_format.value.max_bwd

            # Allocate scales and amaxes
            self.init_fp8_meta_tensors()
            self.fp8_initialized = True
        else:
            # If fp8 isn't enabled, turn off and return.
            self.fp8_initialized = False
            return

    @contextmanager
    def prepare_forward(
        self,
        inp: torch.Tensor,
673
        is_first_microbatch: Union[bool, None],  # pylint: disable=unused-argument
674
        num_gemms: int = 1,
675
        allow_non_contiguous: bool = False,
Jan Bielak's avatar
Jan Bielak committed
676
    ) -> Generator[torch.Tensor, None, None]:
677
678
679
680
681
682
683
684
        """Checks and prep for FWD.
        The context manager is needed because there isn't a way for a module to know
        if it's the last FP8 module in the forward autocast. It is useful
        to setup the forward aggregated amax reduction for every module
        just in case. The autocast exit will pick up the most recent one.
        """
        # Activation recomputation is used and this is the second forward phase.
        if self.fp8 and in_fp8_activation_recompute_phase():
685
            FP8GlobalStateManager.get_old_fp8_meta_tensors_for_recompute(self.fp8_meta)
686
687
688
689
690
691
692
        else:
            assert inp.is_cuda, "TransformerEngine needs CUDA."

            if self.tp_size > 1:
                assert self.tp_group_initialized, "TP group not initialized."

            self.set_activation_dtype(inp)
693
            self.init_fp8_metadata(num_gemms=num_gemms)
694

695
            if self.fp8 and self.sequence_parallel:
696
697
698
699
                assert self.fp8_meta["recipe"].reduce_amax, (
                    "Amax reduction across tensor parallel group is "
                    "necessary when using sequence parallelism with FP8."
                )
700

701
702
            if self.fp8 and not FP8GlobalStateManager.fp8_graph_capturing():
                FP8GlobalStateManager.add_fp8_tensors_to_global_buffer(
703
704
                    self.fp8_meta, fp8_weights=self._get_fp8_params()
                )
705
706

            # Activation recomputation is used and this is the first forward phase.
707
            if self.fp8 and self.training and is_fp8_activation_recompute_enabled():
708
                FP8GlobalStateManager.copy_forward_fp8_meta_tensors_for_recompute(self.fp8_meta)
709
710

        with torch.cuda.nvtx.range(self.__class__.__name__ + " forward"):
711
712
713
714
            if not allow_non_contiguous:
                yield inp.contiguous()
            else:
                yield inp
715
716

        if self.fp8 and in_fp8_activation_recompute_phase():
717
            FP8GlobalStateManager.restore_fp8_meta_tensors(self.fp8_meta)
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
            return

    def set_nccl_overlap_warning_if_tp(self) -> None:
        """When using TP, the NCCL communication needs to be scheduled
        before the GEMM for there to be a guaranteed overlap. From the
        host side in TE, the comm calls are always launched first, but
        to ensure that the GEMM isn't scheduled first, the environment
        variable `CUDA_DEVICE_MAX_CONNECTIONS` needs to be set to 1 to
        force a single channel.
        """
        if self.tp_size == 1:
            return
        num_cuda_work_queues = int(os.getenv("CUDA_DEVICE_MAX_CONNECTIONS", "0"))
        if num_cuda_work_queues != 1:
            warnings.warn(
                "To guarantee overlapping TP and SP collectives with the backward"
                "GEMMs, set environment variable CUDA_DEVICE_MAX_CONNECTIONS = 1"
            )

    @staticmethod
    def grad_output_preprocess(
        ctx, grad_output: torch.Tensor, row_parallel_mode: bool
    ) -> Tuple[Union[torch.Tensor, None], ...]:
        """Utility function for backward.
        Returns tuple in order (all optional/None based on training precion/recipe):
            R1: gathered `grad_output` in higher precision.
            R2: gathered `grad_output` in FP8.
            R3: R2 transposed.
            R4: bias gradient on R1.

        """
749
750
751
752
753
        if isinstance(grad_output, Float8Tensor):
            grad_output._data = grad_output._data.contiguous()
        else:
            grad_output = grad_output.contiguous()
        grad_output_mat = grad_output.view(-1, grad_output.shape[-1])
754
755
756
757
758
        gather_grad_output = row_parallel_mode and ctx.sequence_parallel

        # No-FP8 case: bgrad is fused with wgrad for this case.
        if not ctx.fp8:
            if gather_grad_output:
759
                if not ctx.ub_overlap_ag:
760
                    grad_output_mat, _ = gather_along_first_dim(grad_output_mat, ctx.tp_group)
761
762
763
764
765
                else:
                    ctx.ub_obj_gradout.copy_input_to_ubuf(grad_output, True)
                    grad_output_mat = ctx.ub_obj_gradout.get_ubuf_output(1)
            return grad_output_mat, None, None, None

766
        fp8_dtype_backward = get_fp8_te_dtype(ctx.fp8_meta["recipe"], fprop_tensor=False)
767
768

        # FP8 case with non-FP8 wgrad
769
        if gather_grad_output and ctx.fp8_meta["recipe"].override_linear_precision.wgrad:
770
            assert (
771
                not ctx.ub_overlap_ag
772
            ), "override_linear_precision.wgrad not supported with UB AG overlap"
773
774
775
776
777
778
779
            grad_output_mat, _ = gather_along_first_dim(grad_output_mat, ctx.tp_group)
        # FP8 case with gather: unfused bgrad, cast, transpose for efficient gather
        elif gather_grad_output:
            if ctx.use_bias:
                grad_bias = grad_output_mat.sum(dim=0)
            else:
                grad_bias = None
780
            if ctx.ub_overlap_ag:
781
782
783
                grad_output_c = ctx.ub_obj_gradout.get_ubuf_output(0)
            else:
                grad_output_c = torch.empty_like(grad_output_mat, dtype=torch.uint8)
784
785
786
787
788
789
790
791
792
            if not isinstance(grad_output_mat, Float8Tensor):
                cast_to_fp8(
                    grad_output_mat,
                    ctx.fp8_meta["scaling_bwd"],
                    tex.FP8BwdTensors.GRAD_OUTPUT1,
                    fp8_dtype_backward,
                    out=grad_output_c,
                )
            else:
793
                grad_output_c = grad_output_mat
794
            if not ctx.ub_overlap_ag:
795
                grad_output_c, _ = gather_along_first_dim(grad_output_c, ctx.tp_group)
796
797
798
799
                if not isinstance(grad_output_c, Float8Tensor):
                    grad_output_t = tex.fp8_transpose(grad_output_c, fp8_dtype_backward)
                else:
                    grad_output_t = grad_output_c.transpose_2d()
800
801
802
803
804
805
806
807
            else:
                grad_output_c = ctx.ub_obj_gradout.get_ubuf_output(1)
                grad_output_t = None

            return grad_output_mat, grad_output_c, grad_output_t, grad_bias

        # FP8 case without gather: cast, transpose, bgrad fused
        if ctx.use_bias:
808
809
810
            grad_output_mat_no_fp8 = grad_output_mat
            if isinstance(grad_output_mat, Float8Tensor):
                grad_output_mat_no_fp8 = grad_output_mat.from_float8(grad_output_mat.dtype)
811
            grad_bias, grad_output_c, grad_output_t = fp8_cast_transpose_bgrad_fused(
812
                grad_output_mat_no_fp8,
813
814
815
816
817
818
                ctx.fp8_meta["scaling_bwd"],
                tex.FP8BwdTensors.GRAD_OUTPUT1,
                fp8_dtype_backward,
            )
        else:
            if not ctx.fp8_meta["recipe"].override_linear_precision.wgrad:
819
820
821
822
823
824
825
826
827
828
                if isinstance(grad_output_mat, Float8Tensor):
                    grad_output_c = grad_output_mat
                    grad_output_t = grad_output_c.transpose_2d()
                else:
                    grad_output_c, grad_output_t = fp8_cast_transpose_fused(
                        grad_output_mat,
                        ctx.fp8_meta["scaling_bwd"],
                        tex.FP8BwdTensors.GRAD_OUTPUT1,
                        fp8_dtype_backward,
                    )
829
830
            else:
                grad_output_t = None
831
832
833
834
835
836
837
838
839
                if not isinstance(grad_output_mat, Float8Tensor):
                    grad_output_c = cast_to_fp8(
                        grad_output_mat,
                        ctx.fp8_meta["scaling_bwd"],
                        tex.FP8BwdTensors.GRAD_OUTPUT1,
                        fp8_dtype_backward,
                    )
                else:
                    grad_output_c = grad_output_mat
840
841
842
843
            grad_bias = None

        return grad_output_mat, grad_output_c, grad_output_t, grad_bias

844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
    def register_parameter(self, name, param, **kwargs):
        """
        Thin wrapper around PyTorch parameter registration to stash additional parameter
        metedata used in deferred initialization.
        """
        super().register_parameter(name, param)
        self.param_init_meta[name] = _ParameterInitMeta(**kwargs)

    def reset_parameters(self, defer_init: Optional[bool] = False) -> None:
        """
        Reset all module parameters to initial values. Unless deferred initialization
        is specified, all parameters on a 'meta' device are also materialized on a real cuda
        device before the values are reset to initial.
        """
        if defer_init:
            return

        for name, param in self.named_parameters(recurse=False):
            # Ensure parameter is on a real device
863
864
            if param.device == torch.device("meta"):
                param = torch.empty_like(param, device="cuda")
865
866
867
868
869
870
871

            # Initialize the parameter values on device
            init_fn = self.param_init_meta[name].init_fn
            get_rng_state_tracker = self.param_init_meta[name].get_rng_state_tracker
            if get_rng_state_tracker is None:
                init_fn(param)
            else:
872
873
874
875
876
877
                if hasattr(self, "rng_tracker_name") and self.rng_tracker_name:
                    with get_rng_state_tracker().fork(self.rng_tracker_name):
                        init_fn(param)
                else:
                    with get_rng_state_tracker().fork():
                        init_fn(param)
878
879
880
881

            # If primary weights are in fp8, wrap the parameter as Float8Tensor
            fp8_meta_index = self.param_init_meta[name].fp8_meta_index
            if self.primary_weights_in_fp8 and fp8_meta_index is not None:
882
883
884
885
886
                dummy_amax = torch.empty(
                    (1, 1),
                    dtype=torch.float32,
                    device=param.device,
                )  # Dummy buffer to avoid overwriting amax history
887
888
889
                param = Float8Tensor.to_float8(
                    param,
                    fp8_meta=self.fp8_meta,
890
                    fp8_meta_index=fp8_meta_index,
891
892
                    amax=dummy_amax,
                    with_transpose_cache=torch.is_grad_enabled(),
893
894
895
896
897
898
899
900
                )

            # Redo parameter wrap in case we broke it above
            # NOTE: Currently this can only be broken when primary weights are in Fp8 but
            #       re-applying the nn.Parameter() wrap is a no-op when the input is already
            #       a parameter so we always re-apply it just for extra safety.
            setattr(self, name, torch.nn.Parameter(param))

901
902
903
    @abstractmethod
    def forward(self):
        """Needs override."""
904

905
    def get_fp8_workspace(
906
        self,
907
908
909
910
911
912
913
        *,
        tensor: Optional[torch.Tensor] = None,
        fp8_meta_forward: Optional[bool] = None,
        fp8_meta_index: Optional[int] = None,
        cache_name: Optional[str] = None,
        update_workspace: bool = True,
        skip_update_flag: Optional[torch.Tensor] = None,
914
        fsdp_group: dist_group_type = None,
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
    ) -> Float8Tensor:
        """Get FP8 workspace buffer and maybe update its values

        The workspace buffer may be cached for future function calls.

        Parameters
        ----------
        tensor : torch.Tensor, optional
            Values to copy into workspace. Required if the workspace
            is being constructed or updated.
        fp8_meta_forward: bool, optional
            Whether to access FP8 meta tensors for the forward pass or
            backward pass. Required if the workspace is being
            constructed.
        fp8_meta_index: int, optional
            Index to access in FP8 meta tensors. Required if the
            workspace is being constructed.
        cache_name: str, optional
            Key for caching.
        update_workspace: bool, default = `True`
            Update workspace with values from `tensor`.
        skip_update_flag: torch.Tensor, optional
            GPU flag to skip updating the workspace. Take precedence
            over `update_workspace` if provided.
939
940
        fsdp_group: bool, default = None
            FSDP process group that the weights are distributed over.
941
942
        """

943
        # Try getting workspace from cache
944
945
946
        out = None
        if cache_name is not None:
            out = self._fp8_workspaces.get(cache_name, None)
947

948
949
950
951
952
953
954
955
956
957
958
959
        # Gather cached Fp8 workspace if it's distributed
        # NOTE: FSDP sharding is supported only for Fp8 buffers and will not work
        #       for models initialized with Fp8 primary weights.
        if (
            out is not None
            and not isinstance(out, Float8Tensor)
            and fsdp_group is not None
            and out._data.shape != tensor.data.shape
        ):
            _fsdp_gather_tensors(fsdp_group, [tensor.data.shape], out)

        # Construct workspace if needed
960
        if out is None:
961
962

            # FP8 data
963
            if tensor is None or fp8_meta_forward is None or fp8_meta_index is None:
964
965
966
967
968
969
970
971
                raise ValueError(
                    "tensor, fp8_meta_forward, and fp8_meta_index kwargs "
                    "must be provided to construct FP8 workspace"
                )
            fp8_dtype = get_fp8_te_dtype(
                self.fp8_meta["recipe"],
                fprop_tensor=fp8_meta_forward,
            )
972
            data = torch.empty_like(tensor, dtype=torch.uint8)
973
            scale_inv = torch.empty([1], dtype=torch.float32, device=tensor.device)
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991

            # Transpose cache
            with_transpose_cache = torch.is_grad_enabled()
            if (
                not with_transpose_cache
                and is_fp8_activation_recompute_enabled()
                and not in_fp8_activation_recompute_phase()
            ):
                with_transpose_cache = True
            data_transpose = None
            if with_transpose_cache:
                data_transpose = torch.empty(
                    (tensor.size(-1), tensor.numel() // tensor.size(-1)),
                    dtype=torch.uint8,
                    device=tensor.device,
                )

            # Construct FP8 tensor
992
            out = Float8Tensor(
993
                data=data,
994
995
996
997
998
999
                fp8_meta=self.fp8_meta,
                fp8_meta_forward=fp8_meta_forward,
                fp8_meta_index=fp8_meta_index,
                fp8_dtype=fp8_dtype,
                fp8_scale_inv=scale_inv,
                dtype=tensor.dtype,
1000
                data_transpose=data_transpose,
1001
            )
1002
1003

            # Update cache
1004
1005
1006
1007
1008
1009
1010
1011
1012
1013
            if cache_name is not None:
                self._fp8_workspaces[cache_name] = out
            update_workspace = True
            skip_update_flag = None

        # Update workspace if needed
        if skip_update_flag is not None:
            update_workspace = True
        if update_workspace:
            if tensor is None:
1014
                raise ValueError("tensor kwarg must be provided to update FP8 workspace")
1015
1016
1017
1018
1019
1020
1021
1022
1023
            if is_in_onnx_export_mode():
                # ONNX export does not support fused cast-transpose
                # kernel and requires that FP8 scales can be
                # represented with constant ops.
                transpose_cache = out._transpose
                out._transpose = None
                out.quantize_(tensor)
                out._scale_inv.fill_(out._scale_inv.item())
                out._transpose = transpose_cache
1024
            else:
1025
                out.quantize_(tensor, noop_flag=skip_update_flag)
1026
1027

        return out
1028

1029
1030
1031
    def _load_from_state_dict(
        self, state_dict, prefix, local_metadata, strict, missing_keys, unexpected_keys, error_msgs
    ):
1032
1033
1034
1035
1036
1037
1038
1039
1040
1041
1042
1043
1044
1045
        """
        This function loads tensors and extra state including fp8 metadata.
        This metadata is essential for copying fp8 tensors, as the copy_ function
        uses the scale_inv parameter from fp8_meta to set the correct scaling factor
        for the new tensor.
        Hence, this extra state must be loaded before the tensor copying process,
        not after, as is typically done in _load_from_state_dict.
        Tensors are copied into fp8 tensors only when self.primary_weights_in_fp8=True,
        otherwise, this behavior is not required.
        """
        if self.primary_weights_in_fp8:
            extra_state_key = prefix + torch.nn.modules.module._EXTRA_STATE_KEY_SUFFIX
            if extra_state_key in state_dict:
                self.set_extra_state(state_dict[extra_state_key])
1046
1047
1048
        super()._load_from_state_dict(
            state_dict, prefix, local_metadata, strict, missing_keys, unexpected_keys, error_msgs
        )