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

DSCF-1224

% [reference]
% https://tex.stackexchange.com/questions/34424/how-do-i-calculate-n-modulo-3-in-latex

\RequirePackage{plautopatch}
\documentclass[uplatex]{jsarticle}

\makeatletter

\newcommand{\calcQuotient}[2]{{%
	\newcount\@quotient%

DSCF-1224

FCFLAGS = -ffree-line-length-none -O2 -pedantic -std=f2008 -Wall -Wextra # -fbacktrace -fbounds-check 
OBJS    = ./main.o

main.exe: $(OBJS)
	gfortran $(FCFLAGS) -o ./main.exe $(OBJS)

%.o: %.f08
	gfortran $(FCFLAGS) -c $<

clean:

DSCF-1224

FCFLAGS = -ffree-line-length-none -O2 -pedantic -std=f2008 -Wall -Wextra # -fbacktrace -fbounds-check 
OBJS    = ./mod_utility.o ./mod_eqn_func.o ./mod_solver.o ./main.o

main.exe: $(OBJS)
	gfortran $(FCFLAGS) -o ./main.exe $(OBJS)

%.o: %.f08
	gfortran $(FCFLAGS) -c $<

clean:

DSCF-1224

参考文献

• 日本ニューメリカルアルゴリズムズグループ株式会社（日本 NAG）
  • Fortran Tip集: 環境変数の値を取得する方法
  • Fortran 2003入門: 環境変数へのアクセス

DSCF-1224

FCFLAGS = -ffree-line-length-none -O2 -pedantic -std=f2008 -Wall -Wextra # -fbacktrace -fbounds-check 
OBJS    = ./mod_is_even.o ./mod_test.o ./main.o

main.exe: $(OBJS)
	gfortran $(FCFLAGS) -o ./main.exe $(OBJS)

%.o: %.f08
	gfortran $(FCFLAGS) -c $<

clean:

DSCF-1224

FCFLAGS= -fbacktrace -fbounds-check -ffree-line-length-none -O2 -pedantic -std=f2008 -Wall -Wextra

main.exe: ./main.f08
	gfortran $(FCFLAGS) -o ./main.exe $<

clean:
	rm ./*.mod ./*.o ./*.dat

DSCF-1224

# -*- coding: utf-8 -*-
import numpy

def read_int32(file_stream):
    return numpy.fromfile(file= file_stream, dtype= numpy.int32, count=1)[0]

def read_real64(file_stream):
    return numpy.fromfile(file= file_stream, dtype= numpy.float64, count=1)[0]

DSCF-1224

参考文献

• 数値計算のつぼ | 共立出版
• カハンの加算アルゴリズム - Wikipedia
• Kahan summation algorithm - Wikipedia

出力

 268435456 < length_data

DSCF-1224

program memorandum

    use, intrinsic :: iso_fortran_env

    implicit none

    integer(INT32) :: input_front = -17_INT32
    integer(INT32) :: input_rear  = +32_INT32
    integer(INT64) :: retval