From nobody Sat May 10 05:37:16 2025
Delivered-To: importer@patchew.org
Authentication-Results: mx.zohomail.com;
	dkim=pass;
	spf=pass (zohomail.com: domain of gnu.org designates 209.51.188.17 as
 permitted sender)
  smtp.mailfrom=qemu-devel-bounces+importer=patchew.org@nongnu.org;
	dmarc=pass(p=none dis=none)  header.from=linaro.org
ARC-Seal: i=1; a=rsa-sha256; t=1657551034; cv=none;
	d=zohomail.com; s=zohoarc;
	b=XqpAOG8xuzU8v1XkRbpmLHGJ3e5Yqs4CEaN27LJABbNDis5J37k9ROhTckAhZiHGTWsKJW2kpFkxGaw8BgHp/w1lUp/vGPIMHyvQe3e5A1GZ6v4CgEYh0KszlEY3NgHierdokchjQ/iT3gVMYaKKyobpJMAWHfoTDSzz1aLrlL0=
ARC-Message-Signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=zohomail.com;
 s=zohoarc;
	t=1657551034;
 h=Content-Transfer-Encoding:Date:From:In-Reply-To:List-Subscribe:List-Post:List-Id:List-Archive:List-Help:List-Unsubscribe:MIME-Version:Message-ID:References:Sender:Subject:To;
	bh=cR0PvnQKG0XfbYb9xus9rRJ6s/0aSnCTWRG0QGmXr4w=;
	b=gxjKO+WjXysGvcmivHX9BGdu2K2fNMu7/1OgMbEAN6+Rp3C4d+hNRh4BxxG15ctg1PUAQfO0pST3r3HtH31aNpd+PI1VnxnMNmo4LpRrE2jf0+boXIP9OR2h3M6gOgUaXh1YgQ1/CueeVghelw7ZQI5Xmx3UGki8XkkBwg/ikgI=
ARC-Authentication-Results: i=1; mx.zohomail.com;
	dkim=pass;
	spf=pass (zohomail.com: domain of gnu.org designates 209.51.188.17 as
 permitted sender)
  smtp.mailfrom=qemu-devel-bounces+importer=patchew.org@nongnu.org;
	dmarc=pass header.from=<peter.maydell@linaro.org> (p=none dis=none)
Return-Path: <qemu-devel-bounces+importer=patchew.org@nongnu.org>
Received: from lists.gnu.org (lists.gnu.org [209.51.188.17]) by
 mx.zohomail.com
	with SMTPS id 1657551034318366.3350817805863;
 Mon, 11 Jul 2022 07:50:34 -0700 (PDT)
Received: from localhost ([::1]:48316 helo=lists1p.gnu.org)
	by lists.gnu.org with esmtp (Exim 4.90_1)
	(envelope-from <qemu-devel-bounces+importer=patchew.org@nongnu.org>)
	id 1oAuk4-0006TV-Ur
	for importer@patchew.org; Mon, 11 Jul 2022 10:50:32 -0400
Received: from eggs.gnu.org ([2001:470:142:3::10]:42780)
 by lists.gnu.org with esmtps (TLS1.2:ECDHE_RSA_AES_256_GCM_SHA384:256)
 (Exim 4.90_1) (envelope-from <peter.maydell@linaro.org>)
 id 1oAtvW-0006Zq-GM
 for qemu-devel@nongnu.org; Mon, 11 Jul 2022 09:58:19 -0400
Received: from mail-wm1-x32e.google.com ([2a00:1450:4864:20::32e]:53118)
 by eggs.gnu.org with esmtps (TLS1.2:ECDHE_RSA_AES_128_GCM_SHA256:128)
 (Exim 4.90_1) (envelope-from <peter.maydell@linaro.org>)
 id 1oAtvT-0002rT-J9
 for qemu-devel@nongnu.org; Mon, 11 Jul 2022 09:58:17 -0400
Received: by mail-wm1-x32e.google.com with SMTP id o8so3086158wms.2
 for <qemu-devel@nongnu.org>; Mon, 11 Jul 2022 06:58:13 -0700 (PDT)
Received: from orth.archaic.org.uk (orth.archaic.org.uk. [2001:8b0:1d0::2])
 by smtp.gmail.com with ESMTPSA id
 q123-20020a1c4381000000b00397402ae674sm7595596wma.11.2022.07.11.06.58.11
 for <qemu-devel@nongnu.org>
 (version=TLS1_3 cipher=TLS_AES_256_GCM_SHA384 bits=256/256);
 Mon, 11 Jul 2022 06:58:12 -0700 (PDT)
DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=linaro.org; s=google;
 h=from:to:subject:date:message-id:in-reply-to:references:mime-version
 :content-transfer-encoding;
 bh=cR0PvnQKG0XfbYb9xus9rRJ6s/0aSnCTWRG0QGmXr4w=;
 b=iqQk8F2xsaV4i1VVN7W8K8HT+bnzR5juzxYOwX7mbjXOcFONGHwSXCK2YoXACpUV+o
 y61yBqFqOT9fv9aL6Tsr0ofRC/LFnZCTujyi8i75I4WEN07hMkL8taxB80H57Bw8sDwJ
 czEJLY59sbyeW6aeXWVGN5Osruu29bIGXs0SgAMQMBmUbEkZJ+4aHpnoPxVx5pUKo3zt
 7sfU9uz1aRbUHMlAveD6YJHqqGWJNGhtD4BDHc0DKNknIUpEty3F56F/byghLsgTB/mR
 MzCP/04LYmvX/N+UPHiB+0mVki8KmFXrPzN91cJM4LKFeFR+8od3JJOU0RZjLPOHrwiq
 2KiQ==
X-Google-DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed;
 d=1e100.net; s=20210112;
 h=x-gm-message-state:from:to:subject:date:message-id:in-reply-to
 :references:mime-version:content-transfer-encoding;
 bh=cR0PvnQKG0XfbYb9xus9rRJ6s/0aSnCTWRG0QGmXr4w=;
 b=2fFQ6M0lSh1dM1Nn4bI+TNFQ2M1VEtvxxjpfjhHnLM6hQWbtpcLiS1YL24Jqs1FbmD
 QXqNfORd1utri6xnLf7hynYWIoHNwR9JNhSLHT4TAzg9lU1au6G4HbKQFAMPB5NDe7rN
 svUH80SH+I6Y4NoJmOIORj8+D3Fgf5gGW91avZZ7NZFihDD0rBxyg+MpEkULAZE66UEk
 ui8o97rEqwlAb6NqWcCQKQEoaLtK34UqFUrMTJtuiIN5GqCUng79gzJL3ZPZMWjwWgWg
 nR3hikwhPLWpwoHsaxkUdHuA10+Via/wqGFWAx13MWkbDwGDwc/SSWXkB42dZ9swjltp
 9R4g==
X-Gm-Message-State: AJIora+wxFnjac754fSolzwv6Sc5AeQVIfA5p14FiSf66ritd6T1jOQd
 2q8FLqdu3kIjjqaj10mohU3NEg96B4dN0Q==
X-Google-Smtp-Source: 
 AGRyM1sczBnig7+K3WSbMT1shm40oSw5Mx6AnUamTfMDd4PLHq0WLczMnauukdi0dS7FvG8JQzF4Dg==
X-Received: by 2002:a05:600c:42d3:b0:3a2:e7b8:8430 with SMTP id
 j19-20020a05600c42d300b003a2e7b88430mr6371997wme.42.1657547892517;
 Mon, 11 Jul 2022 06:58:12 -0700 (PDT)
From: Peter Maydell <peter.maydell@linaro.org>
To: qemu-devel@nongnu.org
Subject: [PULL 24/45] target/arm: Implement FMOPA, FMOPS (non-widening)
Date: Mon, 11 Jul 2022 14:57:29 +0100
Message-Id: <20220711135750.765803-25-peter.maydell@linaro.org>
X-Mailer: git-send-email 2.25.1
In-Reply-To: <20220711135750.765803-1-peter.maydell@linaro.org>
References: <20220711135750.765803-1-peter.maydell@linaro.org>
MIME-Version: 1.0
Content-Transfer-Encoding: quoted-printable
Received-SPF: pass (zohomail.com: domain of gnu.org designates 209.51.188.17
 as permitted sender) client-ip=209.51.188.17;
 envelope-from=qemu-devel-bounces+importer=patchew.org@nongnu.org;
 helo=lists.gnu.org;
Received-SPF: pass client-ip=2a00:1450:4864:20::32e;
 envelope-from=peter.maydell@linaro.org; helo=mail-wm1-x32e.google.com
X-Spam_score_int: -20
X-Spam_score: -2.1
X-Spam_bar: --
X-Spam_report: (-2.1 / 5.0 requ) BAYES_00=-1.9, DKIM_SIGNED=0.1,
 DKIM_VALID=-0.1, DKIM_VALID_AU=-0.1, DKIM_VALID_EF=-0.1,
 RCVD_IN_DNSWL_NONE=-0.0001, SPF_HELO_NONE=0.001, SPF_PASS=-0.001,
 T_SCC_BODY_TEXT_LINE=-0.01 autolearn=ham autolearn_force=no
X-Spam_action: no action
X-BeenThere: qemu-devel@nongnu.org
X-Mailman-Version: 2.1.29
Precedence: list
List-Id: <qemu-devel.nongnu.org>
List-Unsubscribe: <https://lists.nongnu.org/mailman/options/qemu-devel>,
 <mailto:qemu-devel-request@nongnu.org?subject=unsubscribe>
List-Archive: <https://lists.nongnu.org/archive/html/qemu-devel>
List-Post: <mailto:qemu-devel@nongnu.org>
List-Help: <mailto:qemu-devel-request@nongnu.org?subject=help>
List-Subscribe: <https://lists.nongnu.org/mailman/listinfo/qemu-devel>,
 <mailto:qemu-devel-request@nongnu.org?subject=subscribe>
Errors-To: qemu-devel-bounces+importer=patchew.org@nongnu.org
Sender: "Qemu-devel" <qemu-devel-bounces+importer=patchew.org@nongnu.org>
X-ZohoMail-DKIM: pass (identity @linaro.org)
X-ZM-MESSAGEID: 1657551036044100001
Content-Type: text/plain; charset="utf-8"

From: Richard Henderson <richard.henderson@linaro.org>

Signed-off-by: Richard Henderson <richard.henderson@linaro.org>
Message-id: 20220708151540.18136-25-richard.henderson@linaro.org
Signed-off-by: Peter Maydell <peter.maydell@linaro.org>
Reviewed-by: Peter Maydell <peter.maydell@linaro.org>
---
 target/arm/helper-sme.h    |  5 +++
 target/arm/sme.decode      |  9 +++++
 target/arm/sme_helper.c    | 69 ++++++++++++++++++++++++++++++++++++++
 target/arm/translate-sme.c | 32 ++++++++++++++++++
 4 files changed, 115 insertions(+)

diff --git a/target/arm/helper-sme.h b/target/arm/helper-sme.h
index 753e9e624cc..f50d0fe1d62 100644
--- a/target/arm/helper-sme.h
+++ b/target/arm/helper-sme.h
@@ -120,3 +120,8 @@ DEF_HELPER_FLAGS_5(sme_addha_s, TCG_CALL_NO_RWG, void, =
ptr, ptr, ptr, ptr, i32)
 DEF_HELPER_FLAGS_5(sme_addva_s, TCG_CALL_NO_RWG, void, ptr, ptr, ptr, ptr,=
 i32)
 DEF_HELPER_FLAGS_5(sme_addha_d, TCG_CALL_NO_RWG, void, ptr, ptr, ptr, ptr,=
 i32)
 DEF_HELPER_FLAGS_5(sme_addva_d, TCG_CALL_NO_RWG, void, ptr, ptr, ptr, ptr,=
 i32)
+
+DEF_HELPER_FLAGS_7(sme_fmopa_s, TCG_CALL_NO_RWG,
+                   void, ptr, ptr, ptr, ptr, ptr, ptr, i32)
+DEF_HELPER_FLAGS_7(sme_fmopa_d, TCG_CALL_NO_RWG,
+                   void, ptr, ptr, ptr, ptr, ptr, ptr, i32)
diff --git a/target/arm/sme.decode b/target/arm/sme.decode
index 8cb6c4053c3..ba4774d1746 100644
--- a/target/arm/sme.decode
+++ b/target/arm/sme.decode
@@ -64,3 +64,12 @@ ADDHA_s         11000000 10 01000 0 ... ... ..... 000 ..=
        @adda_32
 ADDVA_s         11000000 10 01000 1 ... ... ..... 000 ..        @adda_32
 ADDHA_d         11000000 11 01000 0 ... ... ..... 00 ...        @adda_64
 ADDVA_d         11000000 11 01000 1 ... ... ..... 00 ...        @adda_64
+
+### SME Outer Product
+
+&op             zad zn zm pm pn sub:bool
+@op_32          ........ ... zm:5 pm:3 pn:3 zn:5 sub:1 .. zad:2 &op
+@op_64          ........ ... zm:5 pm:3 pn:3 zn:5 sub:1 .  zad:3 &op
+
+FMOPA_s         10000000 100 ..... ... ... ..... . 00 ..        @op_32
+FMOPA_d         10000000 110 ..... ... ... ..... . 0 ...        @op_64
diff --git a/target/arm/sme_helper.c b/target/arm/sme_helper.c
index f1e924db748..7dc76b6a1c3 100644
--- a/target/arm/sme_helper.c
+++ b/target/arm/sme_helper.c
@@ -25,6 +25,7 @@
 #include "exec/cpu_ldst.h"
 #include "exec/exec-all.h"
 #include "qemu/int128.h"
+#include "fpu/softfloat.h"
 #include "vec_internal.h"
 #include "sve_ldst_internal.h"
=20
@@ -918,3 +919,71 @@ void HELPER(sme_addva_d)(void *vzda, void *vzn, void *=
vpn,
         }
     }
 }
+
+void HELPER(sme_fmopa_s)(void *vza, void *vzn, void *vzm, void *vpn,
+                         void *vpm, void *vst, uint32_t desc)
+{
+    intptr_t row, col, oprsz =3D simd_maxsz(desc);
+    uint32_t neg =3D simd_data(desc) << 31;
+    uint16_t *pn =3D vpn, *pm =3D vpm;
+    float_status fpst;
+
+    /*
+     * Make a copy of float_status because this operation does not
+     * update the cumulative fp exception status.  It also produces
+     * default nans.
+     */
+    fpst =3D *(float_status *)vst;
+    set_default_nan_mode(true, &fpst);
+
+    for (row =3D 0; row < oprsz; ) {
+        uint16_t pa =3D pn[H2(row >> 4)];
+        do {
+            if (pa & 1) {
+                void *vza_row =3D vza + tile_vslice_offset(row);
+                uint32_t n =3D *(uint32_t *)(vzn + H1_4(row)) ^ neg;
+
+                for (col =3D 0; col < oprsz; ) {
+                    uint16_t pb =3D pm[H2(col >> 4)];
+                    do {
+                        if (pb & 1) {
+                            uint32_t *a =3D vza_row + H1_4(col);
+                            uint32_t *m =3D vzm + H1_4(col);
+                            *a =3D float32_muladd(n, *m, *a, 0, vst);
+                        }
+                        col +=3D 4;
+                        pb >>=3D 4;
+                    } while (col & 15);
+                }
+            }
+            row +=3D 4;
+            pa >>=3D 4;
+        } while (row & 15);
+    }
+}
+
+void HELPER(sme_fmopa_d)(void *vza, void *vzn, void *vzm, void *vpn,
+                         void *vpm, void *vst, uint32_t desc)
+{
+    intptr_t row, col, oprsz =3D simd_oprsz(desc) / 8;
+    uint64_t neg =3D (uint64_t)simd_data(desc) << 63;
+    uint64_t *za =3D vza, *zn =3D vzn, *zm =3D vzm;
+    uint8_t *pn =3D vpn, *pm =3D vpm;
+    float_status fpst =3D *(float_status *)vst;
+
+    set_default_nan_mode(true, &fpst);
+
+    for (row =3D 0; row < oprsz; ++row) {
+        if (pn[H1(row)] & 1) {
+            uint64_t *za_row =3D &za[tile_vslice_index(row)];
+            uint64_t n =3D zn[row] ^ neg;
+
+            for (col =3D 0; col < oprsz; ++col) {
+                if (pm[H1(col)] & 1) {
+                    uint64_t *a =3D &za_row[col];
+                    *a =3D float64_muladd(n, zm[col], *a, 0, &fpst);
+                }
+            }
+        }
+    }
+}
diff --git a/target/arm/translate-sme.c b/target/arm/translate-sme.c
index d3b9cdd5c4b..fa8f343a7d6 100644
--- a/target/arm/translate-sme.c
+++ b/target/arm/translate-sme.c
@@ -298,3 +298,35 @@ TRANS_FEAT(ADDHA_s, aa64_sme, do_adda, a, MO_32, gen_h=
elper_sme_addha_s)
 TRANS_FEAT(ADDVA_s, aa64_sme, do_adda, a, MO_32, gen_helper_sme_addva_s)
 TRANS_FEAT(ADDHA_d, aa64_sme_i16i64, do_adda, a, MO_64, gen_helper_sme_add=
ha_d)
 TRANS_FEAT(ADDVA_d, aa64_sme_i16i64, do_adda, a, MO_64, gen_helper_sme_add=
va_d)
+
+static bool do_outprod_fpst(DisasContext *s, arg_op *a, MemOp esz,
+                            gen_helper_gvec_5_ptr *fn)
+{
+    int svl =3D streaming_vec_reg_size(s);
+    uint32_t desc =3D simd_desc(svl, svl, a->sub);
+    TCGv_ptr za, zn, zm, pn, pm, fpst;
+
+    if (!sme_smza_enabled_check(s)) {
+        return true;
+    }
+
+    /* Sum XZR+zad to find ZAd. */
+    za =3D get_tile_rowcol(s, esz, 31, a->zad, false);
+    zn =3D vec_full_reg_ptr(s, a->zn);
+    zm =3D vec_full_reg_ptr(s, a->zm);
+    pn =3D pred_full_reg_ptr(s, a->pn);
+    pm =3D pred_full_reg_ptr(s, a->pm);
+    fpst =3D fpstatus_ptr(FPST_FPCR);
+
+    fn(za, zn, zm, pn, pm, fpst, tcg_constant_i32(desc));
+
+    tcg_temp_free_ptr(za);
+    tcg_temp_free_ptr(zn);
+    tcg_temp_free_ptr(pn);
+    tcg_temp_free_ptr(pm);
+    tcg_temp_free_ptr(fpst);
+    return true;
+}
+
+TRANS_FEAT(FMOPA_s, aa64_sme, do_outprod_fpst, a, MO_32, gen_helper_sme_fm=
opa_s)
+TRANS_FEAT(FMOPA_d, aa64_sme_f64f64, do_outprod_fpst, a, MO_64, gen_helper=
_sme_fmopa_d)
--=20
2.25.1