Subversion Repositories HelenOS-historic

Rev

Rev 652 | Rev 688 | 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==b)||(((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;
  62.     };
  63.     a.parts.sign^=a.parts.sign;
  64.     b.parts.sign^=b.parts.sign;
  65.     return (a.binary<b.binary);
  66.            
  67. }
  68.  
  69. /**
  70.  * @return 1, if a>b - but NaNs are not recognized
  71.  */
  72. inline int isFloat32gt(float32 a, float32 b)
  73. {
  74.     if (((a.binary| b.binary)&0x7FFFFFFF)==0) {
  75.         return 0;
  76.     };
  77.     a.parts.sign^=a.parts.sign;
  78.     b.parts.sign^=b.parts.sign;
  79.     return (a.binary>b.binary);
  80.            
  81. }
  82.  
  83.  
  84.