xref: /freebsd/lib/libc/arm/aeabi/aeabi_double.c (revision f39bffc62c1395bde25d152c7f68fdf7cbaab414)
1 /*-
2  * SPDX-License-Identifier: BSD-2-Clause-FreeBSD
3  *
4  * Copyright (C) 2012 Andrew Turner
5  * All rights reserved.
6  *
7  * Redistribution and use in source and binary forms, with or without
8  * modification, are permitted provided that the following conditions
9  * are met:
10  * 1. Redistributions of source code must retain the above copyright
11  *    notice, this list of conditions and the following disclaimer.
12  * 2. Redistributions in binary form must reproduce the above copyright
13  *    notice, this list of conditions and the following disclaimer in the
14  *    documentation and/or other materials provided with the distribution.
15  *
16  * THIS SOFTWARE IS PROVIDED BY THE AUTHOR AND CONTRIBUTORS ``AS IS'' AND
17  * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
18  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
19  * ARE DISCLAIMED.  IN NO EVENT SHALL THE AUTHOR OR CONTRIBUTORS BE LIABLE
20  * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
21  * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
22  * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
23  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
24  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
25  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
26  * SUCH DAMAGE.
27  *
28  */
29 
30 #include <sys/cdefs.h>
31 __FBSDID("$FreeBSD$");
32 
33 #include "softfloat-for-gcc.h"
34 #include "milieu.h"
35 #include "softfloat.h"
36 
37 #include "aeabi_vfp.h"
38 
39 extern int _libc_arm_fpu_present;
40 
41 flag __unorddf2(float64, float64);
42 
43 /* These are written in asm and are only called from this file */
44 int __aeabi_dcmpeq_vfp(float64, float64);
45 int __aeabi_dcmplt_vfp(float64, float64);
46 int __aeabi_dcmple_vfp(float64, float64);
47 int __aeabi_dcmpgt_vfp(float64, float64);
48 int __aeabi_dcmpge_vfp(float64, float64);
49 int __aeabi_dcmpun_vfp(float64, float64);
50 int __aeabi_d2iz_vfp(float64);
51 float32 __aeabi_d2f_vfp(float64);
52 float64 __aeabi_i2d_vfp(int);
53 float64 __aeabi_dadd_vfp(float64, float64);
54 float64 __aeabi_ddiv_vfp(float64, float64);
55 float64 __aeabi_dmul_vfp(float64, float64);
56 float64 __aeabi_dsub_vfp(float64, float64);
57 
58 /*
59  * Depending on the target these will:
60  *  On armv6 with a vfp call the above function, or
61  *  Call the softfloat function in the 3rd argument.
62  */
63 int AEABI_FUNC2(dcmpeq, float64, float64_eq)
64 int AEABI_FUNC2(dcmplt, float64, float64_lt)
65 int AEABI_FUNC2(dcmple, float64, float64_le)
66 int AEABI_FUNC2_REV(dcmpge, float64, float64_le)
67 int AEABI_FUNC2_REV(dcmpgt, float64, float64_lt)
68 int AEABI_FUNC2(dcmpun, float64, __unorddf2)
69 
70 int AEABI_FUNC(d2iz, float64, float64_to_int32_round_to_zero)
71 float32 AEABI_FUNC(d2f, float64, float64_to_float32)
72 float64 AEABI_FUNC(i2d, int, int32_to_float64)
73 
74 float64 AEABI_FUNC2(dadd, float64, float64_add)
75 float64 AEABI_FUNC2(ddiv, float64, float64_div)
76 float64 AEABI_FUNC2(dmul, float64, float64_mul)
77 float64 AEABI_FUNC2(dsub, float64, float64_sub)
78 
79 int
80 __aeabi_cdcmpeq_helper(float64 a, float64 b)
81 {
82 	int quiet = 0;
83 
84 	/* Check if a is a NaN */
85 	if ((a << 1) > 0xffe0000000000000ull) {
86 		/* If it's a signalling NaN we will always signal */
87 		if ((a & 0x0008000000000000ull) == 0)
88 			return (0);
89 
90 		quiet = 1;
91 	}
92 
93 	/* Check if b is a NaN */
94 	if ((b << 1) > 0xffe0000000000000ull) {
95 		/* If it's a signalling NaN we will always signal */
96 		if ((b & 0x0008000000000000ull) == 0)
97 			return (0);
98 
99 		quiet = 1;
100 	}
101 
102 	return (quiet);
103 }
104