1.1 --- a/.hgignore	Sat Mar 03 10:54:39 2012 -0600
     1.2 +++ b/.hgignore	Sat Mar 03 11:03:44 2012 -0600
     1.3 @@ -8,4 +8,5 @@
     1.4  depcomp
     1.5  autoscan.log
     1.6  configure.scan
     1.7 -*Makefile.in*
     1.8 \ No newline at end of file
     1.9 +*Makefile.in*
    1.10 +config.h.in
    1.11 \ No newline at end of file

     2.1 --- a/configure.ac	Sat Mar 03 10:54:39 2012 -0600
     2.2 +++ b/configure.ac	Sat Mar 03 11:03:44 2012 -0600
     2.3 @@ -1,5 +1,26 @@
     2.4 +AC_PREREQ([2.68])
     2.5 +
     2.6  AC_INIT([vba-rlm], [1.0])
     2.7  
     2.8 +dnl TODO: change this to gba.cpp or something
     2.9 +AC_CONFIG_SRCDIR([src/lua/src/lopcodes.c])
    2.10 +AC_CONFIG_HEADERS([config.h])
    2.11 +AC_PROG_RANLIB
    2.12 +
    2.13 +# Checks for header files.
    2.14 +AC_CHECK_HEADERS([limits.h locale.h stddef.h stdlib.h string.h unistd.h])
    2.15 +
    2.16 +# Checks for typedefs, structures, and compiler characteristics.
    2.17 +AC_TYPE_SIZE_T
    2.18 +AC_CHECK_TYPES([ptrdiff_t])
    2.19 +
    2.20 +# Checks for library functions.
    2.21 +AC_FUNC_ERROR_AT_LINE
    2.22 +AC_FUNC_MKTIME
    2.23 +AC_FUNC_REALLOC
    2.24 +AC_FUNC_STRCOLL
    2.25 +AC_CHECK_FUNCS([floor localeconv memchr modf pow setlocale sqrt strchr strcspn strerror strpbrk strrchr strstr strtoul])
    2.26 +
    2.27  AC_CONFIG_FILES([Makefile 
    2.28  	         src/Makefile
    2.29                   src/lua/Makefile])

     3.1 --- a/src/AutoBuild.h	Sat Mar 03 10:54:39 2012 -0600
     3.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
     3.3 @@ -1,30 +0,0 @@
     3.4 -// -*- C++ -*-
     3.5 -// VisualBoyAdvance - Nintendo Gameboy/GameboyAdvance (TM) emulator.
     3.6 -// Copyright (C) 1999-2003 Forgotten
     3.7 -// Copyright (C) 2004 Forgotten and the VBA development team
     3.8 -
     3.9 -// This program is free software; you can redistribute it and/or modify
    3.10 -// it under the terms of the GNU General Public License as published by
    3.11 -// the Free Software Foundation; either version 2, or(at your option)
    3.12 -// any later version.
    3.13 -//
    3.14 -// This program is distributed in the hope that it will be useful,
    3.15 -// but WITHOUT ANY WARRANTY; without even the implied warranty of
    3.16 -// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    3.17 -// GNU General Public License for more details.
    3.18 -//
    3.19 -// You should have received a copy of the GNU General Public License
    3.20 -// along with this program; if not, write to the Free Software Foundation,
    3.21 -// Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
    3.22 -
    3.23 -#ifndef __AUTOBUILD_H__
    3.24 -#define __AUTOBUILD_H__
    3.25 -#include "version.h"
    3.26 -//change the FALSE to TRUE for autoincrement of build number
    3.27 -#define INCREMENT_VERSION FALSE
    3.28 -#define FILEVER        1,8,0,600
    3.29 -#define PRODUCTVER     1,8,0,600
    3.30 -#define STRFILEVER     "1, 8, 0, 600\0"
    3.31 -#define STRPRODUCTVER  "1, 8, 0, 600\0"
    3.32 -#endif //__AUTOBUILD_H__
    3.33 -

     4.1 --- a/src/NLS.h	Sat Mar 03 10:54:39 2012 -0600
     4.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
     4.3 @@ -1,52 +0,0 @@
     4.4 -#ifndef VBS_NLS_H
     4.5 -#define VBA_NLS_H
     4.6 -
     4.7 -#if _MSC_VER > 1000
     4.8 -#pragma once
     4.9 -#endif // _MSC_VER > 1000
    4.10 -
    4.11 -#define N_(String) (String)
    4.12 -
    4.13 -#define MSG_UNSUPPORTED_VBA_SGM             1
    4.14 -#define MSG_CANNOT_LOAD_SGM                 2
    4.15 -#define MSG_SAVE_GAME_NOT_USING_BIOS        3
    4.16 -#define MSG_SAVE_GAME_USING_BIOS            4
    4.17 -#define MSG_UNSUPPORTED_SAVE_TYPE           5
    4.18 -#define MSG_CANNOT_OPEN_FILE                6
    4.19 -#define MSG_BAD_ZIP_FILE                    7
    4.20 -#define MSG_NO_IMAGE_ON_ZIP                 8
    4.21 -#define MSG_ERROR_OPENING_IMAGE             9
    4.22 -#define MSG_ERROR_READING_IMAGE            10
    4.23 -#define MSG_UNSUPPORTED_BIOS_FUNCTION      11
    4.24 -#define MSG_INVALID_BIOS_FILE_SIZE         12
    4.25 -#define MSG_INVALID_CHEAT_CODE             13
    4.26 -#define MSG_UNKNOWN_ARM_OPCODE             14
    4.27 -#define MSG_UNKNOWN_THUMB_OPCODE           15
    4.28 -#define MSG_ERROR_CREATING_FILE            16
    4.29 -#define MSG_FAILED_TO_READ_SGM             17
    4.30 -#define MSG_FAILED_TO_READ_RTC             18
    4.31 -#define MSG_UNSUPPORTED_VB_SGM             19
    4.32 -#define MSG_CANNOT_LOAD_SGM_FOR            20
    4.33 -#define MSG_ERROR_OPENING_IMAGE_FROM       21
    4.34 -#define MSG_ERROR_READING_IMAGE_FROM       22
    4.35 -#define MSG_UNSUPPORTED_ROM_SIZE           23
    4.36 -#define MSG_UNSUPPORTED_RAM_SIZE           24
    4.37 -#define MSG_UNKNOWN_CARTRIDGE_TYPE         25
    4.38 -#define MSG_MAXIMUM_NUMBER_OF_CHEATS       26
    4.39 -#define MSG_INVALID_GAMESHARK_CODE         27
    4.40 -#define MSG_INVALID_GAMEGENIE_CODE         28
    4.41 -#define MSG_INVALID_CHEAT_TO_REMOVE        29
    4.42 -#define MSG_INVALID_CHEAT_CODE_ADDRESS     30
    4.43 -#define MSG_UNSUPPORTED_CHEAT_LIST_VERSION 31
    4.44 -#define MSG_UNSUPPORTED_CHEAT_LIST_TYPE    32
    4.45 -#define MSG_INVALID_GSA_CODE               33
    4.46 -#define MSG_CANNOT_IMPORT_SNAPSHOT_FOR     34
    4.47 -#define MSG_UNSUPPORTED_SNAPSHOT_FILE      35
    4.48 -#define MSG_UNSUPPORTED_ARM_MODE           36
    4.49 -#define MSG_UNSUPPORTED_CODE_FILE          37
    4.50 -#define MSG_GBA_CODE_WARNING               38
    4.51 -#define MSG_INVALID_CBA_CODE               39
    4.52 -#define MSG_CBA_CODE_WARNING               40
    4.53 -#define MSG_OUT_OF_MEMORY                  41
    4.54 -
    4.55 -#endif // VBA_NLS_H

     5.1 --- a/src/Port.h	Sat Mar 03 10:54:39 2012 -0600
     5.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
     5.3 @@ -1,178 +0,0 @@
     5.4 -#ifndef VBA_PORT_H
     5.5 -#define VBA_PORT_H
     5.6 -
     5.7 -#if _MSC_VER > 1000
     5.8 -#pragma once
     5.9 -#endif // _MSC_VER > 1000
    5.10 -
    5.11 -#include <ctime>
    5.12 -
    5.13 -#ifndef NULL
    5.14 -#define NULL 0
    5.15 -#endif
    5.16 -
    5.17 -typedef unsigned char bool8;
    5.18 -
    5.19 -#ifdef HAVE_STDINT_H
    5.20 -#include <stdint.h>
    5.21 -
    5.22 -typedef int8_t   int8;
    5.23 -typedef uint8_t  uint8;
    5.24 -typedef int16_t  int16;
    5.25 -typedef uint16_t uint16;
    5.26 -typedef int32_t  int32;
    5.27 -typedef uint32_t uint32;
    5.28 -typedef int64_t  int64;
    5.29 -typedef uint64_t uint64;
    5.30 -typedef intptr_t pint;
    5.31 -
    5.32 -#else /* Don't have stdint.h */
    5.33 -
    5.34 -#ifdef PTR_NOT_INT
    5.35 -typedef long pint;
    5.36 -#else /* pointer is int */
    5.37 -typedef int pint;
    5.38 -#endif /* PTR_NOT_INT */
    5.39 -
    5.40 -/* FIXME: Refactor this by moving out the BORLAND part and unifying typedefs */
    5.41 -#ifndef WIN32
    5.42 -typedef unsigned char  uint8;
    5.43 -typedef unsigned short uint16;
    5.44 -typedef signed char    int8;
    5.45 -typedef short          int16;
    5.46 -typedef int            int32;
    5.47 -typedef unsigned int   uint32;
    5.48 -# ifdef __GNUC__  /* long long is not part of ISO C++ */
    5.49 -__extension__ typedef long long          int64;
    5.50 -__extension__ typedef unsigned long long uint64;
    5.51 -# else
    5.52 -typedef long long          int64;
    5.53 -typedef unsigned long long uint64;
    5.54 -# endif
    5.55 -#else /* WIN32 */
    5.56 -
    5.57 -# ifdef __BORLANDC__
    5.58 -#   include <systypes.h>
    5.59 -# else
    5.60 -
    5.61 -typedef unsigned char  uint8;
    5.62 -typedef unsigned short uint16;
    5.63 -typedef signed char    int8;
    5.64 -typedef short          int16;
    5.65 -
    5.66 -# ifndef WSAAPI
    5.67 -/* winsock2.h typedefs int32 as well. */
    5.68 -typedef long int32;
    5.69 -# endif
    5.70 -
    5.71 -typedef unsigned int uint32;
    5.72 -
    5.73 -# endif /* __BORLANDC__ */
    5.74 -
    5.75 -typedef __int64          int64;
    5.76 -typedef unsigned __int64 uint64;
    5.77 -
    5.78 -#endif /* WIN32 */
    5.79 -#endif /* HAVE_STDINT_H */
    5.80 -
    5.81 -#ifndef WIN32
    5.82 -
    5.83 -#ifndef PATH_MAX
    5.84 -#define PATH_MAX 1024
    5.85 -#endif
    5.86 -
    5.87 -#define _MAX_DIR PATH_MAX
    5.88 -#define _MAX_DRIVE 1
    5.89 -#define _MAX_FNAME PATH_MAX
    5.90 -#define _MAX_EXT PATH_MAX
    5.91 -#define _MAX_PATH PATH_MAX
    5.92 -
    5.93 -#define ZeroMemory(a, b) memset((a), 0, (b))
    5.94 -
    5.95 -void _makepath(char *path, const char *drive, const char *dir,
    5.96 -               const char *fname, const char *ext);
    5.97 -void _splitpath(const char *path, char *drive, char *dir, char *fname,
    5.98 -                char *ext);
    5.99 -#else /* WIN32 */
   5.100 -#define strcasecmp stricmp
   5.101 -#define strncasecmp strnicmp
   5.102 -#endif
   5.103 -
   5.104 -typedef uint8  u8;
   5.105 -typedef uint16 u16;
   5.106 -typedef uint32 u32;
   5.107 -typedef uint64 u64;
   5.108 -typedef int8   s8;
   5.109 -typedef int16  s16;
   5.110 -typedef int32  s32;
   5.111 -typedef int64  s64;
   5.112 -
   5.113 -// for consistency
   5.114 -static inline u8 swap8(u8 v)
   5.115 -{
   5.116 -	return v;
   5.117 -}
   5.118 -
   5.119 -// swaps a 16-bit value
   5.120 -static inline u16 swap16(u16 v)
   5.121 -{
   5.122 -	return (v<<8)|(v>>8);
   5.123 -}
   5.124 -
   5.125 -// swaps a 32-bit value
   5.126 -static inline u32 swap32(u32 v)
   5.127 -{
   5.128 -	return (v<<24)|((v<<8)&0xff0000)|((v>>8)&0xff00)|(v>>24);
   5.129 -}
   5.130 -
   5.131 -#define READ8LE(x) \
   5.132 -    *((u8 *)x)
   5.133 -
   5.134 -#define WRITE8LE(x, v) \
   5.135 -    *((u8 *)x) = (v)
   5.136 -
   5.137 -#ifdef WORDS_BIGENDIAN
   5.138 -#if defined(__GNUC__) && defined(__ppc__)
   5.139 -
   5.140 -#define READ16LE(base) \
   5.141 -    ({ unsigned short lhbrxResult;       \
   5.142 -       __asm__("lhbrx %0, 0, %1" : "=r" (lhbrxResult) : "r" (base) : "memory"); \
   5.143 -       lhbrxResult; })
   5.144 -
   5.145 -#define READ32LE(base) \
   5.146 -    ({ unsigned long lwbrxResult; \
   5.147 -       __asm__("lwbrx %0, 0, %1" : "=r" (lwbrxResult) : "r" (base) : "memory"); \
   5.148 -       lwbrxResult; })
   5.149 -
   5.150 -#define WRITE16LE(base, value) \
   5.151 -    __asm__("sthbrx %0, 0, %1" : : "r" (value), "r" (base) : "memory")
   5.152 -
   5.153 -#define WRITE32LE(base, value) \
   5.154 -    __asm__("stwbrx %0, 0, %1" : : "r" (value), "r" (base) : "memory")
   5.155 -
   5.156 -#else
   5.157 -#define READ16LE(x) \
   5.158 -    swap16(*((u16 *)(x)))
   5.159 -#define READ32LE(x) \
   5.160 -    swap32(*((u32 *)(x)))
   5.161 -#define WRITE16LE(x, v) \
   5.162 -    *((u16 *)x) = swap16((v))
   5.163 -#define WRITE32LE(x, v) \
   5.164 -    *((u32 *)x) = swap32((v))
   5.165 -#endif
   5.166 -#else
   5.167 -#define READ16LE(x) \
   5.168 -    *((u16 *)x)
   5.169 -#define READ32LE(x) \
   5.170 -    *((u32 *)x)
   5.171 -#define WRITE16LE(x, v) \
   5.172 -    *((u16 *)x) = (v)
   5.173 -#define WRITE32LE(x, v) \
   5.174 -    *((u32 *)x) = (v)
   5.175 -#endif
   5.176 -
   5.177 -#ifndef CTASSERT
   5.178 -#define CTASSERT(x)  typedef char __assert ## y[(x) ? 1 : -1];
   5.179 -#endif
   5.180 -
   5.181 -#endif // VBA_PORT_H

     6.1 --- a/src/SFMT/CHANGE-LOG.txt	Sat Mar 03 10:54:39 2012 -0600
     6.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
     6.3 @@ -1,55 +0,0 @@
     6.4 -ver 1.3.3
     6.5 --------
     6.6 -change condition compile of do_recursion in SFMT.c
     6.7 -
     6.8 -ver 1.3.2
     6.9 --------
    6.10 -bug fix to_res53_mix and genrand_res53_mix.
    6.11 -
    6.12 -ver 1.3.1
    6.13 --------
    6.14 -gcc compile option changed form -O9 to -O3.
    6.15 -add functions genrand_res53_mix and to_res53_mix.
    6.16 -bug fix about definition of ALWAYS_INLINE.
    6.17 -add new definition PRE_ALWAYS for MSC.
    6.18 -
    6.19 -ver 1.3
    6.20 --------
    6.21 -bug fixed: -DONLY64 without -DBIG_ENIAN64 had been generating
    6.22 -wrong sequence.
    6.23 -bug fixed: There is no documentation about BIG_ENDIAN64.
    6.24 -add automatic endian check by __BIG_ENDIAN__ predefined macro.
    6.25 -bug fixed: change == in check.sh to =
    6.26 -add SFMT-params216091.h
    6.27 -add AltiVec parameter format for systems which are not osx.
    6.28 -change Makefile for systems which are not osx and support AltiVec.
    6.29 -change sample2 of howto-compile for Free BSD.
    6.30 -change source files for BORLANDC and Visual Studio.
    6.31 -change period certification code more smart.
    6.32 -add params directory.
    6.33 -
    6.34 -ver 1.2.1
    6.35 --------
    6.36 -Fix typo in SFMT-alti.c SFMT-sse2.c
    6.37 -marge SFMT-alti.c and SFMT-alti.h into SFMT-alti.h
    6.38 -marge SFMT-sse2.c and SFMT-sse2.h into SFMT-sse2.h
    6.39 -This version is not released.
    6.40 -
    6.41 -ver 1.2
    6.42 --------
    6.43 -Support many periods: 2^{607}, 2^{1279}, 2^{2281}, 2^{4253}, 2^{11213},
    6.44 -2^{19937}, 2^{44497}, 2^{86243}, 2^{132049}
    6.45 -Fix typo in LICENSE.txt.
    6.46 -Add cast to vec_perm for SFMT-alti.c, SFMT-alti64.c.
    6.47 -combine source codes.
    6.48 -
    6.49 -ver 1.1
    6.50 --------
    6.51 -The period certification method is changed from constant to function.
    6.52 -The convert functions from 32-bit and 64-bit integer to double are added.
    6.53 -The documentation is changed.
    6.54 -Sample programs are added.
    6.55 -
    6.56 -ver 1.0
    6.57 --------
    6.58 -The first version.

     7.1 --- a/src/SFMT/LICENSE.txt	Sat Mar 03 10:54:39 2012 -0600
     7.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
     7.3 @@ -1,29 +0,0 @@
     7.4 -Copyright (c) 2006,2007 Mutsuo Saito, Makoto Matsumoto and Hiroshima
     7.5 -University. All rights reserved.
     7.6 -
     7.7 -Redistribution and use in source and binary forms, with or without
     7.8 -modification, are permitted provided that the following conditions are
     7.9 -met:
    7.10 -
    7.11 -    * Redistributions of source code must retain the above copyright
    7.12 -      notice, this list of conditions and the following disclaimer.
    7.13 -    * Redistributions in binary form must reproduce the above
    7.14 -      copyright notice, this list of conditions and the following
    7.15 -      disclaimer in the documentation and/or other materials provided
    7.16 -      with the distribution.
    7.17 -    * Neither the name of the Hiroshima University nor the names of
    7.18 -      its contributors may be used to endorse or promote products
    7.19 -      derived from this software without specific prior written
    7.20 -      permission.
    7.21 -
    7.22 -THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
    7.23 -"AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
    7.24 -LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
    7.25 -A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
    7.26 -OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
    7.27 -SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
    7.28 -LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
    7.29 -DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
    7.30 -THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
    7.31 -(INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
    7.32 -OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.

     8.1 --- a/src/SFMT/README.txt	Sat Mar 03 10:54:39 2012 -0600
     8.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
     8.3 @@ -1,22 +0,0 @@
     8.4 - =================================================================
     8.5 - SFMT ver. 1.3.3
     8.6 - SIMD oriented Fast Mersenne Twister(SFMT)
     8.7 -
     8.8 - Mutsuo Saito (Hiroshima University) and
     8.9 - Makoto Matsumoto (Hiroshima University)
    8.10 -
    8.11 - Copyright (C) 2006, 2007 Mutsuo Saito, Makoto Matsumoto and Hiroshima
    8.12 - University. All rights reserved.
    8.13 -
    8.14 - The (modified) BSD License is applied to this software, see LICENSE.txt
    8.15 - =================================================================
    8.16 -
    8.17 - To see documents, see html/index.html.
    8.18 -
    8.19 - To make test program, see html/howto-compile.html
    8.20 -
    8.21 - If you want to redistribute and/or change source files, see LICENSE.txt.
    8.22 -
    8.23 - When you change these files and redistribute them, PLEASE write your
    8.24 - e-mail address in redistribution and write to contact YOU first if
    8.25 - users of your changed source encounter troubles.

     9.1 --- a/src/SFMT/SFMT-alti.h	Sat Mar 03 10:54:39 2012 -0600
     9.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
     9.3 @@ -1,156 +0,0 @@
     9.4 -/** 
     9.5 - * @file SFMT-alti.h 
     9.6 - *
     9.7 - * @brief SIMD oriented Fast Mersenne Twister(SFMT)
     9.8 - * pseudorandom number generator
     9.9 - *
    9.10 - * @author Mutsuo Saito (Hiroshima University)
    9.11 - * @author Makoto Matsumoto (Hiroshima University)
    9.12 - *
    9.13 - * Copyright (C) 2007 Mutsuo Saito, Makoto Matsumoto and Hiroshima
    9.14 - * University. All rights reserved.
    9.15 - *
    9.16 - * The new BSD License is applied to this software.
    9.17 - * see LICENSE.txt
    9.18 - */
    9.19 -
    9.20 -#ifndef SFMT_ALTI_H
    9.21 -#define SFMT_ALTI_H
    9.22 -
    9.23 -inline static vector unsigned int vec_recursion(vector unsigned int a,
    9.24 -						vector unsigned int b,
    9.25 -						vector unsigned int c,
    9.26 -						vector unsigned int d)
    9.27 -    ALWAYSINLINE;
    9.28 -
    9.29 -/**
    9.30 - * This function represents the recursion formula in AltiVec and BIG ENDIAN.
    9.31 - * @param a a 128-bit part of the interal state array
    9.32 - * @param b a 128-bit part of the interal state array
    9.33 - * @param c a 128-bit part of the interal state array
    9.34 - * @param d a 128-bit part of the interal state array
    9.35 - * @return output
    9.36 - */
    9.37 -inline static vector unsigned int vec_recursion(vector unsigned int a,
    9.38 -						vector unsigned int b,
    9.39 -						vector unsigned int c,
    9.40 -						vector unsigned int d) {
    9.41 -
    9.42 -    const vector unsigned int sl1 = ALTI_SL1;
    9.43 -    const vector unsigned int sr1 = ALTI_SR1;
    9.44 -#ifdef ONLY64
    9.45 -    const vector unsigned int mask = ALTI_MSK64;
    9.46 -    const vector unsigned char perm_sl = ALTI_SL2_PERM64;
    9.47 -    const vector unsigned char perm_sr = ALTI_SR2_PERM64;
    9.48 -#else
    9.49 -    const vector unsigned int mask = ALTI_MSK;
    9.50 -    const vector unsigned char perm_sl = ALTI_SL2_PERM;
    9.51 -    const vector unsigned char perm_sr = ALTI_SR2_PERM;
    9.52 -#endif
    9.53 -    vector unsigned int v, w, x, y, z;
    9.54 -    x = vec_perm(a, (vector unsigned int)perm_sl, perm_sl);
    9.55 -    v = a;
    9.56 -    y = vec_sr(b, sr1);
    9.57 -    z = vec_perm(c, (vector unsigned int)perm_sr, perm_sr);
    9.58 -    w = vec_sl(d, sl1);
    9.59 -    z = vec_xor(z, w);
    9.60 -    y = vec_and(y, mask);
    9.61 -    v = vec_xor(v, x);
    9.62 -    z = vec_xor(z, y);
    9.63 -    z = vec_xor(z, v);
    9.64 -    return z;
    9.65 -}
    9.66 -
    9.67 -/**
    9.68 - * This function fills the internal state array with pseudorandom
    9.69 - * integers.
    9.70 - */
    9.71 -inline static void gen_rand_all(void) {
    9.72 -    int i;
    9.73 -    vector unsigned int r, r1, r2;
    9.74 -
    9.75 -    r1 = sfmt[N - 2].s;
    9.76 -    r2 = sfmt[N - 1].s;
    9.77 -    for (i = 0; i < N - POS1; i++) {
    9.78 -	r = vec_recursion(sfmt[i].s, sfmt[i + POS1].s, r1, r2);
    9.79 -	sfmt[i].s = r;
    9.80 -	r1 = r2;
    9.81 -	r2 = r;
    9.82 -    }
    9.83 -    for (; i < N; i++) {
    9.84 -	r = vec_recursion(sfmt[i].s, sfmt[i + POS1 - N].s, r1, r2);
    9.85 -	sfmt[i].s = r;
    9.86 -	r1 = r2;
    9.87 -	r2 = r;
    9.88 -    }
    9.89 -}
    9.90 -
    9.91 -/**
    9.92 - * This function fills the user-specified array with pseudorandom
    9.93 - * integers.
    9.94 - *
    9.95 - * @param array an 128-bit array to be filled by pseudorandom numbers.  
    9.96 - * @param size number of 128-bit pesudorandom numbers to be generated.
    9.97 - */
    9.98 -inline static void gen_rand_array(w128_t *array, int size) {
    9.99 -    int i, j;
   9.100 -    vector unsigned int r, r1, r2;
   9.101 -
   9.102 -    r1 = sfmt[N - 2].s;
   9.103 -    r2 = sfmt[N - 1].s;
   9.104 -    for (i = 0; i < N - POS1; i++) {
   9.105 -	r = vec_recursion(sfmt[i].s, sfmt[i + POS1].s, r1, r2);
   9.106 -	array[i].s = r;
   9.107 -	r1 = r2;
   9.108 -	r2 = r;
   9.109 -    }
   9.110 -    for (; i < N; i++) {
   9.111 -	r = vec_recursion(sfmt[i].s, array[i + POS1 - N].s, r1, r2);
   9.112 -	array[i].s = r;
   9.113 -	r1 = r2;
   9.114 -	r2 = r;
   9.115 -    }
   9.116 -    /* main loop */
   9.117 -    for (; i < size - N; i++) {
   9.118 -	r = vec_recursion(array[i - N].s, array[i + POS1 - N].s, r1, r2);
   9.119 -	array[i].s = r;
   9.120 -	r1 = r2;
   9.121 -	r2 = r;
   9.122 -    }
   9.123 -    for (j = 0; j < 2 * N - size; j++) {
   9.124 -	sfmt[j].s = array[j + size - N].s;
   9.125 -    }
   9.126 -    for (; i < size; i++) {
   9.127 -	r = vec_recursion(array[i - N].s, array[i + POS1 - N].s, r1, r2);
   9.128 -	array[i].s = r;
   9.129 -	sfmt[j++].s = r;
   9.130 -	r1 = r2;
   9.131 -	r2 = r;
   9.132 -    }
   9.133 -}
   9.134 -
   9.135 -#ifndef ONLY64
   9.136 -#if defined(__APPLE__)
   9.137 -#define ALTI_SWAP (vector unsigned char) \
   9.138 -	(4, 5, 6, 7, 0, 1, 2, 3, 12, 13, 14, 15, 8, 9, 10, 11)
   9.139 -#else
   9.140 -#define ALTI_SWAP {4, 5, 6, 7, 0, 1, 2, 3, 12, 13, 14, 15, 8, 9, 10, 11}
   9.141 -#endif
   9.142 -/**
   9.143 - * This function swaps high and low 32-bit of 64-bit integers in user
   9.144 - * specified array.
   9.145 - *
   9.146 - * @param array an 128-bit array to be swaped.
   9.147 - * @param size size of 128-bit array.
   9.148 - */
   9.149 -inline static void swap(w128_t *array, int size) {
   9.150 -    int i;
   9.151 -    const vector unsigned char perm = ALTI_SWAP;
   9.152 -
   9.153 -    for (i = 0; i < size; i++) {
   9.154 -	array[i].s = vec_perm(array[i].s, (vector unsigned int)perm, perm);
   9.155 -    }
   9.156 -}
   9.157 -#endif
   9.158 -
   9.159 -#endif

    10.1 --- a/src/SFMT/SFMT-params.h	Sat Mar 03 10:54:39 2012 -0600
    10.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
    10.3 @@ -1,97 +0,0 @@
    10.4 -#ifndef SFMT_PARAMS_H
    10.5 -#define SFMT_PARAMS_H
    10.6 -
    10.7 -#if !defined(MEXP)
    10.8 -#ifdef __GNUC__
    10.9 -  #warning "MEXP is not defined. I assume MEXP is 19937."
   10.10 -#endif
   10.11 -  #define MEXP 19937
   10.12 -#endif
   10.13 -/*-----------------
   10.14 -  BASIC DEFINITIONS
   10.15 -  -----------------*/
   10.16 -/** Mersenne Exponent. The period of the sequence 
   10.17 - *  is a multiple of 2^MEXP-1.
   10.18 - * #define MEXP 19937 */
   10.19 -/** SFMT generator has an internal state array of 128-bit integers,
   10.20 - * and N is its size. */
   10.21 -#define N (MEXP / 128 + 1)
   10.22 -/** N32 is the size of internal state array when regarded as an array
   10.23 - * of 32-bit integers.*/
   10.24 -#define N32 (N * 4)
   10.25 -/** N64 is the size of internal state array when regarded as an array
   10.26 - * of 64-bit integers.*/
   10.27 -#define N64 (N * 2)
   10.28 -
   10.29 -/*----------------------
   10.30 -  the parameters of SFMT
   10.31 -  following definitions are in paramsXXXX.h file.
   10.32 -  ----------------------*/
   10.33 -/** the pick up position of the array.
   10.34 -#define POS1 122 
   10.35 -*/
   10.36 -
   10.37 -/** the parameter of shift left as four 32-bit registers.
   10.38 -#define SL1 18
   10.39 - */
   10.40 -
   10.41 -/** the parameter of shift left as one 128-bit register. 
   10.42 - * The 128-bit integer is shifted by (SL2 * 8) bits. 
   10.43 -#define SL2 1 
   10.44 -*/
   10.45 -
   10.46 -/** the parameter of shift right as four 32-bit registers.
   10.47 -#define SR1 11
   10.48 -*/
   10.49 -
   10.50 -/** the parameter of shift right as one 128-bit register. 
   10.51 - * The 128-bit integer is shifted by (SL2 * 8) bits. 
   10.52 -#define SR2 1 
   10.53 -*/
   10.54 -
   10.55 -/** A bitmask, used in the recursion.  These parameters are introduced
   10.56 - * to break symmetry of SIMD.
   10.57 -#define MSK1 0xdfffffefU
   10.58 -#define MSK2 0xddfecb7fU
   10.59 -#define MSK3 0xbffaffffU
   10.60 -#define MSK4 0xbffffff6U 
   10.61 -*/
   10.62 -
   10.63 -/** These definitions are part of a 128-bit period certification vector.
   10.64 -#define PARITY1	0x00000001U
   10.65 -#define PARITY2	0x00000000U
   10.66 -#define PARITY3	0x00000000U
   10.67 -#define PARITY4	0xc98e126aU
   10.68 -*/
   10.69 -
   10.70 -#if MEXP == 607
   10.71 -  #include "SFMT-params607.h"
   10.72 -#elif MEXP == 1279
   10.73 -  #include "SFMT-params1279.h"
   10.74 -#elif MEXP == 2281
   10.75 -  #include "SFMT-params2281.h"
   10.76 -#elif MEXP == 4253
   10.77 -  #include "SFMT-params4253.h"
   10.78 -#elif MEXP == 11213
   10.79 -  #include "SFMT-params11213.h"
   10.80 -#elif MEXP == 19937
   10.81 -  #include "SFMT-params19937.h"
   10.82 -#elif MEXP == 44497
   10.83 -  #include "SFMT-params44497.h"
   10.84 -#elif MEXP == 86243
   10.85 -  #include "SFMT-params86243.h"
   10.86 -#elif MEXP == 132049
   10.87 -  #include "SFMT-params132049.h"
   10.88 -#elif MEXP == 216091
   10.89 -  #include "SFMT-params216091.h"
   10.90 -#else
   10.91 -#ifdef __GNUC__
   10.92 -  #error "MEXP is not valid."
   10.93 -  #undef MEXP
   10.94 -#else
   10.95 -  #undef MEXP
   10.96 -#endif
   10.97 -
   10.98 -#endif
   10.99 -
  10.100 -#endif /* SFMT_PARAMS_H */

    11.1 --- a/src/SFMT/SFMT-params11213.h	Sat Mar 03 10:54:39 2012 -0600
    11.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
    11.3 @@ -1,46 +0,0 @@
    11.4 -#ifndef SFMT_PARAMS11213_H
    11.5 -#define SFMT_PARAMS11213_H
    11.6 -
    11.7 -#define POS1	68
    11.8 -#define SL1	14
    11.9 -#define SL2	3
   11.10 -#define SR1	7
   11.11 -#define SR2	3
   11.12 -#define MSK1	0xeffff7fbU
   11.13 -#define MSK2	0xffffffefU
   11.14 -#define MSK3	0xdfdfbfffU
   11.15 -#define MSK4	0x7fffdbfdU
   11.16 -#define PARITY1	0x00000001U
   11.17 -#define PARITY2	0x00000000U
   11.18 -#define PARITY3	0xe8148000U
   11.19 -#define PARITY4	0xd0c7afa3U
   11.20 -
   11.21 -
   11.22 -/* PARAMETERS FOR ALTIVEC */
   11.23 -#if defined(__APPLE__)	/* For OSX */
   11.24 -    #define ALTI_SL1	(vector unsigned int)(SL1, SL1, SL1, SL1)
   11.25 -    #define ALTI_SR1	(vector unsigned int)(SR1, SR1, SR1, SR1)
   11.26 -    #define ALTI_MSK	(vector unsigned int)(MSK1, MSK2, MSK3, MSK4)
   11.27 -    #define ALTI_MSK64 \
   11.28 -	(vector unsigned int)(MSK2, MSK1, MSK4, MSK3)
   11.29 -    #define ALTI_SL2_PERM \
   11.30 -	(vector unsigned char)(3,21,21,21,7,0,1,2,11,4,5,6,15,8,9,10)
   11.31 -    #define ALTI_SL2_PERM64 \
   11.32 -	(vector unsigned char)(3,4,5,6,7,29,29,29,11,12,13,14,15,0,1,2)
   11.33 -    #define ALTI_SR2_PERM \
   11.34 -	(vector unsigned char)(5,6,7,0,9,10,11,4,13,14,15,8,19,19,19,12)
   11.35 -    #define ALTI_SR2_PERM64 \
   11.36 -	(vector unsigned char)(13,14,15,0,1,2,3,4,19,19,19,8,9,10,11,12)
   11.37 -#else	/* For OTHER OSs(Linux?) */
   11.38 -    #define ALTI_SL1	{SL1, SL1, SL1, SL1}
   11.39 -    #define ALTI_SR1	{SR1, SR1, SR1, SR1}
   11.40 -    #define ALTI_MSK	{MSK1, MSK2, MSK3, MSK4}
   11.41 -    #define ALTI_MSK64	{MSK2, MSK1, MSK4, MSK3}
   11.42 -    #define ALTI_SL2_PERM	{3,21,21,21,7,0,1,2,11,4,5,6,15,8,9,10}
   11.43 -    #define ALTI_SL2_PERM64	{3,4,5,6,7,29,29,29,11,12,13,14,15,0,1,2}
   11.44 -    #define ALTI_SR2_PERM	{5,6,7,0,9,10,11,4,13,14,15,8,19,19,19,12}
   11.45 -    #define ALTI_SR2_PERM64	{13,14,15,0,1,2,3,4,19,19,19,8,9,10,11,12}
   11.46 -#endif	/* For OSX */
   11.47 -#define IDSTR	"SFMT-11213:68-14-3-7-3:effff7fb-ffffffef-dfdfbfff-7fffdbfd"
   11.48 -
   11.49 -#endif /* SFMT_PARAMS11213_H */

    12.1 --- a/src/SFMT/SFMT-params1279.h	Sat Mar 03 10:54:39 2012 -0600
    12.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
    12.3 @@ -1,46 +0,0 @@
    12.4 -#ifndef SFMT_PARAMS1279_H
    12.5 -#define SFMT_PARAMS1279_H
    12.6 -
    12.7 -#define POS1	7
    12.8 -#define SL1	14
    12.9 -#define SL2	3
   12.10 -#define SR1	5
   12.11 -#define SR2	1
   12.12 -#define MSK1	0xf7fefffdU
   12.13 -#define MSK2	0x7fefcfffU
   12.14 -#define MSK3	0xaff3ef3fU
   12.15 -#define MSK4	0xb5ffff7fU
   12.16 -#define PARITY1	0x00000001U
   12.17 -#define PARITY2	0x00000000U
   12.18 -#define PARITY3	0x00000000U
   12.19 -#define PARITY4	0x20000000U
   12.20 -
   12.21 -
   12.22 -/* PARAMETERS FOR ALTIVEC */
   12.23 -#if defined(__APPLE__)	/* For OSX */
   12.24 -    #define ALTI_SL1	(vector unsigned int)(SL1, SL1, SL1, SL1)
   12.25 -    #define ALTI_SR1	(vector unsigned int)(SR1, SR1, SR1, SR1)
   12.26 -    #define ALTI_MSK	(vector unsigned int)(MSK1, MSK2, MSK3, MSK4)
   12.27 -    #define ALTI_MSK64 \
   12.28 -	(vector unsigned int)(MSK2, MSK1, MSK4, MSK3)
   12.29 -    #define ALTI_SL2_PERM \
   12.30 -	(vector unsigned char)(3,21,21,21,7,0,1,2,11,4,5,6,15,8,9,10)
   12.31 -    #define ALTI_SL2_PERM64 \
   12.32 -	(vector unsigned char)(3,4,5,6,7,29,29,29,11,12,13,14,15,0,1,2)
   12.33 -    #define ALTI_SR2_PERM \
   12.34 -	(vector unsigned char)(7,0,1,2,11,4,5,6,15,8,9,10,17,12,13,14)
   12.35 -    #define ALTI_SR2_PERM64 \
   12.36 -	(vector unsigned char)(15,0,1,2,3,4,5,6,17,8,9,10,11,12,13,14)
   12.37 -#else	/* For OTHER OSs(Linux?) */
   12.38 -    #define ALTI_SL1	{SL1, SL1, SL1, SL1}
   12.39 -    #define ALTI_SR1	{SR1, SR1, SR1, SR1}
   12.40 -    #define ALTI_MSK	{MSK1, MSK2, MSK3, MSK4}
   12.41 -    #define ALTI_MSK64	{MSK2, MSK1, MSK4, MSK3}
   12.42 -    #define ALTI_SL2_PERM	{3,21,21,21,7,0,1,2,11,4,5,6,15,8,9,10}
   12.43 -    #define ALTI_SL2_PERM64	{3,4,5,6,7,29,29,29,11,12,13,14,15,0,1,2}
   12.44 -    #define ALTI_SR2_PERM	{7,0,1,2,11,4,5,6,15,8,9,10,17,12,13,14}
   12.45 -    #define ALTI_SR2_PERM64	{15,0,1,2,3,4,5,6,17,8,9,10,11,12,13,14}
   12.46 -#endif	/* For OSX */
   12.47 -#define IDSTR	"SFMT-1279:7-14-3-5-1:f7fefffd-7fefcfff-aff3ef3f-b5ffff7f"
   12.48 -
   12.49 -#endif /* SFMT_PARAMS1279_H */

    13.1 --- a/src/SFMT/SFMT-params132049.h	Sat Mar 03 10:54:39 2012 -0600
    13.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
    13.3 @@ -1,46 +0,0 @@
    13.4 -#ifndef SFMT_PARAMS132049_H
    13.5 -#define SFMT_PARAMS132049_H
    13.6 -
    13.7 -#define POS1	110
    13.8 -#define SL1	19
    13.9 -#define SL2	1
   13.10 -#define SR1	21
   13.11 -#define SR2	1
   13.12 -#define MSK1	0xffffbb5fU
   13.13 -#define MSK2	0xfb6ebf95U
   13.14 -#define MSK3	0xfffefffaU
   13.15 -#define MSK4	0xcff77fffU
   13.16 -#define PARITY1	0x00000001U
   13.17 -#define PARITY2	0x00000000U
   13.18 -#define PARITY3	0xcb520000U
   13.19 -#define PARITY4	0xc7e91c7dU
   13.20 -
   13.21 -
   13.22 -/* PARAMETERS FOR ALTIVEC */
   13.23 -#if defined(__APPLE__)	/* For OSX */
   13.24 -    #define ALTI_SL1	(vector unsigned int)(SL1, SL1, SL1, SL1)
   13.25 -    #define ALTI_SR1	(vector unsigned int)(SR1, SR1, SR1, SR1)
   13.26 -    #define ALTI_MSK	(vector unsigned int)(MSK1, MSK2, MSK3, MSK4)
   13.27 -    #define ALTI_MSK64 \
   13.28 -	(vector unsigned int)(MSK2, MSK1, MSK4, MSK3)
   13.29 -    #define ALTI_SL2_PERM \
   13.30 -	(vector unsigned char)(1,2,3,23,5,6,7,0,9,10,11,4,13,14,15,8)
   13.31 -    #define ALTI_SL2_PERM64 \
   13.32 -	(vector unsigned char)(1,2,3,4,5,6,7,31,9,10,11,12,13,14,15,0)
   13.33 -    #define ALTI_SR2_PERM \
   13.34 -	(vector unsigned char)(7,0,1,2,11,4,5,6,15,8,9,10,17,12,13,14)
   13.35 -    #define ALTI_SR2_PERM64 \
   13.36 -	(vector unsigned char)(15,0,1,2,3,4,5,6,17,8,9,10,11,12,13,14)
   13.37 -#else	/* For OTHER OSs(Linux?) */
   13.38 -    #define ALTI_SL1	{SL1, SL1, SL1, SL1}
   13.39 -    #define ALTI_SR1	{SR1, SR1, SR1, SR1}
   13.40 -    #define ALTI_MSK	{MSK1, MSK2, MSK3, MSK4}
   13.41 -    #define ALTI_MSK64	{MSK2, MSK1, MSK4, MSK3}
   13.42 -    #define ALTI_SL2_PERM	{1,2,3,23,5,6,7,0,9,10,11,4,13,14,15,8}
   13.43 -    #define ALTI_SL2_PERM64	{1,2,3,4,5,6,7,31,9,10,11,12,13,14,15,0}
   13.44 -    #define ALTI_SR2_PERM	{7,0,1,2,11,4,5,6,15,8,9,10,17,12,13,14}
   13.45 -    #define ALTI_SR2_PERM64	{15,0,1,2,3,4,5,6,17,8,9,10,11,12,13,14}
   13.46 -#endif	/* For OSX */
   13.47 -#define IDSTR	"SFMT-132049:110-19-1-21-1:ffffbb5f-fb6ebf95-fffefffa-cff77fff"
   13.48 -
   13.49 -#endif /* SFMT_PARAMS132049_H */

    14.1 --- a/src/SFMT/SFMT-params19937.h	Sat Mar 03 10:54:39 2012 -0600
    14.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
    14.3 @@ -1,46 +0,0 @@
    14.4 -#ifndef SFMT_PARAMS19937_H
    14.5 -#define SFMT_PARAMS19937_H
    14.6 -
    14.7 -#define POS1	122
    14.8 -#define SL1	18
    14.9 -#define SL2	1
   14.10 -#define SR1	11
   14.11 -#define SR2	1
   14.12 -#define MSK1	0xdfffffefU
   14.13 -#define MSK2	0xddfecb7fU
   14.14 -#define MSK3	0xbffaffffU
   14.15 -#define MSK4	0xbffffff6U
   14.16 -#define PARITY1	0x00000001U
   14.17 -#define PARITY2	0x00000000U
   14.18 -#define PARITY3	0x00000000U
   14.19 -#define PARITY4	0x13c9e684U
   14.20 -
   14.21 -
   14.22 -/* PARAMETERS FOR ALTIVEC */
   14.23 -#if defined(__APPLE__)	/* For OSX */
   14.24 -    #define ALTI_SL1	(vector unsigned int)(SL1, SL1, SL1, SL1)
   14.25 -    #define ALTI_SR1	(vector unsigned int)(SR1, SR1, SR1, SR1)
   14.26 -    #define ALTI_MSK	(vector unsigned int)(MSK1, MSK2, MSK3, MSK4)
   14.27 -    #define ALTI_MSK64 \
   14.28 -	(vector unsigned int)(MSK2, MSK1, MSK4, MSK3)
   14.29 -    #define ALTI_SL2_PERM \
   14.30 -	(vector unsigned char)(1,2,3,23,5,6,7,0,9,10,11,4,13,14,15,8)
   14.31 -    #define ALTI_SL2_PERM64 \
   14.32 -	(vector unsigned char)(1,2,3,4,5,6,7,31,9,10,11,12,13,14,15,0)
   14.33 -    #define ALTI_SR2_PERM \
   14.34 -	(vector unsigned char)(7,0,1,2,11,4,5,6,15,8,9,10,17,12,13,14)
   14.35 -    #define ALTI_SR2_PERM64 \
   14.36 -	(vector unsigned char)(15,0,1,2,3,4,5,6,17,8,9,10,11,12,13,14)
   14.37 -#else	/* For OTHER OSs(Linux?) */
   14.38 -    #define ALTI_SL1	{SL1, SL1, SL1, SL1}
   14.39 -    #define ALTI_SR1	{SR1, SR1, SR1, SR1}
   14.40 -    #define ALTI_MSK	{MSK1, MSK2, MSK3, MSK4}
   14.41 -    #define ALTI_MSK64	{MSK2, MSK1, MSK4, MSK3}
   14.42 -    #define ALTI_SL2_PERM	{1,2,3,23,5,6,7,0,9,10,11,4,13,14,15,8}
   14.43 -    #define ALTI_SL2_PERM64	{1,2,3,4,5,6,7,31,9,10,11,12,13,14,15,0}
   14.44 -    #define ALTI_SR2_PERM	{7,0,1,2,11,4,5,6,15,8,9,10,17,12,13,14}
   14.45 -    #define ALTI_SR2_PERM64	{15,0,1,2,3,4,5,6,17,8,9,10,11,12,13,14}
   14.46 -#endif	/* For OSX */
   14.47 -#define IDSTR	"SFMT-19937:122-18-1-11-1:dfffffef-ddfecb7f-bffaffff-bffffff6"
   14.48 -
   14.49 -#endif /* SFMT_PARAMS19937_H */

    15.1 --- a/src/SFMT/SFMT-params216091.h	Sat Mar 03 10:54:39 2012 -0600
    15.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
    15.3 @@ -1,46 +0,0 @@
    15.4 -#ifndef SFMT_PARAMS216091_H
    15.5 -#define SFMT_PARAMS216091_H
    15.6 -
    15.7 -#define POS1	627
    15.8 -#define SL1	11
    15.9 -#define SL2	3
   15.10 -#define SR1	10
   15.11 -#define SR2	1
   15.12 -#define MSK1	0xbff7bff7U
   15.13 -#define MSK2	0xbfffffffU
   15.14 -#define MSK3	0xbffffa7fU
   15.15 -#define MSK4	0xffddfbfbU
   15.16 -#define PARITY1	0xf8000001U
   15.17 -#define PARITY2	0x89e80709U
   15.18 -#define PARITY3	0x3bd2b64bU
   15.19 -#define PARITY4	0x0c64b1e4U
   15.20 -
   15.21 -
   15.22 -/* PARAMETERS FOR ALTIVEC */
   15.23 -#if defined(__APPLE__)	/* For OSX */
   15.24 -    #define ALTI_SL1	(vector unsigned int)(SL1, SL1, SL1, SL1)
   15.25 -    #define ALTI_SR1	(vector unsigned int)(SR1, SR1, SR1, SR1)
   15.26 -    #define ALTI_MSK	(vector unsigned int)(MSK1, MSK2, MSK3, MSK4)
   15.27 -    #define ALTI_MSK64 \
   15.28 -	(vector unsigned int)(MSK2, MSK1, MSK4, MSK3)
   15.29 -    #define ALTI_SL2_PERM \
   15.30 -	(vector unsigned char)(3,21,21,21,7,0,1,2,11,4,5,6,15,8,9,10)
   15.31 -    #define ALTI_SL2_PERM64 \
   15.32 -	(vector unsigned char)(3,4,5,6,7,29,29,29,11,12,13,14,15,0,1,2)
   15.33 -    #define ALTI_SR2_PERM \
   15.34 -	(vector unsigned char)(7,0,1,2,11,4,5,6,15,8,9,10,17,12,13,14)
   15.35 -    #define ALTI_SR2_PERM64 \
   15.36 -	(vector unsigned char)(15,0,1,2,3,4,5,6,17,8,9,10,11,12,13,14)
   15.37 -#else	/* For OTHER OSs(Linux?) */
   15.38 -    #define ALTI_SL1	{SL1, SL1, SL1, SL1}
   15.39 -    #define ALTI_SR1	{SR1, SR1, SR1, SR1}
   15.40 -    #define ALTI_MSK	{MSK1, MSK2, MSK3, MSK4}
   15.41 -    #define ALTI_MSK64	{MSK2, MSK1, MSK4, MSK3}
   15.42 -    #define ALTI_SL2_PERM	{3,21,21,21,7,0,1,2,11,4,5,6,15,8,9,10}
   15.43 -    #define ALTI_SL2_PERM64	{3,4,5,6,7,29,29,29,11,12,13,14,15,0,1,2}
   15.44 -    #define ALTI_SR2_PERM	{7,0,1,2,11,4,5,6,15,8,9,10,17,12,13,14}
   15.45 -    #define ALTI_SR2_PERM64	{15,0,1,2,3,4,5,6,17,8,9,10,11,12,13,14}
   15.46 -#endif	/* For OSX */
   15.47 -#define IDSTR	"SFMT-216091:627-11-3-10-1:bff7bff7-bfffffff-bffffa7f-ffddfbfb"
   15.48 -
   15.49 -#endif /* SFMT_PARAMS216091_H */

    16.1 --- a/src/SFMT/SFMT-params2281.h	Sat Mar 03 10:54:39 2012 -0600
    16.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
    16.3 @@ -1,46 +0,0 @@
    16.4 -#ifndef SFMT_PARAMS2281_H
    16.5 -#define SFMT_PARAMS2281_H
    16.6 -
    16.7 -#define POS1	12
    16.8 -#define SL1	19
    16.9 -#define SL2	1
   16.10 -#define SR1	5
   16.11 -#define SR2	1
   16.12 -#define MSK1	0xbff7ffbfU
   16.13 -#define MSK2	0xfdfffffeU
   16.14 -#define MSK3	0xf7ffef7fU
   16.15 -#define MSK4	0xf2f7cbbfU
   16.16 -#define PARITY1	0x00000001U
   16.17 -#define PARITY2	0x00000000U
   16.18 -#define PARITY3	0x00000000U
   16.19 -#define PARITY4	0x41dfa600U
   16.20 -
   16.21 -
   16.22 -/* PARAMETERS FOR ALTIVEC */
   16.23 -#if defined(__APPLE__)	/* For OSX */
   16.24 -    #define ALTI_SL1	(vector unsigned int)(SL1, SL1, SL1, SL1)
   16.25 -    #define ALTI_SR1	(vector unsigned int)(SR1, SR1, SR1, SR1)
   16.26 -    #define ALTI_MSK	(vector unsigned int)(MSK1, MSK2, MSK3, MSK4)
   16.27 -    #define ALTI_MSK64 \
   16.28 -	(vector unsigned int)(MSK2, MSK1, MSK4, MSK3)
   16.29 -    #define ALTI_SL2_PERM \
   16.30 -	(vector unsigned char)(1,2,3,23,5,6,7,0,9,10,11,4,13,14,15,8)
   16.31 -    #define ALTI_SL2_PERM64 \
   16.32 -	(vector unsigned char)(1,2,3,4,5,6,7,31,9,10,11,12,13,14,15,0)
   16.33 -    #define ALTI_SR2_PERM \
   16.34 -	(vector unsigned char)(7,0,1,2,11,4,5,6,15,8,9,10,17,12,13,14)
   16.35 -    #define ALTI_SR2_PERM64 \
   16.36 -	(vector unsigned char)(15,0,1,2,3,4,5,6,17,8,9,10,11,12,13,14)
   16.37 -#else	/* For OTHER OSs(Linux?) */
   16.38 -    #define ALTI_SL1	{SL1, SL1, SL1, SL1}
   16.39 -    #define ALTI_SR1	{SR1, SR1, SR1, SR1}
   16.40 -    #define ALTI_MSK	{MSK1, MSK2, MSK3, MSK4}
   16.41 -    #define ALTI_MSK64	{MSK2, MSK1, MSK4, MSK3}
   16.42 -    #define ALTI_SL2_PERM	{1,2,3,23,5,6,7,0,9,10,11,4,13,14,15,8}
   16.43 -    #define ALTI_SL2_PERM64	{1,2,3,4,5,6,7,31,9,10,11,12,13,14,15,0}
   16.44 -    #define ALTI_SR2_PERM	{7,0,1,2,11,4,5,6,15,8,9,10,17,12,13,14}
   16.45 -    #define ALTI_SR2_PERM64	{15,0,1,2,3,4,5,6,17,8,9,10,11,12,13,14}
   16.46 -#endif	/* For OSX */
   16.47 -#define IDSTR	"SFMT-2281:12-19-1-5-1:bff7ffbf-fdfffffe-f7ffef7f-f2f7cbbf"
   16.48 -
   16.49 -#endif /* SFMT_PARAMS2281_H */

    17.1 --- a/src/SFMT/SFMT-params4253.h	Sat Mar 03 10:54:39 2012 -0600
    17.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
    17.3 @@ -1,46 +0,0 @@
    17.4 -#ifndef SFMT_PARAMS4253_H
    17.5 -#define SFMT_PARAMS4253_H
    17.6 -
    17.7 -#define POS1	17
    17.8 -#define SL1	20
    17.9 -#define SL2	1
   17.10 -#define SR1	7
   17.11 -#define SR2	1
   17.12 -#define MSK1	0x9f7bffffU
   17.13 -#define MSK2	0x9fffff5fU
   17.14 -#define MSK3	0x3efffffbU
   17.15 -#define MSK4	0xfffff7bbU
   17.16 -#define PARITY1	0xa8000001U
   17.17 -#define PARITY2	0xaf5390a3U
   17.18 -#define PARITY3	0xb740b3f8U
   17.19 -#define PARITY4	0x6c11486dU
   17.20 -
   17.21 -
   17.22 -/* PARAMETERS FOR ALTIVEC */
   17.23 -#if defined(__APPLE__)	/* For OSX */
   17.24 -    #define ALTI_SL1	(vector unsigned int)(SL1, SL1, SL1, SL1)
   17.25 -    #define ALTI_SR1	(vector unsigned int)(SR1, SR1, SR1, SR1)
   17.26 -    #define ALTI_MSK	(vector unsigned int)(MSK1, MSK2, MSK3, MSK4)
   17.27 -    #define ALTI_MSK64 \
   17.28 -	(vector unsigned int)(MSK2, MSK1, MSK4, MSK3)
   17.29 -    #define ALTI_SL2_PERM \
   17.30 -	(vector unsigned char)(1,2,3,23,5,6,7,0,9,10,11,4,13,14,15,8)
   17.31 -    #define ALTI_SL2_PERM64 \
   17.32 -	(vector unsigned char)(1,2,3,4,5,6,7,31,9,10,11,12,13,14,15,0)
   17.33 -    #define ALTI_SR2_PERM \
   17.34 -	(vector unsigned char)(7,0,1,2,11,4,5,6,15,8,9,10,17,12,13,14)
   17.35 -    #define ALTI_SR2_PERM64 \
   17.36 -	(vector unsigned char)(15,0,1,2,3,4,5,6,17,8,9,10,11,12,13,14)
   17.37 -#else	/* For OTHER OSs(Linux?) */
   17.38 -    #define ALTI_SL1	{SL1, SL1, SL1, SL1}
   17.39 -    #define ALTI_SR1	{SR1, SR1, SR1, SR1}
   17.40 -    #define ALTI_MSK	{MSK1, MSK2, MSK3, MSK4}
   17.41 -    #define ALTI_MSK64	{MSK2, MSK1, MSK4, MSK3}
   17.42 -    #define ALTI_SL2_PERM	{1,2,3,23,5,6,7,0,9,10,11,4,13,14,15,8}
   17.43 -    #define ALTI_SL2_PERM64	{1,2,3,4,5,6,7,31,9,10,11,12,13,14,15,0}
   17.44 -    #define ALTI_SR2_PERM	{7,0,1,2,11,4,5,6,15,8,9,10,17,12,13,14}
   17.45 -    #define ALTI_SR2_PERM64	{15,0,1,2,3,4,5,6,17,8,9,10,11,12,13,14}
   17.46 -#endif	/* For OSX */
   17.47 -#define IDSTR	"SFMT-4253:17-20-1-7-1:9f7bffff-9fffff5f-3efffffb-fffff7bb"
   17.48 -
   17.49 -#endif /* SFMT_PARAMS4253_H */

    18.1 --- a/src/SFMT/SFMT-params44497.h	Sat Mar 03 10:54:39 2012 -0600
    18.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
    18.3 @@ -1,46 +0,0 @@
    18.4 -#ifndef SFMT_PARAMS44497_H
    18.5 -#define SFMT_PARAMS44497_H
    18.6 -
    18.7 -#define POS1	330
    18.8 -#define SL1	5
    18.9 -#define SL2	3
   18.10 -#define SR1	9
   18.11 -#define SR2	3
   18.12 -#define MSK1	0xeffffffbU
   18.13 -#define MSK2	0xdfbebfffU
   18.14 -#define MSK3	0xbfbf7befU
   18.15 -#define MSK4	0x9ffd7bffU
   18.16 -#define PARITY1	0x00000001U
   18.17 -#define PARITY2	0x00000000U
   18.18 -#define PARITY3	0xa3ac4000U
   18.19 -#define PARITY4	0xecc1327aU
   18.20 -
   18.21 -
   18.22 -/* PARAMETERS FOR ALTIVEC */
   18.23 -#if defined(__APPLE__)	/* For OSX */
   18.24 -    #define ALTI_SL1	(vector unsigned int)(SL1, SL1, SL1, SL1)
   18.25 -    #define ALTI_SR1	(vector unsigned int)(SR1, SR1, SR1, SR1)
   18.26 -    #define ALTI_MSK	(vector unsigned int)(MSK1, MSK2, MSK3, MSK4)
   18.27 -    #define ALTI_MSK64 \
   18.28 -	(vector unsigned int)(MSK2, MSK1, MSK4, MSK3)
   18.29 -    #define ALTI_SL2_PERM \
   18.30 -	(vector unsigned char)(3,21,21,21,7,0,1,2,11,4,5,6,15,8,9,10)
   18.31 -    #define ALTI_SL2_PERM64 \
   18.32 -	(vector unsigned char)(3,4,5,6,7,29,29,29,11,12,13,14,15,0,1,2)
   18.33 -    #define ALTI_SR2_PERM \
   18.34 -	(vector unsigned char)(5,6,7,0,9,10,11,4,13,14,15,8,19,19,19,12)
   18.35 -    #define ALTI_SR2_PERM64 \
   18.36 -	(vector unsigned char)(13,14,15,0,1,2,3,4,19,19,19,8,9,10,11,12)
   18.37 -#else	/* For OTHER OSs(Linux?) */
   18.38 -    #define ALTI_SL1	{SL1, SL1, SL1, SL1}
   18.39 -    #define ALTI_SR1	{SR1, SR1, SR1, SR1}
   18.40 -    #define ALTI_MSK	{MSK1, MSK2, MSK3, MSK4}
   18.41 -    #define ALTI_MSK64	{MSK2, MSK1, MSK4, MSK3}
   18.42 -    #define ALTI_SL2_PERM	{3,21,21,21,7,0,1,2,11,4,5,6,15,8,9,10}
   18.43 -    #define ALTI_SL2_PERM64	{3,4,5,6,7,29,29,29,11,12,13,14,15,0,1,2}
   18.44 -    #define ALTI_SR2_PERM	{5,6,7,0,9,10,11,4,13,14,15,8,19,19,19,12}
   18.45 -    #define ALTI_SR2_PERM64	{13,14,15,0,1,2,3,4,19,19,19,8,9,10,11,12}
   18.46 -#endif	/* For OSX */
   18.47 -#define IDSTR	"SFMT-44497:330-5-3-9-3:effffffb-dfbebfff-bfbf7bef-9ffd7bff"
   18.48 -
   18.49 -#endif /* SFMT_PARAMS44497_H */

    19.1 --- a/src/SFMT/SFMT-params607.h	Sat Mar 03 10:54:39 2012 -0600
    19.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
    19.3 @@ -1,46 +0,0 @@
    19.4 -#ifndef SFMT_PARAMS607_H
    19.5 -#define SFMT_PARAMS607_H
    19.6 -
    19.7 -#define POS1	2
    19.8 -#define SL1	15
    19.9 -#define SL2	3
   19.10 -#define SR1	13
   19.11 -#define SR2	3
   19.12 -#define MSK1	0xfdff37ffU
   19.13 -#define MSK2	0xef7f3f7dU
   19.14 -#define MSK3	0xff777b7dU
   19.15 -#define MSK4	0x7ff7fb2fU
   19.16 -#define PARITY1	0x00000001U
   19.17 -#define PARITY2	0x00000000U
   19.18 -#define PARITY3	0x00000000U
   19.19 -#define PARITY4	0x5986f054U
   19.20 -
   19.21 -
   19.22 -/* PARAMETERS FOR ALTIVEC */
   19.23 -#if defined(__APPLE__)	/* For OSX */
   19.24 -    #define ALTI_SL1	(vector unsigned int)(SL1, SL1, SL1, SL1)
   19.25 -    #define ALTI_SR1	(vector unsigned int)(SR1, SR1, SR1, SR1)
   19.26 -    #define ALTI_MSK	(vector unsigned int)(MSK1, MSK2, MSK3, MSK4)
   19.27 -    #define ALTI_MSK64 \
   19.28 -	(vector unsigned int)(MSK2, MSK1, MSK4, MSK3)
   19.29 -    #define ALTI_SL2_PERM \
   19.30 -	(vector unsigned char)(3,21,21,21,7,0,1,2,11,4,5,6,15,8,9,10)
   19.31 -    #define ALTI_SL2_PERM64 \
   19.32 -	(vector unsigned char)(3,4,5,6,7,29,29,29,11,12,13,14,15,0,1,2)
   19.33 -    #define ALTI_SR2_PERM \
   19.34 -	(vector unsigned char)(5,6,7,0,9,10,11,4,13,14,15,8,19,19,19,12)
   19.35 -    #define ALTI_SR2_PERM64 \
   19.36 -	(vector unsigned char)(13,14,15,0,1,2,3,4,19,19,19,8,9,10,11,12)
   19.37 -#else	/* For OTHER OSs(Linux?) */
   19.38 -    #define ALTI_SL1	{SL1, SL1, SL1, SL1}
   19.39 -    #define ALTI_SR1	{SR1, SR1, SR1, SR1}
   19.40 -    #define ALTI_MSK	{MSK1, MSK2, MSK3, MSK4}
   19.41 -    #define ALTI_MSK64	{MSK2, MSK1, MSK4, MSK3}
   19.42 -    #define ALTI_SL2_PERM	{3,21,21,21,7,0,1,2,11,4,5,6,15,8,9,10}
   19.43 -    #define ALTI_SL2_PERM64	{3,4,5,6,7,29,29,29,11,12,13,14,15,0,1,2}
   19.44 -    #define ALTI_SR2_PERM	{5,6,7,0,9,10,11,4,13,14,15,8,19,19,19,12}
   19.45 -    #define ALTI_SR2_PERM64	{13,14,15,0,1,2,3,4,19,19,19,8,9,10,11,12}
   19.46 -#endif	/* For OSX */
   19.47 -#define IDSTR	"SFMT-607:2-15-3-13-3:fdff37ff-ef7f3f7d-ff777b7d-7ff7fb2f"
   19.48 -
   19.49 -#endif /* SFMT_PARAMS607_H */

    20.1 --- a/src/SFMT/SFMT-params86243.h	Sat Mar 03 10:54:39 2012 -0600
    20.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
    20.3 @@ -1,46 +0,0 @@
    20.4 -#ifndef SFMT_PARAMS86243_H
    20.5 -#define SFMT_PARAMS86243_H
    20.6 -
    20.7 -#define POS1	366
    20.8 -#define SL1	6
    20.9 -#define SL2	7
   20.10 -#define SR1	19
   20.11 -#define SR2	1
   20.12 -#define MSK1	0xfdbffbffU
   20.13 -#define MSK2	0xbff7ff3fU
   20.14 -#define MSK3	0xfd77efffU
   20.15 -#define MSK4	0xbf9ff3ffU
   20.16 -#define PARITY1	0x00000001U
   20.17 -#define PARITY2	0x00000000U
   20.18 -#define PARITY3	0x00000000U
   20.19 -#define PARITY4	0xe9528d85U
   20.20 -
   20.21 -
   20.22 -/* PARAMETERS FOR ALTIVEC */
   20.23 -#if defined(__APPLE__)	/* For OSX */
   20.24 -    #define ALTI_SL1	(vector unsigned int)(SL1, SL1, SL1, SL1)
   20.25 -    #define ALTI_SR1	(vector unsigned int)(SR1, SR1, SR1, SR1)
   20.26 -    #define ALTI_MSK	(vector unsigned int)(MSK1, MSK2, MSK3, MSK4)
   20.27 -    #define ALTI_MSK64 \
   20.28 -	(vector unsigned int)(MSK2, MSK1, MSK4, MSK3)
   20.29 -    #define ALTI_SL2_PERM \
   20.30 -	(vector unsigned char)(25,25,25,25,3,25,25,25,7,0,1,2,11,4,5,6)
   20.31 -    #define ALTI_SL2_PERM64 \
   20.32 -	(vector unsigned char)(7,25,25,25,25,25,25,25,15,0,1,2,3,4,5,6)
   20.33 -    #define ALTI_SR2_PERM \
   20.34 -	(vector unsigned char)(7,0,1,2,11,4,5,6,15,8,9,10,17,12,13,14)
   20.35 -    #define ALTI_SR2_PERM64 \
   20.36 -	(vector unsigned char)(15,0,1,2,3,4,5,6,17,8,9,10,11,12,13,14)
   20.37 -#else	/* For OTHER OSs(Linux?) */
   20.38 -    #define ALTI_SL1	{SL1, SL1, SL1, SL1}
   20.39 -    #define ALTI_SR1	{SR1, SR1, SR1, SR1}
   20.40 -    #define ALTI_MSK	{MSK1, MSK2, MSK3, MSK4}
   20.41 -    #define ALTI_MSK64	{MSK2, MSK1, MSK4, MSK3}
   20.42 -    #define ALTI_SL2_PERM	{25,25,25,25,3,25,25,25,7,0,1,2,11,4,5,6}
   20.43 -    #define ALTI_SL2_PERM64	{7,25,25,25,25,25,25,25,15,0,1,2,3,4,5,6}
   20.44 -    #define ALTI_SR2_PERM	{7,0,1,2,11,4,5,6,15,8,9,10,17,12,13,14}
   20.45 -    #define ALTI_SR2_PERM64	{15,0,1,2,3,4,5,6,17,8,9,10,11,12,13,14}
   20.46 -#endif	/* For OSX */
   20.47 -#define IDSTR	"SFMT-86243:366-6-7-19-1:fdbffbff-bff7ff3f-fd77efff-bf9ff3ff"
   20.48 -
   20.49 -#endif /* SFMT_PARAMS86243_H */

    21.1 --- a/src/SFMT/SFMT-sse2.h	Sat Mar 03 10:54:39 2012 -0600
    21.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
    21.3 @@ -1,121 +0,0 @@
    21.4 -/** 
    21.5 - * @file  SFMT-sse2.h
    21.6 - * @brief SIMD oriented Fast Mersenne Twister(SFMT) for Intel SSE2
    21.7 - *
    21.8 - * @author Mutsuo Saito (Hiroshima University)
    21.9 - * @author Makoto Matsumoto (Hiroshima University)
   21.10 - *
   21.11 - * @note We assume LITTLE ENDIAN in this file
   21.12 - *
   21.13 - * Copyright (C) 2006, 2007 Mutsuo Saito, Makoto Matsumoto and Hiroshima
   21.14 - * University. All rights reserved.
   21.15 - *
   21.16 - * The new BSD License is applied to this software, see LICENSE.txt
   21.17 - */
   21.18 -
   21.19 -#ifndef SFMT_SSE2_H
   21.20 -#define SFMT_SSE2_H
   21.21 -
   21.22 -PRE_ALWAYS static __m128i mm_recursion(__m128i *a, __m128i *b, __m128i c,
   21.23 -				   __m128i d, __m128i mask) ALWAYSINLINE;
   21.24 -
   21.25 -/**
   21.26 - * This function represents the recursion formula.
   21.27 - * @param a a 128-bit part of the interal state array
   21.28 - * @param b a 128-bit part of the interal state array
   21.29 - * @param c a 128-bit part of the interal state array
   21.30 - * @param d a 128-bit part of the interal state array
   21.31 - * @param mask 128-bit mask
   21.32 - * @return output
   21.33 - */
   21.34 -PRE_ALWAYS static __m128i mm_recursion(__m128i *a, __m128i *b, 
   21.35 -				   __m128i c, __m128i d, __m128i mask) {
   21.36 -    __m128i v, x, y, z;
   21.37 -    
   21.38 -    x = _mm_load_si128(a);
   21.39 -    y = _mm_srli_epi32(*b, SR1);
   21.40 -    z = _mm_srli_si128(c, SR2);
   21.41 -    v = _mm_slli_epi32(d, SL1);
   21.42 -    z = _mm_xor_si128(z, x);
   21.43 -    z = _mm_xor_si128(z, v);
   21.44 -    x = _mm_slli_si128(x, SL2);
   21.45 -    y = _mm_and_si128(y, mask);
   21.46 -    z = _mm_xor_si128(z, x);
   21.47 -    z = _mm_xor_si128(z, y);
   21.48 -    return z;
   21.49 -}
   21.50 -
   21.51 -/**
   21.52 - * This function fills the internal state array with pseudorandom
   21.53 - * integers.
   21.54 - */
   21.55 -inline static void gen_rand_all(void) {
   21.56 -    int i;
   21.57 -    __m128i r, r1, r2, mask;
   21.58 -    mask = _mm_set_epi32(MSK4, MSK3, MSK2, MSK1);
   21.59 -
   21.60 -    r1 = _mm_load_si128(&sfmt[N - 2].si);
   21.61 -    r2 = _mm_load_si128(&sfmt[N - 1].si);
   21.62 -    for (i = 0; i < N - POS1; i++) {
   21.63 -	r = mm_recursion(&sfmt[i].si, &sfmt[i + POS1].si, r1, r2, mask);
   21.64 -	_mm_store_si128(&sfmt[i].si, r);
   21.65 -	r1 = r2;
   21.66 -	r2 = r;
   21.67 -    }
   21.68 -    for (; i < N; i++) {
   21.69 -	r = mm_recursion(&sfmt[i].si, &sfmt[i + POS1 - N].si, r1, r2, mask);
   21.70 -	_mm_store_si128(&sfmt[i].si, r);
   21.71 -	r1 = r2;
   21.72 -	r2 = r;
   21.73 -    }
   21.74 -}
   21.75 -
   21.76 -/**
   21.77 - * This function fills the user-specified array with pseudorandom
   21.78 - * integers.
   21.79 - *
   21.80 - * @param array an 128-bit array to be filled by pseudorandom numbers.  
   21.81 - * @param size number of 128-bit pesudorandom numbers to be generated.
   21.82 - */
   21.83 -inline static void gen_rand_array(w128_t *array, int size) {
   21.84 -    int i, j;
   21.85 -    __m128i r, r1, r2, mask;
   21.86 -    mask = _mm_set_epi32(MSK4, MSK3, MSK2, MSK1);
   21.87 -
   21.88 -    r1 = _mm_load_si128(&sfmt[N - 2].si);
   21.89 -    r2 = _mm_load_si128(&sfmt[N - 1].si);
   21.90 -    for (i = 0; i < N - POS1; i++) {
   21.91 -	r = mm_recursion(&sfmt[i].si, &sfmt[i + POS1].si, r1, r2, mask);
   21.92 -	_mm_store_si128(&array[i].si, r);
   21.93 -	r1 = r2;
   21.94 -	r2 = r;
   21.95 -    }
   21.96 -    for (; i < N; i++) {
   21.97 -	r = mm_recursion(&sfmt[i].si, &array[i + POS1 - N].si, r1, r2, mask);
   21.98 -	_mm_store_si128(&array[i].si, r);
   21.99 -	r1 = r2;
  21.100 -	r2 = r;
  21.101 -    }
  21.102 -    /* main loop */
  21.103 -    for (; i < size - N; i++) {
  21.104 -	r = mm_recursion(&array[i - N].si, &array[i + POS1 - N].si, r1, r2,
  21.105 -			 mask);
  21.106 -	_mm_store_si128(&array[i].si, r);
  21.107 -	r1 = r2;
  21.108 -	r2 = r;
  21.109 -    }
  21.110 -    for (j = 0; j < 2 * N - size; j++) {
  21.111 -	r = _mm_load_si128(&array[j + size - N].si);
  21.112 -	_mm_store_si128(&sfmt[j].si, r);
  21.113 -    }
  21.114 -    for (; i < size; i++) {
  21.115 -	r = mm_recursion(&array[i - N].si, &array[i + POS1 - N].si, r1, r2,
  21.116 -			 mask);
  21.117 -	_mm_store_si128(&array[i].si, r);
  21.118 -	_mm_store_si128(&sfmt[j++].si, r);
  21.119 -	r1 = r2;
  21.120 -	r2 = r;
  21.121 -    }
  21.122 -}
  21.123 -
  21.124 -#endif

    22.1 --- a/src/SFMT/SFMT.c	Sat Mar 03 10:54:39 2012 -0600
    22.2 +++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
    22.3 @@ -1,620 +0,0 @@
    22.4 -/** 
    22.5 - * @file  SFMT.c
    22.6 - * @brief SIMD oriented Fast Mersenne Twister(SFMT)
    22.7 - *
    22.8 - * @author Mutsuo Saito (Hiroshima University)
    22.9 - * @author Makoto Matsumoto (Hiroshima University)
   22.10 - *
   22.11 - * Copyright (C) 2006,2007 Mutsuo Saito, Makoto Matsumoto and Hiroshima
   22.12 - * University. All rights reserved.
   22.13 - *
   22.14 - * The new BSD License is applied to this software, see LICENSE.txt
   22.15 - */
   22.16 -#include <string.h>
   22.17 -#include <assert.h>
   22.18 -#include "SFMT.h"
   22.19 -#include "SFMT-params.h"
   22.20 -
   22.21 -#if defined(__BIG_ENDIAN__) && !defined(__amd64) && !defined(BIG_ENDIAN64)
   22.22 -#define BIG_ENDIAN64 1
   22.23 -#endif
   22.24 -#if defined(HAVE_ALTIVEC) && !defined(BIG_ENDIAN64)
   22.25 -#define BIG_ENDIAN64 1
   22.26 -#endif
   22.27 -#if defined(ONLY64) && !defined(BIG_ENDIAN64)
   22.28 -  #if defined(__GNUC__)
   22.29 -    #error "-DONLY64 must be specified with -DBIG_ENDIAN64"
   22.30 -  #endif
   22.31 -#undef ONLY64
   22.32 -#endif
   22.33 -/*------------------------------------------------------
   22.34 -  128-bit SIMD data type for Altivec, SSE2 or standard C
   22.35 -  ------------------------------------------------------*/
   22.36 -#if defined(HAVE_ALTIVEC)
   22.37 -  #if !defined(__APPLE__)
   22.38 -    #include <altivec.h>
   22.39 -  #endif
   22.40 -/** 128-bit data structure */
   22.41 -union W128_T {
   22.42 -    vector unsigned int s;
   22.43 -    uint32_t u[4];
   22.44 -};
   22.45 -/** 128-bit data type */
   22.46 -typedef union W128_T w128_t;
   22.47 -
   22.48 -#elif defined(HAVE_SSE2)
   22.49 -  #include <emmintrin.h>
   22.50 -
   22.51 -/** 128-bit data structure */
   22.52 -union W128_T {
   22.53 -    __m128i si;
   22.54 -    uint32_t u[4];
   22.55 -};
   22.56 -/** 128-bit data type */
   22.57 -typedef union W128_T w128_t;
   22.58 -
   22.59 -#else
   22.60 -
   22.61 -/** 128-bit data structure */
   22.62 -struct W128_T {
   22.63 -    uint32_t u[4];
   22.64 -};
   22.65 -/** 128-bit data type */
   22.66 -typedef struct W128_T w128_t;
   22.67 -
   22.68 -#endif
   22.69 -
   22.70 -/*--------------------------------------
   22.71 -  FILE GLOBAL VARIABLES
   22.72 -  internal state, index counter and flag 
   22.73 -  --------------------------------------*/
   22.74 -/** the 128-bit internal state array */
   22.75 -static w128_t sfmt[N];
   22.76 -/** the 32bit integer pointer to the 128-bit internal state array */
   22.77 -static uint32_t *psfmt32 = &sfmt[0].u[0];
   22.78 -#if !defined(BIG_ENDIAN64) || defined(ONLY64)
   22.79 -/** the 64bit integer pointer to the 128-bit internal state array */
   22.80 -static uint64_t *psfmt64 = (uint64_t *)&sfmt[0].u[0];
   22.81 -#endif
   22.82 -/** index counter to the 32-bit internal state array */
   22.83 -static int idx;
   22.84 -/** a flag: it is 0 if and only if the internal state is not yet
   22.85 - * initialized. */
   22.86 -static int initialized = 0;
   22.87 -/** a parity check vector which certificate the period of 2^{MEXP} */
   22.88 -static uint32_t parity[4] = {PARITY1, PARITY2, PARITY3, PARITY4};
   22.89 -
   22.90 -/*----------------
   22.91 -  STATIC FUNCTIONS
   22.92 -  ----------------*/
   22.93 -inline static int idxof(int i);
   22.94 -inline static void rshift128(w128_t *out,  w128_t const *in, int shift);
   22.95 -inline static void lshift128(w128_t *out,  w128_t const *in, int shift);
   22.96 -inline static void gen_rand_all(void);
   22.97 -inline static void gen_rand_array(w128_t *array, int size);
   22.98 -inline static uint32_t func1(uint32_t x);
   22.99 -inline static uint32_t func2(uint32_t x);
  22.100 -static void period_certification(void);
  22.101 -#if defined(BIG_ENDIAN64) && !defined(ONLY64)
  22.102 -inline static void swap(w128_t *array, int size);
  22.103 -#endif
  22.104 -
  22.105 -#if defined(HAVE_ALTIVEC)
  22.106 -  #include "SFMT-alti.h"
  22.107 -#elif defined(HAVE_SSE2)
  22.108 -  #include "SFMT-sse2.h"
  22.109 -#endif
  22.110 -
  22.111 -/**
  22.112 - * This function simulate a 64-bit index of LITTLE ENDIAN 
  22.113 - * in BIG ENDIAN machine.
  22.114 - */
  22.115 -#ifdef ONLY64
  22.116 -inline static int idxof(int i) {
  22.117 -    return i ^ 1;
  22.118 -}
  22.119 -#else
  22.120 -inline static int idxof(int i) {
  22.121 -    return i;
  22.122 -}
  22.123 -#endif
  22.124 -/**
  22.125 - * This function simulates SIMD 128-bit right shift by the standard C.
  22.126 - * The 128-bit integer given in in is shifted by (shift * 8) bits.
  22.127 - * This function simulates the LITTLE ENDIAN SIMD.
  22.128 - * @param out the output of this function
  22.129 - * @param in the 128-bit data to be shifted
  22.130 - * @param shift the shift value
  22.131 - */
  22.132 -#ifdef ONLY64
  22.133 -inline static void rshift128(w128_t *out, w128_t const *in, int shift) {
  22.134 -    uint64_t th, tl, oh, ol;
  22.135 -
  22.136 -    th = ((uint64_t)in->u[2] << 32) | ((uint64_t)in->u[3]);
  22.137 -    tl = ((uint64_t)in->u[0] << 32) | ((uint64_t)in->u[1]);
  22.138 -
  22.139 -    oh = th >> (shift * 8);
  22.140 -    ol = tl >> (shift * 8);
  22.141 -    ol |= th << (64 - shift * 8);
  22.142 -    out->u[0] = (uint32_t)(ol >> 32);
  22.143 -    out->u[1] = (uint32_t)(ol & 0xffffffff);
  22.144 -    out->u[2] = (uint32_t)(oh >> 32);
  22.145 -    out->u[3] = (uint32_t)(oh & 0xffffffff);
  22.146 -}
  22.147 -#else
  22.148 -inline static void rshift128(w128_t *out, w128_t const *in, int shift) {
  22.149 -    uint64_t th, tl, oh, ol;
  22.150 -
  22.151 -    th = ((uint64_t)in->u[3] << 32) | ((uint64_t)in->u[2]);
  22.152 -    tl = ((uint64_t)in->u[1] << 32) | ((uint64_t)in->u[0]);
  22.153 -
  22.154 -    oh = th >> (shift * 8);
  22.155 -    ol = tl >> (shift * 8);
  22.156 -    ol |= th << (64 - shift * 8);
  22.157 -    out->u[1] = (uint32_t)(ol >> 32);
  22.158 -    out->u[0] = (uint32_t)(ol & 0xffffffff);
  22.159 -    out->u[3] = (uint32_t)(oh >> 32);
  22.160 -    out->u[2] = (uint32_t)(oh & 0xffffffff);
  22.161 -}
  22.162 -#endif
  22.163 -/**
  22.164 - * This function simulates SIMD 128-bit left shift by the standard C.
  22.165 - * The 128-bit integer given in in is shifted by (shift * 8) bits.
  22.166 - * This function simulates the LITTLE ENDIAN SIMD.
  22.167 - * @param out the output of this function
  22.168 - * @param in the 128-bit data to be shifted
  22.169 - * @param shift the shift value
  22.170 - */
  22.171 -#ifdef ONLY64
  22.172 -inline static void lshift128(w128_t *out, w128_t const *in, int shift) {
  22.173 -    uint64_t th, tl, oh, ol;
  22.174 -
  22.175 -    th = ((uint64_t)in->u[2] << 32) | ((uint64_t)in->u[3]);
  22.176 -    tl = ((uint64_t)in->u[0] << 32) | ((uint64_t)in->u[1]);
  22.177 -
  22.178 -    oh = th << (shift * 8);
  22.179 -    ol = tl << (shift * 8);
  22.180 -    oh |= tl >> (64 - shift * 8);
  22.181 -    out->u[0] = (uint32_t)(ol >> 32);
  22.182 -    out->u[1] = (uint32_t)(ol & 0xffffffff);
  22.183 -    out->u[2] = (uint32_t)(oh >> 32);
  22.184 -    out->u[3] = (uint32_t)(oh & 0xffffffff);
  22.185 -}
  22.186 -#else
  22.187 -inline static void lshift128(w128_t *out, w128_t const *in, int shift) {
  22.188 -    uint64_t th, tl, oh, ol;
  22.189 -
  22.190 -    th = ((uint64_t)in->u[3] << 32) | ((uint64_t)in->u[2]);
  22.191 -    tl = ((uint64_t)in->u[1] << 32) | ((uint64_t)in->u[0]);
  22.192 -
  22.193 -    oh = th << (shift * 8);
  22.194 -    ol = tl << (shift * 8);
  22.195 -    oh |= tl >> (64 - shift * 8);
  22.196 -    out->u[1] = (uint32_t)(ol >> 32);
  22.197 -    out->u[0] = (uint32_t)(ol & 0xffffffff);
  22.198 -    out->u[3] = (uint32_t)(oh >> 32);
  22.199 -    out->u[2] = (uint32_t)(oh & 0xffffffff);
  22.200 -}
  22.201 -#endif
  22.202 -
  22.203 -/**
  22.204 - * This function represents the recursion formula.
  22.205 - * @param r output
  22.206 - * @param a a 128-bit part of the internal state array
  22.207 - * @param b a 128-bit part of the internal state array
  22.208 - * @param c a 128-bit part of the internal state array
  22.209 - * @param d a 128-bit part of the internal state array
  22.210 - */
  22.211 -#if (!defined(HAVE_ALTIVEC)) && (!defined(HAVE_SSE2))
  22.212 -#ifdef ONLY64
  22.213 -inline static void do_recursion(w128_t *r, w128_t *a, w128_t *b, w128_t *c,
  22.214 -				w128_t *d) {
  22.215 -    w128_t x;
  22.216 -    w128_t y;
  22.217 -
  22.218 -    lshift128(&x, a, SL2);
  22.219 -    rshift128(&y, c, SR2);
  22.220 -    r->u[0] = a->u[0] ^ x.u[0] ^ ((b->u[0] >> SR1) & MSK2) ^ y.u[0] 
  22.221 -	^ (d->u[0] << SL1);
  22.222 -    r->u[1] = a->u[1] ^ x.u[1] ^ ((b->u[1] >> SR1) & MSK1) ^ y.u[1] 
  22.223 -	^ (d->u[1] << SL1);
  22.224 -    r->u[2] = a->u[2] ^ x.u[2] ^ ((b->u[2] >> SR1) & MSK4) ^ y.u[2] 
  22.225 -	^ (d->u[2] << SL1);
  22.226 -    r->u[3] = a->u[3] ^ x.u[3] ^ ((b->u[3] >> SR1) & MSK3) ^ y.u[3] 
  22.227 -	^ (d->u[3] << SL1);
  22.228 -}
  22.229 -#else
  22.230 -inline static void do_recursion(w128_t *r, w128_t *a, w128_t *b, w128_t *c,
  22.231 -				w128_t *d) {
  22.232 -    w128_t x;
  22.233 -    w128_t y;
  22.234 -
  22.235 -    lshift128(&x, a, SL2);
  22.236 -    rshift128(&y, c, SR2);
  22.237 -    r->u[0] = a->u[0] ^ x.u[0] ^ ((b->u[0] >> SR1) & MSK1) ^ y.u[0] 
  22.238 -	^ (d->u[0] << SL1);
  22.239 -    r->u[1] = a->u[1] ^ x.u[1] ^ ((b->u[1] >> SR1) & MSK2) ^ y.u[1] 
  22.240 -	^ (d->u[1] << SL1);
  22.241 -    r->u[2] = a->u[2] ^ x.u[2] ^ ((b->u[2] >> SR1) & MSK3) ^ y.u[2] 
  22.242 -	^ (d->u[2] << SL1);
  22.243 -    r->u[3] = a->u[3] ^ x.u[3] ^ ((b->u[3] >> SR1) & MSK4) ^ y.u[3] 
  22.244 -	^ (d->u[3] << SL1);
  22.245 -}
  22.246 -#endif
  22.247 -#endif
  22.248 -
  22.249 -#if (!defined(HAVE_ALTIVEC)) && (!defined(HAVE_SSE2))
  22.250 -/**
  22.251 - * This function fills the internal state array with pseudorandom
  22.252 - * integers.
  22.253 - */
  22.254 -inline static void gen_rand_all(void) {
  22.255 -    int i;
  22.256 -    w128_t *r1, *r2;
  22.257 -
  22.258 -    r1 = &sfmt[N - 2];
  22.259 -    r2 = &sfmt[N - 1];
  22.260 -    for (i = 0; i < N - POS1; i++) {
  22.261 -	do_recursion(&sfmt[i], &sfmt[i], &sfmt[i + POS1], r1, r2);
  22.262 -	r1 = r2;
  22.263 -	r2 = &sfmt[i];
  22.264 -    }
  22.265 -    for (; i < N; i++) {
  22.266 -	do_recursion(&sfmt[i], &sfmt[i], &sfmt[i + POS1 - N], r1, r2);
  22.267 -	r1 = r2;
  22.268 -	r2 = &sfmt[i];
  22.269 -    }
  22.270 -}
  22.271 -
  22.272 -/**
  22.273 - * This function fills the user-specified array with pseudorandom
  22.274 - * integers.
  22.275 - *
  22.276 - * @param array an 128-bit array to be filled by pseudorandom numbers.  
  22.277 - * @param size number of 128-bit pseudorandom numbers to be generated.
  22.278 - */
  22.279 -inline static void gen_rand_array(w128_t *array, int size) {
  22.280 -    int i, j;
  22.281 -    w128_t *r1, *r2;
  22.282 -
  22.283 -    r1 = &sfmt[N - 2];
  22.284 -    r2 = &sfmt[N - 1];
  22.285 -    for (i = 0; i < N - POS1; i++) {
  22.286 -	do_recursion(&array[i], &sfmt[i], &sfmt[i + POS1], r1, r2);
  22.287 -	r1 = r2;
  22.288 -	r2 = &array[i];
  22.289 -    }
  22.290 -    for (; i < N; i++) {
  22.291 -	do_recursion(&array[i], &sfmt[i], &array[i + POS1 - N], r1, r2);
  22.292 -	r1 = r2;
  22.293 -	r2 = &array[i];
  22.294 -    }
  22.295 -    for (; i < size - N; i++) {
  22.296 -	do_recursion(&array[i], &array[i - N], &array[i + POS1 - N], r1, r2);
  22.297 -	r1 = r2;
  22.298 -	r2 = &array[i];
  22.299 -    }
  22.300 -    for (j = 0; j < 2 * N - size; j++) {
  22.301 -	sfmt[j] = array[j + size - N];
  22.302 -    }
  22.303 -    for (; i < size; i++, j++) {
  22.304 -	do_recursion(&array[i], &array[i - N], &array[i + POS1 - N], r1, r2);
  22.305 -	r1 = r2;
  22.306 -	r2 = &array[i];
  22.307 -	sfmt[j] = array[i];
  22.308 -    }
  22.309 -}
  22.310 -#endif
  22.311 -
  22.312 -#if defined(BIG_ENDIAN64) && !defined(ONLY64) && !defined(HAVE_ALTIVEC)
  22.313 -inline static void swap(w128_t *array, int size) {
  22.314 -    int i;
  22.315 -    uint32_t x, y;
  22.316 -
  22.317 -    for (i = 0; i < size; i++) {
  22.318 -	x = array[i].u[0];
  22.319 -	y = array[i].u[2];
  22.320 -	array[i].u[0] = array[i].u[1];
  22.321 -	array[i].u[2] = array[i].u[3];
  22.322 -	array[i].u[1] = x;
  22.323 -	array[i].u[3] = y;
  22.324 -    }
  22.325 -}
  22.326 -#endif
  22.327 -/**
  22.328 - * This function represents a function used in the initialization
  22.329 - * by init_by_array
  22.330 - * @param x 32-bit integer
  22.331 - * @return 32-bit integer
  22.332 - */
  22.333 -static uint32_t func1(uint32_t x) {
  22.334 -    return (x ^ (x >> 27)) * (uint32_t)1664525UL;
  22.335 -}
  22.336 -
  22.337 -/**
  22.338 - * This function represents a function used in the initialization
  22.339 - * by init_by_array
  22.340 - * @param x 32-bit integer
  22.341 - * @return 32-bit integer
  22.342 - */
  22.343 -static uint32_t func2(uint32_t x) {
  22.344 -    return (x ^ (x >> 27)) * (uint32_t)1566083941UL;
  22.345 -}
  22.346 -
  22.347 -/**
  22.348 - * This function certificate the period of 2^{MEXP}
  22.349 - */
  22.350 -static void period_certification(void) {
  22.351 -    int inner = 0;
  22.352 -    int i, j;
  22.353 -    uint32_t work;
  22.354 -
  22.355 -    for (i = 0; i < 4; i++)
  22.356 -	inner ^= psfmt32[idxof(i)] & parity[i];
  22.357 -    for (i = 16; i > 0; i >>= 1)
  22.358 -	inner ^= inner >> i;
  22.359 -    inner &= 1;
  22.360 -    /* check OK */
  22.361 -    if (inner == 1) {
  22.362 -	return;
  22.363 -    }
  22.364 -    /* check NG, and modification */
  22.365 -    for (i = 0; i < 4; i++) {
  22.366 -	work = 1;
  22.367 -	for (j = 0; j < 32; j++) {
  22.368 -	    if ((work & parity[i]) != 0) {
  22.369 -		psfmt32[idxof(i)] ^= work;
  22.370 -		return;
  22.371 -	    }
  22.372 -	    work = work << 1;
  22.373 -	}
  22.374 -    }
  22.375 -}
  22.376 -
  22.377 -/*----------------
  22.378 -  PUBLIC FUNCTIONS
  22.379 -  ----------------*/
  22.380 -/**
  22.381 - * This function returns the identification string.
  22.382 - * The string shows the word size, the Mersenne exponent,
  22.383 - * and all parameters of this generator.
  22.384 - */
  22.385 -const char *get_idstring(void) {
  22.386 -    return IDSTR;
  22.387 -}
  22.388 -
  22.389 -/**
  22.390 - * This function returns the minimum size of array used for \b
  22.391 - * fill_array32() function.
  22.392 - * @return minimum size of array used for fill_array32() function.
  22.393 - */
  22.394 -int get_min_array_size32(void) {
  22.395 -    return N32;
  22.396 -}
  22.397 -
  22.398 -/**
  22.399 - * This function returns the minimum size of array used for \b
  22.400 - * fill_array64() function.
  22.401 - * @return minimum size of array used for fill_array64() function.
  22.402 - */
  22.403 -int get_min_array_size64(void) {
  22.404 -    return N64;
  22.405 -}
  22.406 -
  22.407 -#ifndef ONLY64
  22.408 -/**
  22.409 - * This function generates and returns 32-bit pseudorandom number.
  22.410 - * init_gen_rand or init_by_array must be called before this function.
  22.411 - * @return 32-bit pseudorandom number
  22.412 - */
  22.413 -uint32_t gen_rand32(void) {
  22.414 -    uint32_t r;
  22.415 -
  22.416 -    assert(initialized);
  22.417 -    if (idx >= N32) {
  22.418 -	gen_rand_all();
  22.419 -	idx = 0;
  22.420 -    }
  22.421 -    r = psfmt32[idx++];
  22.422 -    return r;
  22.423 -}
  22.424 -#endif
  22.425 -/**
  22.426 - * This function generates and returns 64-bit pseudorandom number.
  22.427 - * init_gen_rand or init_by_array must be called before this function.
  22.428 - * The function gen_rand64 should not be called after gen_rand32,
  22.429 - * unless an initialization is again executed. 
  22.430 - * @return 64-bit pseudorandom number
  22.431 - */
  22.432 -uint64_t gen_rand64(void) {
  22.433 -#if defined(BIG_ENDIAN64) && !defined(ONLY64)
  22.434 -    uint32_t r1, r2;
  22.435 -#else
  22.436 -    uint64_t r;
  22.437 -#endif
  22.438 -
  22.439 -    assert(initialized);
  22.440 -    assert(idx % 2 == 0);
  22.441 -
  22.442 -    if (idx >= N32) {
  22.443 -	gen_rand_all();
  22.444 -	idx = 0;
  22.445 -    }
  22.446 -#if defined(BIG_ENDIAN64) && !defined(ONLY64)
  22.447 -    r1 = psfmt32[idx];
  22.448 -    r2 = psfmt32[idx + 1];
  22.449 -    idx += 2;
  22.450 -    return ((uint64_t)r2 << 32) | r1;
  22.451 -#else
  22.452 -    r = psfmt64[idx / 2];
  22.453 -    idx += 2;
  22.454 -    return r;
  22.455 -#endif
  22.456 -}
  22.457 -
  22.458 -#ifndef ONLY64
  22.459 -/**
  22.460 - * This function generates pseudorandom 32-bit integers in the
  22.461 - * specified array[] by one call. The number of pseudorandom integers
  22.462 - * is specified by the argument size, which must be at least 624 and a
  22.463 - * multiple of four.  The generation by this function is much faster
  22.464 - * than the following gen_rand function.
  22.465 - *
  22.466 - * For initialization, init_gen_rand or init_by_array must be called
  22.467 - * before the first call of this function. This function can not be
  22.468 - * used after calling gen_rand function, without initialization.
  22.469 - *
  22.470 - * @param array an array where pseudorandom 32-bit integers are filled
  22.471 - * by this function.  The pointer to the array must be \b "aligned"
  22.472 - * (namely, must be a multiple of 16) in the SIMD version, since it
  22.473 - * refers to the address of a 128-bit integer.  In the standard C
  22.474 - * version, the pointer is arbitrary.
  22.475 - *
  22.476 - * @param size the number of 32-bit pseudorandom integers to be
  22.477 - * generated.  size must be a multiple of 4, and greater than or equal
  22.478 - * to (MEXP / 128 + 1) * 4.
  22.479 - *
  22.480 - * @note \b memalign or \b posix_memalign is available to get aligned
  22.481 - * memory. Mac OSX doesn't have these functions, but \b malloc of OSX
  22.482 - * returns the pointer to the aligned memory block.
  22.483 - */
  22.484 -void fill_array32(uint32_t *array, int size) {
  22.485 -    assert(initialized);
  22.486 -    assert(idx == N32);
  22.487 -    assert(size % 4 == 0);
  22.488 -    assert(size >= N32);
  22.489 -
  22.490 -    gen_rand_array((w128_t *)array, size / 4);
  22.491 -    idx = N32;
  22.492 -}
  22.493 -#endif
  22.494 -
  22.495 -/**
  22.496 - * This function generates pseudorandom 64-bit integers in the
  22.497 - * specified array[] by one call. The number of pseudorandom integers
  22.498 - * is specified by the argument size, which must be at least 312 and a
  22.499 - * multiple of two.  The generation by this function is much faster
  22.500 - * than the following gen_rand function.
  22.501 - *
  22.502 - * For initialization, init_gen_rand or init_by_array must be called
  22.503 - * before the first call of this function. This function can not be
  22.504 - * used after calling gen_rand function, without initialization.
  22.505 - *
  22.506 - * @param array an array where pseudorandom 64-bit integers are filled
  22.507 - * by this function.  The pointer to the array must be "aligned"
  22.508 - * (namely, must be a multiple of 16) in the SIMD version, since it
  22.509 - * refers to the address of a 128-bit integer.  In the standard C
  22.510 - * version, the pointer is arbitrary.
  22.511 - *
  22.512 - * @param size the number of 64-bit pseudorandom integers to be
  22.513 - * generated.  size must be a multiple of 2, and greater than or equal
  22.514 - * to (MEXP / 128 + 1) * 2
  22.515 - *
  22.516 - * @note \b memalign or \b posix_memalign is available to get aligned
  22.517 - * memory. Mac OSX doesn't have these functions, but \b malloc of OSX
  22.518 - * returns the pointer to the aligned memory block.
  22.519 - */
  22.520 -void fill_array64(uint64_t *array, int size) {
  22.521 -    assert(initialized);
  22.522 -    assert(idx == N32);
  22.523 -    assert(size % 2 == 0);
  22.524 -    assert(size >= N64);
  22.525 -
  22.526 -    gen_rand_array((w128_t *)array, size / 2);
  22.527 -    idx = N32;
  22.528 -
  22.529 -#if defined(BIG_ENDIAN64) && !defined(ONLY64)
  22.530 -    swap((w128_t *)array, size /2);
  22.531 -#endif
  22.532 -}
  22.533 -
  22.534 -/**
  22.535 - * This function initializes the internal state array with a 32-bit
  22.536 - * integer seed.
  22.537 - *
  22.538 - * @param seed a 32-bit integer used as the seed.
  22.539 - */
  22.540 -void init_gen_rand(uint32_t seed) {
  22.541 -    int i;
  22.542 -
  22.543 -    psfmt32[idxof(0)] = seed;
  22.544 -    for (i = 1; i < N32; i++) {
  22.545 -	psfmt32[idxof(i)] = 1812433253UL * (psfmt32[idxof(i - 1)] 
  22.546 -					    ^ (psfmt32[idxof(i - 1)] >> 30))
  22.547 -	    + i;
  22.548 -    }
  22.549 -    idx = N32;
  22.550 -    period_certification();
  22.551 -    initialized = 1;
  22.552 -}
  22.553 -
  22.554 -/**
  22.555 - * This function initializes the internal state array,
  22.556 - * with an array of 32-bit integers used as the seeds
  22.557 - * @param init_key the array of 32-bit integers, used as a seed.
  22.558 - * @param key_length the length of init_key.
  22.559 - */
  22.560 -void init_by_array(uint32_t *init_key, int key_length) {
  22.561 -    int i, j, count;
  22.562 -    uint32_t r;
  22.563 -    int lag;
  22.564 -    int mid;
  22.565 -    int size = N * 4;
  22.566 -
  22.567 -    if (size >= 623) {
  22.568 -	lag = 11;
  22.569 -    } else if (size >= 68) {
  22.570 -	lag = 7;
  22.571 -    } else if (size >= 39) {
  22.572 -	lag = 5;
  22.573 -    } else {
  22.574 -	lag = 3;
  22.575 -    }
  22.576 -    mid = (size - lag) / 2;
  22.577 -
  22.578 -    memset(sfmt, 0x8b, sizeof(sfmt));
  22.579 -    if (key_length + 1 > N32) {
  22.580 -	count = key_length + 1;
  22.581 -    } else {
  22.582 -	count = N32;
  22.583 -    }
  22.584 -    r = func1(psfmt32[idxof(0)] ^ psfmt32[idxof(mid)] 
  22.585 -	      ^ psfmt32[idxof(N32 - 1)]);
  22.586 -    psfmt32[idxof(mid)] += r;
  22.587 -    r += key_length;
  22.588 -    psfmt32[idxof(mid + lag)] += r;
  22.589 -    psfmt32[idxof(0)] = r;
  22.590 -
  22.591 -    count--;
  22.592 -    for (i = 1, j = 0; (j < count) && (j < key_length); j++) {
  22.593 -	r = func1(psfmt32[idxof(i)] ^ psfmt32[idxof((i + mid) % N32)] 
  22.594 -		  ^ psfmt32[idxof((i + N32 - 1) % N32)]);
  22.595 -	psfmt32[idxof((i + mid) % N32)] += r;
  22.596 -	r += init_key[j] + i;
  22.597 -	psfmt32[idxof((i + mid + lag) % N32)] += r;
  22.598 -	psfmt32[idxof(i)] = r;
  22.599 -	i = (i + 1) % N32;
  22.600 -    }
  22.601 -    for (; j < count; j++) {
  22.602 -	r = func1(psfmt32[idxof(i)] ^ psfmt32[idxof((i + mid) % N32)] 
  22.603 -		  ^ psfmt32[idxof((i + N32 - 1) % N32)]);
  22.604 -	psfmt32[idxof((i + mid) % N32)] += r;
  22.605 -	r += i;
  22.606 -	psfmt32[idxof((i + mid + lag) % N32)] += r;
  22.607 -	psfmt32[idxof(i)] = r;
  22.608 -	i = (i + 1) % N32;
  22.609 -    }
  22.610 -    for (j = 0; j < N32; j++) {
  22.611 -	r = func2(psfmt32[idxof(i)] + psfmt32[idxof((i + mid) % N32)] 
  22.612 -		  + psfmt32[idxof((i + N32 - 1) % N32)]);
  22.613 -	psfmt32[idxof((i + mid) % N32)] ^= r;
  22.614 -	r -= i;
  22.615 -	psfmt32[idxof((i + mid + lag) % N32)] ^= r;
  22.616 -	psfmt32[idxof(i)] = r;
  22.617 -	i = (i + 1) % N32;
  22.618 -    }
  22.619 -
  22.620 -    idx = N32;
  22.621 -    period_certification();
  22.622 -    initialized = 1;
  22.623 -}