aclnnPromptFlashAttentionV3

📄 查看源码

产品支持情况

产品 是否支持
Ascend 950PR/Ascend 950DT
Atlas A3 训练系列产品/Atlas A3 推理系列产品
Atlas A2 训练系列产品/Atlas A2 推理系列产品
Atlas 200I/500 A2 推理产品 ×
Atlas 推理系列产品
Atlas 训练系列产品 ×

功能说明

  • 接口功能:全量推理场景的FlashAttention算子,支持sparse优化、actualSeqLengthsKv优化、int8量化功能、innerPrecise参数(用于支持高精度或者高性能模式选择)。

  • 计算公式:

    self-attention(自注意力)利用输入样本自身的关系构建了一种注意力模型。其原理是假设有一个长度为nn的输入样本序列xxxx的每个元素都是一个dd维向量,可以将每个dd维向量看作一个token embedding,将这样一条序列经过3个权重矩阵变换得到3个维度为n∗dn*d的矩阵。

    self-attention的计算公式一般定义如下,其中QQKKVV为输入样本的重要属性元素,是输入样本经过空间变换得到,且可以统一到一个特征空间中。公式及算子名称中的"Attention"为"self-attention"的简写。

    Attention(Q,K,V)=Score(Q,K)VAttention(Q,K,V)=Score(Q,K)V

    本算子中Score函数采用Softmax函数,self-attention计算公式为:

    Attention(Q,K,V)=Softmax(QKTd)VAttention(Q,K,V)=Softmax(\frac{QK^T}{\sqrt{d}})V

    其中:QQKTK^T的乘积代表输入xx的注意力,为避免该值变得过大,通常除以dd的开根号进行缩放,并对每行进行softmax归一化,与VV相乘后得到一个n∗dn*d的矩阵。

函数原型

算子执行接口为两段式接口,必须先调用“aclnnPromptFlashAttentionV3GetWorkspaceSize”接口获取计算所需workspace大小以及包含了算子计算流程的执行器,再调用“aclnnPromptFlashAttentionV3”接口执行计算。

aclnnStatus aclnnPromptFlashAttentionV3GetWorkspaceSize(
    const aclTensor   *query,
    const aclTensor   *key,
    const aclTensor   *value,
    const aclTensor   *pseShift,
    const aclTensor   *attenMask,
    const aclIntArray *actualSeqLengths,
    const aclIntArray *actualSeqLengthsKv,
    const aclTensor   *deqScale1,
    const aclTensor   *quantScale1,
    const aclTensor   *deqScale2,
    const aclTensor   *quantScale2,
    const aclTensor   *quantOffset2,
    int64_t            numHeads,
    double             scaleValue,
    int64_t            preTokens,
    int64_t            nextTokens,
    char              *inputLayout,
    int64_t            numKeyValueHeads,
    int64_t            sparseMode,
    int64_t            innerPrecise,
    const aclTensor   *attentionOut,
    uint64_t          *workspaceSize,
    aclOpExecutor     **executor)
aclnnStatus aclnnPromptFlashAttentionV3(
     void              *workspace,
     uint64_t           workspaceSize,
     aclOpExecutor     *executor,
     const aclrtStream  stream)

aclnnPromptFlashAttentionV3GetWorkspaceSize

  • 参数说明

    参数名 输入/输出 描述 使用说明 数据类型 数据格式 维度(shape) 非连续Tensor
    query 输入 公式中的输入Q。 保持与key、value的数据类型一致。 FLOAT16、BFLOAT16、INT8 ND 3-4 ×
    key 输入 公式中的输入K。 保持与query、value的数据类型一致。 FLOAT16、BFLOAT16、INT8 ND 3-4 ×
    value 输入 公式中的输入V。 保持与query、key的数据类型一致。 FLOAT16、BFLOAT16、INT8 ND 3-4 ×
    pseShift 输入 位置编码
    • 不使用该功能时可传入nullptr。
    FLOAT16、BFLOAT16 ND 4 ×
    attenMask 输入 mask矩阵
    • 不使用该功能可传入nullptr。
    • 综合约束请见约束说明
    BOOL、INT8、UINT8 ND 2-4 ×
    actualSeqLengths 输入 不同Batch中query的有效序列长度。
    • 不指定序列长度可传入nullptr。
    • 综合约束请见约束说明
      INT64 ND 1 -
      actualSeqLengthsKv 输入 不同Batch中key/value的有效序列长度。
      • 不指定序列长度可传入nullptr。
      INT64 ND 1 -
      deqScale1 输入 BMM1后面的反量化因子。
      • 支持per-tensor。
      • 不使用该功能时可传入nullptr。
      UINT64、FLOAT32 ND 1 -
      quantScale1 输入 BMM2前面的量化因子。
      • 支持per-tensor。
      • 不使用该功能时可传入nullptr。
      FLOAT32 ND 1 -
      deqScale2 输入 BMM2后面的反量化因子。
      • 支持per-tensor。
      • 不使用该功能时可传入nullptr。
      UINT64、FLOAT32 ND 1 -
      quantScale2 输入 输出的量化因子。
      • 支持per-tensor,per-channel。
      • 不使用该功能时可传入nullptr。
      FLOAT32、BFLOAT16 ND 1 -
      quantOffset2 输入 输出的量化偏移。
      • 支持per-tensor,per-channel。
      • 不使用该功能时可传入nullptr。
      FLOAT32、BFLOAT16 ND 1 -
      numHeads 输入 query的head个数。 - INT64 ND 1 -
      scaleValue 输入 公式中d开根号的倒数。
      • 数据类型与query的数据类型需满足数据类型推导规则。
      • 用户不特意指定时建议传入1.0。
      DOUBLE - 1 -
      preTokens 输入 表示attention需要和前几个Token计算关联。 不特意指定时建议传入2147483647。 INT64 - 1 -
      nextTokens 输入 表示attention需要和后几个Token计算关联。 不特意指定时建议传入0。 INT64 - 1 -
      inputLayout 输入 标识输入query、key、value的数据排布格式。
      • 不特意指定时建议传入"BSH"。
      CHAR - - -
      numKeyValueHeads 输入 key、value中head个数。
      • 用户不特意指定时建议传入0,表示key/value和query的head个数相等。
      INT64 - - -
      sparseMode 输入 sparse的模式。 综合约束请见约束说明 INT64 - 1 -
      innerPrecise 输入 表示高精度或者高性能选择。 综合约束请见约束说明 INT8 - 1 -
      attentionOut 输出 公式中的输出。 - FLOAT16、BFLOAT16、INT8 ND 3-4 -
      workspaceSize 输出 返回用户需要在Device侧申请的workspace大小。 - - - 1 -
      executor 输出 返回op执行器,包含了算子计算流程。 - - - 1 -
    • 返回值

      返回aclnnStatus状态码,具体参见aclnn返回码

      第一段接口完成入参校验,若出现以下错误码,则对应原因为:

      返回值 错误码 描述
      ACLNN_ERR_PARAM_NULLPTR 161001 如果传入参数是必选输入,输出或者必选属性,且是空指针,则返回161001。
      ACLNN_ERR_PARAM_INVALID 161002 query、key、value、pseShift、attenMask、attentionOut的数据类型和数据格式不在支持的范围内。
      ACLNN_ERR_RUNTIME_ERROR 361001 API内存调用npu runtime的接口异常。

    aclnnPromptFlashAttentionV3

    • 参数说明

      参数名 输入/输出 描述
      workspace 输入 在Device侧申请的workspace内存地址。
      workspaceSize 输入 在Device侧申请的workspace大小,由第一段接口aclnnPromptFlashAttentionV3GetWorkspaceSize获取。
      executor 输入 op执行器,包含了算子计算流程。
      stream 输入 指定执行任务的Stream。
    • 返回值

      返回aclnnStatus状态码,具体参见aclnn返回码

    约束说明

    • 确定性计算:

      • aclnnPromptFlashAttentionV3默认确定性实现。
    • 该接口与PyTorch配合使用时,需要保证CANN相关包与PyTorch相关包的版本匹配。

    • 入参为空的处理:算子内部需要判断参数query是否为空,如果是空则直接返回。参数query不为空Tensor,参数key、value为空tensor(即S2为0),则attentionOut填充为全零。attentionOut为空Tensor时,AscendCLNN框架会处理。其余在上述参数说明中标注了“可传入nullptr”的入参为空指针时,不进行处理。

    • query、key、value数据排布格式支持从多种维度解读,其中B(Batch)表示输入样本批量大小、S(Seq-Length)表示输入样本序列长度、H(Head-Size)表示隐藏层的大小、N(Head-Num)表示多头数、D(Head-Dim)表示隐藏层最小的单元尺寸,且满足D=H/N,T表示所有Batch输入样本序列长度的累加和。

    • query,key,value输入,功能使用限制如下:

      • Atlas A2 训练系列产品/Atlas A2 推理系列产品、Ascend 950PR/Ascend 950DT:

        • 支持B轴小于等于65536(64k),输入类型包含INT8时D轴非32对齐或输入类型为FLOAT16或BFLOAT16时D轴非16对齐时,B轴仅支持到128。

        • 支持N轴小于等于256。

        • S支持小于等于20971520(20M)。部分长序列场景下,如果计算量过大可能会导致pfa算子执行超时(aicore error类型报错,errorStr为:timeout or trap error),此场景下建议做S切分处理,注:这里计算量会受B、S、N、D等的影响,值越大计算量越大。典型的会超时的长序列(即B、S、N、D的乘积较大)场景包括但不限于:

          B Q_N Q_S D KV_N KV_S
          1 20 2097152 256 1 2097152
          1 2 20971520 256 2 20971520
          20 1 2097152 256 1 2097152
          1 10 2097152 512 1 2097152
        • 支持D轴小于等于512。inputLayout为BSH或者BSND时,要求N*D小于65535。

      • Atlas A2 训练系列产品/Atlas A2 推理系列产品: TND场景下query,key,value输入的综合限制:

        • T小于等于65536。
        • N等于8/16/32/64/128,且Q_N、K_N、V_N相等。
        • Q_D、K_D等于192,V_D等于128/192。
        • 数据类型仅支持BFLOAT16。
        • sparse模式仅支持sparse=0且不传mask,或sparse=3且传入mask。
        • 当sparse=3时,要求每个batch单独的actualSeqLengths < actualSeqLengthsKv。
      • Atlas 推理系列加速卡产品:

        • 支持B轴小于等于128;支持N轴小于等于256;支持S轴小于等于65535(64k), Q_S或KV_S非128对齐,Q_S和KV_S不等长的场景不支持配置atten_mask;支持D轴小于等于512。
      • Atlas A2 训练系列产品/Atlas A2 推理系列产品、Ascend 950PR/Ascend 950DT:数据类型支持FLOAT16、BFLOAT16、INT8。

      • Atlas 推理系列加速卡产品:数据类型仅支持FLOAT16。

    • pseShift功能使用限制如下:

      • 输入shape类型需为(B, N, Q_S, KV_S)或(1, N, Q_S, KV_S),其中Q_S为query的shape中的S,KV_S为key和value的shape中的S。

      • Q_S需大于等于query的S长度,KV_S需大于等于key的S长度。

      • 对于pseShift的KV_S为非32对齐的场景,建议padding到32字节来提高性能,多余部分的填充值不做要求。如不使用该功能时可传入nullptr。

      • Atlas A2 训练系列产品/Atlas A2 推理系列产品、Ascend 950PR/Ascend 950DT:数据类型支持FLOAT16、BFLOAT16,且在pseShift为FLOAT16类型时,要求此时的query为FLOAT16或INT8类型,而在pseShift为BFLOAT16类型时,要求此时的query为BFLOAT16类型。在query、key、value为FLOAT16且pseShift存在的情况下,默认走高精度模式,对应的限制继承自高精度模式的限制。

      • Atlas 推理系列加速卡产品:仅支持nullptr。

    • attenMask功能使用限制如下:

      • 对于attenMask的KV_S为非32对齐的场景,建议padding到32对齐来提高性能,多余部分填充成1。
      • 通常建议shape输入Q_S, KV_S; B, Q_S, KV_S; 1, Q_S, KV_S; B, 1, Q_S, KV_S; 1, 1, Q_S, KV_S,其中Q_S为query的shape中的S,KV_S为key和value的shape中的S。
      • Atlas A2 训练系列产品/Atlas A2 推理系列产品、Ascend 950PR/Ascend 950DT:数据类型支持BOOL、INT8和UINT8。
      • Atlas 推理系列加速卡产品:仅支持BOOL。
    • actualSeqLengths,actualSeqLengthsKv输入,功能使用限制如下:

      • 对于actualSeqLengths,如果不指定序列长度,可以传入nullptr,这表示有效序列长度与query的shape中的S长度相同。需要注意的是,该参数中每个batch的有效序列长度不应超过query中对应batch的序列长度。
      • 对于actualSeqLengthsKv,如果不指定序列长度,可以传入nullptr,这表示有效序列长度与key/value的shape中的S长度相同。需要注意的是,该参数中每个batch的有效序列长度不应超过key/value中对应batch的序列长度。
      • 关于seqlen的传入长度有以下规则:当传入长度为1时,所有Batch将使用相同的seqlen;当传入长度大于或等于Batch数量时,将取seqlen的前Batch个数值;其他长度的传入将不被支持。
      • 当query的inputLayout为TND时,该入参必须传入,且以该入参元素的数量作为Batch值。该入参中每个元素的值表示当前Batch与之前所有Batch的Sequence Length和,因此后一个元素的值必须大于等于前一个元素的值,且不能出现负值。
      • Atlas A2 训练系列产品/Atlas A2 推理系列产品:数据类型支持INT64,支持TND格式。
      • Atlas 推理系列加速卡产品:数据类型支持INT64。
      • Ascend 950PR/Ascend 950DT:数据类型支持INT64。
    • deqScale1,deqScale2输入,功能使用限制如下:

      • Atlas A2 训练系列产品/Atlas A2 推理系列产品、Ascend 950PR/Ascend 950DT:数据类型支持UINT64、FLOAT32。
      • Atlas 推理系列加速卡产品:仅支持nullptr。
    • quantScale1输入,功能使用限制如下:

      • Atlas A2 训练系列产品/Atlas A2 推理系列产品、Ascend 950PR/Ascend 950DT:数据类型支持FLOAT32。
      • Atlas 推理系列加速卡产品:仅支持nullptr。
    • quantScale2,quantOffset2输入,功能使用限制如下:

      • Atlas A2 训练系列产品/Atlas A2 推理系列产品、Ascend 950PR/Ascend 950DT:数据类型支持FLOAT32和BFLOAT16。
      • Atlas 推理系列加速卡产品:仅支持nullptr。
    • preTokens输入,功能使用限制如下:

      • Atlas A2 训练系列产品/Atlas A2 推理系列产品、Ascend 950PR/Ascend 950DT:数据类型支持INT64。
      • Atlas 推理系列加速卡产品:仅支持取值2147483647。
    • nextTokens输入,功能使用限制如下:

      • Atlas A2 训练系列产品/Atlas A2 推理系列产品、Ascend 950PR/Ascend 950DT:数据类型支持INT64。
      • Atlas 推理系列加速卡产品:仅支持取值0和2147483647。
    • inputLayout输入,功能使用限制如下:

      • 当前支持BSH、BSND、BNSD、BNSD_BSND(输入为BNSD时,输出格式为BSND)。用户不特意指定时建议传入"BSH"。
      • Atlas A2 训练系列产品/Atlas A2 推理系列产品:除了上述格式,还支持TND(不支持pse、全量化、后量化)。
    • numKeyValueHeads输入,功能使用限制如下:

      • 需要满足numHeads整除numKeyValueHeads,numHeads与numKeyValueHeads的比值不能大于64,且在BSND、BNSD、BNSD_BSND场景下,需要与shape中的key/value的N轴shape值相同,否则报错。
      • Atlas A2 训练系列产品/Atlas A2 推理系列产品、Ascend 950PR/Ascend 950DT:数据类型支持INT64。
      • Atlas 推理系列加速卡产品:仅支持取值0。
    • sparseMode输入,功能使用限制如下:

      • Atlas A2 训练系列产品/Atlas A2 推理系列产品、Ascend 950PR/Ascend 950DT:
        • sparseMode为0时,代表defaultMask模式,如果attenmask未传入则不做mask操作,忽略preTokens和nextTokens(内部赋值为INT_MAX);如果传入,则需要传入完整的attenmask矩阵(S1 * S2),表示preTokens和nextTokens之间的部分需要计算。
        • sparseMode为1时,代表allMask,必须传入完整的attenmask矩阵(S1 * S2)。
        • sparseMode为2时,代表leftUpCausal模式的mask,需要传入优化后的attenmask矩阵(2048*2048)。
        • sparseMode为3时,代表rightDownCausal模式的mask,对应以右顶点为划分的下三角场景,需要传入优化后的attenmask矩阵(2048*2048)。
        • sparseMode为4时,代表band模式的mask,需要传入优化后的attenmask矩阵(2048*2048)。
        • sparseMode为5、6、7、8时,分别代表prefix、global、dilated、block_local,均暂不支持。用户不特意指定时建议传入0。
      • Atlas A2 训练系列产品/Atlas A2 推理系列产品:当inputLayout为TND时,sparseMode仅支持取值0、3。
      • Atlas 推理系列加速卡产品:仅支持取值0。
    • innerPrecise输入,功能使用限制如下:

      • 一共4种模式:0、1、2、3。一共两位bit位,第0位(bit0)表示高精度或者高性能选择,第1位(bit1)表示是否做行无效修正。

        innerPrecise 模式 行无效修正
        0 高精度 ×
        1 高性能 ×
        2 高精度
        3 高性能
      • Q_S>1时,sparse_mode为0或1,并传入用户自定义mask的情况下,建议开启行无效。

      • BFLOAT16和INT8不区分高精度和高性能,行无效修正对FLOAT16、BFLOAT16和INT8均生效。 当前0、1为保留配置值,当计算过程中“参与计算的mask部分”存在某整行全为1的情况时,精度可能会有损失。此时可以尝试将该参数配置为2或3来使能行无效功能以提升精度,但是该配置会导致性能下降。 如果算子可判断出存在无效行场景,会自动使能无效行计算,例如sparse_mode为3,Sq > Skv场景。

    • attentionOut输出,功能使用限制如下:

      • 当inputLayout为BNSD_BSND时,输入query的shape是BNSD,输出shape为BSND;其余情况该入参的shape需要与入参query的shape保持一致。
      • Atlas A2 训练系列产品/Atlas A2 推理系列产品、Ascend 950PR/Ascend 950DT:数据类型支持FLOAT16、BFLOAT16、INT8。
      • Atlas 推理系列加速卡产品:仅支持FLOAT16。
    • int8量化相关入参数量与输入、输出数据格式的综合限制:

      • 输入为INT8,输出为INT8的场景:入参deqScale1、quantScale1、deqScale2、quantScale2需要同时存在,quantOffset2可选,不传时默认为0。
      • 输入为INT8,输出为FLOAT16的场景:入参deqScale1、quantScale1、deqScale2需要同时存在,若存在入参quantOffset2 或 quantScale2(即不为nullptr),则报错并返回。
      • 输入为FLOAT16或BFLOAT16,输出为INT8的场景:入参quantScale2需存在,quantOffset2可选,不传时默认为0,若存在入参deqScale1 或 quantScale1 或 deqScale2(即不为nullptr),则报错并返回。
      • 入参 quantScale2 和 quantOffset2 支持 per-tensor/per-channel 两种格式和 FLOAT32/BFLOAT16 两种数据类型。若传入 quantOffset2 ,需保证其类型和shape信息与 quantScale2 一致。当输入为BFLOAT16时,同时支持 FLOAT32和BFLOAT16 ,否则仅支持 FLOAT32 。per-channel 格式,当输出layout为BSH时,要求 quantScale2 所有维度的乘积等于H;其他layout要求乘积等于N*D。(建议输出layout为BSH时,quantScale2 shape传入[1,1,H]或[H];输出为BNSD时,建议传入[1,N,1,D]或[N,D];输出为BSND时,建议传入[1,1,N,D]或[N,D])
      • 输出为int8,quantScale2 和 quantOffset2 为 per-channel时,暂不支持左padding、Ring Attention或者D非32Byte对齐的场景。
      • 输出为int8时,暂不支持sparse为band且preTokens/nextTokens为负数。
    • 当输出为INT8,入参quantOffset2传入非空指针和非空tensor值,并且sparseMode、preTokens和nextTokens满足以下条件,矩阵会存在某几行不参与计算的情况,导致计算结果误差,该场景会拦截(解决方案:如果希望该场景不被拦截,需要在PFA接口外部做后量化操作,不在PFA接口内部使能):

      • parseMode = 0,attenMask如果非空指针,每个batch actualSeqLengths - actualSeqLengthsKV - preTokens > 0 或 nextTokens < 0 时,满足拦截条件。
      • sparseMode = 1 或 2,不会出现满足拦截条件的情况。
      • sparseMode = 3,每个batch actualSeqLengthsKV - actualSeqLengths < 0,满足拦截条件。
      • sparseMode = 4,preTokens < 0 或 每个batch nextTokens + actualSeqLengthsKV - actualSeqLengths < 0 时,满足拦截条件。

    调用示例

    示例代码如下,仅供参考,具体编译和执行过程请参考编译与运行样例

    #include <iostream>
    #include <vector>
    #include <cmath>
    #include <cstring>
    #include "acl/acl.h"
    #include "aclnn/opdev/fp16_t.h"
    #include "aclnnop/aclnn_prompt_flash_attention_v3.h"
    #include "securec.h"
     
    using namespace std;
    
    namespace {
     
    #define CHECK_RET(cond) ((cond) ? true :(false))
     
    #define LOG_PRINT(message, ...)     \
      do {                              \
        (void)printf(message, ##__VA_ARGS__); \
      } while (0)
     
    int64_t GetShapeSize(const std::vector<int64_t>& shape) {
      int64_t shapeSize = 1;
      for (auto i : shape) {
        shapeSize *= i;
      }
      return shapeSize;
    }
     
    int Init(int32_t deviceId, aclrtStream* stream) {
      auto ret = aclInit(nullptr);
      if (!CHECK_RET(ret == ACL_SUCCESS)) { 
        LOG_PRINT("aclInit failed. ERROR: %d\n", ret); 
        return ret;
      }
      ret = aclrtSetDevice(deviceId);
      if (!CHECK_RET(ret == ACL_SUCCESS)) { 
        LOG_PRINT("aclrtSetDevice failed. ERROR: %d\n", ret); 
        return ret;
      }
      ret = aclrtCreateStream(stream);
      if (!CHECK_RET(ret == ACL_SUCCESS)) { 
        LOG_PRINT("aclrtCreateStream failed. ERROR: %d\n", ret); 
        return ret;
      }
      return 0;
    }
     
    template <typename T>
    int CreateAclTensor(const std::vector<T>& hostData, const std::vector<int64_t>& shape, void** deviceAddr,
                        aclDataType dataType, aclTensor** tensor) {
      auto size = GetShapeSize(shape) * sizeof(T);
      auto ret = aclrtMalloc(deviceAddr, size, ACL_MEM_MALLOC_HUGE_FIRST);
      if (!CHECK_RET(ret == ACL_SUCCESS)) { 
        LOG_PRINT("aclrtMalloc failed. ERROR: %d\n", ret); 
        return ret;
      }
      
      ret = aclrtMemcpy(*deviceAddr, size, hostData.data(), size, ACL_MEMCPY_HOST_TO_DEVICE);
      if (!CHECK_RET(ret == ACL_SUCCESS)) { 
        LOG_PRINT("aclrtMemcpy failed. ERROR: %d\n", ret); 
        return ret;
      }
     
      std::vector<int64_t> strides(shape.size(), 1);
      for (int64_t i = shape.size() - 2; i >= 0; i--) {
        strides[i] = shape[i + 1] * strides[i + 1];
      }
     
      *tensor = aclCreateTensor(shape.data(), shape.size(), dataType, strides.data(), 0, aclFormat::ACL_FORMAT_ND,
                                shape.data(), shape.size(), *deviceAddr);
      return 0;
    }
    
    struct TensorResources {
        void* queryDeviceAddr = nullptr;
        void* keyDeviceAddr = nullptr;
        void* valueDeviceAddr = nullptr;
        void* attenDeviceAddr = nullptr;
        void* outDeviceAddr = nullptr;
        aclTensor* queryTensor = nullptr;
        aclTensor* keyTensor = nullptr;
        aclTensor* valueTensor = nullptr;
        aclTensor* attenTensor = nullptr;
        aclTensor* outTensor = nullptr;
        aclIntArray* actualSeqLengths = nullptr;
    };
    
    int InitializeTensors(TensorResources& resources) {
        std::vector<int64_t> queryShape = {1, 2, 1, 16};
        std::vector<int64_t> keyShape = {1, 2, 2, 16};
        std::vector<int64_t> valueShape = {1, 2, 2, 16};
        std::vector<int64_t> attenShape = {1, 1, 1, 2};
        std::vector<int64_t> outShape = {1, 2, 1, 16};
        
        int64_t queryShapeSize = GetShapeSize(queryShape);
        int64_t keyShapeSize = GetShapeSize(keyShape);
        int64_t valueShapeSize = GetShapeSize(valueShape);
        int64_t attenShapeSize = GetShapeSize(attenShape);
        int64_t outShapeSize = GetShapeSize(outShape);
        
        std::vector<float> queryHostData(queryShapeSize, 1);
        std::vector<float> keyHostData(keyShapeSize, 1);
        std::vector<float> valueHostData(valueShapeSize, 1);
        std::vector<float> attenHostData(attenShapeSize, 1);
        std::vector<float> outHostData(outShapeSize, 1);
    
        int ret = CreateAclTensor(queryHostData, queryShape, &resources.queryDeviceAddr, 
                                 aclDataType::ACL_FLOAT16, &resources.queryTensor);
        if (!CHECK_RET(ret == ACL_SUCCESS)) {
          return ret;
        }
    
        ret = CreateAclTensor(keyHostData, keyShape, &resources.keyDeviceAddr, 
                             aclDataType::ACL_FLOAT16, &resources.keyTensor);
        if (!CHECK_RET(ret == ACL_SUCCESS)) {
          return ret;
        }
    
        ret = CreateAclTensor(valueHostData, valueShape, &resources.valueDeviceAddr, 
                             aclDataType::ACL_FLOAT16, &resources.valueTensor);
        if (!CHECK_RET(ret == ACL_SUCCESS)) {
          return ret;
        }
    
        ret = CreateAclTensor(attenHostData, attenShape, &resources.attenDeviceAddr, 
                             aclDataType::ACL_BOOL, &resources.attenTensor);
        if (!CHECK_RET(ret == ACL_SUCCESS)) {
          return ret;
        }
    
        ret = CreateAclTensor(outHostData, outShape, &resources.outDeviceAddr, 
                             aclDataType::ACL_FLOAT16, &resources.outTensor);
        if (!CHECK_RET(ret == ACL_SUCCESS)) {
          return ret;
        }
    
        std::vector<int64_t> actualSeqlenVector = {2};
        resources.actualSeqLengths = aclCreateIntArray(actualSeqlenVector.data(), 
                                                      actualSeqlenVector.size());
    
        return ACL_SUCCESS;
    }
    
    int ExecutePromptFlashAttention(TensorResources& resources, aclrtStream stream, 
                                   void** workspaceAddr, uint64_t* workspaceSize) {
        int64_t numHeads = 2;
        int64_t numKeyValueHeads = numHeads;
        double scaleValue = 1 / sqrt(2);
        int64_t preTokens = 65535;
        int64_t nextTokens = 65535;
        int64_t sparseMode = 0;
        int64_t innerPrecise = 1;
        
        constexpr const char LAYER_OUT_STR[] = "BNSD";
        constexpr size_t LAYER_OUT_LEN = sizeof(LAYER_OUT_STR);  
        char layerOut[LAYER_OUT_LEN];
        memcpy(layerOut, LAYER_OUT_STR, LAYER_OUT_LEN);
    
        aclOpExecutor* executor;
        int ret = aclnnPromptFlashAttentionV3GetWorkspaceSize(
            resources.queryTensor, resources.keyTensor, resources.valueTensor, 
            nullptr, nullptr, nullptr, nullptr, nullptr, nullptr, nullptr, nullptr, nullptr, 
            numHeads, scaleValue, preTokens, nextTokens, layerOut, numKeyValueHeads, 
            sparseMode, innerPrecise, resources.outTensor, workspaceSize, &executor);
            
        if (!CHECK_RET(ret == ACL_SUCCESS)) {
            LOG_PRINT("aclnnPromptFlashAttentionV3GetWorkspaceSize failed. ERROR: %d\n", ret);
            return ret;
        }
    
        if (*workspaceSize > 0ULL) {
            ret = aclrtMalloc(workspaceAddr, *workspaceSize, ACL_MEM_MALLOC_HUGE_FIRST);
            if (!CHECK_RET(ret == ACL_SUCCESS)) {
                LOG_PRINT("allocate workspace failed. ERROR: %d\n", ret);
                return ret;
            }
        }
    
        ret = aclnnPromptFlashAttentionV3(*workspaceAddr, *workspaceSize, executor, stream);
        if (!CHECK_RET(ret == ACL_SUCCESS)) {
            LOG_PRINT("aclnnPromptFlashAttentionV3 failed. ERROR: %d\n", ret);
            return ret;
        }
    
        return ACL_SUCCESS;
    }
    
    int ProcessResults(TensorResources& resources, const std::vector<int64_t>& outShape) {
        auto size = GetShapeSize(outShape);
        std::vector<op::fp16_t> resultData(size, 0);
        
        int ret = aclrtMemcpy(resultData.data(), resultData.size() * sizeof(resultData[0]), 
                             resources.outDeviceAddr, size * sizeof(resultData[0]), 
                             ACL_MEMCPY_DEVICE_TO_HOST);
        if (!CHECK_RET(ret == ACL_SUCCESS)) {
            LOG_PRINT("copy result from device to host failed. ERROR: %d\n", ret);
            return ret;
        }
        
        for (int64_t i = 0; i < size; i++) {
            std::cout << "index: " << i << ": " << static_cast<float>(resultData[i]) << std::endl;
        }
        
        return ACL_SUCCESS;
    }
    
    void CleanupResources(TensorResources& resources, void* workspaceAddr, 
                         aclrtStream stream, int32_t deviceId) {
        if (resources.queryTensor) {
          aclDestroyTensor(resources.queryTensor);
        }
        if (resources.keyTensor) {
          aclDestroyTensor(resources.keyTensor);
        }
        if (resources.valueTensor) {
          aclDestroyTensor(resources.valueTensor);
        }
        if (resources.attenTensor) {
          aclDestroyTensor(resources.attenTensor);
        }
        if (resources.outTensor) {
          aclDestroyTensor(resources.outTensor);
        }
        if (resources.actualSeqLengths) {
          aclDestroyIntArray(resources.actualSeqLengths);
        }
        
        if (resources.queryDeviceAddr) {
          aclrtFree(resources.queryDeviceAddr);
        }
        if (resources.keyDeviceAddr) {
          aclrtFree(resources.keyDeviceAddr);
        }
        if (resources.valueDeviceAddr) {
          aclrtFree(resources.valueDeviceAddr);
        }
        if (resources.attenDeviceAddr) {
          aclrtFree(resources.attenDeviceAddr);
        }
        if (resources.outDeviceAddr) {
          aclrtFree(resources.outDeviceAddr);
        }
        
        if (workspaceAddr) {
          aclrtFree(workspaceAddr);
        }
        if (stream) {
          aclrtDestroyStream(stream);
        }
        
        aclrtResetDevice(deviceId);
        aclFinalize();
    }
    
    } // namespace
     
    int main() {
        int32_t deviceId = 0;
        aclrtStream stream = nullptr;
        TensorResources resources = {};
        void* workspaceAddr = nullptr;
        uint64_t workspaceSize = 0;
        std::vector<int64_t> outShape = {1, 2, 1, 16};
        int ret = ACL_SUCCESS;
    
        // 1. Initialize device and stream
        ret = Init(deviceId, &stream);
        if (!CHECK_RET(ret == ACL_SUCCESS)) {
            LOG_PRINT("Init acl failed. ERROR: %d\n", ret);
            return ret;
        }
    
        // 2. Initialize tensors
        ret = InitializeTensors(resources);
        if (!CHECK_RET(ret == ACL_SUCCESS)) {
            CleanupResources(resources, workspaceAddr, stream, deviceId);
            return ret;
        }
    
        // 3. Execute the operation
        ret = ExecutePromptFlashAttention(resources, stream, &workspaceAddr, &workspaceSize);
        if (!CHECK_RET(ret == ACL_SUCCESS)) {
            CleanupResources(resources, workspaceAddr, stream, deviceId);
            return ret;
        }
    
        // 4. Synchronize stream
        ret = aclrtSynchronizeStream(stream);
        if (!CHECK_RET(ret == ACL_SUCCESS)) {
            LOG_PRINT("aclrtSynchronizeStream failed. ERROR: %d\n", ret);
            CleanupResources(resources, workspaceAddr, stream, deviceId);
            return ret;
        }
    
        // 5. Process results
        ret = ProcessResults(resources, outShape);
        if (!CHECK_RET(ret == ACL_SUCCESS)) {
            CleanupResources(resources, workspaceAddr, stream, deviceId);
            return ret;
        }
    
        // 6. Cleanup resources
        CleanupResources(resources, workspaceAddr, stream, deviceId);
        return 0;
    }