mas01mc@292: /* 
mas01mc@292:    A C-program for MT19937, with initialization improved 2002/1/26.
mas01mc@292:    Coded by Takuji Nishimura and Makoto Matsumoto.
mas01mc@292: 
mas01mc@292:    Before using, initialize the state by using init_genrand(seed)  
mas01mc@292:    or init_by_array(init_key, key_length).
mas01mc@292: 
mas01mc@292:    Copyright (C) 1997 - 2002, Makoto Matsumoto and Takuji Nishimura,
mas01mc@292:    All rights reserved.
mas01mc@292:    Copyright (C) 2005, Mutsuo Saito,
mas01mc@292:    All rights reserved.
mas01mc@292: 
mas01mc@292:    Redistribution and use in source and binary forms, with or without
mas01mc@292:    modification, are permitted provided that the following conditions
mas01mc@292:    are met:
mas01mc@292: 
mas01mc@292:      1. Redistributions of source code must retain the above copyright
mas01mc@292:         notice, this list of conditions and the following disclaimer.
mas01mc@292: 
mas01mc@292:      2. Redistributions in binary form must reproduce the above copyright
mas01mc@292:         notice, this list of conditions and the following disclaimer in the
mas01mc@292:         documentation and/or other materials provided with the distribution.
mas01mc@292: 
mas01mc@292:      3. The names of its contributors may not be used to endorse or promote 
mas01mc@292:         products derived from this software without specific prior written 
mas01mc@292:         permission.
mas01mc@292: 
mas01mc@292:    THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
mas01mc@292:    "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
mas01mc@292:    LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
mas01mc@292:    A PARTICULAR PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE COPYRIGHT OWNER OR
mas01mc@292:    CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
mas01mc@292:    EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
mas01mc@292:    PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
mas01mc@292:    PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
mas01mc@292:    LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
mas01mc@292:    NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
mas01mc@292:    SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
mas01mc@292: 
mas01mc@292: 
mas01mc@292:    Any feedback is very welcome.
mas01mc@292:    http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/emt.html
mas01mc@292:    email: m-mat @ math.sci.hiroshima-u.ac.jp (remove space)
mas01mc@292: */
mas01mc@292: 
mas01mc@292: #include <stdio.h>
mas01mc@292: #include "mt19937ar.h"
mas01mc@292: 
mas01mc@292: int main(void)
mas01mc@292: {
mas01mc@292:     int i;
mas01mc@292:     unsigned long init[4]={0x123, 0x234, 0x345, 0x456}, length=4;
mas01mc@292:     init_by_array(init, length);
mas01mc@292:     printf("1000 outputs of genrand_int32()\n");
mas01mc@292:     for (i=0; i<1000; i++) {
mas01mc@292:       printf("%10lu ", genrand_int32());
mas01mc@292:       if (i%5==4) printf("\n");
mas01mc@292:     }
mas01mc@292:     printf("\n1000 outputs of genrand_real2()\n");
mas01mc@292:     for (i=0; i<1000; i++) {
mas01mc@292:       printf("%10.8f ", genrand_real2());
mas01mc@292:       if (i%5==4) printf("\n");
mas01mc@292:     }
mas01mc@292:     return 0;
mas01mc@292: }