triton.language.extra.cann.libdevice.logb

目录

triton.language.extra.cann.libdevice.logb#

triton.language.extra.cann.libdevice.logb(arg0: None, _semantic: None = 'None')#

Extracts the exponent value of a floating-point number.

参数:

arg0 (scalar or tl.tensor) -- x. Supported dtype(s): float32.

返回:

The exponent value of the input parameter.

返回类型:

float32

示例

import os

os.environ.setdefault("TRITON_ENABLE_LIBDEVICE_SIMT", "1")

import pytest
import math
import triton
import triton.language as tl
import triton.language.extra.cann.libdevice as libdevice
import torch
from triton.backends.ascend.utils import triton_enable_libdevice_simt

_SIMT_SKIP_MSG = ("SIMT libdevice ops require an Ascend 950 target "
                  "with TRITON_ENABLE_LIBDEVICE_SIMT=1; skipping.")


def torch_logb_reference(x0):
    assert x0.device.type == "cpu"
    assert x0.dtype == torch.float32
    original_shape = x0.shape
    exponents = [math.frexp(value.item())[1] - 1 if value.item() != 0 else 0 for value in x0.flatten()]
    result = torch.tensor(exponents, dtype=torch.int32).reshape(original_shape)
    return result


@triton.jit
def triton_kernel(input0, output, n_elements, XBLOCK: tl.constexpr, XBLOCK_SUB: tl.constexpr):
    offset = tl.program_id(0) * XBLOCK
    base = tl.arange(0, XBLOCK_SUB)
    loops: tl.constexpr = XBLOCK // XBLOCK_SUB
    for loop in range(loops):
        x0 = offset + (loop * XBLOCK_SUB) + base
        mask = x0 < n_elements
        tmp0 = tl.load(input0 + (x0), mask=mask)
        tmp1 = libdevice.logb(tmp0)
        tl.store(output + (x0), tmp1, mask=mask)


@pytest.mark.skipif(not triton_enable_libdevice_simt(), reason=_SIMT_SKIP_MSG)
def test_logb():
    x0 = (torch.rand((8, )) + 0.1).to(torch.float32)
    expected = (torch_logb_reference(x0).to(torch.float32)).npu()
    x0 = x0.npu()
    output = torch.empty(8, dtype=torch.float32, device='npu')
    triton_kernel[(1, )](x0, output, 8, XBLOCK=8, XBLOCK_SUB=8, compile_mode='simt_only')
    output = output.cpu()
    expected = expected.cpu()
    torch.testing.assert_close(output, expected, rtol=1e-03, atol=1e-03, equal_nan=True)


if __name__ == "__main__":
    if not triton_enable_libdevice_simt():
        print(_SIMT_SKIP_MSG)
    else:
        test_logb()

特殊说明

Platform and compilation modes:

  • Ascend 950 supports SIMD, SIMT.

目录