Subversion Repositories HelenOS-historic

Rev

Rev 804 | Go to most recent revision | 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. #include<sftypes.h>
  30. #include<comparison.h>
  31.  
  32. inline int isFloat32NaN(float32 f)
  33. {   /* NaN : exp = 0xff and nonzero fraction */
  34.     return ((f.parts.exp==0xFF)&&(f.parts.fraction));
  35. }
  36.  
  37. inline int isFloat64NaN(float64 d)
  38. {   /* NaN : exp = 0x7ff and nonzero fraction */
  39.     return ((d.parts.exp==0x7FF)&&(d.parts.fraction));
  40. }
  41.  
  42. inline int isFloat32SigNaN(float32 f)
  43. {   /* SigNaN : exp = 0xff fraction = 0xxxxx..x (binary), where at least one x is nonzero */
  44.     return ((f.parts.exp==0xFF)&&(f.parts.fraction<0x400000)&&(f.parts.fraction));
  45. }
  46.  
  47. inline int isFloat64SigNaN(float64 d)
  48. {   /* SigNaN : exp = 0x7ff fraction = 0xxxxx..x (binary), where at least one x is nonzero */
  49.     return ((d.parts.exp==0x7FF)&&(d.parts.fraction)&&(d.parts.fraction<0x8000000000000ll));
  50. }
  51.  
  52. inline int isFloat32Infinity(float32 f)
  53. {
  54.     return ((f.parts.exp==0xFF)&&(f.parts.fraction==0x0));
  55. }
  56.  
  57. inline int isFloat64Infinity(float64 d)
  58. {
  59.     return ((d.parts.exp==0x7FF)&&(d.parts.fraction==0x0));
  60. }
  61.  
  62. inline int isFloat32Zero(float32 f)
  63. {
  64.     return (((f.binary) & 0x7FFFFFFF) == 0);
  65. }
  66.  
  67. inline int isFloat64Zero(float64 d)
  68. {
  69.     return (((d.binary) & 0x7FFFFFFFFFFFFFFFll) == 0);
  70. }
  71.  
  72. /**
  73.  * @return 1, if both floats are equal - but NaNs are not recognized
  74.  */
  75. inline int isFloat32eq(float32 a, float32 b)
  76. {
  77.     return ((a.binary==b.binary)||(((a.binary| b.binary)&0x7FFFFFFF)==0)); /* a equals to b or both are zeros (with any sign) */
  78. }
  79.  
  80. /**
  81.  * @return 1, if a<b - but NaNs are not recognized
  82.  */
  83. inline int isFloat32lt(float32 a, float32 b)
  84. {
  85.     if (((a.binary| b.binary)&0x7FFFFFFF)==0) {
  86.         return 0; /* +- zeroes */
  87.     };
  88.    
  89.     if ((a.parts.sign)&&(b.parts.sign)) {
  90.         /*if both are negative, smaller is that with greater binary value*/
  91.         return (a.binary>b.binary);
  92.         };
  93.    
  94.     /* lets negate signs - now will be positive numbers allways bigger than negative (first bit will be set for unsigned integer comparison)*/
  95.     a.parts.sign=!a.parts.sign;
  96.     b.parts.sign=!b.parts.sign;
  97.     return (a.binary<b.binary);
  98.            
  99. }
  100.  
  101. /**
  102.  * @return 1, if a>b - but NaNs are not recognized
  103.  */
  104. inline int isFloat32gt(float32 a, float32 b)
  105. {
  106.     if (((a.binary| b.binary)&0x7FFFFFFF)==0) {
  107.         return 0; /* zeroes are equal with any sign */
  108.     };
  109.    
  110.     if ((a.parts.sign)&&(b.parts.sign)) {
  111.         /*if both are negative, greater is that with smaller binary value*/
  112.         return (a.binary<b.binary);
  113.         };
  114.    
  115.     /* lets negate signs - now will be positive numbers allways bigger than negative (first bit will be set for unsigned integer comparison)*/
  116.     a.parts.sign=!a.parts.sign;
  117.     b.parts.sign=!b.parts.sign;
  118.     return (a.binary>b.binary);
  119.            
  120. }
  121.  
  122.  
  123.  
  124.