← Back to branch summary

~ubuntu-branches/ubuntu/raring/avr-libc/raring-proposed

« back to all changes in this revision

Viewing changes to libm/fplib/fmax.S

  • Committer: Bazaar Package Importer
  • Author(s): Hakan Ardo
  • Date: 2008-04-04 17:05:32 UTC
  • mto: This revision was merged to the branch mainline in revision 8.
  • Revision ID: james.westby@ubuntu.com-20080404170532-t84huz4qk2928pt3
Tags: upstream-1.6.2
Import upstream version 1.6.2

Show diffs side-by-side

added added

removed removed

Lines of Context:
libm/fplib/fmax.S
 
1
/* Copyright (c) 2002  Michael Stumpf  <mistumpf@de.pepperl-fuchs.com>
 
2
   Copyright (c) 2006  Dmitry Xmelkov
 
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 are met:
 
7
 
 
8
   * Redistributions of source code must retain the above copyright
 
9
     notice, this list of conditions and the following disclaimer.
 
10
   * Redistributions in binary form must reproduce the above copyright
 
11
     notice, this list of conditions and the following disclaimer in
 
12
     the documentation and/or other materials provided with the
 
13
     distribution.
 
14
   * Neither the name of the copyright holders nor the names of
 
15
     contributors may be used to endorse or promote products derived
 
16
     from this software without specific prior written permission.
 
17
 
 
18
   THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
 
19
   AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
 
20
   IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
 
21
   ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
 
22
   LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
 
23
   CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
 
24
   SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
 
25
   INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
 
26
   CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
 
27
   ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
 
28
   POSSIBILITY OF SUCH DAMAGE. */
 
29
 
 
30
/* $Id: fmax.S,v 1.1 2007/01/14 14:51:04 dmix Exp $ */
 
31
 
 
32
#include "fp32def.h"
 
33
#include "asmdef.h"
 
34
 
 
35
/* float fmax (float x, float y);
 
36
     The fmax() function returns the greater of the two values x and y.
 
37
     If an argument is NaN, the other argument is returned. If both
 
38
     arguments are NaN, NaN is returned.
 
39
 */
 
40
 
 
41
ENTRY fmax
 
42
        lsl     rA3
 
43
        sbc     rAE, rAE        ; rAE = (A < 0) ? -1 : 0
 
44
        lsl     rB3
 
45
        sbc     rBE, rBE        ; rBE = (B < 0) ? -1 : 0
 
46
  ; isnan(A) ?
 
47
        ldi     ZL, 0x80        ; NaN: 0x{f/7}f800001..0x{f/7}fffffff
 
48
        ldi     ZH, 0xfe
 
49
        cp      r1, rA0
 
50
        cpc     r1, rA1
 
51
        cpc     ZL, rA2
 
52
        cpc     ZH, rA3
 
53
        brlo    .L_retB         ; A is NaN
 
54
  ; isnan(B) ?
 
55
        cp      r1, rB0
 
56
        cpc     r1, rB1
 
57
        cpc     ZL, rB2
 
58
        cpc     ZH, rB3
 
59
        brlo    .L_retA
 
60
  ; compare signs
 
61
        cp      rAE, rBE
 
62
        brlt    .L_retB         ; A < 0, B >= 0
 
63
        brne    .L_retA         ; A >= 0, B < 0
 
64
  ; compare absolute values
 
65
        cp      rA0, rB0
 
66
        cpc     rA1, rB1
 
67
        cpc     rA2, rB2
 
68
        cpc     rA3, rB3        ; C is set, if A < B
 
69
  ; analize
 
70
        ror     rBE             ; V = C ^ sign
 
71
        brvc    .L_retA
 
72
.L_retB:
 
73
        X_movw  rA0, rB0
 
74
        X_movw  rA2, rB2
 
75
        mov     rAE, rBE
 
76
.L_retA:
 
77
        lsr     rAE
 
78
        ror     rA3
 
79
        ret
 
80
ENDFUNC