xref: /freebsd/lib/libc/arm/aeabi/aeabi_double.c (revision 6f1af0d7d2af54b339b5212434cd6d4fda628d80)
1 /*-
2  * SPDX-License-Identifier: BSD-2-Clause
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 #include "softfloat-for-gcc.h"
32 #include "milieu.h"
33 #include "softfloat.h"
34 
35 #include "aeabi_vfp.h"
36 
37 extern int _libc_arm_fpu_present;
38 
39 flag __unorddf2(float64, float64);
40 
41 /* These are written in asm and are only called from this file */
42 int __aeabi_dcmpeq_vfp(float64, float64);
43 int __aeabi_dcmplt_vfp(float64, float64);
44 int __aeabi_dcmple_vfp(float64, float64);
45 int __aeabi_dcmpgt_vfp(float64, float64);
46 int __aeabi_dcmpge_vfp(float64, float64);
47 int __aeabi_dcmpun_vfp(float64, float64);
48 int __aeabi_d2iz_vfp(float64);
49 float32 __aeabi_d2f_vfp(float64);
50 float64 __aeabi_i2d_vfp(int);
51 float64 __aeabi_dadd_vfp(float64, float64);
52 float64 __aeabi_ddiv_vfp(float64, float64);
53 float64 __aeabi_dmul_vfp(float64, float64);
54 float64 __aeabi_dsub_vfp(float64, float64);
55 
56 /*
57  * Depending on the target these will:
58  *  On armv6 with a vfp call the above function, or
59  *  Call the softfloat function in the 3rd argument.
60  */
61 int AEABI_FUNC2(dcmpeq, float64, float64_eq)
62 int AEABI_FUNC2(dcmplt, float64, float64_lt)
63 int AEABI_FUNC2(dcmple, float64, float64_le)
64 int AEABI_FUNC2_REV(dcmpge, float64, float64_le)
65 int AEABI_FUNC2_REV(dcmpgt, float64, float64_lt)
66 int AEABI_FUNC2(dcmpun, float64, __unorddf2)
67 
68 int AEABI_FUNC(d2iz, float64, float64_to_int32_round_to_zero)
69 float32 AEABI_FUNC(d2f, float64, float64_to_float32)
70 float64 AEABI_FUNC(i2d, int, int32_to_float64)
71 
72 float64 AEABI_FUNC2(dadd, float64, float64_add)
73 float64 AEABI_FUNC2(ddiv, float64, float64_div)
74 float64 AEABI_FUNC2(dmul, float64, float64_mul)
75 float64 AEABI_FUNC2(dsub, float64, float64_sub)
76 
77 int
78 __aeabi_cdcmpeq_helper(float64 a, float64 b)
79 {
80 	int quiet = 0;
81 
82 	/* Check if a is a NaN */
83 	if ((a << 1) > 0xffe0000000000000ull) {
84 		/* If it's a signalling NaN we will always signal */
85 		if ((a & 0x0008000000000000ull) == 0)
86 			return (0);
87 
88 		quiet = 1;
89 	}
90 
91 	/* Check if b is a NaN */
92 	if ((b << 1) > 0xffe0000000000000ull) {
93 		/* If it's a signalling NaN we will always signal */
94 		if ((b & 0x0008000000000000ull) == 0)
95 			return (0);
96 
97 		quiet = 1;
98 	}
99 
100 	return (quiet);
101 }
102