C# Mersenne Twister 随机整数生成器实现（SFMT）蒙特卡罗模拟

到目前为止我一直在使用C#梅森扭转者 http://en.wikipedia.org/wiki/Mersenne_twister在这里找到生成随机数：

http://www.centerspace.net/resources.php http://www.centerspace.net/resources.php

我刚刚发现SFMT http://en.wikipedia.org/wiki/Mersenne_twister#SFMT这里的速度应该是两倍：

http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/SFMT/ http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/SFMT/

谁能告诉我 SFMT 的 C# 实现?

我的要求是生成一个介于（含）0 和 2^20 (1048576) 之间的整数。

我需要这样做每天数万亿次对于在 24 小时时钟上运行的模拟，因此我准备花几天时间将其调整到完美。

目前，我已经通过添加新方法来调整 Center Space Mersenne Twister 来满足我的要求：

public uint Next20()
{            
    return (uint)(genrand_int32() >> 12);
}

使用方法genrand_int32()我想制作自己的版本，genrand_int20()，生成 0 和 2^20 之间（含）之间的整数以保存上方投射并移动但我不懂数学。我到底该怎么做呢？

还正在使用一个uint 会比 int 更快，或者只是可寻址数字的问题？因为我只需要最多 1048576，所以我只关心速度。

这也将运行在带 .NET 2 的 Windows Server 2003 R2 SP2（32 位）盒。处理器是AMD 皓龙 275 http://en.wikipedia.org/wiki/List_of_AMD_Opteron_microprocessors#Opteron_200-series_.22Italy.22_.28E6.2C_90_nm.29(4核).

你能做的就是下载来源来自链接 http://www.codeproject.com/KB/DLL/SFMT_dll.aspx?msg=3130186您在 Code Project 上发现的。解压缩它，在 Visual Studio 中加载解决方案并编译它。这将为您提供源代码、非托管 c dll 和 .lib 文件。

您可以 P/Invoke 此 dll 中的函数（仅导出 5 个简单函数，您只需要其中两个）或者您可以使用此 dll、lib 和 SFMT 头文件来创建您可以使用的托管包装器 dll在没有 P/Invoke 的 C# 中。我刚刚尝试了这个方法，而且做起来非常简单。没有涉及明确的编组。

就是这样。下载并编译后source（除了 dll 之外，您还需要创建的头文件和 lib 文件）创建一个新的 C++ CLR 类库项目。称之为 WrapSFMT 或其他名称。转到项目属性。在 C++/预编译头下，更改为“不使用预编译头”。在链接器/常规/附加库目录下，输入 SFMT.lib 的路径。在链接器/输入/附加依赖项下，添加 SFMT.lib。关闭属性页。将 SFMT.h 复制到您的项目文件夹并将其包含在项目中。

编辑 WrapSFMT.h 如下：

#pragma once
#include "SFMT.H"

using namespace System;

namespace WrapSFMT {

public ref class SRandom
{
public:SRandom(UInt32);
public:UInt32 Rand32(void);
};
}

它们声明了您的类中的方法。现在编辑 WrapSFMT.cpp 如下：

#include "WrapSFMT.h"

namespace WrapSFMT {

SRandom::SRandom(UInt32 seed)
{
    init_gen_rand(seed);
}

UInt32 SRandom::Rand32()
{
    return gen_rand32();
}
}

它们实现了您在头文件中声明的方法。您所做的只是从 SFMT.dll 调用函数，C++/CLI 会自动处理从非托管到托管的转换。现在您应该能够构建 WrapSFMT.dll 并在 C# 项目中引用它。确保 SFMT.dll 在路径中，应该不会有任何问题。