1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
|
* $NetBSD: MONADIC.R3V6,v 1.2 1994/10/26 07:48:44 cgd Exp $
* MOTOROLA MICROPROCESSOR & MEMORY TECHNOLOGY GROUP
* M68000 Hi-Performance Microprocessor Division
* M68040 Software Package
*
* M68040 Software Package Copyright (c) 1993, 1994 Motorola Inc.
* All rights reserved.
*
* THE SOFTWARE is provided on an "AS IS" basis and without warranty.
* To the maximum extent permitted by applicable law,
* MOTOROLA DISCLAIMS ALL WARRANTIES WHETHER EXPRESS OR IMPLIED,
* INCLUDING IMPLIED WARRANTIES OF MERCHANTABILITY OR FITNESS FOR A
* PARTICULAR PURPOSE and any warranty against infringement with
* regard to the SOFTWARE (INCLUDING ANY MODIFIED VERSIONS THEREOF)
* and any accompanying written materials.
*
* To the maximum extent permitted by applicable law,
* IN NO EVENT SHALL MOTOROLA BE LIABLE FOR ANY DAMAGES WHATSOEVER
* (INCLUDING WITHOUT LIMITATION, DAMAGES FOR LOSS OF BUSINESS
* PROFITS, BUSINESS INTERRUPTION, LOSS OF BUSINESS INFORMATION, OR
* OTHER PECUNIARY LOSS) ARISING OF THE USE OR INABILITY TO USE THE
* SOFTWARE. Motorola assumes no responsibility for the maintenance
* and support of the SOFTWARE.
*
* You are hereby granted a copyright license to use, modify, and
* distribute the SOFTWARE so long as this entire notice is retained
* without alteration in any modified and/or redistributed versions,
* and that such modified versions are clearly identified as such.
* No licenses are granted by implication, estoppel or otherwise
* under any patents or trademarks of Motorola, Inc.
*
* MONADIC.R3V6 1.3 4/30/91
*
* MONADIC.R3V6 --- MONADIC template for MCD R3V6 native C compiler
*
* The MCD compiler is old. It returns float and double values
* as a double stored in d0/d1. There is no support for single or extended
* precision operations. It's not clear whether the float registers
* should be preserved, so for speed, they're not.
*
xref tag
xref _OPr_
xref _OPz_
xref _OPi_
xref _OPn_
xref _OPm_
xdef _OPd_
_OPd_:
link a6,#-LOCAL_SIZE
fmove.l fpcr,d1 ; user's rounding mode/precision
fmove.l #0,fpcr ; force rounding mode/prec to extended,rn
*
* copy, convert and tag input argument
*
fmove.d 8(a6),fp0
fmove.x fp0,ETEMP(a6)
lea ETEMP(a6),a0
bsr tag
move.b d0,STAG(a6)
tst.b d0
bne.b _TMP_2
bsr _OPr_ ; normalized (regular) number
bra.b _TMP_6
_TMP_2:
cmp.b #$20,d0 ; zero?
bne.b _TMP_3
bsr _OPz_
bra.b _TMP_6
_TMP_3:
cmp.b #$40,d0 ; infinity?
bne.b _TMP_4
bsr _OPi_
bra.b _TMP_6
_TMP_4:
cmp.b #$60,d0 ; NaN?
bne.b _TMP_5
bsr _OPn_
bra.b _TMP_6
_TMP_5:
bsr _OPm_ ; assuming a denorm...
_TMP_6:
fmove.d fp0,USER_D0(a6) ; result goes into d0/d1 pair
movem.l USER_D0(a6),d0-d1
unlk a6
rts
|
