From mboxrd@z Thu Jan 1 00:00:00 1970 Return-Path: Received: by sourceware.org (Postfix, from userid 2155) id 8A86F3857C45; Mon, 4 Dec 2023 09:33:41 +0000 (GMT) DKIM-Filter: OpenDKIM Filter v2.11.0 sourceware.org 8A86F3857C45 DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=sourceware.org; s=default; t=1701682421; bh=xdYlHJVjopyJmuj8LrwV4hduPMVL1RJhYUG97lIiXLU=; h=From:To:Subject:Date:From; b=U9lm2EnJTXwb/muPiH3Pkq3U+CaOhAaybqb83LCEzUqQ/PfjRJsiZTO3H2eqgSrVL X0tiK9lsmxGp7up7cScYGQhh7L5etjXeF4VwVHHiPUS4fMQvbHwkYLeMTHUYKKexic iRFrYNbekLvreb2a0qxJMf453YSB14szb6tRXa6M= Content-Type: text/plain; charset="utf-8" MIME-Version: 1.0 Content-Transfer-Encoding: quoted-printable From: Corinna Vinschen To: newlib-cvs@sourceware.org Subject: [newlib-cygwin/main] RISC-V: Support long double math X-Act-Checkin: newlib-cygwin X-Git-Author: Kito Cheng X-Git-Refname: refs/heads/main X-Git-Oldrev: bcd1883560b79b67f06fa1df6dea60d18fdd00c7 X-Git-Newrev: 04798b7bb69571452d2cfc7e0b052a9bbd3b619d Message-Id: <20231204093341.8A86F3857C45@sourceware.org> Date: Mon, 4 Dec 2023 09:33:41 +0000 (GMT) List-Id: https://sourceware.org/git/gitweb.cgi?p=3Dnewlib-cygwin.git;h=3D04798b7bb69= 571452d2cfc7e0b052a9bbd3b619d commit 04798b7bb69571452d2cfc7e0b052a9bbd3b619d Author: Kito Cheng AuthorDate: Mon Dec 4 15:41:39 2023 +0800 Commit: Corinna Vinschen CommitDate: Mon Dec 4 10:26:41 2023 +0100 RISC-V: Support long double math =20 Long double for RISC-V is using 128 bit IEEE 754 format like Aarch64, so we reference AArch64 to support that. Diff: --- newlib/Makefile.in | 63 +++++++++++++++++++------= --- newlib/libc/include/math.h | 3 +- newlib/libc/machine/riscv/machine/_fpmath.h | 64 +++++++++++++++++++++++++= ++++ newlib/libm/Makefile.inc | 1 + 4 files changed, 110 insertions(+), 21 deletions(-) diff --git a/newlib/Makefile.in b/newlib/Makefile.in index 9cbb4254ec23..3cba156d3f9b 100644 --- a/newlib/Makefile.in +++ b/newlib/Makefile.in @@ -947,14 +947,17 @@ check_PROGRAMS =3D @HAVE_LIBM_MACHINE_SPARC_TRUE@am__append_159 =3D $(libm_machine_sparc_src) @HAVE_LIBM_MACHINE_SPU_TRUE@am__append_160 =3D $(libm_machine_spu_src) @HAVE_LIBM_MACHINE_RISCV_TRUE@am__append_161 =3D $(libm_machine_riscv_src) -@HAVE_LIBM_MACHINE_X86_64_TRUE@am__append_162 =3D $(libm_machine_x86_64_sr= c) -@HAVE_LIBM_MACHINE_X86_64_TRUE@@HAVE_LONG_DOUBLE_TRUE@am__append_163 =3D $= (libm_ld80_lsrc) -@HAVE_LIBM_MACHINE_X86_64_TRUE@am__append_164 =3D=20 -@HAVE_LIBM_MACHINE_X86_64_TRUE@am__append_165 =3D=20 -@HAVE_LIBM_MACHINE_XTENSA_TRUE@@XTENSA_XCHAL_HAVE_FP_SQRT_TRUE@am__append_= 166 =3D \ +@HAVE_LIBM_MACHINE_RISCV_TRUE@@HAVE_LONG_DOUBLE_TRUE@am__append_162 =3D $(= libm_ld128_lsrc) +@HAVE_LIBM_MACHINE_RISCV_TRUE@am__append_163 =3D +@HAVE_LIBM_MACHINE_RISCV_TRUE@am__append_164 =3D +@HAVE_LIBM_MACHINE_X86_64_TRUE@am__append_165 =3D $(libm_machine_x86_64_sr= c) +@HAVE_LIBM_MACHINE_X86_64_TRUE@@HAVE_LONG_DOUBLE_TRUE@am__append_166 =3D $= (libm_ld80_lsrc) +@HAVE_LIBM_MACHINE_X86_64_TRUE@am__append_167 =3D +@HAVE_LIBM_MACHINE_X86_64_TRUE@am__append_168 =3D +@HAVE_LIBM_MACHINE_XTENSA_TRUE@@XTENSA_XCHAL_HAVE_FP_SQRT_TRUE@am__append_= 169 =3D \ @HAVE_LIBM_MACHINE_XTENSA_TRUE@@XTENSA_XCHAL_HAVE_FP_SQRT_TRUE@ libm/machi= ne/xtensa/ef_sqrt.c =20 -@HAVE_LIBM_MACHINE_XTENSA_TRUE@am__append_167 =3D $(libm_machine_xtensa_sr= c) +@HAVE_LIBM_MACHINE_XTENSA_TRUE@am__append_170 =3D $(libm_machine_xtensa_sr= c) subdir =3D . ACLOCAL_M4 =3D $(top_srcdir)/aclocal.m4 am__aclocal_m4_deps =3D $(top_srcdir)/../config/depstand.m4 \ @@ -3198,6 +3201,16 @@ am__objects_159 =3D libm/fenv/libm_a-feclearexcept.$= (OBJEXT) \ @HAVE_LIBM_MACHINE_AARCH64_TRUE@ libm/machine/aarch64/libm_a-fetestexcept.= $(OBJEXT) \ @HAVE_LIBM_MACHINE_AARCH64_TRUE@ libm/machine/aarch64/libm_a-feupdateenv.$= (OBJEXT) @HAVE_LIBM_MACHINE_AARCH64_TRUE@am__objects_161 =3D $(am__objects_160) +@HAVE_LIBM_MACHINE_AARCH64_FALSE@@HAVE_LIBM_MACHINE_RISCV_TRUE@am__objects= _162 =3D libm/ld128/libm_a-e_powl.$(OBJEXT) \ +@HAVE_LIBM_MACHINE_AARCH64_FALSE@@HAVE_LIBM_MACHINE_RISCV_TRUE@ libm/ld128= /libm_a-s_erfl.$(OBJEXT) \ +@HAVE_LIBM_MACHINE_AARCH64_FALSE@@HAVE_LIBM_MACHINE_RISCV_TRUE@ libm/ld128= /libm_a-s_exp2l.$(OBJEXT) \ +@HAVE_LIBM_MACHINE_AARCH64_FALSE@@HAVE_LIBM_MACHINE_RISCV_TRUE@ libm/ld128= /libm_a-s_expl.$(OBJEXT) \ +@HAVE_LIBM_MACHINE_AARCH64_FALSE@@HAVE_LIBM_MACHINE_RISCV_TRUE@ libm/ld128= /libm_a-s_logl.$(OBJEXT) \ +@HAVE_LIBM_MACHINE_AARCH64_FALSE@@HAVE_LIBM_MACHINE_RISCV_TRUE@ libm/ld128= /libm_a-b_tgammal.$(OBJEXT) \ +@HAVE_LIBM_MACHINE_AARCH64_FALSE@@HAVE_LIBM_MACHINE_RISCV_TRUE@ libm/ld128= /libm_a-invtrig.$(OBJEXT) \ +@HAVE_LIBM_MACHINE_AARCH64_FALSE@@HAVE_LIBM_MACHINE_RISCV_TRUE@ libm/ld128= /libm_a-e_lgammal_r.$(OBJEXT) \ +@HAVE_LIBM_MACHINE_AARCH64_FALSE@@HAVE_LIBM_MACHINE_RISCV_TRUE@ libm/ld128= /libm_a-k_cosl.$(OBJEXT) \ +@HAVE_LIBM_MACHINE_AARCH64_FALSE@@HAVE_LIBM_MACHINE_RISCV_TRUE@ libm/ld128= /libm_a-k_sinl.$(OBJEXT) @HAVE_LIBM_MACHINE_AARCH64_TRUE@am__objects_162 =3D libm/ld128/libm_a-e_po= wl.$(OBJEXT) \ @HAVE_LIBM_MACHINE_AARCH64_TRUE@ libm/ld128/libm_a-s_erfl.$(OBJEXT) \ @HAVE_LIBM_MACHINE_AARCH64_TRUE@ libm/ld128/libm_a-s_exp2l.$(OBJEXT) \ @@ -3615,7 +3628,8 @@ am__objects_159 =3D libm/fenv/libm_a-feclearexcept.$(= OBJEXT) \ @HAVE_LIBM_MACHINE_RISCV_TRUE@ libm/machine/riscv/libm_a-s_llround.$(OBJEX= T) \ @HAVE_LIBM_MACHINE_RISCV_TRUE@ libm/machine/riscv/libm_a-sf_llround.$(OBJE= XT) @HAVE_LIBM_MACHINE_RISCV_TRUE@am__objects_187 =3D $(am__objects_186) -@HAVE_LIBM_MACHINE_X86_64_TRUE@am__objects_188 =3D libm/machine/x86_64/lib= m_a-feclearexcept.$(OBJEXT) \ +@HAVE_LIBM_MACHINE_RISCV_TRUE@@HAVE_LONG_DOUBLE_TRUE@am__objects_188 =3D $= (am__objects_162) +@HAVE_LIBM_MACHINE_X86_64_TRUE@am__objects_189 =3D libm/machine/x86_64/lib= m_a-feclearexcept.$(OBJEXT) \ @HAVE_LIBM_MACHINE_X86_64_TRUE@ libm/machine/x86_64/libm_a-fegetenv.$(OBJE= XT) \ @HAVE_LIBM_MACHINE_X86_64_TRUE@ libm/machine/x86_64/libm_a-fegetexceptflag= .$(OBJEXT) \ @HAVE_LIBM_MACHINE_X86_64_TRUE@ libm/machine/x86_64/libm_a-fegetround.$(OB= JEXT) \ @@ -3627,10 +3641,10 @@ am__objects_159 =3D libm/fenv/libm_a-feclearexcept.= $(OBJEXT) \ @HAVE_LIBM_MACHINE_X86_64_TRUE@ libm/machine/x86_64/libm_a-fesetround.$(OB= JEXT) \ @HAVE_LIBM_MACHINE_X86_64_TRUE@ libm/machine/x86_64/libm_a-fetestexcept.$(= OBJEXT) \ @HAVE_LIBM_MACHINE_X86_64_TRUE@ libm/machine/x86_64/libm_a-feupdateenv.$(O= BJEXT) -@HAVE_LIBM_MACHINE_X86_64_TRUE@am__objects_189 =3D $(am__objects_188) -@HAVE_LIBM_MACHINE_X86_64_TRUE@@HAVE_LONG_DOUBLE_TRUE@am__objects_190 =3D = $(am__objects_170) -@HAVE_LIBM_MACHINE_XTENSA_TRUE@@XTENSA_XCHAL_HAVE_FP_SQRT_TRUE@am__objects= _191 =3D libm/machine/xtensa/libm_a-ef_sqrt.$(OBJEXT) -@HAVE_LIBM_MACHINE_XTENSA_TRUE@am__objects_192 =3D libm/machine/xtensa/lib= m_a-feclearexcept.$(OBJEXT) \ +@HAVE_LIBM_MACHINE_X86_64_TRUE@am__objects_190 =3D $(am__objects_189) +@HAVE_LIBM_MACHINE_X86_64_TRUE@@HAVE_LONG_DOUBLE_TRUE@am__objects_191 =3D = $(am__objects_170) +@HAVE_LIBM_MACHINE_XTENSA_TRUE@@XTENSA_XCHAL_HAVE_FP_SQRT_TRUE@am__objects= _192 =3D libm/machine/xtensa/libm_a-ef_sqrt.$(OBJEXT) +@HAVE_LIBM_MACHINE_XTENSA_TRUE@am__objects_193 =3D libm/machine/xtensa/lib= m_a-feclearexcept.$(OBJEXT) \ @HAVE_LIBM_MACHINE_XTENSA_TRUE@ libm/machine/xtensa/libm_a-fegetenv.$(OBJE= XT) \ @HAVE_LIBM_MACHINE_XTENSA_TRUE@ libm/machine/xtensa/libm_a-fegetexcept.$(O= BJEXT) \ @HAVE_LIBM_MACHINE_XTENSA_TRUE@ libm/machine/xtensa/libm_a-fegetexceptflag= .$(OBJEXT) \ @@ -3639,8 +3653,8 @@ am__objects_159 =3D libm/fenv/libm_a-feclearexcept.$(= OBJEXT) \ @HAVE_LIBM_MACHINE_XTENSA_TRUE@ libm/machine/xtensa/libm_a-feraiseexcept.$= (OBJEXT) \ @HAVE_LIBM_MACHINE_XTENSA_TRUE@ libm/machine/xtensa/libm_a-fetestexcept.$(= OBJEXT) \ @HAVE_LIBM_MACHINE_XTENSA_TRUE@ libm/machine/xtensa/libm_a-feupdateenv.$(O= BJEXT) \ -@HAVE_LIBM_MACHINE_XTENSA_TRUE@ $(am__objects_191) -@HAVE_LIBM_MACHINE_XTENSA_TRUE@am__objects_193 =3D $(am__objects_192) +@HAVE_LIBM_MACHINE_XTENSA_TRUE@ $(am__objects_192) +@HAVE_LIBM_MACHINE_XTENSA_TRUE@am__objects_194 =3D $(am__objects_193) am_libm_a_OBJECTS =3D $(am__objects_145) $(am__objects_149) \ $(am__objects_150) $(am__objects_151) $(am__objects_153) \ $(am__objects_155) $(am__objects_156) $(am__objects_157) \ @@ -3649,7 +3663,8 @@ am_libm_a_OBJECTS =3D $(am__objects_145) $(am__object= s_149) \ $(am__objects_169) $(am__objects_171) $(am__objects_173) \ $(am__objects_177) $(am__objects_179) $(am__objects_181) \ $(am__objects_183) $(am__objects_185) $(am__objects_187) \ - $(am__objects_189) $(am__objects_190) $(am__objects_193) + $(am__objects_188) $(am__objects_190) $(am__objects_191) \ + $(am__objects_194) libm_a_OBJECTS =3D $(am_libm_a_OBJECTS) am_libm_test_test_OBJECTS =3D libm/test/test.$(OBJEXT) \ libm/test/string.$(OBJEXT) libm/test/convert.$(OBJEXT) \ @@ -4209,7 +4224,8 @@ libm_a_SOURCES =3D $(am__append_133) $(am__append_136= ) \ $(am__append_149) $(am__append_150) $(am__append_153) \ $(am__append_156) $(am__append_157) $(am__append_158) \ $(am__append_159) $(am__append_160) $(am__append_161) \ - $(am__append_162) $(am__append_163) $(am__append_167) + $(am__append_162) $(am__append_165) $(am__append_166) \ + $(am__append_170) libm_a_CFLAGS =3D $(AM_CFLAGS) $(libm_a_CFLAGS_$(subst /,_,$(@D))) $(libm_= a_CFLAGS_$(subst /,_,$(@D)_$( + * Copyright (2) 2014 The FreeBSD Foundation + * All rights reserved. + * + * Redistribution and use in source and binary forms, with or without + * modification, are permitted provided that the following conditions + * are met: + * 1. Redistributions of source code must retain the above copyright + * notice, this list of conditions and the following disclaimer. + * 2. Redistributions in binary form must reproduce the above copyright + * notice, this list of conditions and the following disclaimer in the + * documentation and/or other materials provided with the distribution. + * + * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND + * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE + * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURP= OSE + * ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE + * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENT= IAL + * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS + * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) + * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STR= ICT + * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY W= AY + * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF + * SUCH DAMAGE. + * + * $FreeBSD$ + */ + +/* + * Change unsigned int/long used by FreeBSD to fixed width types because + * ilp32 has a different size for unsigned long. --joel (20 Aug 2022) + */ +#include + +union IEEEl2bits { + long double e; + struct { + uint64_t manl :64; + uint64_t manh :48; + uint32_t exp :15; + uint32_t sign :1; + } bits; + /* TODO andrew: Check the packing here */ + struct { + uint64_t manl :64; + uint64_t manh :48; + uint32_t expsign :16; + } xbits; +}; + +#define LDBL_NBIT 0 +#define LDBL_IMPLICIT_NBIT +#define mask_nbit_l(u) ((void)0) + +#define LDBL_MANH_SIZE 48 +#define LDBL_MANL_SIZE 64 + +#define LDBL_TO_ARRAY32(u, a) do { \ + (a)[0] =3D (uint32_t)(u).bits.manl; \ + (a)[1] =3D (uint32_t)((u).bits.manl >> 32); \ + (a)[2] =3D (uint32_t)(u).bits.manh; \ + (a)[3] =3D (uint32_t)((u).bits.manh >> 32); \ +} while(0) diff --git a/newlib/libm/Makefile.inc b/newlib/libm/Makefile.inc index 28b84bfe2ede..bf31b1be788c 100644 --- a/newlib/libm/Makefile.inc +++ b/newlib/libm/Makefile.inc @@ -87,6 +87,7 @@ include %D%/machine/spu/Makefile.inc endif if HAVE_LIBM_MACHINE_RISCV include %D%/machine/riscv/Makefile.inc +include %D%/ld128/Makefile.inc endif if HAVE_LIBM_MACHINE_X86_64 include %D%/machine/x86_64/Makefile.inc