Home | History | Annotate | Line # | Download | only in spmath
divsir.c revision 1.1.2.2
      1  1.1.2.2  jdolecek /*	$NetBSD: divsir.c,v 1.1.2.2 2002/06/23 17:37:13 jdolecek Exp $	*/
      2  1.1.2.2  jdolecek 
      3  1.1.2.2  jdolecek /*	$OpenBSD: divsir.c,v 1.5 2001/03/29 03:58:18 mickey Exp $	*/
      4  1.1.2.2  jdolecek 
      5  1.1.2.2  jdolecek /*
      6  1.1.2.2  jdolecek  * Copyright 1996 1995 by Open Software Foundation, Inc.
      7  1.1.2.2  jdolecek  *              All Rights Reserved
      8  1.1.2.2  jdolecek  *
      9  1.1.2.2  jdolecek  * Permission to use, copy, modify, and distribute this software and
     10  1.1.2.2  jdolecek  * its documentation for any purpose and without fee is hereby granted,
     11  1.1.2.2  jdolecek  * provided that the above copyright notice appears in all copies and
     12  1.1.2.2  jdolecek  * that both the copyright notice and this permission notice appear in
     13  1.1.2.2  jdolecek  * supporting documentation.
     14  1.1.2.2  jdolecek  *
     15  1.1.2.2  jdolecek  * OSF DISCLAIMS ALL WARRANTIES WITH REGARD TO THIS SOFTWARE
     16  1.1.2.2  jdolecek  * INCLUDING ALL IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
     17  1.1.2.2  jdolecek  * FOR A PARTICULAR PURPOSE.
     18  1.1.2.2  jdolecek  *
     19  1.1.2.2  jdolecek  * IN NO EVENT SHALL OSF BE LIABLE FOR ANY SPECIAL, INDIRECT, OR
     20  1.1.2.2  jdolecek  * CONSEQUENTIAL DAMAGES OR ANY DAMAGES WHATSOEVER RESULTING FROM
     21  1.1.2.2  jdolecek  * LOSS OF USE, DATA OR PROFITS, WHETHER IN ACTION OF CONTRACT,
     22  1.1.2.2  jdolecek  * NEGLIGENCE, OR OTHER TORTIOUS ACTION, ARISING OUT OF OR IN CONNECTION
     23  1.1.2.2  jdolecek  * WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
     24  1.1.2.2  jdolecek  *
     25  1.1.2.2  jdolecek  */
     26  1.1.2.2  jdolecek /*
     27  1.1.2.2  jdolecek  * pmk1.1
     28  1.1.2.2  jdolecek  */
     29  1.1.2.2  jdolecek /*
     30  1.1.2.2  jdolecek  * (c) Copyright 1986 HEWLETT-PACKARD COMPANY
     31  1.1.2.2  jdolecek  *
     32  1.1.2.2  jdolecek  * To anyone who acknowledges that this file is provided "AS IS"
     33  1.1.2.2  jdolecek  * without any express or implied warranty:
     34  1.1.2.2  jdolecek  *     permission to use, copy, modify, and distribute this file
     35  1.1.2.2  jdolecek  * for any purpose is hereby granted without fee, provided that
     36  1.1.2.2  jdolecek  * the above copyright notice and this notice appears in all
     37  1.1.2.2  jdolecek  * copies, and that the name of Hewlett-Packard Company not be
     38  1.1.2.2  jdolecek  * used in advertising or publicity pertaining to distribution
     39  1.1.2.2  jdolecek  * of the software without specific, written prior permission.
     40  1.1.2.2  jdolecek  * Hewlett-Packard Company makes no representations about the
     41  1.1.2.2  jdolecek  * suitability of this software for any purpose.
     42  1.1.2.2  jdolecek  */
     43  1.1.2.2  jdolecek 
     44  1.1.2.2  jdolecek 
     45  1.1.2.2  jdolecek #include "md.h"
     46  1.1.2.2  jdolecek 
     47  1.1.2.2  jdolecek void
     48  1.1.2.2  jdolecek divsir(opnd1,opnd2,result)
     49  1.1.2.2  jdolecek 
     50  1.1.2.2  jdolecek int opnd1, opnd2;
     51  1.1.2.2  jdolecek struct mdsfu_register *result;
     52  1.1.2.2  jdolecek {
     53  1.1.2.2  jdolecek 	int sign, op1_sign;
     54  1.1.2.2  jdolecek 
     55  1.1.2.2  jdolecek 	/* check divisor for zero */
     56  1.1.2.2  jdolecek 	if (opnd2 == 0) {
     57  1.1.2.2  jdolecek 		overflow = TRUE;
     58  1.1.2.2  jdolecek 		return;
     59  1.1.2.2  jdolecek 	}
     60  1.1.2.2  jdolecek 
     61  1.1.2.2  jdolecek 	/* get sign of result */
     62  1.1.2.2  jdolecek 	sign = opnd1 ^ opnd2;
     63  1.1.2.2  jdolecek 
     64  1.1.2.2  jdolecek 	/* get absolute value of operands */
     65  1.1.2.2  jdolecek 	if (opnd1 < 0) {
     66  1.1.2.2  jdolecek 		opnd1 = -opnd1;
     67  1.1.2.2  jdolecek 		op1_sign = TRUE;
     68  1.1.2.2  jdolecek 	}
     69  1.1.2.2  jdolecek 	else op1_sign = FALSE;
     70  1.1.2.2  jdolecek 	if (opnd2 < 0) opnd2 = -opnd2;
     71  1.1.2.2  jdolecek 
     72  1.1.2.2  jdolecek 	/* check for opnd2 = -2**31 */
     73  1.1.2.2  jdolecek 	if (opnd2 + opnd2 == 0) {
     74  1.1.2.2  jdolecek 		if (opnd1 == opnd2) {
     75  1.1.2.2  jdolecek 			result_hi = 0;    /* remainder = 0 */
     76  1.1.2.2  jdolecek 			result_lo = 1;
     77  1.1.2.2  jdolecek 		}
     78  1.1.2.2  jdolecek 		else {
     79  1.1.2.2  jdolecek 			result_hi = opnd1;  /* remainder = opnd1 */
     80  1.1.2.2  jdolecek 			result_lo = 0;
     81  1.1.2.2  jdolecek 		}
     82  1.1.2.2  jdolecek 	}
     83  1.1.2.2  jdolecek 	else {
     84  1.1.2.2  jdolecek 		/* do the divide */
     85  1.1.2.2  jdolecek 		divu(0,opnd1,opnd2,result);
     86  1.1.2.2  jdolecek 
     87  1.1.2.2  jdolecek 		/*
     88  1.1.2.2  jdolecek 		 * check for overflow
     89  1.1.2.2  jdolecek 		 *
     90  1.1.2.2  jdolecek 		 * at this point, the only way we can get overflow
     91  1.1.2.2  jdolecek 		 * is with opnd1 = -2**31 and opnd2 = -1
     92  1.1.2.2  jdolecek 		 */
     93  1.1.2.2  jdolecek 		if (sign>0 && result_lo<0) {
     94  1.1.2.2  jdolecek 			overflow = TRUE;
     95  1.1.2.2  jdolecek 			return;
     96  1.1.2.2  jdolecek 		}
     97  1.1.2.2  jdolecek 	}
     98  1.1.2.2  jdolecek 	overflow = FALSE;
     99  1.1.2.2  jdolecek 
    100  1.1.2.2  jdolecek 	/* return appropriately signed remainder and result */
    101  1.1.2.2  jdolecek 	if (op1_sign) result_hi = -result_hi;
    102  1.1.2.2  jdolecek 	if (sign<0) result_lo = -result_lo;
    103  1.1.2.2  jdolecek 	return;
    104  1.1.2.2  jdolecek }
    105