usr/src/lib/libc/port/gen/ldivide.c

   1 /*
   2  * CDDL HEADER START
   3  *
   4  * The contents of this file are subject to the terms of the
   5  * Common Development and Distribution License (the "License").
   6  * You may not use this file except in compliance with the License.
   7  *
   8  * You can obtain a copy of the license at usr/src/OPENSOLARIS.LICENSE
   9  * or http://www.opensolaris.org/os/licensing.
  10  * See the License for the specific language governing permissions
  11  * and limitations under the License.
  12  *
  13  * When distributing Covered Code, include this CDDL HEADER in each
  14  * file and include the License file at usr/src/OPENSOLARIS.LICENSE.
  15  * If applicable, add the following below this CDDL HEADER, with the
  16  * fields enclosed by brackets "[]" replaced with your own identifying
  17  * information: Portions Copyright [yyyy] [name of copyright owner]
  18  *
  19  * CDDL HEADER END
  20  */
  21
  22 /*
  23  * Copyright 2008 Sun Microsystems, Inc.  All rights reserved.
  24  * Use is subject to license terms.
  25  */
  26
  27 /*      Copyright (c) 1988 AT&T */
  28 /*        All Rights Reserved   */
  29
  30 #pragma ident   "%Z%%M% %I%     %E% SMI"
  31
  32 #include        "lint.h"
  33 #include        <sys/types.h>
  34 #include        <sys/dl.h>
  35 #include        "libc.h"
  36
  37 dl_t
  38 ldivide(dl_t lop, dl_t rop)
  39 {
  40         int             cnt;
  41         dl_t            ans;
  42         dl_t            tmp;
  43         dl_t            div;
  44
  45         if (lsign(lop))
  46                 lop = lsub(lzero, lop);
  47         if (lsign(rop))
  48                 rop = lsub(lzero, rop);
  49
  50         ans = lzero;
  51         div = lzero;
  52
  53         for (cnt = 0; cnt < 63; cnt++) {
  54                 div = lshiftl(div, 1);
  55                 lop = lshiftl(lop, 1);
  56                 if (lsign(lop))
  57                         div.dl_lop |= 1;
  58                 tmp = lsub(div, rop);
  59                 ans = lshiftl(ans, 1);
  60                 if (lsign(tmp) == 0) {
  61                         ans.dl_lop |= 1;
  62                         div = tmp;
  63                 }
  64         }
  65
  66         return (ans);
  67 }