modules/src/flt_arith/b64_sft.c

   1 /*
   2   (c) copyright 1988 by the Vrije Universiteit, Amsterdam, The Netherlands.
   3   See the copyright notice in the ACK home directory, in the file "Copyright".
   4 */
   5
   6 /* $Id: b64_sft.c,v 1.7 1994/06/24 11:15:30 ceriel Exp $ */
   7
   8 #include "flt_misc.h"
   9
  10 void
  11 flt_b64_sft(e,n)
  12         register struct flt_mantissa *e;
  13         register int n;
  14 {
  15         if (n > 63 || n < -63) {
  16                 e->flt_l_32 = 0;
  17                 e->flt_h_32 = 0;
  18                 return;
  19         }
  20         if (n >= 32) {
  21                 e->flt_l_32 = e->flt_h_32;
  22                 e->flt_h_32 = 0;
  23                 n -= 32;
  24         }
  25         if (n > 0) {
  26                 e->flt_l_32 = (e->flt_l_32 >> 1) & 0x7FFFFFFF;
  27                 e->flt_l_32 >>= (n - 1);
  28                 if (e->flt_h_32 != 0) {
  29                         e->flt_l_32 |= (e->flt_h_32 << (32 - n)) & 0xFFFFFFFF;
  30                         e->flt_h_32 = (e->flt_h_32 >> 1) & 0x7FFFFFFF;
  31                         e->flt_h_32 >>= (n - 1);
  32                 }
  33         }
  34         n = -n;
  35         if (n >= 32) {
  36                 e->flt_h_32 = e->flt_l_32;
  37                 e->flt_l_32 = 0;
  38                 n -= 32;
  39         }
  40         if (n > 0) {
  41                 e->flt_h_32 = (e->flt_h_32 << n) & 0xFFFFFFFF;
  42                 if (e->flt_l_32 != 0) {
  43                         long l = (e->flt_l_32 >> 1) & 0x7FFFFFFF;
  44                         e->flt_h_32 |= (l >> (31 - n));
  45                         e->flt_l_32 = (e->flt_l_32 << n) & 0xFFFFFFFF;
  46                 }
  47         }
  48 }