Add version 2.

Author: phcoding

fortran/build.cmd

@@ -0,0 +1 @@
+gfortran util.f90 ising.f90 ising_main.f90 -o ising

fortran/README.md renamed to fortran/v1/README.md

@@ -0,0 +1 @@
+gfortran util.f90 ising.f90 handle.f90 main.f90 -o ising

fortran/v1/build.cmd

@@ -0,0 +1,92 @@
+module handle
+    use ising
+    implicit none
+    real, parameter :: esp = 0.000001, J1 = 1
+contains
+    real function get_delta_energy(spin, row, col) result(delta_energy)
+        implicit none
+        integer, intent(in)     :: row, col
+        type(spins), intent(in) :: spin
+        integer(kind=1)         :: sr, sl, st, sb
+
+        if (row < 2) then
+            st = spin%data(spin%rows, col)
+            sb = spin%data(row + 1  , col)
+        else if(row > spin%rows-1) then
+            st = spin%data(row - 1  , col)
+            sb = spin%data(1        , col)
+        else
+            st = spin%data(row - 1  , col)
+            sb = spin%data(row + 1  , col)
+        end if
+
+        if (col < 2) then
+            sl = spin%data(row, spin%cols)
+            sr = spin%data(row  , col + 1)
+        else if(col > spin%cols-1) then
+            sl = spin%data(row  , col - 1)
+            sr = spin%data(row        , 1)
+        else
+            sl = spin%data(row  , col - 1)
+            sr = spin%data(row  , col + 1)
+        end if
+        delta_energy = 2*spin%data(row, col)*(J1*(sl+sr+st+sb)+spin%magn)
+    end function get_delta_energy
+
+    real function get_total_energy(spin) result(total_energy)
+        implicit none
+        type(spins)     :: spin
+        integer         :: row, col
+        integer(kind=1) :: sr, sb
+        
+        total_energy = 0
+        do col = 1, spin%cols
+            do row = 1, spin%rows
+                if (row > spin%rows-1) then
+                    sb = spin%data(1, col)
+                else
+                    sb = spin%data(row+1, col)
+                end if
+                if (col > spin%cols-1) then
+                    sr = spin%data(row, 1)
+                else
+                    sr = spin%data(row, col+1)
+                end if
+                total_energy = total_energy - spin%data(row, col)*(J1*(sr+sb)+spin%magn)
+            end do
+        end do
+    end function get_total_energy
+
+    real function ising_average(spin, value) result(value_average)
+        ! 
+        ! Calculate average value of ising quantity by given value.
+        ! 
+        !   spin    [type(spins)]   spin data storage of ising system.
+        !   value   [real]          value of ising quantity.
+        !   return  [real]          average value of ising quantity.
+        ! 
+        implicit none
+        real :: value
+        type(spins) :: spin
+        value_average = value/(spin%rows*spin%cols)
+    end function ising_average
+
+    ! update function to update ising system.
+    subroutine update(spin, step, time, row, col)
+        implicit none
+        integer(kind=4) :: step, time
+        integer         :: row, col
+        type(spins)     :: spin
+        real            :: delta_energy, u
+
+        call random_number(u)
+        delta_energy = get_delta_energy(spin, row, col)
+        if (delta_energy < esp .or. exp(-delta_energy/spin%temp) > u) then
+            spin%data(row, col) = -spin%data(row, col)
+        endif
+
+        if (step .eq. 1) then
+            write(*,"(T4,I0,T12,F0.6)") time, ising_average(spin,get_total_energy(spin))
+        endif
+    end subroutine update
+end module handle

fortran/img/demo_cmd.png renamed to fortran/v1/img/demo_cmd.png

@@ -0,0 +1,73 @@
+module ising
+    ! 
+    ! Ising 2d simulation module.
+    ! 
+    ! Author: ph
+    ! Email: huipan.hnu@qq.com
+    ! Date: 2017.07.02
+    ! 
+    use util
+    implicit none
+    type :: spins
+        ! spins type for storage of spin data
+        integer :: rows, cols
+        real    :: temp, magn
+        integer(kind=1), pointer :: data(:, :)
+    end type spins
+contains
+    subroutine init_spins(spin, rows, cols, temp, magn)
+        ! 
+        ! Init spins according to given spin and other pramaters.
+        ! 
+        !   spin    [type(spins)]     spin data storage of ising system.
+        !   rows    [integer]         2d square ising system rows.
+        !   cols    [integer]         2d square ising system cols.
+        !   temp    [real]            temperature of ising system.
+        !   magn    [real]            magnet filed of ising system.
+        !   
+        implicit none
+        integer         ::  rows, cols
+        real            ::  temp, magn, rdata(rows, cols)
+        type(spins)     ::  spin
+        integer         ::  row , col
+        
+        call init_random_seed()
+        call random_number(rdata)
+        allocate(spin%data(rows, cols))
+        where(rdata < 0.5)
+            spin%data = -1
+        elsewhere
+            spin%data = 1
+        end where
+        spin%rows = rows
+        spin%cols = cols
+        spin%temp = temp
+        spin%magn = magn
+    end subroutine init_spins
+
+    subroutine ising_simulate(spin, steps, times, handle)
+        ! 
+        ! Ising module simulation based on Monte-Carl method.
+        ! 
+        !   spin    [type(spins)]   spin data storage of ising system.
+        !   steps   [integer]       steps of inner loop, also sampling steps.
+        !   times   [integer]       outter loop iteration, also sampling times.
+        !   handle  [function]      supply update function for Monte-Carl method.
+        ! 
+        implicit none
+        type(spins)         :: spin
+        integer(kind=4)     :: steps, times, step, time
+        integer             :: row, col
+        real                :: position(2)
+        external            :: handle
+
+        do time = 1, times
+            do step = 1, steps
+                call random_number(position)
+                row = int(position(1)*spin%rows+1)
+                col = int(position(2)*spin%cols+1)
+                call handle(spin, step, time, row, col)
+            end do
+        end do
+    end subroutine ising_simulate
+end module ising

fortran/img/demo_runtime.png renamed to fortran/v1/img/demo_runtime.png

@@ -0,0 +1,45 @@
+#!/usr/bin/env python3
+
+import numpy as np
+import pylab as plb
+import os
+import argparse
+
+def main():
+    """
+    Draw spin file to image and display or save to image file.
+
+    Usage: `spin_plot  [-l] [fdata] [-o fsave]`
+    Params:
+        --fdata     data file from ising simulating.
+        --fsave     image file to save in.
+    """
+    parser = argparse.ArgumentParser()
+    parser.add_argument("-l", "--log", help="plot log data.", action="store_true")
+    parser.add_argument("fdata", help="data file from ising simulation.")
+    parser.add_argument("-o", "--fsave", help="image file to save in.")
+    args = parser.parse_args()
+    if not os.path.isfile(args.fdata):
+        print('Error: data file not exist !')
+    else:
+        # load data file.
+        data = np.loadtxt(args.fdata)
+
+        # plot data to image.
+        if args.log:
+            plb.plot(data[:,0], data[:,1])
+        else:
+            plb.imshow(data)
+            plb.set_cmap('gray')
+
+        # save image file
+        if args.fsave is not None:
+            if args.log:
+                plb.savefig(args.fsave)
+            else:
+                plb.imsave(args.fsave, data)
+        else:
+            plb.show()
+
+if __name__ == '__main__':
+    main()

fortran/img/demo_spin.png renamed to fortran/v1/img/demo_spin.png

@@ -0,0 +1,103 @@
+subroutine show_usage()
+    ! 
+    ! show help message of this pogram.
+    !
+    implicit none
+    write(*,"(A,/)") 'Usage: ising.exe [rows] [cols] [temp] [magn] [steps] [times] [outfile]'
+    write(*,"(A,/)") 'Params:'
+    write(*,"(T4,A,/)") "--[rows]     lattic system rows."
+    write(*,"(T4,A,/)") "--[cols]     lattic system cols."
+    write(*,"(T4,A,/)") "--[temp]     temperature of system."
+    write(*,"(T4,A,/)") "--[magn]     magnet filed of system."
+    write(*,"(T4,A,/)") "--[steps]    steps of sampling in Monte Carlo."
+    write(*,"(T4,A,/)") "--[times]    sampling times."
+    write(*,"(T4,A,/)") "--[outfile]  file to save simulation result of spin data."
+end subroutine show_usage
+
+program main
+    ! 
+    ! ising module 2d simulation main function.
+    ! 
+    ! Usage: `ising [rows] [cols] [temp] [magn] [steps] [times] [outfile]`
+    ! Params:
+    !   --[rows]     lattic system rows.
+    !   --[cols]     lattic system cols.
+    !   --[temp]     temperature of system.
+    !   --[magn]     magnet filed of system.
+    !   --[steps]    steps of sampling in Monte Carlo.
+    !   --[times]    sampling times.
+    !   --[outfile]  optional, file to save simulation result of spin data.
+    ! 
+    ! Example: 
+    !   cmd:
+    !       `ising 64 64 2.27 0. 1000 10 spin.csv`
+    !   out:
+    !          1       -1.975586
+    !          2       -1.974609
+    !          3       -1.974609
+    !          4       -1.979004
+    !          5       -1.977539
+    !          6       -1.976562
+    !          7       -1.981445
+    !          8       -1.986328
+    !          9       -1.985840
+    !          10      -1.986816
+    ! 
+    ! Author: ph
+    ! Email: huipan.hnu@qq.com
+    ! Date: 2017.06.20
+    ! 
+    use util
+    use ising
+    use handle
+    implicit none
+
+    type(spins) :: spin
+    real :: temp, magn
+    character(len=32) :: arg
+    integer :: rows, cols, steps, times, i, narg
+
+    narg = iargc()
+
+    ! only if number of paramters is 5 or 6, program can be continued.
+    if (narg /= 6 .and. narg /= 7) then
+        write(*,"(A,/)") 'Error: number of args not fit !'
+        call show_usage()
+        stop
+    endif
+
+    ! get input command args.
+    do i = 1, 6
+        call getarg(i,arg)
+        select case(i)
+        case(1)
+            read(arg, *) rows
+        case(2)
+            read(arg, *) cols
+        case(3)
+            read(arg, *) temp
+        case(4)
+            read(arg, *) magn
+        case(5)
+            read(arg, *) steps
+        case(6)
+            read(arg, *) times
+        end select
+    end do
+
+    ! init spin site data.
+    call init_spins(spin, rows, cols, temp, magn)
+
+    ! start simulating.
+    call ising_simulate(spin, steps, times, update)
+
+    ! save result spin site data to file.
+    if (narg .eq. 7) then
+        call getarg(i,arg)
+        open(unit=12, file=arg)
+    endif
+    call show_matrix(int(spin%data), 12)
+
+    ! deallocate spin data.
+    deallocate(spin%data)
+end program main

fortran/ising.f90 renamed to fortran/v1/ising.f90

@@ -0,0 +1,24 @@
+program main
+    ! 
+    ! Execute program and count it's running time.
+    ! 
+    ! Usage: runtime [command]
+    ! 
+    ! Author: ph
+    ! Email: huipan.hnu@qq.com
+    ! Date: 2017.06.20
+    ! 
+    implicit none
+    character(len=512) :: cmd_line
+    character(len=32) :: args(32)
+    integer :: i, narg, time(0:1)
+    narg = iargc()
+    do i = 1, narg
+        call getarg(i,args(i))
+    end do
+    write(cmd_line,"(TR1,32A)") args(1:narg)
+    call system_clock(time(0))
+    call system(cmd_line)
+    call system_clock(time(1))
+    write(*,"('[Finished in ',(F0.3),'s]')") real(time(1)-time(0))/1000
+end program main

fortran/ising_main.f90 renamed to fortran/v1/ising_main.f90

@@ -0,0 +1,4 @@
+program main
+    implicit none
+    write(*,"(A,' ',A)") "123","ph"
+end program main