Home | History | Annotate | Line # | Download | only in intrinsics
      1      1.1  mrg /* Implementation of ishftc intrinsic.
      2  1.1.1.3  mrg    Copyright (C) 2002-2022 Free Software Foundation, Inc.
      3      1.1  mrg    Contributed by Paul Brook <paul (at) nowt.org>
      4      1.1  mrg 
      5      1.1  mrg This file is part of the GNU Fortran 95 runtime library (libgfortran).
      6      1.1  mrg 
      7      1.1  mrg Libgfortran is free software; you can redistribute it and/or
      8      1.1  mrg modify it under the terms of the GNU General Public
      9      1.1  mrg License as published by the Free Software Foundation; either
     10      1.1  mrg version 3 of the License, or (at your option) any later version.
     11      1.1  mrg 
     12      1.1  mrg Libgfortran is distributed in the hope that it will be useful,
     13      1.1  mrg but WITHOUT ANY WARRANTY; without even the implied warranty of
     14      1.1  mrg MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
     15      1.1  mrg GNU General Public License for more details.
     16      1.1  mrg 
     17      1.1  mrg Under Section 7 of GPL version 3, you are granted additional
     18      1.1  mrg permissions described in the GCC Runtime Library Exception, version
     19      1.1  mrg 3.1, as published by the Free Software Foundation.
     20      1.1  mrg 
     21      1.1  mrg You should have received a copy of the GNU General Public License and
     22      1.1  mrg a copy of the GCC Runtime Library Exception along with this program;
     23      1.1  mrg see the files COPYING3 and COPYING.RUNTIME respectively.  If not, see
     24      1.1  mrg <http://www.gnu.org/licenses/>.  */
     25      1.1  mrg 
     26      1.1  mrg #include "libgfortran.h"
     27      1.1  mrg 
     28      1.1  mrg extern GFC_INTEGER_4 ishftc4 (GFC_INTEGER_4, GFC_INTEGER_4, GFC_INTEGER_4);
     29      1.1  mrg export_proto(ishftc4);
     30      1.1  mrg 
     31      1.1  mrg GFC_INTEGER_4
     32      1.1  mrg ishftc4 (GFC_INTEGER_4 i, GFC_INTEGER_4 shift, GFC_INTEGER_4 size)
     33      1.1  mrg {
     34      1.1  mrg   GFC_UINTEGER_4 mask, bits;
     35      1.1  mrg 
     36      1.1  mrg   if (shift < 0)
     37      1.1  mrg     shift = shift + size;
     38      1.1  mrg 
     39      1.1  mrg   if (shift == 0 || shift == size)
     40      1.1  mrg     return i;
     41      1.1  mrg 
     42      1.1  mrg   /* In C, the result of the shift operator is undefined if the right operand
     43      1.1  mrg      is greater than or equal to the number of bits in the left operand. So we
     44      1.1  mrg      have to special case it for fortran.  */
     45      1.1  mrg   mask = ~((size == 32) ? (GFC_UINTEGER_4)0 : (~(GFC_UINTEGER_4)0 << size));
     46      1.1  mrg 
     47      1.1  mrg   bits = i & mask;
     48      1.1  mrg 
     49      1.1  mrg   return (i & ~mask) | ((bits << shift) & mask) | (bits >> (size - shift));
     50      1.1  mrg }
     51      1.1  mrg 
     52      1.1  mrg extern GFC_INTEGER_8 ishftc8 (GFC_INTEGER_8, GFC_INTEGER_4, GFC_INTEGER_4);
     53      1.1  mrg export_proto(ishftc8);
     54      1.1  mrg 
     55      1.1  mrg GFC_INTEGER_8
     56      1.1  mrg ishftc8 (GFC_INTEGER_8 i, GFC_INTEGER_4 shift, GFC_INTEGER_4 size)
     57      1.1  mrg {
     58      1.1  mrg   GFC_UINTEGER_8 mask, bits;
     59      1.1  mrg 
     60      1.1  mrg   if (shift < 0)
     61      1.1  mrg     shift = shift + size;
     62      1.1  mrg 
     63      1.1  mrg   if (shift == 0 || shift == size)
     64      1.1  mrg     return i;
     65      1.1  mrg 
     66      1.1  mrg   /* In C, the result of the shift operator is undefined if the right operand
     67      1.1  mrg      is greater than or equal to the number of bits in the left operand. So we
     68      1.1  mrg      have to special case it for fortran.  */
     69      1.1  mrg   mask = ~((size == 64) ? (GFC_UINTEGER_8)0 : (~(GFC_UINTEGER_8)0 << size));
     70      1.1  mrg 
     71      1.1  mrg   bits = i & mask;
     72      1.1  mrg 
     73      1.1  mrg   return (i & ~mask) | ((bits << shift) & mask) | (bits >> (size - shift));
     74      1.1  mrg }
     75      1.1  mrg 
     76      1.1  mrg #ifdef HAVE_GFC_INTEGER_16
     77      1.1  mrg extern GFC_INTEGER_16 ishftc16 (GFC_INTEGER_16, GFC_INTEGER_4, GFC_INTEGER_4);
     78      1.1  mrg export_proto(ishftc16);
     79      1.1  mrg 
     80      1.1  mrg GFC_INTEGER_16
     81      1.1  mrg ishftc16 (GFC_INTEGER_16 i, GFC_INTEGER_4 shift, GFC_INTEGER_4 size)
     82      1.1  mrg {
     83      1.1  mrg   GFC_UINTEGER_16 mask, bits;
     84      1.1  mrg 
     85      1.1  mrg   if (shift < 0)
     86      1.1  mrg     shift = shift + size;
     87      1.1  mrg 
     88      1.1  mrg   if (shift == 0 || shift == size)
     89      1.1  mrg     return i;
     90      1.1  mrg 
     91      1.1  mrg   /* In C, the result of the shift operator is undefined if the right operand
     92      1.1  mrg      is greater than or equal to the number of bits in the left operand. So we
     93      1.1  mrg      have to special case it for fortran.  */
     94      1.1  mrg   mask = ~((size == 128) ? (GFC_UINTEGER_16)0 : (~(GFC_UINTEGER_16)0 << size));
     95      1.1  mrg 
     96      1.1  mrg   bits = i & mask;
     97      1.1  mrg 
     98      1.1  mrg   return (i & ~mask) | ((bits << shift) & mask) | (bits >> (size - shift));
     99      1.1  mrg }
    100      1.1  mrg #endif
    101