xref: /freebsd/lib/libc/arm/aeabi/aeabi_double.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 __unorddf2(float64, float64);
3901298368SAndrew Turner 
40a06918a3SAndrew Turner /* These are written in asm and are only called from this file */
41a06918a3SAndrew Turner int __aeabi_dcmpeq_vfp(float64, float64);
42a06918a3SAndrew Turner int __aeabi_dcmplt_vfp(float64, float64);
43a06918a3SAndrew Turner int __aeabi_dcmple_vfp(float64, float64);
44a06918a3SAndrew Turner int __aeabi_dcmpgt_vfp(float64, float64);
45a06918a3SAndrew Turner int __aeabi_dcmpge_vfp(float64, float64);
46a06918a3SAndrew Turner int __aeabi_dcmpun_vfp(float64, float64);
47a06918a3SAndrew Turner int __aeabi_d2iz_vfp(float64);
48a06918a3SAndrew Turner float32 __aeabi_d2f_vfp(float64);
49a06918a3SAndrew Turner float64 __aeabi_i2d_vfp(int);
50a06918a3SAndrew Turner float64 __aeabi_dadd_vfp(float64, float64);
51a06918a3SAndrew Turner float64 __aeabi_ddiv_vfp(float64, float64);
52a06918a3SAndrew Turner float64 __aeabi_dmul_vfp(float64, float64);
53a06918a3SAndrew Turner float64 __aeabi_dsub_vfp(float64, float64);
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(dcmpeq,float64,float64_eq)60a06918a3SAndrew Turner int AEABI_FUNC2(dcmpeq, float64, float64_eq)
61a06918a3SAndrew Turner int AEABI_FUNC2(dcmplt, float64, float64_lt)
62a06918a3SAndrew Turner int AEABI_FUNC2(dcmple, float64, float64_le)
63a06918a3SAndrew Turner int AEABI_FUNC2_REV(dcmpge, float64, float64_le)
64a06918a3SAndrew Turner int AEABI_FUNC2_REV(dcmpgt, float64, float64_lt)
65a06918a3SAndrew Turner int AEABI_FUNC2(dcmpun, float64, __unorddf2)
6601298368SAndrew Turner 
67a06918a3SAndrew Turner int AEABI_FUNC(d2iz, float64, float64_to_int32_round_to_zero)
68a06918a3SAndrew Turner float32 AEABI_FUNC(d2f, float64, float64_to_float32)
69a06918a3SAndrew Turner float64 AEABI_FUNC(i2d, int, int32_to_float64)
7001298368SAndrew Turner 
71a06918a3SAndrew Turner float64 AEABI_FUNC2(dadd, float64, float64_add)
72a06918a3SAndrew Turner float64 AEABI_FUNC2(ddiv, float64, float64_div)
73a06918a3SAndrew Turner float64 AEABI_FUNC2(dmul, float64, float64_mul)
74a06918a3SAndrew Turner float64 AEABI_FUNC2(dsub, float64, float64_sub)
7501298368SAndrew Turner 
76cb8bcafaSAndrew Turner int
77cb8bcafaSAndrew Turner __aeabi_cdcmpeq_helper(float64 a, float64 b)
78cb8bcafaSAndrew Turner {
79cb8bcafaSAndrew Turner 	int quiet = 0;
80cb8bcafaSAndrew Turner 
81cb8bcafaSAndrew Turner 	/* Check if a is a NaN */
82cb8bcafaSAndrew Turner 	if ((a << 1) > 0xffe0000000000000ull) {
83cb8bcafaSAndrew Turner 		/* If it's a signalling NaN we will always signal */
84cb8bcafaSAndrew Turner 		if ((a & 0x0008000000000000ull) == 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) > 0xffe0000000000000ull) {
92cb8bcafaSAndrew Turner 		/* If it's a signalling NaN we will always signal */
93cb8bcafaSAndrew Turner 		if ((b & 0x0008000000000000ull) == 0)
94cb8bcafaSAndrew Turner 			return (0);
95cb8bcafaSAndrew Turner 
96cb8bcafaSAndrew Turner 		quiet = 1;
97cb8bcafaSAndrew Turner 	}
98cb8bcafaSAndrew Turner 
99cb8bcafaSAndrew Turner 	return (quiet);
100cb8bcafaSAndrew Turner }
101