Subversion Repositories HelenOS

Rev

Rev 2479 | Blame | Compare with Previous | Last modification | View Log | Download | RSS feed

  1. /*
  2.  * Copyright (c) 2005 Josef Cejka
  3.  * All rights reserved.
  4.  *
  5.  * Redistribution and use in source and binary forms, with or without
  6.  * modification, are permitted provided that the following conditions
  7.  * are met:
  8.  *
  9.  * - Redistributions of source code must retain the above copyright
  10.  *   notice, this list of conditions and the following disclaimer.
  11.  * - Redistributions in binary form must reproduce the above copyright
  12.  *   notice, this list of conditions and the following disclaimer in the
  13.  *   documentation and/or other materials provided with the distribution.
  14.  * - The name of the author may not be used to endorse or promote products
  15.  *   derived from this software without specific prior written permission.
  16.  *
  17.  * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR
  18.  * IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
  19.  * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED.
  20.  * IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, INDIRECT,
  21.  * INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT
  22.  * NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
  23.  * DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
  24.  * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
  25.  * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF
  26.  * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
  27.  */
  28.  
  29. /** @addtogroup softfloat  
  30.  * @{
  31.  */
  32. /** @file
  33.  */
  34.  
  35. #include<sftypes.h>
  36. #include<comparison.h>
  37.  
  38. inline int isFloat32NaN(float32 f)
  39. {   /* NaN : exp = 0xff and nonzero fraction */
  40.     return ((f.parts.exp==0xFF)&&(f.parts.fraction));
  41. }
  42.  
  43. inline int isFloat64NaN(float64 d)
  44. {   /* NaN : exp = 0x7ff and nonzero fraction */
  45.     return ((d.parts.exp==0x7FF)&&(d.parts.fraction));
  46. }
  47.  
  48. inline int isFloat32SigNaN(float32 f)
  49. {   /* SigNaN : exp = 0xff fraction = 0xxxxx..x (binary), where at least one x is nonzero */
  50.     return ((f.parts.exp==0xFF)&&(f.parts.fraction<0x400000)&&(f.parts.fraction));
  51. }
  52.  
  53. inline int isFloat64SigNaN(float64 d)
  54. {   /* SigNaN : exp = 0x7ff fraction = 0xxxxx..x (binary), where at least one x is nonzero */
  55.     return ((d.parts.exp==0x7FF)&&(d.parts.fraction)&&(d.parts.fraction<0x8000000000000ll));
  56. }
  57.  
  58. inline int isFloat32Infinity(float32 f)
  59. {
  60.     return ((f.parts.exp==0xFF)&&(f.parts.fraction==0x0));
  61. }
  62.  
  63. inline int isFloat64Infinity(float64 d)
  64. {
  65.     return ((d.parts.exp==0x7FF)&&(d.parts.fraction==0x0));
  66. }
  67.  
  68. inline int isFloat32Zero(float32 f)
  69. {
  70.     return (((f.binary) & 0x7FFFFFFF) == 0);
  71. }
  72.  
  73. inline int isFloat64Zero(float64 d)
  74. {
  75.     return (((d.binary) & 0x7FFFFFFFFFFFFFFFll) == 0);
  76. }
  77.  
  78. /**
  79.  * @return 1, if both floats are equal - but NaNs are not recognized
  80.  */
  81. inline int isFloat32eq(float32 a, float32 b)
  82. {
  83.     return ((a.binary==b.binary)||(((a.binary| b.binary)&0x7FFFFFFF)==0)); /* a equals to b or both are zeros (with any sign) */
  84. }
  85.  
  86. /**
  87.  * @return 1, if a<b - but NaNs are not recognized
  88.  */
  89. inline int isFloat32lt(float32 a, float32 b)
  90. {
  91.     if (((a.binary| b.binary)&0x7FFFFFFF)==0) {
  92.         return 0; /* +- zeroes */
  93.     };
  94.    
  95.     if ((a.parts.sign)&&(b.parts.sign)) {
  96.         /*if both are negative, smaller is that with greater binary value*/
  97.         return (a.binary>b.binary);
  98.         };
  99.    
  100.     /* lets negate signs - now will be positive numbers allways bigger than negative (first bit will be set for unsigned integer comparison)*/
  101.     a.parts.sign=!a.parts.sign;
  102.     b.parts.sign=!b.parts.sign;
  103.     return (a.binary<b.binary);
  104.            
  105. }
  106.  
  107. /**
  108.  * @return 1, if a>b - but NaNs are not recognized
  109.  */
  110. inline int isFloat32gt(float32 a, float32 b)
  111. {
  112.     if (((a.binary| b.binary)&0x7FFFFFFF)==0) {
  113.         return 0; /* zeroes are equal with any sign */
  114.     };
  115.    
  116.     if ((a.parts.sign)&&(b.parts.sign)) {
  117.         /*if both are negative, greater is that with smaller binary value*/
  118.         return (a.binary<b.binary);
  119.         };
  120.    
  121.     /* lets negate signs - now will be positive numbers allways bigger than negative (first bit will be set for unsigned integer comparison)*/
  122.     a.parts.sign=!a.parts.sign;
  123.     b.parts.sign=!b.parts.sign;
  124.     return (a.binary>b.binary);
  125.            
  126. }
  127.  
  128. /** @}
  129.  */
  130.