lib/libc/arm/aeabi/aeabi_double.c

   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 }