Xorshift128Plus (empty) → 0.1.0.0

raw patch · 4 files changed

+261/−0 lines, 4 filesdep +basesetup-changed

Dependencies added: base

Files

• LICENSE +117/−0
• Setup.hs +2/−0
• System/Random/Xorshift128Plus.hs +104/−0
• Xorshift128Plus.cabal +38/−0

@@ -0,0 +1,117 @@+CC0 1.0 Universal++Statement of Purpose++The laws of most jurisdictions throughout the world automatically confer+exclusive Copyright and Related Rights (defined below) upon the creator and+subsequent owner(s) (each and all, an "owner") of an original work of+authorship and/or a database (each, a "Work").++Certain owners wish to permanently relinquish those rights to a Work for the+purpose of contributing to a commons of creative, cultural and scientific+works ("Commons") that the public can reliably and without fear of later+claims of infringement build upon, modify, incorporate in other works, reuse+and redistribute as freely as possible in any form whatsoever and for any+purposes, including without limitation commercial purposes. These owners may+contribute to the Commons to promote the ideal of a free culture and the+further production of creative, cultural and scientific works, or to gain+reputation or greater distribution for their Work in part through the use and+efforts of others.++For these and/or other purposes and motivations, and without any expectation+of additional consideration or compensation, the person associating CC0 with a+Work (the "Affirmer"), to the extent that he or she is an owner of Copyright+and Related Rights in the Work, voluntarily elects to apply CC0 to the Work+and publicly distribute the Work under its terms, with knowledge of his or her+Copyright and Related Rights in the Work and the meaning and intended legal+effect of CC0 on those rights.++1. Copyright and Related Rights. A Work made available under CC0 may be+protected by copyright and related or neighboring rights ("Copyright and+Related Rights"). Copyright and Related Rights include, but are not limited+to, the following:++  i. the right to reproduce, adapt, distribute, perform, display, communicate,+  and translate a Work;++  ii. moral rights retained by the original author(s) and/or performer(s);++  iii. publicity and privacy rights pertaining to a person's image or likeness+  depicted in a Work;++  iv. rights protecting against unfair competition in regards to a Work,+  subject to the limitations in paragraph 4(a), below;++  v. rights protecting the extraction, dissemination, use and reuse of data in+  a Work;++  vi. database rights (such as those arising under Directive 96/9/EC of the+  European Parliament and of the Council of 11 March 1996 on the legal+  protection of databases, and under any national implementation thereof,+  including any amended or successor version of such directive); and++  vii. other similar, equivalent or corresponding rights throughout the world+  based on applicable law or treaty, and any national implementations thereof.++2. Waiver. To the greatest extent permitted by, but not in contravention of,+applicable law, Affirmer hereby overtly, fully, permanently, irrevocably and+unconditionally waives, abandons, and surrenders all of Affirmer's Copyright+and Related Rights and associated claims and causes of action, whether now+known or unknown (including existing as well as future claims and causes of+action), in the Work (i) in all territories worldwide, (ii) for the maximum+duration provided by applicable law or treaty (including future time+extensions), (iii) in any current or future medium and for any number of+copies, and (iv) for any purpose whatsoever, including without limitation+commercial, advertising or promotional purposes (the "Waiver"). Affirmer makes+the Waiver for the benefit of each member of the public at large and to the+detriment of Affirmer's heirs and successors, fully intending that such Waiver+shall not be subject to revocation, rescission, cancellation, termination, or+any other legal or equitable action to disrupt the quiet enjoyment of the Work+by the public as contemplated by Affirmer's express Statement of Purpose.++3. Public License Fallback. Should any part of the Waiver for any reason be+judged legally invalid or ineffective under applicable law, then the Waiver+shall be preserved to the maximum extent permitted taking into account+Affirmer's express Statement of Purpose. In addition, to the extent the Waiver+is so judged Affirmer hereby grants to each affected person a royalty-free,+non transferable, non sublicensable, non exclusive, irrevocable and+unconditional license to exercise Affirmer's Copyright and Related Rights in+the Work (i) in all territories worldwide, (ii) for the maximum duration+provided by applicable law or treaty (including future time extensions), (iii)+in any current or future medium and for any number of copies, and (iv) for any+purpose whatsoever, including without limitation commercial, advertising or+promotional purposes (the "License"). The License shall be deemed effective as+of the date CC0 was applied by Affirmer to the Work. Should any part of the+License for any reason be judged legally invalid or ineffective under+applicable law, such partial invalidity or ineffectiveness shall not+invalidate the remainder of the License, and in such case Affirmer hereby+affirms that he or she will not (i) exercise any of his or her remaining+Copyright and Related Rights in the Work or (ii) assert any associated claims+and causes of action with respect to the Work, in either case contrary to+Affirmer's express Statement of Purpose.++4. Limitations and Disclaimers.++  a. No trademark or patent rights held by Affirmer are waived, abandoned,+  surrendered, licensed or otherwise affected by this document.++  b. Affirmer offers the Work as-is and makes no representations or warranties+  of any kind concerning the Work, express, implied, statutory or otherwise,+  including without limitation warranties of title, merchantability, fitness+  for a particular purpose, non infringement, or the absence of latent or+  other defects, accuracy, or the present or absence of errors, whether or not+  discoverable, all to the greatest extent permissible under applicable law.++  c. Affirmer disclaims responsibility for clearing rights of other persons+  that may apply to the Work or any use thereof, including without limitation+  any person's Copyright and Related Rights in the Work. Further, Affirmer+  disclaims responsibility for obtaining any necessary consents, permissions+  or other rights required for any use of the Work.++  d. Affirmer understands and acknowledges that Creative Commons is not a+  party to this document and has no duty or obligation with respect to this+  CC0 or use of the Work.++For more information, please see+<http://creativecommons.org/publicdomain/zero/1.0/>+

@@ -0,0 +1,2 @@+import Distribution.Simple+main = defaultMain

@@ -0,0 +1,104 @@+{-# LANGUAGE MagicHash, UnboxedTuples, BangPatterns #-}++{-+Written in 2015 by Kanai Hiroki <kanai.hiroki12@gmail.com>+Originaly written by Sebastiano Vigna.++To the extent possible under law, the author has dedicated all copyright+and related and neighboring rights to this software to the public domain+worldwide. This software is distributed without any warranty.++See <http://creativecommons.org/publicdomain/zero/1.0/>.+-}++-- |+-- Module    : System.Random.Xorshift128Plus+-- License   : Public Commons+--+-- Maintainer  : kanai.hiroki12@gmail.com+-- Stability   : experimental+-- Portability : portable+--+-- This module is implementation of xorshift128+ random number generator.+-- Read <http://xorshift.di.unimi.it/xorshift128plus.c original implementation and description>+-- is strictly recommended.+--+-- The generator state is stored in the 'Gen' data type and+-- it is created by 'initialize' function with seed value or+-- calling 'Gen' data constructor with generator state.+--+-- To generate random values, first 'initialize' a random number+-- generator by seed value.+--+-- @+--   let gen = 'initialize' 3748374974327+-- @+--+-- Then generate a random number by 'next' or 'next01'.+--+-- @+--   let (v, gen') = 'next' gen+-- @+--+-- @+--   let (v, gen') = 'next01' gen+-- @+module System.Random.Xorshift128Plus(+  Gen(..),+  initialize,+  next,+  next01,+  next#+  ) where++import Data.Word(Word64(..))+import Data.Bits++-- | Random number generator.+data Gen = Gen {-# UNPACK #-} !Word64 {-# UNPACK #-} !Word64+         deriving (Eq,Show)++-- | Create a generator by seed value. Please do not supply 0.+initialize :: Word64 -> Gen+initialize s = let s0 = 1812433253 * (s `xor` (s `shiftR` 30)) + 1+                   s1 = 1812433253 * (s0 `xor` (s0 `shiftR` 30)) + 2+               in  Gen s0 s1+{-# INLINE initialize #-}++-- | Generate a 64bit random value and next generator state.+next :: Gen -> (Word64, Gen)+next g = let (# n, g' #) = next# g+         in  (n, g')+{-# INLINE next #-}++-- | Generate a random value between 0 and 1. And next generator state.+next01 :: Gen -> (Double, Gen)+next01 g = let (# n, g' #) = next# g+           in (fromIntegral n / word64max, g')+{-# INLINE next01 #-}++-- | Same as 'next', but return values with unboxed tuple.+next# :: Gen -> (# Word64, Gen #)+next# g = let g'@(Gen s0 s1) = step g+          in  (# (s0 + s1), g' #)+{-# INLINE next# #-}++step :: Gen -> Gen+step (Gen s0 s1) =+  let s1' = s0 `xor` (s0 `shiftL` 23)+  in  Gen s1 (s1' `xor` s1 `xor` (s1' `shiftR` 17) `xor` (s1 `shiftR` 26))+{-# INLINE step #-}++word64max :: Double+word64max = fromIntegral (maxBound :: Word64)+{-# INLINE word64max #-}++-- TODO: unit test+-- test+-- main :: IO ()+-- main = do+--   let g = initialize 3748374974327+--       (# d1, g1 #) = next01# g+--       (# d2, g2 #) = next01# g1+--       (# d3, g3 #) = next01# g2+--   print (d1, d2, d3) -- => (0.4837001127332096,8.873189590555798e-2,0.7393015945448153)

@@ -0,0 +1,38 @@+-- Initial Xorshift128Plus.cabal generated by cabal init.  For further+-- documentation, see http://haskell.org/cabal/users-guide/++name:                Xorshift128Plus+version:             0.1.0.0+synopsis:            Pure haskell implementation of xorshift128plus random number generator.+description:+            This package contains pure haskell implementation of+            xorshift128+ random number generator which is extremely fast.+            .+            The random number generator's inner state is just 128bit memory block+            so that which is not encapusulated for simplicity.+            .+            Please see <http://xorshift.di.unimi.it/xorshift128plus.c original C implementation>+            and <http://xorshift.di.unimi.it/ comparison of otther algorithms>.++license:             PublicDomain+license-file:        LICENSE+homepage:            https://github.com/kanaihiroki/Xorshift128Plus+bug-reports:         https://github.com/kanaihiroki/Xorshift128Plus/issues+author:              Kanai Hiroki <kanai.hiroki12@gmail.com>+maintainer:          Kanai Hiroki <kanai.hiroki12@gmail.com>+category:            Math+build-type:          Simple+cabal-version:       >=1.10++library+  exposed-modules:     System.Random.Xorshift128Plus+  -- other-modules:+  -- other-extensions:+  build-depends:       base >=4.8 && <4.9+  -- hs-source-dirs:+  default-language:    Haskell2010+  ghc-options:         -O3 -fexcess-precision++source-repository head+  type:     git+  location: https://github.com/kanaihiroki/Xorshift128Plus.git