#include "cpu/aarch64/brgemm/brgemm.hpp"
#include "cpu/aarch64/brgemm/brgemm_types.hpp"
#include "cpu/aarch64/brgemm/brgemm_utils.hpp"
#include "common/c_types_map.hpp"
#include "common/nstl.hpp"
#include "common/type_helpers.hpp"
#include "common/utils.hpp"
#include "cpu/aarch64/cpu_barrier.hpp"
#include "cpu/aarch64/injectors/jit_uni_postops_injector.hpp"
#include "cpu/platform.hpp"
namespace dnnl {
namespace impl {
namespace cpu {
namespace aarch64 {
using namespace dnnl::impl::status;
using namespace dnnl::impl::utils;
using namespace prop_kind;
using namespace data_type;
using namespace brgemm_utils;
void brgemm_kernel_execute(const brgemm_kernel_t *brg_kernel, int bs,
const brgemm_batch_element_t *batch, void *ptr_C, void *scratch) {
brgemm_kernel_params_t brgemm_p;
brgemm_p.batch = batch;
brgemm_p.ptr_A = nullptr;
brgemm_p.ptr_B = nullptr;
brgemm_p.ptr_C = ptr_C;
brgemm_p.ptr_D = ptr_C;
brgemm_p.ptr_buf = scratch;
brgemm_p.ptr_bias = nullptr;
brgemm_p.do_post_ops = 0;
brgemm_p.do_apply_comp = 0;
brgemm_p.skip_accm = 0;
brgemm_p.BS = bs;
assert(brg_kernel);
(*brg_kernel)(&brgemm_p);
}
void brgemm_kernel_execute(const brgemm_kernel_t *brg_kernel, int bs,
const void *addr_A, const void *addr_B,
const brgemm_batch_element_t *batch, void *ptr_C, void *scratch) {
brgemm_kernel_params_t brgemm_p;
brgemm_p.batch = batch;
brgemm_p.ptr_A = addr_A;
brgemm_p.ptr_B = addr_B;
brgemm_p.ptr_C = ptr_C;
brgemm_p.ptr_D = ptr_C;
brgemm_p.ptr_buf = scratch;
brgemm_p.ptr_bias = nullptr;
brgemm_p.do_post_ops = 0;
brgemm_p.do_apply_comp = 0;
brgemm_p.skip_accm = 0;
brgemm_p.BS = bs;
assert(brg_kernel);
(*brg_kernel)(&brgemm_p);
}
void brgemm_kernel_execute_postops(const brgemm_kernel_t *brg_kernel, int bs,
const brgemm_batch_element_t *batch, void *ptr_C, void *ptr_D,
const brgemm_post_ops_data_t &post_ops_data, void *scratch) {
brgemm_kernel_params_t brgemm_p;
brgemm_p.batch = batch;
brgemm_p.ptr_A = nullptr;
brgemm_p.ptr_B = nullptr;
brgemm_p.ptr_C = ptr_C;
brgemm_p.ptr_D = ptr_D;
brgemm_p.ptr_buf = scratch;
brgemm_p.ptr_bias = post_ops_data.bias;
brgemm_p.ptr_scales = post_ops_data.scales;
brgemm_p.do_post_ops
= post_ops_data.do_only_comp || post_ops_data.do_only_zp_a_val ? 0
: 1;
brgemm_p.do_apply_comp = post_ops_data.do_only_zp_a_val ? 0 : 1;
brgemm_p.skip_accm = post_ops_data.skip_accumulation ? 1 : 0;
brgemm_p.BS = bs;
brgemm_p.zp_a_val = post_ops_data.zp_a_val;
brgemm_p.post_ops_binary_rhs_arg_vec = post_ops_data.binary_post_ops_rhs;
brgemm_p.oc_logical_off = post_ops_data.oc_logical_off;
brgemm_p.dst_row_logical_off = post_ops_data.dst_row_logical_off;
brgemm_p.data_C_ptr_ = post_ops_data.data_C_ptr_;
brgemm_p.first_mb_matrix_addr_off = post_ops_data.first_mb_matrix_addr_off;
brgemm_p.a_zp_compensations = post_ops_data.a_zp_compensations;
brgemm_p.b_zp_compensations = post_ops_data.b_zp_compensations;
brgemm_p.c_zp_values = post_ops_data.c_zp_values;
brgemm_p.ptr_dst_scales = post_ops_data.dst_scales;
assert(brg_kernel);
(*brg_kernel)(&brgemm_p);
}
void brgemm_kernel_execute_postops(const brgemm_kernel_t *brg_kernel, int bs,
const void *addr_A, const void *addr_B,
const brgemm_batch_element_t *batch, void *ptr_C, void *ptr_D,
const brgemm_post_ops_data_t &post_ops_data, void *scratch) {
brgemm_kernel_params_t brgemm_p;
brgemm_p.batch = batch;
brgemm_p.ptr_A = addr_A;
brgemm_p.ptr_B = addr_B;
brgemm_p.ptr_C = ptr_C;
brgemm_p.ptr_D = ptr_D;
brgemm_p.ptr_buf = scratch;
brgemm_p.ptr_bias = post_ops_data.bias;
brgemm_p.ptr_scales = post_ops_data.scales;
brgemm_p.do_post_ops
= post_ops_data.do_only_comp || post_ops_data.do_only_zp_a_val ? 0
: 1;
brgemm_p.do_apply_comp = post_ops_data.do_only_zp_a_val ? 0 : 1;
brgemm_p.skip_accm = post_ops_data.skip_accumulation ? 1 : 0;
brgemm_p.BS = bs;
brgemm_p.zp_a_val = post_ops_data.zp_a_val;
brgemm_p.post_ops_binary_rhs_arg_vec = post_ops_data.binary_post_ops_rhs;
brgemm_p.oc_logical_off = post_ops_data.oc_logical_off;
brgemm_p.dst_row_logical_off = post_ops_data.dst_row_logical_off;
brgemm_p.data_C_ptr_ = post_ops_data.data_C_ptr_;
brgemm_p.first_mb_matrix_addr_off = post_ops_data.first_mb_matrix_addr_off;
brgemm_p.a_zp_compensations = post_ops_data.a_zp_compensations;
brgemm_p.b_zp_compensations = post_ops_data.b_zp_compensations;
brgemm_p.c_zp_values = post_ops_data.c_zp_values;
brgemm_p.ptr_dst_scales = post_ops_data.dst_scales;
assert(brg_kernel);
(*brg_kernel)(&brgemm_p);
}
status_t brgemm_desc_init(brgemm_desc_t *brg, cpu_isa_t isa,
brgemm_batch_kind_t type, impl::data_type_t dt_a,
impl::data_type_t dt_b, bool transA, bool transB,
brgemm_layout_t layout, float alpha, float beta, dim_t LDA, dim_t LDB,
dim_t LDC, dim_t M, dim_t N, dim_t K, const brgemm_strides_t *strides) {
if (brg == nullptr) return status::invalid_arguments;
if (transA || transB) return status::unimplemented;
CHECK(brgemm_utils::init_brgemm_conf(brg, isa, type, dt_a, dt_b, layout,
alpha, beta, LDA, LDB, LDC, M, N, K, strides));
if (M <= 0 || N <= 0 || K <= 0) return status::invalid_arguments;
bool ldx_check = (brg->is_gemv || brg->is_row_major())
? (LDA < K)
: (LDA < M || LDB < K || LDC < M);
if (ldx_check) return status::invalid_arguments;
if (utils::everyone_is(
false, brg->is_int8, brg->is_bf16, brg->is_f32, brg->is_f16))
return status::unimplemented;
CHECK(brgemm_blocking(brg));
return status::success;
}
status_t brdgmm_desc_init(brgemm_desc_t *brg, cpu_isa_t isa,
brgemm_batch_kind_t type, impl::data_type_t dt_a,
impl::data_type_t dt_b, bool transA, brgemm_layout_t layout,
float alpha, float beta, dim_t LDA, dim_t LDC, dim_t M, dim_t N,
const brgemm_strides_t *strides) {
if (brg == nullptr) return status::invalid_arguments;
if (transA || layout != brgemm_row_major || alpha != 1.0f || beta != 0.f)
return status::unimplemented;
CHECK(brgemm_utils::init_brdgmm_conf(brg, isa, type, dt_a, dt_b, layout,
alpha, beta, LDA, LDC, M, N, strides));
const bool ldx_check = (LDA < N || LDC < N);
if (ldx_check) return status::invalid_arguments;
if (utils::everyone_is(
false, brg->is_int8, brg->is_bf16, brg->is_f32, brg->is_f16))
return status::unimplemented;
CHECK(brdgmm_blocking(brg));
return status::success;
}
status_t brgemm_desc_set_postops(brgemm_desc_t *brg,
const primitive_attr_t *attr, const memory_desc_t *dst_md, int LDD,
impl::data_type_t dt_bias) {
if (!brg || !dst_md) return status::invalid_arguments;
brg->attr = attr;
brg->dst_md = dst_md;
brg->with_bias = (dt_bias == data_type::undef) ? false : true;
brg->dt_bias = dt_bias;
brg->typesize_bias = (dt_bias == data_type::undef)
? 0
: types::data_type_size(brg->dt_bias);
brg->LDD = LDD;
const auto dt_d = dst_md->data_type;
if (brg->is_int8) {
if ((brg->dt_a == data_type::s8 && brg->dt_b == data_type::u8)
|| (!one_of(dt_bias, data_type::undef, data_type::s8,
data_type::u8, data_type::f32, data_type::s32))
|| (!one_of(dt_d, data_type::s8, data_type::u8, data_type::s32,
data_type::f32, data_type::bf16)))
return status::unimplemented;
}
if ((brg->dt_a == data_type::bf16 && brg->dt_b == data_type::bf16)
&& ((!one_of(dt_d, data_type::bf16, data_type::f32))
|| (!one_of(dt_bias, data_type::undef, data_type::bf16,
data_type::f32, data_type::s32))))
return status::unimplemented;
if ((brg->dt_a == data_type::f32 && brg->dt_b == data_type::f32)
&& (!one_of(dt_d, data_type::f32))
&& (!one_of(dt_bias, data_type::undef, data_type::s8, data_type::u8,
data_type::f32, data_type::s32)))
return status::unimplemented;
brg->dt_d = dt_d;
brg->typesize_D = types::data_type_size(brg->dt_d);
if (!brg->attr) return status::success;
using namespace injector;
const auto &post_ops = brg->attr->post_ops_;
const memory_desc_wrapper dst_d(dst_md);
const auto binary_ind = post_ops.find(primitive_kind::binary);
const auto prelu_ind = post_ops.find(primitive_kind::prelu);
brg->with_binary = !everyone_is(-1, binary_ind, prelu_ind);
const cpu_isa_t isa = get_max_cpu_isa();
if ((brg->with_binary && !dst_md)
|| !injector::post_ops_ok(
post_ops_ok_args_t(isa, {sum, eltwise, binary}, post_ops,
&dst_d, false ,
false ,
false ,
true ,
{broadcasting_strategy_t::per_oc,
broadcasting_strategy_t::scalar,
broadcasting_strategy_t::per_mb_spatial,
broadcasting_strategy_t::per_mb_w,
broadcasting_strategy_t::per_w,
broadcasting_strategy_t::no_broadcast})))
return status::unimplemented;
const auto sum_idx = post_ops.find(primitive_kind::sum);
const bool with_sum = sum_idx != -1;
brg->with_sum = with_sum;
brg->sum_scale = with_sum ? post_ops.entry_[sum_idx].sum.scale : 0;
brg->sum_zp = with_sum ? post_ops.entry_[sum_idx].sum.zero_point : 0;
const auto sum_dt
= with_sum ? post_ops.entry_[sum_idx].sum.dt : data_type::undef;
brg->sum_dt = sum_dt != data_type::undef ? sum_dt : dt_d;
const auto eltwise_ind = post_ops.find(primitive_kind::eltwise);
brg->with_eltwise = eltwise_ind != -1;
const auto &src_scales = attr->scales_.get(DNNL_ARG_SRC);
const auto &wei_scales = attr->scales_.get(DNNL_ARG_WEIGHTS);
const bool has_src_scales = !src_scales.has_default_values();
const bool has_wei_scales = !wei_scales.has_default_values();
brg->with_scales = has_src_scales || has_wei_scales
|| brg->with_weights_scale_adjust;
if (brg->with_scales) {
brg->is_oc_scale = wei_scales.get_mask() > 0;
}
const auto &dst_scales = attr->scales_.get(DNNL_ARG_DST);
const bool has_dst_scales = !dst_scales.has_default_values();
brg->with_dst_scales = has_dst_scales;
const bool scales_ok
= IMPLICATION(has_src_scales, src_scales.get_mask() == 0)
&& IMPLICATION(has_dst_scales, dst_scales.get_mask() == 0)
&& attr->scales_.has_default_values(
{DNNL_ARG_SRC, DNNL_ARG_WEIGHTS, DNNL_ARG_DST});
if (!scales_ok) return status::unimplemented;
auto init_zp_type
= [&](brgemm_broadcast_t &zp_type, int mem_arg) -> status_t {
const auto &zp = attr->zero_points_;
zp_type = brgemm_broadcast_t::none;
if (!zp.has_default_values(mem_arg)) {
int mask = zp.get_mask(mem_arg);
if (mask == 0) {
zp_type = brgemm_broadcast_t::per_tensor;
} else if (mask == (1 << 1)) {
zp_type = brgemm_broadcast_t::per_n;
} else {
return status::unimplemented;
}
}
return status::success;
};
init_zp_type(brg->zp_type_a, DNNL_ARG_SRC);
init_zp_type(brg->zp_type_b, DNNL_ARG_WEIGHTS);
init_zp_type(brg->zp_type_c, DNNL_ARG_DST);
const bool is_zp_src = brg->zp_type_a != brgemm_broadcast_t::none;
if (brg->is_dgmm) {
if (is_zp_src) CHECK(brdgmm_blocking(brg));
} else if (is_zp_src || brg->is_bf16_emu)
CHECK(brgemm_blocking(brg));
return status::success;
}
status_t brgemm_desc_set_attr(
brgemm_desc_t *brg, const brgemm_attr_t &brgattr) {
if (brg == nullptr) return status::invalid_arguments;
if (brgattr.max_top_vpad < 0 || brgattr.max_bottom_vpad < 0)
return status::unimplemented;
if (!brg->is_dgmm) {
if (brgattr.max_top_vpad > brgemm_desc_t::MAX_VPAD
|| brgattr.max_bottom_vpad > brgemm_desc_t::MAX_VPAD)
return status::unimplemented;
}
if ((brgattr.max_top_vpad > 0 || brgattr.max_bottom_vpad > 0)
&& brg->layout != brgemm_row_major)
return status::unimplemented;
brg->brgattr = brgattr;
if (brgattr.fpmath_mode != fpmath_mode::strict) maybe_try_bf32(brg);
const int max_vpad = nstl::max(brgattr.max_top_vpad,
brgattr.max_bottom_vpad); bool hint_blocking_set
= (brgattr.hint_bd_block != 0 || brgattr.hint_bd_block2 != 0
|| brgattr.hint_ld_block != 0 || brgattr.hint_ld_block2 != 0
|| brgattr.hint_load_nt_A != brgemm_hint_nt_undef
|| brgattr.hint_load_nt_B != brgemm_hint_nt_undef);
if (brgattr.use_uker || hint_blocking_set || brgattr.bd_mask_level
|| brgattr.fpmath_mode != fpmath_mode::strict || max_vpad > 0) {
if (brg->is_dgmm)
CHECK(brdgmm_blocking(brg));
else
CHECK(brgemm_blocking(brg));
}
if (!brg->is_dgmm) {
const int min_bd_block
= brg->bdb_tail > 0 ? brg->bdb_tail : brg->bd_block;
if ((max_vpad > min_bd_block)) return status::unimplemented;
}
brg->LDA2 = (brgattr.LDA2 != 0) ? brgattr.LDA2 : brg->LDA;
brg->LDB2 = (brgattr.LDB2 != 0) ? brgattr.LDB2 : brg->LDB;
brg->LDC2_M = (brgattr.LDC2_M != 0) ? brgattr.LDC2_M : brg->LDC;
brg->LDC2_N = (brgattr.LDC2_N != 0) ? brgattr.LDC2_N : brg->ld_block;
brg->is_blocked = (brg->LDA2 != brg->LDA || brg->LDB2 != brg->LDB
|| brg->LDC2_M != brg->LDC || brg->LDC2_N != brg->ld_block);
if (!IMPLICATION(brg->is_blocked, brg->layout = brgemm_row_major))
return status::invalid_arguments;
brg->prfA = brgattr.hint_prfA;
brg->prfB = brgattr.hint_prfB;
brg->prfC = brgattr.hint_prfC;
if (brgattr.hint_innermost_loop != brgemm_innermost_undef)
brg->innermost_loop = brgattr.hint_innermost_loop;
if (brgattr.hint_prefetching == brgemm_kernel_prefetching_t::brgemm_prf1
&& brg->prfC.dist1 < 0)
brg->prfC.dist1 = 0;
if (brgattr.hint_prefetching == brgemm_kernel_prefetching_t::brgemm_prf2
&& brg->prfC.dist2 < 0)
brg->prfC.dist2 = 0;
return status::success;
}
status_t brgemm_desc_finalize(brgemm_desc_t *brg) {
return status::success;
}
status_t brgemm_kernel_create(
brgemm_kernel_t **brg_kernel, const brgemm_desc_t &brg) {
if (!brg_kernel) return status::invalid_arguments;
*brg_kernel = nullptr;
if (brg.is_dgmm) {
CHECK(safe_ptr_assign<brgemm_kernel_t>(
*brg_kernel, new brdgmm_kernel_t(brg)));
} else {
CHECK(safe_ptr_assign<brgemm_kernel_t>(
*brg_kernel, new brgemm_kernel_common_t(brg)));
}
if (!(*brg_kernel)) return status::unimplemented;
return (*brg_kernel)->create_kernel();
}
status_t brgemm_kernel_destroy(brgemm_kernel_t *brg_kernel) {
delete brg_kernel;
return status::success;
}
status_t brgemm_init_tiles(const brgemm_desc_t &brg, char palette[64]) {
return status::unimplemented;
}
namespace {
template <typename T>
inline int sign(T v) {
return (v > 0) ? 1 : ((v < 0) ? -1 : 0);
}
int brgemm_cmp(const brgemm_desc_t &lhs, const brgemm_desc_t &rhs) {
#define CMP_BRGEMM_FIELD(x) \
if ((lhs.x) != (rhs.x)) return sign((lhs.x) - (rhs.x))
CMP_BRGEMM_FIELD(bcast_dim);
CMP_BRGEMM_FIELD(load_dim);
CMP_BRGEMM_FIELD(reduce_dim);
CMP_BRGEMM_FIELD(LDA);
CMP_BRGEMM_FIELD(LDB);
CMP_BRGEMM_FIELD(LDC);
CMP_BRGEMM_FIELD(LDD);
CMP_BRGEMM_FIELD(isa_user);
CMP_BRGEMM_FIELD(isa_impl);
CMP_BRGEMM_FIELD(alpha);
CMP_BRGEMM_FIELD(beta);
CMP_BRGEMM_FIELD(dt_a);
CMP_BRGEMM_FIELD(dt_b);
CMP_BRGEMM_FIELD(dt_c);
CMP_BRGEMM_FIELD(dt_d);
CMP_BRGEMM_FIELD(dt_bias);
CMP_BRGEMM_FIELD(stride_a);
CMP_BRGEMM_FIELD(stride_b);
CMP_BRGEMM_FIELD(layout);
CMP_BRGEMM_FIELD(type);
CMP_BRGEMM_FIELD(is_dgmm);
CMP_BRGEMM_FIELD(with_sum);
CMP_BRGEMM_FIELD(req_cal_comp_pads);
CMP_BRGEMM_FIELD(sum_scale);
CMP_BRGEMM_FIELD(sum_zp);
CMP_BRGEMM_FIELD(sum_dt);
CMP_BRGEMM_FIELD(with_eltwise);
CMP_BRGEMM_FIELD(with_binary);
CMP_BRGEMM_FIELD(with_scales);
CMP_BRGEMM_FIELD(zp_type_a);
CMP_BRGEMM_FIELD(zp_type_b);
CMP_BRGEMM_FIELD(zp_type_c);
CMP_BRGEMM_FIELD(is_oc_scale);
CMP_BRGEMM_FIELD(with_dst_scales);
CMP_BRGEMM_FIELD(brgattr.max_bs);
CMP_BRGEMM_FIELD(brgattr.max_top_vpad);
CMP_BRGEMM_FIELD(brgattr.max_bottom_vpad);
CMP_BRGEMM_FIELD(brgattr.hint_expected_A_size);
CMP_BRGEMM_FIELD(brgattr.hint_expected_B_size);
CMP_BRGEMM_FIELD(brgattr.hint_expected_C_size);
CMP_BRGEMM_FIELD(brgattr.hint_innermost_loop);
CMP_BRGEMM_FIELD(brgattr.hint_loop_order);
CMP_BRGEMM_FIELD(brgattr.hint_prefetching);
CMP_BRGEMM_FIELD(brgattr.hint_prfA.dist1);
CMP_BRGEMM_FIELD(brgattr.hint_prfA.dist2);
CMP_BRGEMM_FIELD(brgattr.hint_prfB.dist1);
CMP_BRGEMM_FIELD(brgattr.hint_prfB.dist2);
CMP_BRGEMM_FIELD(brgattr.hint_prfC.dist1);
CMP_BRGEMM_FIELD(brgattr.hint_prfC.dist2);
CMP_BRGEMM_FIELD(brgattr.wary_A_k_tail_read);
CMP_BRGEMM_FIELD(brgattr.extendable_k);
CMP_BRGEMM_FIELD(brgattr.generate_skip_accumulation);
CMP_BRGEMM_FIELD(brgattr.bd_mask_level);
CMP_BRGEMM_FIELD(brgattr.use_uker);
CMP_BRGEMM_FIELD(brgattr.use_interleave_stores);
CMP_BRGEMM_FIELD(brgattr.b_is_vnni);
CMP_BRGEMM_FIELD(brgattr.fpmath_mode);
CMP_BRGEMM_FIELD(brgattr.LDA2);
CMP_BRGEMM_FIELD(brgattr.LDB2);
CMP_BRGEMM_FIELD(brgattr.LDC2_M);
CMP_BRGEMM_FIELD(brgattr.LDC2_N);
CMP_BRGEMM_FIELD(brgattr.var_bs);
CMP_BRGEMM_FIELD(brgattr.postops_only);
CMP_BRGEMM_FIELD(brgattr.hint_bd_block);
CMP_BRGEMM_FIELD(brgattr.hint_ld_block);
CMP_BRGEMM_FIELD(brgattr.hint_bd_block2);
CMP_BRGEMM_FIELD(brgattr.hint_ld_block2);
CMP_BRGEMM_FIELD(brgattr.hint_ununroll_bd_loop);
CMP_BRGEMM_FIELD(brgattr.hint_load_nt_A);
CMP_BRGEMM_FIELD(brgattr.hint_load_nt_B);
CMP_BRGEMM_FIELD(brgattr.K_koef);
if (lhs.brgattr.bd_mask_level > 0)
for (int i = 0; i < lhs.bcast_dim; i++) {
CMP_BRGEMM_FIELD(brgattr.bd_mask[i]);
}
if (lhs.type == brgemm_static_offs)
for (int i = 0; i < lhs.brgattr.max_bs; i++) {
CMP_BRGEMM_FIELD(brgattr.static_offsets[i].offset.A);
CMP_BRGEMM_FIELD(brgattr.static_offsets[i].offset.B);
}
#undef CMP_BRGEMM_FIELD
return 0;
}
}
bool brgemm_desc_t::operator==(const brgemm_desc_t &rhs) const {
return (brgemm_cmp(*this, rhs) == 0);
}
bool brgemm_desc_t::operator<(const brgemm_desc_t &rhs) const {
return (brgemm_cmp(*this, rhs) < 0);
}
} } } }