home *** CD-ROM | disk | FTP | other *** search
/ InfoMagic Source Code 1993 July / THE_SOURCE_CODE_CD_ROM.iso / bsd_srcs / lib / libc / tahoe / gen / modf.s < prev    next >
Encoding:
Text File  |  1991-04-12  |  3.1 KB  |  89 lines
  1. /*
  2.  * Copyright (c) 1988 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, with or without
  9.  * modification, are permitted provided that the following conditions
  10.  * are met:
  11.  * 1. Redistributions of source code must retain the above copyright
  12.  *    notice, this list of conditions and the following disclaimer.
  13.  * 2. Redistributions in binary form must reproduce the above copyright
  14.  *    notice, this list of conditions and the following disclaimer in the
  15.  *    documentation and/or other materials provided with the distribution.
  16.  * 3. All advertising materials mentioning features or use of this software
  17.  *    must display the following acknowledgement:
  18.  *    This product includes software developed by the University of
  19.  *    California, Berkeley and its contributors.
  20.  * 4. Neither the name of the University nor the names of its contributors
  21.  *    may be used to endorse or promote products derived from this software
  22.  *    without specific prior written permission.
  23.  *
  24.  * THIS SOFTWARE IS PROVIDED BY THE REGENTS AND CONTRIBUTORS ``AS IS'' AND
  25.  * ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
  26.  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
  27.  * ARE DISCLAIMED.  IN NO EVENT SHALL THE REGENTS OR CONTRIBUTORS BE LIABLE
  28.  * FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
  29.  * DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS
  30.  * OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
  31.  * HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
  32.  * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY
  33.  * OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF
  34.  * SUCH DAMAGE.
  35.  */
  36.  
  37. #if defined(LIBC_SCCS) && !defined(lint)
  38.     .asciz "@(#)modf.s    5.1 (Berkeley) 1/30/91"
  39. #endif /* LIBC_SCCS and not lint */
  40.  
  41. /*
  42.  * double modf (value, iptr)
  43.  * double value, *iptr;
  44.  *
  45.  * Modf returns the fractional part of "value",
  46.  * and stores the integer part indirectly through "iptr".
  47.  *
  48.  * This version uses floating point (look in ../fpe for
  49.  * a much slower integer version).
  50.  */
  51.  
  52. #include "DEFS.h"
  53.  
  54. ENTRY(modf, 0)
  55.     ldd    4(fp)        # value
  56.     cvdl    r2        # integerize
  57.     bvs    1f        # did integer part overflow?
  58.     cvld    r2        # integer part
  59.     std    r0
  60.     std    *12(fp)        # *iptr = r2
  61.     ldd    4(fp)
  62.     subd    r0        # value-(int)value
  63.     std    r0        # return fraction
  64.     ret
  65. 1:
  66.     /*
  67.      * If the integer portion overflowed, mask out the fractional
  68.      * bits in the double word instead of cvdl-ing.
  69.      */
  70.     ldd    4(fp)
  71.     std    r0        # (r0,r1) = value
  72.     shrl    $23,r0,r2    # extract sign,exponent of value
  73.     andl2    $255,r2        # exponent
  74.     subl2    $152,r2        # e-152
  75.     /*
  76.      * If it overflowed then value>=2^31 and e>=160
  77.      * so we mask only r1 (low bits of fraction), not r0
  78.      */
  79.     mnegl    $1,r3
  80.     shrl    r2,r3,r3    # -1>>(e-152) is neg mask to clear fraction
  81.     mcoml    r3,r3        # complement mask
  82.     andl2    r3,r1        # mask off truly fractional bits from fraction
  83.     ldd    r0        # now (r0,r1) = integerized value
  84.     std    *12(fp)        # *iptr = integerized
  85.     ldd    4(fp)
  86.     subd    r0
  87.     std    r0        # return fraction
  88.     ret
  89.