8817e257e3927ef32ecabc0c11d09cc186699cbeb87222ae428ed99e4a0170bb

Download raw

Created: 2023-12-19 03:40:10.302489 UTC

Size: 462

Preview (limited to 512 bytes)

"module Algo.Spectral ( spectral ) where\n\nimport Data.Vector.Unboxed as V\n\nimport Data.Bits\n\nspectral :: Vector Double -> Vector Double\n{-# NOINLINE spectral #-}\nspectral us = us `seq` V.map row (V.enumFromTo 0 (n-1))\n where\n n = V.length us\n\n row i = i `seq` V.sum (V.imap (\\j u -> eval_A i j * u) us)\n\n eval_A i j = 1 / fromIntegral r\n where\n r = u + (i+1)\n u = t `shiftR` 1\n t = n * (n+1)\n n = i+j\n\n"

Casa is a service provided by the Haskell Foundation │ Originally developed by FP Complete