1.file "fmax.s"
2
3
4// Copyright (c) 2001 - 2003, Intel Corporation
5// All rights reserved.
6//
7//
8// Redistribution and use in source and binary forms, with or without
9// modification, are permitted provided that the following conditions are
10// met:
11//
12// * Redistributions of source code must retain the above copyright
13// notice, this list of conditions and the following disclaimer.
14//
15// * Redistributions in binary form must reproduce the above copyright
16// notice, this list of conditions and the following disclaimer in the
17// documentation and/or other materials provided with the distribution.
18//
19// * The name of Intel Corporation may not be used to endorse or promote
20// products derived from this software without specific prior written
21// permission.
22
23// THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
24// "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
25// LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
26// A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL INTEL OR ITS
27// CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
28// EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
29// PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
30// PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY
31// OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY OR TORT (INCLUDING
32// NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
33// SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
34//
35// Intel Corporation is the author of this code, and requests that all
36// problem reports or change requests be submitted to it directly at
37// http://www.intel.com/software/products/opensource/libraries/num.htm.
38//
39// History
40//==============================================================
41// 05/31/01 Initial version
42// 05/20/02 Cleaned up namespace and sf0 syntax
43// 02/10/03 Reordered header: .section, .global, .proc, .align
44//
45// API
46//==============================================================
47// double fmax  (double x, double y)
48//
49// Overview of operation
50//==============================================================
51// returns the algebraic maximum of 2 input values
52//
53// Special cases:
54//    fmax(x, nan) returns x if x is numeric  // Must special case this one
55//    fmax(nan, y) returns y if y is numeric
56//    fmax(nan1, nan2) returns quietized nan2
57//    fmax(+0,+0) returns +0
58//    fmax(-0,+0) returns +0
59//    fmax(-0,-0) returns -0
60//    fmax(+0,-0) returns +0                  // Must special case this one
61//
62// SNaN causes invalid to be set
63
64// floating-point registers used: 2
65// f8,  input x, output
66// f9,  input y
67
68.section .text
69GLOBAL_LIBM_ENTRY(fmax)
70
71{ .mfi
72      nop.m 999
73      fcmp.unord.s0 p6,p7 = f8, f9 // Is x or y a nan? Raise invalid or denormal
74      nop.i 999
75}
76;;
77
78{ .mfi
79      nop.m 999
80(p7)  fclass.m.unc  p8,p9 = f9, 0x06 // If no nan, is y=-0?
81      nop.i 999
82}
83;;
84
85{ .mfi
86      nop.m 999
87(p6)  fclass.m.unc  p10,p0 = f8, 0xc3 // Is x nan?
88      nop.i 999
89}
90;;
91
92{ .mfi
93      nop.m 999
94(p9)  fmax.s0  f8 = f8, f9    // Normal case, no nan and y not -0
95      nop.i 999
96}
97;;
98
99{ .mfi
100      nop.m 999
101(p8)  fmax.s0  f8 = f9, f8    // No nan and y -0
102      nop.i 999
103}
104;;
105
106{ .mfb
107      nop.m 999
108(p10) fmerge.s  f8 = f9, f9   // If x nan, return y, else do nothing (returns x)
109      br.ret.sptk    b0
110}
111;;
112
113GLOBAL_LIBM_END(fmax)
114libm_alias_double_other (fmax, fmax)
115