From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: from mail-ot1-x32a.google.com (mail-ot1-x32a.google.com [IPv6:2607:f8b0:4864:20::32a]) by sourceware.org (Postfix) with ESMTPS id 429653858423 for ; Mon, 13 Mar 2023 17:08:20 +0000 (GMT) DMARC-Filter: OpenDMARC Filter v1.4.2 sourceware.org 429653858423 Authentication-Results: sourceware.org; dmarc=pass (p=none dis=none) header.from=linaro.org Authentication-Results: sourceware.org; spf=pass smtp.mailfrom=linaro.org Received: by mail-ot1-x32a.google.com with SMTP id l16-20020a9d4c10000000b006944b17058cso7050319otf.2 for ; Mon, 13 Mar 2023 10:08:20 -0700 (PDT) DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=linaro.org; s=google; t=1678727299; h=content-transfer-encoding:in-reply-to:organization:from:references :to:content-language:subject:user-agent:mime-version:date:message-id :from:to:cc:subject:date:message-id:reply-to; bh=uzfuKqEmQT16iU26xqm5kECj1+OXRgXI+3Jq9kFXjOI=; b=nisCsxWYULyNOVp/nuwqQKSiSKHpeQjZ4cejliq/xWcYNFSYlmpbDERLAtOvHC3TS4 lDVcF6wyRAe2rUi3Bs+/y7cQnfPCEvTCT6xJSC1ccwx7e+099l9Z/F/ykbmvcRivdrp/ B564+sB2ZssbpIhHojOZ1CF9ySsv06+A/4Ybcu12iW36n0JvMDnWZXEVPZRkdd3iW9tu X0IeCaxx3FFxNr01wpi89Z4IwiziP0Sa6GhcTE3WqyooLCeLYJ5gVQnS1oUWsTkHJEos ayZwg6b+PiQ5awWalqu8IIwcepwoC/+InStvN54TQrhL/rvujYTiPsbTU6IyR42SasBD QyTw== X-Google-DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=1e100.net; s=20210112; t=1678727299; h=content-transfer-encoding:in-reply-to:organization:from:references :to:content-language:subject:user-agent:mime-version:date:message-id :x-gm-message-state:from:to:cc:subject:date:message-id:reply-to; bh=uzfuKqEmQT16iU26xqm5kECj1+OXRgXI+3Jq9kFXjOI=; b=L6bocE6B6AolVKrKbxZ7cbHj6rIsW0uF0sVHPvt5I6x2VqNtpRLySFBOfp8w3R0fU2 +u7F8J+J7R546ziLpIF3BhhsQXYRzZvF3HCVWDS19TOutU+bFMtdlnKabpiR56ztiJJI Se8U2m6/KdS5h5b5bVK+9NDaCLD4L82RkQr4fm6zvGbJ4Kq1rfYl3A33i5xwb1X6yL9e aVQEdCqO+/CXF3YdU/oh3dZLy01dX25zUbL+6klXSZwI+hbDAnjR2rigSp2Ub0qXCJh7 Q/lmOKUuJ471xYtezJ8oEauxyuLfdBwikGDwtOtU2ZuJaII3jhloPODM5HNTNTA8bOPG 8Acg== X-Gm-Message-State: AO0yUKVd2J1arltQkcH3P4tftxdwHA28qjPHjNUM0hc7spSOOnIv3lHn /UNmR7Z5cupFDsxURMnIOyb7GO4LNAwLklaGAMn4ng== X-Google-Smtp-Source: AK7set/Pygbfxkhnb9Qgb0XvQJzVYxihXXcgpwAsURbYl6jX3h+K+pXYffAmGGv/1Z1UfBrC3QmCng== X-Received: by 2002:a05:6830:411b:b0:68b:cb6a:d1c5 with SMTP id w27-20020a056830411b00b0068bcb6ad1c5mr22792585ott.36.1678727299009; Mon, 13 Mar 2023 10:08:19 -0700 (PDT) Received: from ?IPV6:2804:1b3:a7c0:544b:c5ae:2fda:455d:6e31? ([2804:1b3:a7c0:544b:c5ae:2fda:455d:6e31]) by smtp.gmail.com with ESMTPSA id i17-20020a9d6251000000b0069457b86060sm176549otk.47.2023.03.13.10.08.17 (version=TLS1_3 cipher=TLS_AES_128_GCM_SHA256 bits=128/128); Mon, 13 Mar 2023 10:08:18 -0700 (PDT) Message-ID: Date: Mon, 13 Mar 2023 14:08:15 -0300 MIME-Version: 1.0 User-Agent: Mozilla/5.0 (Macintosh; Intel Mac OS X 10.15; rv:102.0) Gecko/20100101 Thunderbird/102.8.0 Subject: Re: [PATCH] Soft-fp: Add brain format support Content-Language: en-US To: Jakub Jelinek , libc-alpha@sourceware.org References: From: Adhemerval Zanella Netto Organization: Linaro In-Reply-To: Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 7bit X-Spam-Status: No, score=-12.5 required=5.0 tests=BAYES_00,DKIM_SIGNED,DKIM_VALID,DKIM_VALID_AU,DKIM_VALID_EF,GIT_PATCH_0,KAM_SHORT,NICE_REPLY_A,RCVD_IN_DNSWL_NONE,SPF_HELO_NONE,SPF_PASS,TXREP autolearn=ham autolearn_force=no version=3.4.6 X-Spam-Checker-Version: SpamAssassin 3.4.6 (2021-04-09) on server2.sourceware.org List-Id: On 09/03/23 14:54, Jakub Jelinek via Libc-alpha wrote: > In https://gcc.gnu.org/r13-3292 I've added brain format support > (std::bfloat16_t) on the GCC side, but as glibc has the master copy > of soft-fp, the following patch adds the files from that commit and > also from > https://gcc.gnu.org/pipermail/gcc-patches/2022-November/606398.html > to glibc. > LGTM, thanks. Reviewed-by: Adhemerval Zanella > diff --git a/soft-fp/brain.h b/soft-fp/brain.h > new file mode 100644 > index 0000000000..15b63a9a34 > --- /dev/null > +++ b/soft-fp/brain.h > @@ -0,0 +1,172 @@ > +/* Software floating-point emulation. > + Definitions for Brain Floating Point format (bfloat16). > + Copyright (C) 1997-2023 Free Software Foundation, Inc. > + This file is part of the GNU C Library. > + > + The GNU C Library is free software; you can redistribute it and/or > + modify it under the terms of the GNU Lesser General Public > + License as published by the Free Software Foundation; either > + version 2.1 of the License, or (at your option) any later version. > + > + In addition to the permissions in the GNU Lesser General Public > + License, the Free Software Foundation gives you unlimited > + permission to link the compiled version of this file into > + combinations with other programs, and to distribute those > + combinations without any restriction coming from the use of this > + file. (The Lesser General Public License restrictions do apply in > + other respects; for example, they cover modification of the file, > + and distribution when not linked into a combine executable.) > + > + The GNU C Library is distributed in the hope that it will be useful, > + but WITHOUT ANY WARRANTY; without even the implied warranty of > + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU > + Lesser General Public License for more details. > + > + You should have received a copy of the GNU Lesser General Public > + License along with the GNU C Library; if not, see > + . */ > + > +#ifndef SOFT_FP_BRAIN_H > +#define SOFT_FP_BRAIN_H 1 > + > +#if _FP_W_TYPE_SIZE < 32 > +# error "Here's a nickel kid. Go buy yourself a real computer." > +#endif > + > +#define _FP_FRACTBITS_B (_FP_W_TYPE_SIZE) > + > +#define _FP_FRACTBITS_DW_B (_FP_W_TYPE_SIZE) > + > +#define _FP_FRACBITS_B 8 > +#define _FP_FRACXBITS_B (_FP_FRACTBITS_B - _FP_FRACBITS_B) > +#define _FP_WFRACBITS_B (_FP_WORKBITS + _FP_FRACBITS_B) > +#define _FP_WFRACXBITS_B (_FP_FRACTBITS_B - _FP_WFRACBITS_B) > +#define _FP_EXPBITS_B 8 > +#define _FP_EXPBIAS_B 127 > +#define _FP_EXPMAX_B 255 > + > +#define _FP_QNANBIT_B ((_FP_W_TYPE) 1 << (_FP_FRACBITS_B-2)) > +#define _FP_QNANBIT_SH_B ((_FP_W_TYPE) 1 << (_FP_FRACBITS_B-2+_FP_WORKBITS)) > +#define _FP_IMPLBIT_B ((_FP_W_TYPE) 1 << (_FP_FRACBITS_B-1)) > +#define _FP_IMPLBIT_SH_B ((_FP_W_TYPE) 1 << (_FP_FRACBITS_B-1+_FP_WORKBITS)) > +#define _FP_OVERFLOW_B ((_FP_W_TYPE) 1 << (_FP_WFRACBITS_B)) > + > +#define _FP_WFRACBITS_DW_B (2 * _FP_WFRACBITS_B) > +#define _FP_WFRACXBITS_DW_B (_FP_FRACTBITS_DW_B - _FP_WFRACBITS_DW_B) > +#define _FP_HIGHBIT_DW_B \ > + ((_FP_W_TYPE) 1 << (_FP_WFRACBITS_DW_B - 1) % _FP_W_TYPE_SIZE) > + > +/* The implementation of _FP_MUL_MEAT_B and _FP_DIV_MEAT_B should be > + chosen by the target machine. */ > + > +typedef float BFtype __attribute__ ((mode (BF))); > + > +union _FP_UNION_B > +{ > + BFtype flt; > + struct _FP_STRUCT_LAYOUT > + { > +#if __BYTE_ORDER == __BIG_ENDIAN > + unsigned sign : 1; > + unsigned exp : _FP_EXPBITS_B; > + unsigned frac : _FP_FRACBITS_B - (_FP_IMPLBIT_B != 0); > +#else > + unsigned frac : _FP_FRACBITS_B - (_FP_IMPLBIT_B != 0); > + unsigned exp : _FP_EXPBITS_B; > + unsigned sign : 1; > +#endif > + } bits; > +}; > + > +#define FP_DECL_B(X) _FP_DECL (1, X) > +#define FP_UNPACK_RAW_B(X, val) _FP_UNPACK_RAW_1 (B, X, (val)) > +#define FP_UNPACK_RAW_BP(X, val) _FP_UNPACK_RAW_1_P (B, X, (val)) > +#define FP_PACK_RAW_B(val, X) _FP_PACK_RAW_1 (B, (val), X) > +#define FP_PACK_RAW_BP(val, X) \ > + do \ > + { \ > + if (!FP_INHIBIT_RESULTS) \ > + _FP_PACK_RAW_1_P (B, (val), X); \ > + } \ > + while (0) > + > +#define FP_UNPACK_B(X, val) \ > + do \ > + { \ > + _FP_UNPACK_RAW_1 (B, X, (val)); \ > + _FP_UNPACK_CANONICAL (B, 1, X); \ > + } \ > + while (0) > + > +#define FP_UNPACK_BP(X, val) \ > + do \ > + { \ > + _FP_UNPACK_RAW_1_P (B, X, (val)); \ > + _FP_UNPACK_CANONICAL (B, 1, X); \ > + } \ > + while (0) > + > +#define FP_UNPACK_SEMIRAW_B(X, val) \ > + do \ > + { \ > + _FP_UNPACK_RAW_1 (B, X, (val)); \ > + _FP_UNPACK_SEMIRAW (B, 1, X); \ > + } \ > + while (0) > + > +#define FP_UNPACK_SEMIRAW_BP(X, val) \ > + do \ > + { \ > + _FP_UNPACK_RAW_1_P (B, X, (val)); \ > + _FP_UNPACK_SEMIRAW (B, 1, X); \ > + } \ > + while (0) > + > +#define FP_PACK_B(val, X) \ > + do \ > + { \ > + _FP_PACK_CANONICAL (B, 1, X); \ > + _FP_PACK_RAW_1 (B, (val), X); \ > + } \ > + while (0) > + > +#define FP_PACK_BP(val, X) \ > + do \ > + { \ > + _FP_PACK_CANONICAL (B, 1, X); \ > + if (!FP_INHIBIT_RESULTS) \ > + _FP_PACK_RAW_1_P (B, (val), X); \ > + } \ > + while (0) > + > +#define FP_PACK_SEMIRAW_B(val, X) \ > + do \ > + { \ > + _FP_PACK_SEMIRAW (B, 1, X); \ > + _FP_PACK_RAW_1 (B, (val), X); \ > + } \ > + while (0) > + > +#define FP_PACK_SEMIRAW_BP(val, X) \ > + do \ > + { \ > + _FP_PACK_SEMIRAW (B, 1, X); \ > + if (!FP_INHIBIT_RESULTS) \ > + _FP_PACK_RAW_1_P (B, (val), X); \ > + } \ > + while (0) > + > +#define FP_TO_INT_B(r, X, rsz, rsg) _FP_TO_INT (B, 1, (r), X, (rsz), (rsg)) > +#define FP_TO_INT_ROUND_B(r, X, rsz, rsg) \ > + _FP_TO_INT_ROUND (B, 1, (r), X, (rsz), (rsg)) > +#define FP_FROM_INT_B(X, r, rs, rt) _FP_FROM_INT (B, 1, X, (r), (rs), rt) > + > +/* BFmode arithmetic is not implemented. */ > + > +#define _FP_FRAC_HIGH_B(X) _FP_FRAC_HIGH_1 (X) > +#define _FP_FRAC_HIGH_RAW_B(X) _FP_FRAC_HIGH_1 (X) > +#define _FP_FRAC_HIGH_DW_B(X) _FP_FRAC_HIGH_1 (X) > + > +#define FP_CMP_EQ_B(r, X, Y, ex) _FP_CMP_EQ (B, 1, (r), X, Y, (ex)) > + > +#endif /* !SOFT_FP_BRAIN_H */ > diff --git a/soft-fp/extendbfsf2.c b/soft-fp/extendbfsf2.c > new file mode 100644 > index 0000000000..3589e57e31 > --- /dev/null > +++ b/soft-fp/extendbfsf2.c > @@ -0,0 +1,49 @@ > +/* Software floating-point emulation. > + Return an bfloat16 converted to IEEE single > + Copyright (C) 2022-2023 Free Software Foundation, Inc. > + This file is part of the GNU C Library. > + > + The GNU C Library is free software; you can redistribute it and/or > + modify it under the terms of the GNU Lesser General Public > + License as published by the Free Software Foundation; either > + version 2.1 of the License, or (at your option) any later version. > + > + In addition to the permissions in the GNU Lesser General Public > + License, the Free Software Foundation gives you unlimited > + permission to link the compiled version of this file into > + combinations with other programs, and to distribute those > + combinations without any restriction coming from the use of this > + file. (The Lesser General Public License restrictions do apply in > + other respects; for example, they cover modification of the file, > + and distribution when not linked into a combine executable.) > + > + The GNU C Library is distributed in the hope that it will be useful, > + but WITHOUT ANY WARRANTY; without even the implied warranty of > + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU > + Lesser General Public License for more details. > + > + You should have received a copy of the GNU Lesser General Public > + License along with the GNU C Library; if not, see > + . */ > + > +#define FP_NO_EXACT_UNDERFLOW > +#include "soft-fp.h" > +#include "brain.h" > +#include "single.h" > + > +SFtype > +__extendbfsf2 (BFtype a) > +{ > + FP_DECL_EX; > + FP_DECL_B (A); > + FP_DECL_S (R); > + SFtype r; > + > + FP_INIT_EXCEPTIONS; > + FP_UNPACK_RAW_B (A, a); > + FP_EXTEND (S, B, 1, 1, R, A); > + FP_PACK_RAW_S (r, R); > + FP_HANDLE_EXCEPTIONS; > + > + return r; > +} > diff --git a/soft-fp/floattibf.c b/soft-fp/floattibf.c > new file mode 100644 > index 0000000000..f3336a95b0 > --- /dev/null > +++ b/soft-fp/floattibf.c > @@ -0,0 +1,45 @@ > +/* Software floating-point emulation. > + Convert a 128bit signed integer to bfloat16 > + Copyright (C) 2007-2023 Free Software Foundation, Inc. > + This file is part of the GNU C Library. > + > + The GNU C Library is free software; you can redistribute it and/or > + modify it under the terms of the GNU Lesser General Public > + License as published by the Free Software Foundation; either > + version 2.1 of the License, or (at your option) any later version. > + > + In addition to the permissions in the GNU Lesser General Public > + License, the Free Software Foundation gives you unlimited > + permission to link the compiled version of this file into > + combinations with other programs, and to distribute those > + combinations without any restriction coming from the use of this > + file. (The Lesser General Public License restrictions do apply in > + other respects; for example, they cover modification of the file, > + and distribution when not linked into a combine executable.) > + > + The GNU C Library is distributed in the hope that it will be useful, > + but WITHOUT ANY WARRANTY; without even the implied warranty of > + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU > + Lesser General Public License for more details. > + > + You should have received a copy of the GNU Lesser General Public > + License along with the GNU C Library; if not, see > + . */ > + > +#include "soft-fp.h" > +#include "brain.h" > + > +BFtype > +__floattibf (TItype i) > +{ > + FP_DECL_EX; > + FP_DECL_B (A); > + BFtype a; > + > + FP_INIT_ROUNDMODE; > + FP_FROM_INT_B (A, i, TI_BITS, UTItype); > + FP_PACK_RAW_B (a, A); > + FP_HANDLE_EXCEPTIONS; > + > + return a; > +} > diff --git a/soft-fp/floatuntibf.c b/soft-fp/floatuntibf.c > new file mode 100644 > index 0000000000..362b0f6c1b > --- /dev/null > +++ b/soft-fp/floatuntibf.c > @@ -0,0 +1,45 @@ > +/* Software floating-point emulation. > + Convert a 128bit unsigned integer to bfloat16 > + Copyright (C) 2007-2023 Free Software Foundation, Inc. > + This file is part of the GNU C Library. > + > + The GNU C Library is free software; you can redistribute it and/or > + modify it under the terms of the GNU Lesser General Public > + License as published by the Free Software Foundation; either > + version 2.1 of the License, or (at your option) any later version. > + > + In addition to the permissions in the GNU Lesser General Public > + License, the Free Software Foundation gives you unlimited > + permission to link the compiled version of this file into > + combinations with other programs, and to distribute those > + combinations without any restriction coming from the use of this > + file. (The Lesser General Public License restrictions do apply in > + other respects; for example, they cover modification of the file, > + and distribution when not linked into a combine executable.) > + > + The GNU C Library is distributed in the hope that it will be useful, > + but WITHOUT ANY WARRANTY; without even the implied warranty of > + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU > + Lesser General Public License for more details. > + > + You should have received a copy of the GNU Lesser General Public > + License along with the GNU C Library; if not, see > + . */ > + > +#include "soft-fp.h" > +#include "brain.h" > + > +BFtype > +__floatuntibf (UTItype i) > +{ > + FP_DECL_EX; > + FP_DECL_B (A); > + BFtype a; > + > + FP_INIT_ROUNDMODE; > + FP_FROM_INT_B (A, i, TI_BITS, UTItype); > + FP_PACK_RAW_B (a, A); > + FP_HANDLE_EXCEPTIONS; > + > + return a; > +} > diff --git a/soft-fp/truncbfhf2.c b/soft-fp/truncbfhf2.c > new file mode 100644 > index 0000000000..655c5f45ef > --- /dev/null > +++ b/soft-fp/truncbfhf2.c > @@ -0,0 +1,75 @@ > +/* Software floating-point emulation. > + Truncate bfloat16 into IEEE half. > + Copyright (C) 2022-2023 Free Software Foundation, Inc. > + This file is part of the GNU C Library. > + > + The GNU C Library is free software; you can redistribute it and/or > + modify it under the terms of the GNU Lesser General Public > + License as published by the Free Software Foundation; either > + version 2.1 of the License, or (at your option) any later version. > + > + In addition to the permissions in the GNU Lesser General Public > + License, the Free Software Foundation gives you unlimited > + permission to link the compiled version of this file into > + combinations with other programs, and to distribute those > + combinations without any restriction coming from the use of this > + file. (The Lesser General Public License restrictions do apply in > + other respects; for example, they cover modification of the file, > + and distribution when not linked into a combine executable.) > + > + The GNU C Library is distributed in the hope that it will be useful, > + but WITHOUT ANY WARRANTY; without even the implied warranty of > + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU > + Lesser General Public License for more details. > + > + You should have received a copy of the GNU Lesser General Public > + License along with the GNU C Library; if not, see > + . */ > + > +#include "soft-fp.h" > +#include "half.h" > +#include "brain.h" > +#include "single.h" > + > +/* BFtype and HFtype are unordered, neither is a superset or subset > + of each other. Convert BFtype to SFtype (lossless) and then > + truncate to HFtype. */ > + > +HFtype > +__truncbfhf2 (BFtype a) > +{ > + FP_DECL_EX; > + FP_DECL_H (A); > + FP_DECL_S (B); > + FP_DECL_B (R); > + SFtype b; > + HFtype r; > + > + FP_INIT_ROUNDMODE; > + /* Optimize BFtype to SFtype conversion to simple left shift > + by 16 if possible, we don't need to raise exceptions on sNaN > + here as the SFtype to HFtype truncation should do that too. */ > + if (sizeof (BFtype) == 2 > + && sizeof (unsigned short) == 2 > + && sizeof (SFtype) == 4 > + && sizeof (unsigned int) == 4) > + { > + union { BFtype a; unsigned short b; } u1; > + union { SFtype a; unsigned int b; } u2; > + u1.a = a; > + u2.b = (u1.b << 8) << 8; > + b = u2.a; > + } > + else > + { > + FP_UNPACK_RAW_B (A, a); > + FP_EXTEND (S, B, 1, 1, B, A); > + FP_PACK_RAW_S (b, B); > + } > + FP_UNPACK_SEMIRAW_S (B, b); > + FP_TRUNC (H, S, 1, 1, R, B); > + FP_PACK_SEMIRAW_H (r, R); > + FP_HANDLE_EXCEPTIONS; > + > + return r; > +} > diff --git a/soft-fp/truncdfbf2.c b/soft-fp/truncdfbf2.c > new file mode 100644 > index 0000000000..fd3ff525a1 > --- /dev/null > +++ b/soft-fp/truncdfbf2.c > @@ -0,0 +1,52 @@ > +/* Software floating-point emulation. > + Truncate IEEE double into bfloat16. > + Copyright (C) 2022-2023 Free Software Foundation, Inc. > + This file is part of the GNU C Library. > + > + The GNU C Library is free software; you can redistribute it and/or > + modify it under the terms of the GNU Lesser General Public > + License as published by the Free Software Foundation; either > + version 2.1 of the License, or (at your option) any later version. > + > + In addition to the permissions in the GNU Lesser General Public > + License, the Free Software Foundation gives you unlimited > + permission to link the compiled version of this file into > + combinations with other programs, and to distribute those > + combinations without any restriction coming from the use of this > + file. (The Lesser General Public License restrictions do apply in > + other respects; for example, they cover modification of the file, > + and distribution when not linked into a combine executable.) > + > + The GNU C Library is distributed in the hope that it will be useful, > + but WITHOUT ANY WARRANTY; without even the implied warranty of > + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU > + Lesser General Public License for more details. > + > + You should have received a copy of the GNU Lesser General Public > + License along with the GNU C Library; if not, see > + . */ > + > +#include "soft-fp.h" > +#include "brain.h" > +#include "double.h" > + > +BFtype > +__truncdfbf2 (DFtype a) > +{ > + FP_DECL_EX; > + FP_DECL_D (A); > + FP_DECL_B (R); > + BFtype r; > + > + FP_INIT_ROUNDMODE; > + FP_UNPACK_SEMIRAW_D (A, a); > +#if _FP_W_TYPE_SIZE < _FP_FRACBITS_D > + FP_TRUNC (B, D, 1, 2, R, A); > +#else > + FP_TRUNC (B, D, 1, 1, R, A); > +#endif > + FP_PACK_SEMIRAW_B (r, R); > + FP_HANDLE_EXCEPTIONS; > + > + return r; > +} > diff --git a/soft-fp/trunchfbf2.c b/soft-fp/trunchfbf2.c > new file mode 100644 > index 0000000000..248bed830c > --- /dev/null > +++ b/soft-fp/trunchfbf2.c > @@ -0,0 +1,58 @@ > +/* Software floating-point emulation. > + Truncate IEEE half into bfloat16. > + Copyright (C) 2022-2023 Free Software Foundation, Inc. > + This file is part of the GNU C Library. > + > + The GNU C Library is free software; you can redistribute it and/or > + modify it under the terms of the GNU Lesser General Public > + License as published by the Free Software Foundation; either > + version 2.1 of the License, or (at your option) any later version. > + > + In addition to the permissions in the GNU Lesser General Public > + License, the Free Software Foundation gives you unlimited > + permission to link the compiled version of this file into > + combinations with other programs, and to distribute those > + combinations without any restriction coming from the use of this > + file. (The Lesser General Public License restrictions do apply in > + other respects; for example, they cover modification of the file, > + and distribution when not linked into a combine executable.) > + > + The GNU C Library is distributed in the hope that it will be useful, > + but WITHOUT ANY WARRANTY; without even the implied warranty of > + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU > + Lesser General Public License for more details. > + > + You should have received a copy of the GNU Lesser General Public > + License along with the GNU C Library; if not, see > + . */ > + > +#include "soft-fp.h" > +#include "brain.h" > +#include "half.h" > +#include "single.h" > + > +/* BFtype and HFtype are unordered, neither is a superset or subset > + of each other. Convert HFtype to SFtype (lossless) and then > + truncate to BFtype. */ > + > +BFtype > +__trunchfbf2 (HFtype a) > +{ > + FP_DECL_EX; > + FP_DECL_H (A); > + FP_DECL_S (B); > + FP_DECL_B (R); > + SFtype b; > + BFtype r; > + > + FP_INIT_ROUNDMODE; > + FP_UNPACK_RAW_H (A, a); > + FP_EXTEND (S, H, 1, 1, B, A); > + FP_PACK_RAW_S (b, B); > + FP_UNPACK_SEMIRAW_S (B, b); > + FP_TRUNC (B, S, 1, 1, R, B); > + FP_PACK_SEMIRAW_B (r, R); > + FP_HANDLE_EXCEPTIONS; > + > + return r; > +} > diff --git a/soft-fp/truncsfbf2.c b/soft-fp/truncsfbf2.c > new file mode 100644 > index 0000000000..9a5a9a0648 > --- /dev/null > +++ b/soft-fp/truncsfbf2.c > @@ -0,0 +1,48 @@ > +/* Software floating-point emulation. > + Truncate IEEE single into bfloat16. > + Copyright (C) 2022-2023 Free Software Foundation, Inc. > + This file is part of the GNU C Library. > + > + The GNU C Library is free software; you can redistribute it and/or > + modify it under the terms of the GNU Lesser General Public > + License as published by the Free Software Foundation; either > + version 2.1 of the License, or (at your option) any later version. > + > + In addition to the permissions in the GNU Lesser General Public > + License, the Free Software Foundation gives you unlimited > + permission to link the compiled version of this file into > + combinations with other programs, and to distribute those > + combinations without any restriction coming from the use of this > + file. (The Lesser General Public License restrictions do apply in > + other respects; for example, they cover modification of the file, > + and distribution when not linked into a combine executable.) > + > + The GNU C Library is distributed in the hope that it will be useful, > + but WITHOUT ANY WARRANTY; without even the implied warranty of > + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU > + Lesser General Public License for more details. > + > + You should have received a copy of the GNU Lesser General Public > + License along with the GNU C Library; if not, see > + . */ > + > +#include "soft-fp.h" > +#include "brain.h" > +#include "single.h" > + > +BFtype > +__truncsfbf2 (SFtype a) > +{ > + FP_DECL_EX; > + FP_DECL_S (A); > + FP_DECL_B (R); > + BFtype r; > + > + FP_INIT_ROUNDMODE; > + FP_UNPACK_SEMIRAW_S (A, a); > + FP_TRUNC (B, S, 1, 1, R, A); > + FP_PACK_SEMIRAW_B (r, R); > + FP_HANDLE_EXCEPTIONS; > + > + return r; > +} > diff --git a/soft-fp/trunctfbf2.c b/soft-fp/trunctfbf2.c > new file mode 100644 > index 0000000000..9b6bf93bbb > --- /dev/null > +++ b/soft-fp/trunctfbf2.c > @@ -0,0 +1,52 @@ > +/* Software floating-point emulation. > + Truncate IEEE quad into bfloat16. > + Copyright (C) 2022-2023 Free Software Foundation, Inc. > + This file is part of the GNU C Library. > + > + The GNU C Library is free software; you can redistribute it and/or > + modify it under the terms of the GNU Lesser General Public > + License as published by the Free Software Foundation; either > + version 2.1 of the License, or (at your option) any later version. > + > + In addition to the permissions in the GNU Lesser General Public > + License, the Free Software Foundation gives you unlimited > + permission to link the compiled version of this file into > + combinations with other programs, and to distribute those > + combinations without any restriction coming from the use of this > + file. (The Lesser General Public License restrictions do apply in > + other respects; for example, they cover modification of the file, > + and distribution when not linked into a combine executable.) > + > + The GNU C Library is distributed in the hope that it will be useful, > + but WITHOUT ANY WARRANTY; without even the implied warranty of > + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU > + Lesser General Public License for more details. > + > + You should have received a copy of the GNU Lesser General Public > + License along with the GNU C Library; if not, see > + . */ > + > +#include "soft-fp.h" > +#include "brain.h" > +#include "quad.h" > + > +BFtype > +__trunctfbf2 (TFtype a) > +{ > + FP_DECL_EX; > + FP_DECL_Q (A); > + FP_DECL_B (R); > + BFtype r; > + > + FP_INIT_ROUNDMODE; > + FP_UNPACK_SEMIRAW_Q (A, a); > +#if _FP_W_TYPE_SIZE < 64 > + FP_TRUNC (B, Q, 1, 4, R, A); > +#else > + FP_TRUNC (B, Q, 1, 2, R, A); > +#endif > + FP_PACK_SEMIRAW_B (r, R); > + FP_HANDLE_EXCEPTIONS; > + > + return r; > +} > diff --git a/soft-fp/truncxfbf2.c b/soft-fp/truncxfbf2.c > new file mode 100644 > index 0000000000..2c3f288c92 > --- /dev/null > +++ b/soft-fp/truncxfbf2.c > @@ -0,0 +1,52 @@ > +/* Software floating-point emulation. > + Truncate IEEE extended into bfloat16. > + Copyright (C) 2022-2023 Free Software Foundation, Inc. > + This file is part of the GNU C Library. > + > + The GNU C Library is free software; you can redistribute it and/or > + modify it under the terms of the GNU Lesser General Public > + License as published by the Free Software Foundation; either > + version 2.1 of the License, or (at your option) any later version. > + > + In addition to the permissions in the GNU Lesser General Public > + License, the Free Software Foundation gives you unlimited > + permission to link the compiled version of this file into > + combinations with other programs, and to distribute those > + combinations without any restriction coming from the use of this > + file. (The Lesser General Public License restrictions do apply in > + other respects; for example, they cover modification of the file, > + and distribution when not linked into a combine executable.) > + > + The GNU C Library is distributed in the hope that it will be useful, > + but WITHOUT ANY WARRANTY; without even the implied warranty of > + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU > + Lesser General Public License for more details. > + > + You should have received a copy of the GNU Lesser General Public > + License along with the GNU C Library; if not, see > + . */ > + > +#include "soft-fp.h" > +#include "brain.h" > +#include "extended.h" > + > +BFtype > +__truncxfbf2 (XFtype a) > +{ > + FP_DECL_EX; > + FP_DECL_E (A); > + FP_DECL_B (R); > + BFtype r; > + > + FP_INIT_ROUNDMODE; > + FP_UNPACK_SEMIRAW_E (A, a); > +#if _FP_W_TYPE_SIZE < 64 > + FP_TRUNC (B, E, 1, 4, R, A); > +#else > + FP_TRUNC (B, E, 1, 2, R, A); > +#endif > + FP_PACK_SEMIRAW_B (r, R); > + FP_HANDLE_EXCEPTIONS; > + > + return r; > +} > > Jakub >