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