ising_MC_drv.F90

program ising_drv
  implicit none
  real    :: x1,x2,x3,y1,y2,y3,T
  integer :: niter,i
  !integer, allocatable :: CCC(:,:,:)
  !real, allocatable :: ZZZc(:,:,:)
  integer :: CCC(4,4,2**16),NX,NC
  real :: ZZZc(2**16)

     
  !write(*,*) " x1,x2,x3 , NC , NX"
  !read(*,*) x1,x2,x3,NC,NX
  !write(*,*) x1,x2,x3 ,NC,NX

  NC = 10*(10**6)  !10**7 !3*(10**6) ! 10**6
  NX = 64
  x2=0.
  x3=0.
  !do i=1,15
  do i=4,4
     T=2.28  !0.5*i
     x1=1./T 
     !call rng_sub( x1,x2,x3,y1,y2,y3 )
     !call ising_4x4_partfun( x1,x2,x3,ZZZc,CCC)
     !call ising_MC_partfun( NC,NX,x1,x2,x3)  !,ZZZc,CCC)
     call ising_MC_rlzn( NC,NX,x1,x2,x3)  !,ZZZc,CCC)
   end do

     
end program ising_drv