Subversion Repositories HelenOS-historic

Rev

Rev 653 | Rev 697 | 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 mantisa */
  34.     return ((f.parts.exp==0xFF)&&(f.parts.mantisa));
  35. };
  36.  
  37. inline int isFloat32SigNaN(float32 f)
  38. {   /* SigNaN : exp = 0xff mantisa = 1xxxxx..x (binary), where at least one x is nonzero */
  39.     return ((f.parts.exp==0xFF)&&(f.parts.mantisa>0x400000));
  40. };
  41.  
  42. inline int isFloat32Infinity(float32 f)
  43. {
  44.     return ((f.parts.exp==0xFF)&&(f.parts.mantisa==0x0));
  45. };
  46.  
  47. /**
  48.  * @return 1, if both floats are equal - but NaNs are not recognized
  49.  */
  50. inline int isFloat32eq(float32 a, float32 b)
  51. {
  52.     return ((a.binary==b.binary)||(((a.binary| b.binary)&0x7FFFFFFF)==0)); /* a equals to b or both are zeros (with any sign) */
  53. }
  54.  
  55. /**
  56.  * @return 1, if a<b - but NaNs are not recognized
  57.  */
  58. inline int isFloat32lt(float32 a, float32 b)
  59. {
  60.     if (((a.binary| b.binary)&0x7FFFFFFF)==0) {
  61.         return 0; /* +- zeroes */
  62.     };
  63.    
  64.     if ((a.parts.sign)&&(b.parts.sign)) {
  65.         /*if both are negative, smaller is that with greater binary value*/
  66.         return (a.binary>b.binary);
  67.         };
  68.    
  69.     /* lets negate signs - now will be positive numbers allways bigger than negative (first bit will be set for unsigned integer comparison)*/
  70.     a.parts.sign=!a.parts.sign;
  71.     b.parts.sign=!b.parts.sign;
  72.     return (a.binary<b.binary);
  73.            
  74. }
  75.  
  76. /**
  77.  * @return 1, if a>b - but NaNs are not recognized
  78.  */
  79. inline int isFloat32gt(float32 a, float32 b)
  80. {
  81.     if (((a.binary| b.binary)&0x7FFFFFFF)==0) {
  82.         return 0; /* zeroes are equal with any sign */
  83.     };
  84.    
  85.     if ((a.parts.sign)&&(b.parts.sign)) {
  86.         /*if both are negative, greater is that with smaller binary value*/
  87.         return (a.binary<b.binary);
  88.         };
  89.    
  90.     /* lets negate signs - now will be positive numbers allways bigger than negative (first bit will be set for unsigned integer comparison)*/
  91.     a.parts.sign=!a.parts.sign;
  92.     b.parts.sign=!b.parts.sign;
  93.     return (a.binary>b.binary);
  94.            
  95. }
  96.  
  97.  
  98.