*       $NetBSD: MONADIC.GEN,v 1.3 1994/10/26 07:48:42 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.GEN 1.5 5/18/92
*
*       MONADIC.GEN 1.4 1/16/92
*
*       MONADIC.GEN 1.3 4/30/91
*
*       MONADIC.GEN --- generic MONADIC template
*
*       This version saves all registers that will be used by the emulation
*       routines and restores all but FP0 on exit.  The FPSR is
*       updated to reflect the result of the operation.  Return value
*       is placed in FP0 for single, double and extended results.
*
*       The package subroutines expect the incoming FPCR to be zeroed
*       since they need extended precision to work properly.  The
*       'final' FPCR is expected in d1 so that the calculated result
*       can be properly sized and rounded.  Also, if the incoming FPCR
*       has enabled any exceptions, the exception will be taken on the
*       final fmovem in this template.
*
*       Customizations:
*               1. Remove the movem.l at the entry and exit of
*                  each routine if your compiler treats those
*                  registers as scratch.
*               2. Likewise, don't save FP0/FP1 if they are scratch
*                  registers.
*               3. Delete handling of the fpsr if you only care about
*                  the result.
*               4. Some (most?) C compilers convert all float arguments
*                  to double, and provide no support at all for extended
*                  precision so remove the _OPs_ and _OPx_ entry points.
*               5. Move the result to d0/d1 if the compiler is that old.
*

       xref    tag
       xref    _OPr_
       xref    _OPz_
       xref    _OPi_
       xref    _OPn_
       xref    _OPm_

       xdef    _OPs_
_OPs_:
       link    a6,#-LOCAL_SIZE
       movem.l d0-d1/a0-a1,USER_DA(a6)
       fmovem.x fp0-fp3,USER_FP0(a6)
       fmove.l fpsr,USER_FPSR(a6)
       fmove.l fpcr,USER_FPCR(a6)
       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.s 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.l fpsr,d0         ; update status register
       or.b    FPSR_AEXCEPT(a6),d0     ;add previously accrued exceptions
       fmove.l d0,fpsr
*
*       Result is now in FP0
*
       movem.l USER_DA(a6),d0-d1/a0-a1
       fmovem.x USER_FP1(a6),fp1-fp3   ; note: FP0 not restored
       fmove.l USER_FPCR(a6),fpcr      ; fpcr restored
       unlk    a6
       rts

       xdef    _OPd_
_OPd_:
       link    a6,#-LOCAL_SIZE
       movem.l d0-d1/a0-a1,USER_DA(a6)
       fmovem.x fp0-fp3,USER_FP0(a6)
       fmove.l fpsr,USER_FPSR(a6)
       fmove.l fpcr,USER_FPCR(a6)
       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_7
       bsr     _OPr_           ; normalized (regular) number
       bra.b   _TMP_B
_TMP_7:
       cmp.b   #$20,d0         ; zero?
       bne.b   _TMP_8
       bsr     _OPz_
       bra.b   _TMP_B
_TMP_8:
       cmp.b   #$40,d0         ; infinity?
       bne.b   _TMP_9
       bsr     _OPi_
       bra.b   _TMP_B
_TMP_9:
       cmp.b   #$60,d0         ; NaN?
       bne.b   _TMP_A
       bsr     _OPn_
       bra.b   _TMP_B
_TMP_A:
       bsr     _OPm_           ; assuming a denorm...

_TMP_B:
       fmove.l fpsr,d0         ; update status register
       or.b    FPSR_AEXCEPT(a6),d0     ;add previously accrued exceptions
       fmove.l d0,fpsr
*
*       Result is now in FP0
*
       movem.l USER_DA(a6),d0-d1/a0-a1
       fmovem.x USER_FP1(a6),fp1-fp3   ; note: FP0 not restored
       fmove.l USER_FPCR(a6),fpcr      ; fpcr restored
       unlk    a6
       rts

       xdef    _OPx_
_OPx_:
       link    a6,#-LOCAL_SIZE
       movem.l d0-d1/a0-a1,USER_DA(a6)
       fmovem.x fp0-fp3,USER_FP0(a6)
       fmove.l fpsr,USER_FPSR(a6)
       fmove.l fpcr,USER_FPCR(a6)
       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.x 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_C
       bsr     _OPr_           ; normalized (regular) number
       bra.b   _TMP_G
_TMP_C:
       cmp.b   #$20,d0         ; zero?
       bne.b   _TMP_D
       bsr     _OPz_
       bra.b   _TMP_G
_TMP_D:
       cmp.b   #$40,d0         ; infinity?
       bne.b   _TMP_E
       bsr     _OPi_
       bra.b   _TMP_G
_TMP_E:
       cmp.b   #$60,d0         ; NaN?
       bne.b   _TMP_F
       bsr     _OPn_
       bra.b   _TMP_G
_TMP_F:
       bsr     _OPm_           ; assuming a denorm...

_TMP_G:
       fmove.l fpsr,d0         ; update status register
       or.b    FPSR_AEXCEPT(a6),d0     ;add previously accrued exceptions
       fmove.l d0,fpsr
*
*       Result is now in FP0
*
       movem.l USER_DA(a6),d0-d1/a0-a1
       fmovem.x USER_FP1(a6),fp1-fp3   ; note: FP0 not restored
       fmove.l USER_FPCR(a6),fpcr      ; fpcr restored
       unlk    a6
       rts