* Copyright (c) 1986 Regents of the University of California.
* This code is derived from software contributed to Berkeley by
* Redistribution and use in source and binary forms, with or without
* modification, are permitted provided that the following conditions
* 1. Redistributions of source code must retain the above copyright
* notice, this list of conditions and the following disclaimer.
* 2. Redistributions in binary form must reproduce the above copyright
* notice, this list of conditions and the following disclaimer in the
* documentation and/or other materials provided with the distribution.
* 3. All advertising materials mentioning features or use of this software
* must display the following acknowledgement:
* This product includes software developed by the University of
* California, Berkeley and its contributors.
* 4. Neither the name of the University nor the names of its contributors
* may be used to endorse or promote products derived from this software
* without specific prior written permission.
* THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
* ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
* IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
* ARE DISCLAIMED. IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
* FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
* DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
* OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
* HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
* LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
* OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
#if defined(SYSLIBC_SCCS) && !defined(lint)
.asciz "@(#)addf.s 1.3 (Berkeley) 6/1/90"
#endif /* SYSLIBC_SCCS and not lint */
#include <tahoemath//fp.h>
XENTRY(addf, R2|R3|R4|R5|R6|R7|R8|R9|R10)
* see which operand has a greater exponent
* The greater one will be fetched into r0,r2,r3.
* r0- 'pure' fraction, r2 - exponent, r3 - sign).
* The smaller operand will be fetched into r4,r6,r7.
movl 12(fp),r0 # bigger operand to r0
movl 4(fp),r4 # smaller operand to r4
movl 4(fp),r0 # bigger operand to r0
movl 12(fp),r4 # smaller operand to r4
andl3 $EXPMASK,r0,r2 # r2 will hold the exponent of greater operand.
jeql is_res1 # check for reserved operand.
andl3 $EXPMASK,r4,r6 # r6 will hold the exponent of smaller operand.
jeql is_res2 # check for reserved operand.
ret # return the bigger number.
bbc $31,r0,sign2 # if negative remember it.(R3=1)
bbc $31,r4,frac # if negative remember it.(R7=1)
*compute 'pure' fraction:
# clear the non fraction parts.
andl2 $(0!(EXPMASK | SIGNBIT)),r0
# clear the non fraction parts.
andl2 $(0!(EXPMASK | SIGNBIT)),r4
*shift the smaller operand:
bbs $0,r3,negative # the bigger operand was negative.
movl r4,r0 # return the smaller operand.
ret # we allready have the 'result' in r0,r1.