diff --git a/LICENSE b/LICENSE new file mode 100644 --- /dev/null +++ b/LICENSE @@ -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 + + diff --git a/Setup.hs b/Setup.hs new file mode 100644 --- /dev/null +++ b/Setup.hs @@ -0,0 +1,2 @@ +import Distribution.Simple +main = defaultMain diff --git a/System/Random/Xorshift128Plus.hs b/System/Random/Xorshift128Plus.hs new file mode 100644 --- /dev/null +++ b/System/Random/Xorshift128Plus.hs @@ -0,0 +1,104 @@ +{-# LANGUAGE MagicHash, UnboxedTuples, BangPatterns #-} + +{- +Written in 2015 by Kanai Hiroki +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 . +-} + +-- | +-- 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 +-- 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) diff --git a/Xorshift128Plus.cabal b/Xorshift128Plus.cabal new file mode 100644 --- /dev/null +++ b/Xorshift128Plus.cabal @@ -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 + and . + +license: PublicDomain +license-file: LICENSE +homepage: https://github.com/kanaihiroki/Xorshift128Plus +bug-reports: https://github.com/kanaihiroki/Xorshift128Plus/issues +author: Kanai Hiroki +maintainer: Kanai Hiroki +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