在C中,如果我想將int除以2,x%2
應該運行得像(x%10)% 2
一樣快,因爲好的編譯器只會查看最後一位。但是如何用無限精度的算術語言呢?無限精度整數:除以2
特別是在Haskell中,它會更快(或者它們會是相同的速度):even x
或even (quot x 10)
?
在C中,如果我想將int除以2,x%2
應該運行得像(x%10)% 2
一樣快,因爲好的編譯器只會查看最後一位。但是如何用無限精度的算術語言呢?無限精度整數:除以2
特別是在Haskell中,它會更快(或者它們會是相同的速度):even x
或even (quot x 10)
?
好吧,我會咬:
import Control.DeepSeq
import Criterion.Main
import Data.Bits
import System.Random
lotsOfBigNumbers :: [Integer]
lotsOfBigNumbers = take 10000 $ randomRs (2^128, 2^132) (mkStdGen 42)
evenRem, evenBits :: Integer -> Bool
evenRem x = even (x `rem` 10)
evenBits x = (x .&. 1) == 0
remRem x = ((x `rem` 10) `rem` 2) == 0
main = do
deepseq lotsOfBigNumbers (return())
defaultMain
[ bench "even" $ nf (map even ) lotsOfBigNumbers
, bench "evenRem" $ nf (map evenRem) lotsOfBigNumbers
, bench "evenBits" $ nf (map evenBits) lotsOfBigNumbers
, bench "remRem" $ nf (map remRem ) lotsOfBigNumbers
]
而且結果:
sorghum:~/programming% ghc -O2 test && ./test
[1 of 1] Compiling Main (test.hs, test.o)
Linking test ...
warming up
estimating clock resolution...
mean is 1.920340 us (320001 iterations)
found 51108 outliers among 319999 samples (16.0%)
46741 (14.6%) low severe
4367 (1.4%) high severe
estimating cost of a clock call...
mean is 83.20445 ns (16 iterations)
found 4 outliers among 16 samples (25.0%)
2 (12.5%) low mild
2 (12.5%) high severe
benchmarking even
mean: 1.508542 ms, lb 1.503661 ms, ub 1.514950 ms, ci 0.950
std dev: 28.53574 us, lb 23.35796 us, ub 35.19769 us, ci 0.950
found 18 outliers among 100 samples (18.0%)
17 (17.0%) high severe
variance introduced by outliers: 11.371%
variance is moderately inflated by outliers
benchmarking evenRem
mean: 1.937735 ms, lb 1.930118 ms, ub 1.949699 ms, ci 0.950
std dev: 48.17240 us, lb 34.95334 us, ub 71.22055 us, ci 0.950
found 14 outliers among 100 samples (14.0%)
3 (3.0%) high mild
11 (11.0%) high severe
variance introduced by outliers: 18.989%
variance is moderately inflated by outliers
benchmarking evenBits
mean: 996.3537 us, lb 992.2839 us, ub 1.003864 ms, ci 0.950
std dev: 27.37875 us, lb 17.51923 us, ub 43.98499 us, ci 0.950
found 15 outliers among 100 samples (15.0%)
2 (2.0%) high mild
13 (13.0%) high severe
variance introduced by outliers: 21.905%
variance is moderately inflated by outliers
benchmarking remRem
mean: 1.925495 ms, lb 1.918590 ms, ub 1.935869 ms, ci 0.950
std dev: 43.00092 us, lb 31.67173 us, ub 57.83841 us, ci 0.950
found 13 outliers among 100 samples (13.0%)
13 (13.0%) high severe
variance introduced by outliers: 15.198%
variance is moderately inflated by outliers
結論:一個額外的費用rem
多一點,和.&.
費用少一點。
'evenBits x =(x。|。1)== 0' < - 應該是'。&。'。這使得這個比較快(〜20%)。但實際上,它應該是'evenBits x =(fromInteger x。&。(1 :: Int))== 0',因爲'Integer'上的位操作並不是特別快。 –
@DanielFischer哎呀,這是一個糟糕的錯誤,謝謝指出! –
我注意到另一個,'evenquot x = even(x'''10)''檢查第二個最低有效位是否是偶數。應該是「even(x'rem' 10)''。然後差別幾乎消失了,因爲''x'rem' 10''變成'S#i #',然後'even'使用'Int#'除法(在測試它是'S#'後)而且比GMP的調用速度快得多,所以'整數'除法+'內部'除法與「整數」除法幾乎相同。 –
我不太確定'x%2'與C中的'(x%10)%2'一樣快。如果編譯器認識到'%10'部分不改變結果,那麼它可能會被刪除,但這是一種算術簡化 - 與位擺動和數字表示無關,因此同樣適用於任意精度整數。如果它沒有被移除,那麼它可能確實會首先計算餘數w/10,這是通過屏蔽掉位來實現的。 – delnan
@delnan這是一個很公平的點,謝謝。我想我不是很嚴謹。 –