terasakisatoshi

import gistcheck
gistcheck.commit()

terasakisatoshi

terasakisatoshi

import random 
import linalg 
from math import ln,sqrt,exp
import times 

const N = 100
const beta:float64 = ln(1+sqrt(2.0))/2

randomize() 

terasakisatoshi

package main

import (
	"fmt"
	"github.com/gonum/matrix/mat64"
	"math"
	"math/rand"
	"time"
)

terasakisatoshi

[package]
name = "ising"
version = "0.1.0"
authors = ["SatoshiTeraaki <[email protected]>"]

[dependencies]
ndarray = "0.11.0"
rand = "0.3"
ndarray-rand="0.7.0"

terasakisatoshi

# Reference: http://nbviewer.jupyter.org/gist/genkuroki/4fa46c68c56ee0f3b1a6fc8ec628b9d7
# Author MathSorcerer

from math import log, sqrt, exp
from random import choice, random
from copy import deepcopy
import time
from numba import jit
import numpy as np

terasakisatoshi

# Reference: http://nbviewer.jupyter.org/gist/genkuroki/4fa46c68c56ee0f3b1a6fc8ec628b9d7
# Author MathSorcerer

from math import exp
import array
from cpython cimport array as carray
from random import choice, random
from libc.stdlib cimport rand, RAND_MAX
import time
cimport cython

terasakisatoshi

Ising Model niter=1e9の場合の計測時間表．

PC環境
---------------------------
julia>versioninfo()
Platform Info:
  OS: macOS (x86_64-apple-darwin14.5.0)
  CPU: Intel(R) Core(TM) M-5Y51 CPU @ 1.10GHz
  WORD_SIZE: 64
  BLAS: libopenblas (USE64BITINT DYNAMIC_ARCH NO_AFFINITY Haswell)

terasakisatoshi

terasakisatoshi