source: trunk/minix/lib/float/shifter.c@ 12

/*
  (c) copyright 1988 by the Vrije Universiteit, Amsterdam, The Netherlands.
  See the copyright notice in the ACK home directory, in the file "Copyright".
*/

/* $Header: /cvsup/minix/src/lib/float/shifter.c,v 1.1.1.1 2005/04/21 14:56:11 beng Exp $ */

# include "FP_types.h"

void
b64_sft(e1,n)
B64     *e1;
int     n;
{
        if (n > 0) {
                if (n > 63) {
                        e1->l_32 = 0;
                        e1->h_32 = 0;
                        return;
                }
                if (n >= 32) {
                        e1->l_32 = e1->h_32;
                        e1->h_32 = 0;
                        n -= 32;
                }
                if (n > 0) {
                        e1->l_32 >>= n;
                        if (e1->h_32 != 0) {
                                e1->l_32 |= (e1->h_32 << (32 - n));
                                e1->h_32 >>= n;
                        }
                }
                return;
        }
        n = -n;
        if (n > 0) {
                if (n > 63) {
                        e1->l_32 = 0;
                        e1->h_32 = 0;
                        return;
                }
                if (n >= 32) {
                        e1->h_32 = e1->l_32;
                        e1->l_32 = 0;
                        n -= 32;
                }
                if (n > 0) {
                        e1->h_32 <<= n;
                        if (e1->l_32 != 0) {
                                e1->h_32 |= (e1->l_32 >> (32 - n));
                                e1->l_32 <<= n;
                        }
                }
        }
}

void
b64_lsft(e1)
B64     *e1;
{
        /*      shift left 1 bit */
        e1->h_32 <<= 1;
        if (e1->l_32 & 0x80000000L) e1->h_32 |= 1;
        e1->l_32 <<= 1;
}

void
b64_rsft(e1)
B64     *e1;
{
        /*      shift right 1 bit */
        e1->l_32 >>= 1;
        if (e1->h_32 & 1) e1->l_32 |= 0x80000000L;
        e1->h_32 >>= 1;
}