DSCF-1224

#include <random>
#include <immintrin.h>

#ifndef PRNG_HPP_
#define PRNG_HPP_

namespace prng
{
	class Prng32
	{

DSCF-1224

# The Ziggurat Method for generating random variables - Marsaglia and Tsang
# Paper and reference code: http://www.jstatsoft.org/v05/i08/

module ZigguratTable

    function compute_core( table_size::Int, pdf0::T, r::T, pdf::Core.Function, integral_pdf::Core.Function, inverse_pdf::Core.Function ) where { T <: Core.AbstractFloat }

        v     =    ( r * pdf( r ) ) + integral_pdf( r )
        x_ref =    r
        x_new =    Base.zero( T )

DSCF-1224

! original: https://github.com/vnglst/pong-wars
! gfortran -O3 -Wall -std=f2008 main.f90 && time ./a.out && gnuplot ./score.plt

program pong_wars

    use , intrinsic :: iso_fortran_env

    implicit none

DSCF-1224

参考文献

• 変分モンテカルロ法で遊ぶ
• 変分モンテカルロ(VMC)をやってみる水素原子 - MiyanTarumi’s blog

実行例

Rev1

DSCF-1224

fig_1_1_hello.cpp

#include <iostream>
#include <sycl/sycl.hpp>

const std::string secret_string{ "Ifmmp-!xpsme\"\012J(n!tpssz-!Ebwf/!"    "J(n!bgsbje!J!dbo(u!ep!uibu/!.!IBM\01" };
const auto secret_string_size = secret_string.size();

int main(void)

DSCF-1224

import numpy
import os
import pickle
import requests



def convert_data(targetData, fileName):

    with open(fileName, 'wb') as writingFileStream:

DSCF-1224

# compiler selection
FC = gfortran

# compiler option: Common Options
FFLAGS_COMMON = -ffree-line-length-none -fimplicit-none -pedantic -std=f2008 -Wall -Werror -Wextra

# compiler option: Release mode
FFLAGS_RELEASE = ${FFLAGS_COMMON} -O3

# compiler option: Debug mode

DSCF-1224

# [reference]
# https://qiita.com/tell/items/5209b92d5f525ca7b028
# https://ja.wikipedia.org/wiki/ヒストグラム

reset session

set datafile separator comma

# DATA_FILE_PATH = '10e3.csv'
# DATA_FILE_PATH = '10e4.csv'

DSCF-1224

! gfortran, gcc version 7.4.0
! https://stdlib.fortran-lang.org/sourcefile/stdlib_strings.fypp.html

program test

    use , intrinsic :: iso_fortran_env
    implicit none

    character(len=1) , parameter :: ASCII_HT = achar(int(Z'09'))
    character(len=1) , parameter :: ASCII_LF = achar(int(Z'0A'))

DSCF-1224

README

mt19937-64.out.txt に掲載されている、1000 outputs of genrand64_int64() を、Fortran で読み取れるように、 バイナリデータは保持しつつ、unsigned long から signed long に変換するための C++ と Fortran のソースコード。

参考文献

• Mersenne Twister: A random number generator (since 1997/10)
  • mt19937-64.out.txt