|
|
1.1 root 1: /*
2: * Copyright (c) 1986 Regents of the University of California.
3: * All rights reserved.
4: *
5: * This code is derived from software contributed to Berkeley by
6: * Computer Consoles Inc.
7: *
8: * Redistribution and use in source and binary forms are permitted
9: * provided that: (1) source distributions retain this entire copyright
10: * notice and comment, and (2) distributions including binaries display
11: * the following acknowledgement: ``This product includes software
12: * developed by the University of California, Berkeley and its contributors''
13: * in the documentation or other materials provided with the distribution
14: * and in all advertising materials mentioning features or use of this
15: * software. Neither the name of the University nor the names of its
16: * contributors may be used to endorse or promote products derived
17: * from this software without specific prior written permission.
18: * THIS SOFTWARE IS PROVIDED ``AS IS'' AND WITHOUT ANY EXPRESS OR
19: * IMPLIED WARRANTIES, INCLUDING, WITHOUT LIMITATION, THE IMPLIED
20: * WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE.
21: */
22:
23: #if defined(SYSLIBC_SCCS) && !defined(lint)
24: .asciz "@(#)sfnorm.s 1.3 (Berkeley) 6/1/90"
25: #endif /* SYSLIBC_SCCS and not lint */
26:
27: #include <tahoemath/fp.h>
28: #include "DEFS.h"
29:
30: XENTRY(sfnorm, R2|R3|R4|R5|R6)
31: clrl r1
32: movl r0,r4 /* copy to temporary. */
33: jeql retzero
34: clrl r3 /* r3 - pos of m.s.b */
35: inr00: ffs r4,r6
36: incl r6
37: addl2 r6,r3
38: shrl r6,r4,r4
39: jneq inr00
40:
41: cmpshift:
42: /* compute the shift (r4). */
43: subl3 r3,$HID_POS,r4
44: jlss shiftr /* if less then zero we shift right. */
45: shll r4,r0,r0 /* else we shift left. */
46: subl2 r4,r2 /* uodate exponent. */
47: jleq underflow /* if less then 0 (biased) it is underflow. */
48: jmp combine /* go to combine exponent and fraction. */
49: shiftr:
50: mnegl r4,r4
51: shrl r4,r0,r0 /* shift right. */
52: addl2 r4,r2 /* update exponent */
53: cmpl r2,$256
54: jgeq overflow /* check for overflow. */
55: combine:
56: andl2 $CLEARHID,r0 /* clear the hidden bit. */
57: shal $EXPSHIFT,r2,r2 /* shift the exponent to its proper place. */
58: orl2 r2,r0
59: ret
60:
61: underflow:
62: callf $4,sfpunder
63: ret
64:
65: overflow:
66: callf $4,sfpover
67: ret
68: retzero:
69: clrl r0
70: ret
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.