home *** CD-ROM | disk | FTP | other *** search
/ InfoMagic Source Code 1993 July / THE_SOURCE_CODE_CD_ROM.iso / bsd_srcs / lib / libc / vax / gen / ldexp.s < prev    next >
Encoding:
Text File  |  1991-04-12  |  3.2 KB  |  90 lines
  1. /*
  2.  * Copyright (c) 1983 Regents of the University of California.
  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
  7.  * are met:
  8.  * 1. Redistributions of source code must retain the above copyright
  9.  *    notice, this list of conditions and the following disclaimer.
  10.  * 2. Redistributions in binary form must reproduce the above copyright
  11.  *    notice, this list of conditions and the following disclaimer in the
  12.  *    documentation and/or other materials provided with the distribution.
  13.  * 3. All advertising materials mentioning features or use of this software
  14.  *    must display the following acknowledgement:
  15.  *    This product includes software developed by the University of
  16.  *    California, Berkeley and its contributors.
  17.  * 4. Neither the name of the University nor the names of its contributors
  18.  *    may be used to endorse or promote products derived from this software
  19.  *    without specific prior written permission.
  20.  *
  21.  * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
  22.  * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  23.  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  24.  * ARE DISCLAIMED.  IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
  25.  * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
  26.  * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
  27.  * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
  28.  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
  29.  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
  30.  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  31.  * SUCH DAMAGE.
  32.  */
  33.  
  34. #if defined(LIBC_SCCS) && !defined(lint)
  35.     .asciz "@(#)ldexp.s    5.6 (Berkeley) 6/1/90"
  36. #endif /* LIBC_SCCS and not lint */
  37.  
  38. /*
  39.  * double ldexp (value, exp)
  40.  * double value;
  41.  * int exp;
  42.  *
  43.  * Ldexp returns value*2**exp, if that result is in range.
  44.  * If underflow occurs, it returns zero.  If overflow occurs,
  45.  * it returns a value of appropriate sign and largest
  46.  * possible magnitude.  In case of either overflow or underflow,
  47.  * errno is set to ERANGE.  Note that errno is not modified if
  48.  * no error occurs.
  49.  */
  50.  
  51. #include "DEFS.h"
  52.  
  53. /*
  54.  * don't include errno.h, ANSI C says it defines errno.
  55.  *
  56.  * #include <errno.h>
  57.  */
  58. #define    ERANGE    34
  59.  
  60.     .globl    _errno
  61.  
  62. ENTRY(ldexp, 0)
  63.     movd    4(ap),r0    /* fetch "value" */
  64.     extzv    $7,$8,r0,r2    /* r2 := biased exponent */
  65.     jeql    1f        /* if zero, done */
  66.  
  67.     addl2    12(ap),r2    /* r2 := new biased exponent */
  68.     jleq    2f        /* if <= 0, underflow */
  69.     cmpl    r2,$256        /* otherwise check if too big */
  70.     jgeq    3f        /* jump if overflow */
  71.     insv    r2,$7,$8,r0    /* put exponent back in result */
  72. 1:
  73.     ret
  74. 2:
  75.     clrd    r0
  76.     jbr    1f
  77. 3:
  78.     movd    huge,r0        /* largest possible floating magnitude */
  79.     jbc    $15,4(ap),1f    /* jump if argument was positive */
  80.     mnegd    r0,r0        /* if arg < 0, make result negative */
  81. 1:
  82.     movl    $ERANGE,_errno
  83.     ret
  84.  
  85.     .data
  86. huge:    .word    0x7fff        /* the largest number that can */
  87.     .word    0xffff        /*   be represented in a long floating */
  88.     .word    0xffff        /*   number.  This is given in hex in order */
  89.     .word    0xffff        /*   to avoid floating conversions */
  90.