0

我想为大整数的按位运算实现一个高效的库。我编写了以下覆盖的函数BTEST：

FUNCTION testb_i2b(n,i)
  INTEGER(I8B), DIMENSION(0:), INTENT(IN) :: n
  INTEGER(I2B), INTENT(IN) :: i
  INTEGER(I2B) :: j
  LOGICAL :: testb_i2b
  j = ISHFT(i,-6)
  IF ( j .LE. UBOUND(n,1) ) THEN
    testb_i2b = BTEST(n(j),i-ISHFT(j,6))
  ELSE
    testb_i2b = .FALSE.
  END IF
END FUNCTION testb_i2b

该数组n包含64*(SIZE(n)-1)我的大整数的位。有没有更有效的方法来获得相同的功能？

4

1 回答 1

3

我不知道这是否比你的版本更快，我会让你测试一下，但它涉及的操作更少，并且没有明确的 if 语句，所以可能是这样。对于我运行的几个测试，它给出的结果与您的代码相同。我已经将 bignum 中整数的大小固定为 64 位，如果您愿意，可以将其设为参数。

  LOGICAL FUNCTION btest_bignum(bn,ix)
    IMPLICIT NONE
    INTEGER(int64), DIMENSION(0:), INTENT(in) :: bn
    INTEGER(int16), INTENT(in) :: ix
    INTEGER :: array_ix
    array_ix = ix/64
    btest_bignum = BTEST(bn(array_ix), ix-(array_ix*64))
  END FUNCTION btest_bignum

请注意，我使用了现在标准的种类声明int64和int16

于 2013-03-12T16:18:01.740 回答