forked from OSchip/llvm-project
62 lines
2.0 KiB
C
62 lines
2.0 KiB
C
//===-- lib/gesf2.c - Single-precision comparisons ----------------*- C -*-===//
|
|
//
|
|
// The LLVM Compiler Infrastructure
|
|
//
|
|
// This file is dual licensed under the MIT and the University of Illinois Open
|
|
// Source Licenses. See LICENSE.TXT for details.
|
|
//
|
|
//===----------------------------------------------------------------------===//
|
|
//
|
|
// This file implements the following soft-fp_t comparison routines:
|
|
//
|
|
// __eqsf2 __gesf2 __unordsf2
|
|
// __lesf2 __gtsf2
|
|
// __ltsf2
|
|
// __nesf2
|
|
//
|
|
// The semantics of the routines grouped in each column are identical, so there
|
|
// is a single implementation for each, and wrappers to provide the other names.
|
|
//
|
|
// The main routines behave as follows:
|
|
//
|
|
// __lesf2(a,b) returns -1 if a < b
|
|
// 0 if a == b
|
|
// 1 if a > b
|
|
// 1 if either a or b is NaN
|
|
//
|
|
// __gesf2(a,b) returns -1 if a < b
|
|
// 0 if a == b
|
|
// 1 if a > b
|
|
// -1 if either a or b is NaN
|
|
//
|
|
// __unordsf2(a,b) returns 0 if both a and b are numbers
|
|
// 1 if either a or b is NaN
|
|
//
|
|
// Note that __lesf2( ) and __gesf2( ) are identical except in their handling of
|
|
// NaN values.
|
|
//
|
|
//===----------------------------------------------------------------------===//
|
|
|
|
#define SINGLE_PRECISION
|
|
#include "fp_lib.h"
|
|
|
|
enum GE_RESULT __gesf2(fp_t a, fp_t b) {
|
|
|
|
const srep_t aInt = toRep(a);
|
|
const srep_t bInt = toRep(b);
|
|
const rep_t aAbs = aInt & absMask;
|
|
const rep_t bAbs = bInt & absMask;
|
|
|
|
if (aAbs > infRep || bAbs > infRep) return GE_UNORDERED;
|
|
if ((aAbs | bAbs) == 0) return GE_EQUAL;
|
|
if ((aInt & bInt) >= 0) {
|
|
if (aInt < bInt) return GE_LESS;
|
|
else if (aInt == bInt) return GE_EQUAL;
|
|
else return GE_GREATER;
|
|
} else {
|
|
if (aInt > bInt) return GE_LESS;
|
|
else if (aInt == bInt) return GE_EQUAL;
|
|
else return GE_GREATER;
|
|
}
|
|
}
|