xref: /freebsd/lib/libc/arm/aeabi/aeabi_float.c (revision 6eeedf35c3ecf50cc3def1eaa40a37b2257b85ac)
1d915a14eSPedro F. Giffuni /*-
24d846d26SWarner Losh  * SPDX-License-Identifier: BSD-2-Clause
3d915a14eSPedro F. Giffuni  *
401298368SAndrew Turner  * Copyright (C) 2012 Andrew Turner
501298368SAndrew Turner  * All rights reserved.
601298368SAndrew Turner  *
701298368SAndrew Turner  * Redistribution and use in source and binary forms, with or without
801298368SAndrew Turner  * modification, are permitted provided that the following conditions
901298368SAndrew Turner  * are met:
1001298368SAndrew Turner  * 1. Redistributions of source code must retain the above copyright
1101298368SAndrew Turner  *    notice, this list of conditions and the following disclaimer.
1201298368SAndrew Turner  * 2. Redistributions in binary form must reproduce the above copyright
1301298368SAndrew Turner  *    notice, this list of conditions and the following disclaimer in the
1401298368SAndrew Turner  *    documentation and/or other materials provided with the distribution.
1501298368SAndrew Turner  *
1601298368SAndrew Turner  * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND
1701298368SAndrew Turner  * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
1801298368SAndrew Turner  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
1901298368SAndrew Turner  * ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
2001298368SAndrew Turner  * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
2101298368SAndrew Turner  * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
2201298368SAndrew Turner  * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
2301298368SAndrew Turner  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
2401298368SAndrew Turner  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
2501298368SAndrew Turner  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
2601298368SAndrew Turner  * SUCH DAMAGE.
2701298368SAndrew Turner  *
2801298368SAndrew Turner  */
2901298368SAndrew Turner 
3001298368SAndrew Turner #include "softfloat-for-gcc.h"
3101298368SAndrew Turner #include "milieu.h"
3201298368SAndrew Turner #include "softfloat.h"
3301298368SAndrew Turner 
34a06918a3SAndrew Turner #include "aeabi_vfp.h"
35a06918a3SAndrew Turner 
36a06918a3SAndrew Turner extern int _libc_arm_fpu_present;
37a06918a3SAndrew Turner 
3801298368SAndrew Turner flag __unordsf2(float32, float32);
3901298368SAndrew Turner 
40a06918a3SAndrew Turner /* These are written in asm and are only called from this file */
41a06918a3SAndrew Turner int __aeabi_fcmpeq_vfp(float32, float32);
42a06918a3SAndrew Turner int __aeabi_fcmplt_vfp(float32, float32);
43a06918a3SAndrew Turner int __aeabi_fcmple_vfp(float32, float32);
44a06918a3SAndrew Turner int __aeabi_fcmpgt_vfp(float32, float32);
45a06918a3SAndrew Turner int __aeabi_fcmpge_vfp(float32, float32);
46a06918a3SAndrew Turner int __aeabi_fcmpun_vfp(float32, float32);
47a06918a3SAndrew Turner int __aeabi_f2iz_vfp(float32);
48a06918a3SAndrew Turner float64 __aeabi_f2d_vfp(float32);
49a06918a3SAndrew Turner float32 __aeabi_i2f_vfp(int);
50a06918a3SAndrew Turner float32 __aeabi_fadd_vfp(float32, float32);
51a06918a3SAndrew Turner float32 __aeabi_fdiv_vfp(float32, float32);
52a06918a3SAndrew Turner float32 __aeabi_fmul_vfp(float32, float32);
53a06918a3SAndrew Turner float32 __aeabi_fsub_vfp(float32, float32);
5401298368SAndrew Turner 
55a06918a3SAndrew Turner /*
56a06918a3SAndrew Turner  * Depending on the target these will:
57*6eeedf35SAndrew Turner  *  On armv7 with a vfp call the above function, or
58aac3f7c7SAndrew Turner  *  Call the softfloat function in the 3rd argument.
59a06918a3SAndrew Turner  */
AEABI_FUNC2(fcmpeq,float32,float32_eq)60a06918a3SAndrew Turner int AEABI_FUNC2(fcmpeq, float32, float32_eq)
61a06918a3SAndrew Turner int AEABI_FUNC2(fcmplt, float32, float32_lt)
62a06918a3SAndrew Turner int AEABI_FUNC2(fcmple, float32, float32_le)
63a06918a3SAndrew Turner int AEABI_FUNC2_REV(fcmpge, float32, float32_le)
64a06918a3SAndrew Turner int AEABI_FUNC2_REV(fcmpgt, float32, float32_lt)
65a06918a3SAndrew Turner int AEABI_FUNC2(fcmpun, float32, __unordsf2)
6601298368SAndrew Turner 
67a06918a3SAndrew Turner int AEABI_FUNC(f2iz, float32, float32_to_int32_round_to_zero)
68a06918a3SAndrew Turner float64 AEABI_FUNC(f2d, float32, float32_to_float64)
69a06918a3SAndrew Turner float32 AEABI_FUNC(i2f, int, int32_to_float32)
7001298368SAndrew Turner 
71a06918a3SAndrew Turner float32 AEABI_FUNC2(fadd, float32, float32_add)
72a06918a3SAndrew Turner float32 AEABI_FUNC2(fdiv, float32, float32_div)
73a06918a3SAndrew Turner float32 AEABI_FUNC2(fmul, float32, float32_mul)
74a06918a3SAndrew Turner float32 AEABI_FUNC2(fsub, float32, float32_sub)
7501298368SAndrew Turner 
76cb8bcafaSAndrew Turner int
77cb8bcafaSAndrew Turner __aeabi_cfcmpeq_helper(float32 a, float32 b)
78cb8bcafaSAndrew Turner {
79cb8bcafaSAndrew Turner 	int quiet = 0;
80cb8bcafaSAndrew Turner 
81cb8bcafaSAndrew Turner 	/* Check if a is a NaN */
82cb8bcafaSAndrew Turner 	if ((a << 1) > 0xff000000u) {
83cb8bcafaSAndrew Turner 		/* If it's a signalling NaN we will always signal */
84cb8bcafaSAndrew Turner 		if ((a & 0x00400000u) == 0)
85cb8bcafaSAndrew Turner 			return (0);
86cb8bcafaSAndrew Turner 
87cb8bcafaSAndrew Turner 		quiet = 1;
88cb8bcafaSAndrew Turner 	}
89cb8bcafaSAndrew Turner 
90cb8bcafaSAndrew Turner 	/* Check if b is a NaN */
91cb8bcafaSAndrew Turner 	if ((b << 1) > 0xff000000u) {
92cb8bcafaSAndrew Turner 		/* If it's a signalling NaN we will always signal */
93cb8bcafaSAndrew Turner 		if ((b & 0x00400000u) == 0)
94cb8bcafaSAndrew Turner 			return (0);
95cb8bcafaSAndrew Turner 
96cb8bcafaSAndrew Turner 		quiet = 1;
97cb8bcafaSAndrew Turner 	}
98cb8bcafaSAndrew Turner 
99cb8bcafaSAndrew Turner 	return (quiet);
100cb8bcafaSAndrew Turner }
101