XorShift On Tcl, Tcl

Language: Layout:

#XorShift Module
oo::class create XorShift {
    variable x
    variable y
    variable z
    variable w
    variable seeds
    method seeds {} {return [array get seeds]}
    variable randCount
    method randCount {} {return $randCount}

constructor {
        {_w "[clock clicks -milliseconds]"}
        {_x "$w<<13"}
        {_y "($w>>9)^($x<<6)"}
        {_z "$y>>7"}
    } {
        set randCount 0
        set w [expr $_w]
        set x [expr $_x]
        set y [expr $_y]
        set z [expr $_z]

array set seeds [list x $x y $y z $z w $w]
    }

method rand {} {
        incr randCount
        set t [expr $x^($x<<11&0xFFFFFFFF)]
        set x $y
        set y $z
        set z $w
        return [set w [expr ($w^($w>>19))^($t^($t>>8))]]
    }

method randInt {{min 0} {max 0x7FFFFFFF}} {
        return [expr [my rand]%($max-$min+1)+$min]
    }

method randFloat {{min 0} {max 1}} {
        return [expr [my rand]%0xFFFF*1./0xFFFF*($max-$min)+$min]
    }

method shuffle {arr} {
        for {set i 0} {$i<=[llength $arr]-2} {incr i} {
            set r [my randInt $i [expr [llength $arr]-1]]
            set tmp [lindex $arr $i]
            lset arr $i [lindex $arr $r]
            lset arr $r $tmp
        }
        return $arr
    }

}
namespace eval XorShift {
    variable defaults
    array set defaults {
        x 123456789
        y 362436069
        z 521288629
        w 88675123
    }
    oo::class create defaultSeed {
        superclass XorShift
        constructor {} {
            next $XorShift::defaults(w) \
                 $XorShift::defaults(x) \
                 $XorShift::defaults(y) \
                 $XorShift::defaults(z)
        }
    }
}
#XorShift End

#source XorShift.tcl

apply {{} {
    #XorShift乱数ジェネレータの初期化
    # 論文デフォルトシード
    set r_def [XorShift::defaultSeed new]
    # 固定値シード
    set r_var [XorShift new 100]
    # 時間シード
    set r [XorShift new]

#デフォルトシード値の取得
    puts ">> defaults"
    puts [array get XorShift::defaults]
    #for(var i in XorShift.defaults){
    #   WScript.Echo(i+":"+XorShift.defaults[i]);
    #}

#適用したシード値の取得
    puts ">> seeds in r"
    array set r_seeds [$r seeds]
    puts [array get r_seeds]

#乱数の生データを取得
    puts ">> rand 0 to UInt32Max"
    for {set i 0} {$i<5} {incr i} {
        puts [$r_def rand]
    }

#0-100の乱数(100含む)を整数で取得
    puts ">> randInt 0 to 100"
    for {set i 0} {$i<5} {incr i} {
        puts [$r_var randInt 0 100]
    }

#0-1の乱数を浮遊小数点で取得
    puts ">> randFloat 0 to 1"
    for {set i 0} {$i<5} {incr i} {
        puts [$r randFloat]
    }

#listのシャッフル
    #値渡しとなるので元のlistは破壊されない
    puts ">> shuffle list"
    set a {}
    for {set i 0} {$i<20} {incr i} {set a [concat $a $i]}

puts [$r shuffle $a]
    puts $a

#今の乱数を回した回数
    puts ">> randCount in r"
    puts [$r randCount]
}}

[ + ] Show input
Absolute running time: 0.09 sec, cpu time: 0.04 sec, memory peak: 3 Mb, absolute service time: 0,11 sec