* Copyright (c) 2025 Huawei Technologies Co., Ltd.
* This program is free software, you can redistribute it and/or modify it under the terms and conditions of
* CANN Open Software License Agreement Version 2.0 (the "License").
* Please refer to the License for details. You may not use this file except in compliance with the License.
* THIS SOFTWARE IS PROVIDED ON AN "AS IS" BASIS, WITHOUT WARRANTIES OF ANY KIND, EITHER EXPRESS OR IMPLIED,
* INCLUDING BUT NOT LIMITED TO NON-INFRINGEMENT, MERCHANTABILITY, OR FITNESS FOR A PARTICULAR PURPOSE.
* See LICENSE in the root of the software repository for the full text of the License.
*/
#include "catlass/gemm/kernel/optimized_matmul.hpp"
#include <acl/acl.h>
#include "catlass/arch/arch.hpp"
#include "catlass/catlass.hpp"
#include "catlass/gemm/block/block_mmad.hpp"
#include "catlass/gemm/block/block_swizzle.hpp"
#include "catlass/gemm/device/device_gemm.hpp"
#include "catlass/gemm/dispatch_policy.hpp"
#include "catlass/gemm/gemm_type.hpp"
#include "catlass/layout/layout.hpp"
#include "catlass/status.hpp"
#include "catlass_kernel.h"
#include "common.hpp"
namespace Catlass {
template <
class ArchTag,
class AType,
class BType,
class CType,
class BiasType = void>
struct TileCopyOpt : public Catlass::Gemm::Tile::TileCopy<ArchTag, AType, BType, CType, BiasType> {
using Base = Catlass::Gemm::Tile::TileCopy<ArchTag, AType, BType, CType, BiasType>;
using ElementA = typename Base::ElementA;
using ElementB = typename Base::ElementB;
using ElementAccumulator = typename Base::ElementAccumulator;
using CopyGmToL1A = typename Base::CopyGmToL1A;
using CopyGmToL1B = typename Base::CopyGmToL1B;
using CopyL1ToL0A = typename Base::CopyL1ToL0A;
using CopyL1ToL0B = typename Base::CopyL1ToL0B;
using CopyL0CToGm = typename Base::CopyL0CToGm;
using BiasTypeSelector = typename Base::BiasTypeSelector;
using CopyGmToL1Bias = typename Base::CopyGmToL1Bias;
using CopyL1ToBT = typename Base::CopyL1ToBT;
};
}
namespace CatlassKernel {
using namespace Catlass;
template <class LayoutA, class LayoutB, class LayoutC, class InDType, class OutDType, bool M_GT_N>
void OptimizedMatmulImpl(const uint32_t blockNum, aclrtStream stream, const KernelInfo &kernelInfo) {
using ArchTag = Arch::AtlasA2;
constexpr uint32_t alignByByte = 512;
constexpr uint32_t alignByElement = alignByByte / sizeof(InDType);
constexpr bool ENABLE_UNIT_FLAG = true;
constexpr bool ENABLE_SHUFFLE_K = true;
using ElementA = InDType;
using ElementB = InDType;
using ElementC = OutDType;
using AType = Gemm::GemmType<ElementA, LayoutA>;
using BType = Gemm::GemmType<ElementB, LayoutB>;
using CType = Gemm::GemmType<ElementC, LayoutC>;
using DispatchPolicy = Gemm::MmadAtlasA2Preload<ENABLE_UNIT_FLAG, ENABLE_SHUFFLE_K>;
using PaddingTag = Catlass::Gemm::Kernel::PaddingTag;
constexpr PaddingTag paddingTagA = (std::is_same_v<LayoutA, layout::zN> || std::is_same_v<LayoutA, layout::nZ>)
? PaddingTag::NO_PADDING
: PaddingTag::PADDING_BLOCK_ND;
constexpr PaddingTag paddingTagB = (std::is_same_v<LayoutB, layout::zN> || std::is_same_v<LayoutB, layout::nZ>)
? PaddingTag::NO_PADDING
: PaddingTag::PADDING_BLOCK_ND;
static const uint32_t COMPUTE_LENGTH_A = 96 * 1024 / sizeof(ElementA);
using PaddingBuilderA = Catlass::Gemm::Kernel::PaddingBuilder<
paddingTagA, ArchTag, ElementA, LayoutA, COMPUTE_LENGTH_A>;
using GlobalPaddingA = typename PaddingBuilderA::Padding;
static const uint32_t COMPUTE_LENGTH_B = 96 * 1024 / sizeof(ElementB);
using PaddingBuilderB = Catlass::Gemm::Kernel::PaddingBuilder<
paddingTagB, ArchTag, ElementB, LayoutB, COMPUTE_LENGTH_B>;
using GlobalPaddingB = typename PaddingBuilderB::Padding;
uint32_t fftsLen{0};
uint64_t fftsAddr{0};
RT_CHECK(rtGetC2cCtrlAddr(&fftsAddr, &fftsLen));
GemmCoord problemShape{kernelInfo.m, kernelInfo.n, kernelInfo.k};
LayoutA layoutA{kernelInfo.m, kernelInfo.k};
LayoutB layoutB{kernelInfo.k, kernelInfo.n};
LayoutC layoutC{kernelInfo.m, kernelInfo.n};
uint8_t *deviceA = kernelInfo.inputAddr.at(0);
uint8_t *deviceB = kernelInfo.inputAddr.at(1);
uint8_t *deviceC = kernelInfo.outputAddr.at(0);
using L1TileShape =
std::conditional_t<std::is_same_v<LayoutA, layout::ColumnMajor> && std::is_same_v<LayoutB, layout::ColumnMajor>,
GemmShape<256, 128, 256>, GemmShape<128, 256, 256>>;
using L0TileShape =
std::conditional_t<std::is_same_v<LayoutA, layout::ColumnMajor> && std::is_same_v<LayoutB, layout::ColumnMajor>,
GemmShape<256, 128, 64>, GemmShape<128, 256, 64>>;
using BlockScheduler = std::conditional_t<M_GT_N, Catlass::Gemm::Block::GemmIdentityBlockSwizzle<3, 0>,
Catlass::Gemm::Block::GemmIdentityBlockSwizzle<3, 1>>;
using BlockEpilogue = void;
bool isNeedPaddingA = IsNeedPadding(layoutA, alignByElement);
bool isNeedPaddingB = IsNeedPadding(layoutB, alignByElement);
if (isNeedPaddingA && isNeedPaddingB) {
using LayoutMmadA = typename PaddingBuilderA::LayoutAfterPadding;
using LayoutMmadB = typename PaddingBuilderB::LayoutAfterPadding;
using ATypeMmad = Gemm::GemmType<ElementA, LayoutMmadA>;
using BTypeMmad = Gemm::GemmType<ElementB, LayoutMmadB>;
using TileCopy = TileCopyOpt<ArchTag, ATypeMmad, BTypeMmad, CType>;
using BlockMmadOpt = Gemm::Block::BlockMmad<DispatchPolicy, L1TileShape, L0TileShape, ATypeMmad, BTypeMmad,
CType, void, TileCopy>;
using MatmulKernel =
Gemm::Kernel::OptimizedMatmul<GlobalPaddingA, GlobalPaddingB, BlockMmadOpt, BlockEpilogue, BlockScheduler>;
typename MatmulKernel::Arguments arguments{problemShape, deviceA, deviceB, deviceC};
using MatmulAdapter = Gemm::Device::DeviceGemm<MatmulKernel>;
MatmulAdapter matmulOp;
RunAdapter(matmulOp, arguments, stream, blockNum, fftsAddr);
} else if (isNeedPaddingA) {
using LayoutMmadA = typename PaddingBuilderA::LayoutAfterPadding;
using ATypeMmad = Gemm::GemmType<ElementA, LayoutMmadA>;
using TileCopy = TileCopyOpt<ArchTag, ATypeMmad, BType, CType>;
using BlockMmadOpt =
Gemm::Block::BlockMmad<DispatchPolicy, L1TileShape, L0TileShape, ATypeMmad, BType, CType, void, TileCopy>;
using MatmulKernel =
Gemm::Kernel::OptimizedMatmul<GlobalPaddingA, void, BlockMmadOpt, BlockEpilogue, BlockScheduler>;
typename MatmulKernel::Arguments arguments{problemShape, deviceA, deviceB, deviceC};
using MatmulAdapter = Gemm::Device::DeviceGemm<MatmulKernel>;
MatmulAdapter matmulOp;
RunAdapter(matmulOp, arguments, stream, blockNum, fftsAddr);
} else if (isNeedPaddingB) {
using LayoutMmadB = typename PaddingBuilderB::LayoutAfterPadding;
using BTypeMmad = Gemm::GemmType<ElementB, LayoutMmadB>;
using TileCopy = TileCopyOpt<ArchTag, AType, BTypeMmad, CType>;
using BlockMmadOpt =
Gemm::Block::BlockMmad<DispatchPolicy, L1TileShape, L0TileShape, AType, BTypeMmad, CType, void, TileCopy>;
using MatmulKernel =
Gemm::Kernel::OptimizedMatmul<void, GlobalPaddingB, BlockMmadOpt, BlockEpilogue, BlockScheduler>;
typename MatmulKernel::Arguments arguments{problemShape, deviceA, deviceB, deviceC};
using MatmulAdapter = Gemm::Device::DeviceGemm<MatmulKernel>;
MatmulAdapter matmulOp;
RunAdapter(matmulOp, arguments, stream, blockNum, fftsAddr);
} else {
using TileCopy = TileCopyOpt<ArchTag, AType, BType, CType>;
using BlockMmadOpt =
Gemm::Block::BlockMmad<DispatchPolicy, L1TileShape, L0TileShape, AType, BType, CType, void, TileCopy>;
using MatmulKernel = Gemm::Kernel::OptimizedMatmul<void, void, BlockMmadOpt, BlockEpilogue, BlockScheduler>;
typename MatmulKernel::Arguments arguments{problemShape, deviceA, deviceB, deviceC};
using MatmulAdapter = Gemm::Device::DeviceGemm<MatmulKernel>;
MatmulAdapter matmulOp;
RunAdapter(matmulOp, arguments, stream, blockNum, fftsAddr);
}
}
void OptimizedMatmul(const uint32_t blockNum, aclrtStream stream, const KernelInfo &kernelInfo) {
if (!kernelInfo.transA && !kernelInfo.transB && kernelInfo.inputDataType == ACL_FLOAT16
&& kernelInfo.outputDataType == ACL_FLOAT16) {
if (kernelInfo.m > kernelInfo.n) {
OptimizedMatmulImpl<layout::RowMajor, layout::RowMajor, layout::RowMajor, half, half, true>(blockNum,
stream,
kernelInfo);
} else {
OptimizedMatmulImpl<layout::RowMajor, layout::RowMajor, layout::RowMajor, half, half, false>(blockNum,
stream,
kernelInfo);
}
}
}
}