
################################################################################

                   Computer simulation of physical systems I                    
            Blocking analysis for error estimation on time averages

################################################################################


--------------------------------------------------------------------------------
 Files provided and description
--------------------------------------------------------------------------------

 The following files are provided:

 *DOCUMENTATION
   - README.txt (this file)
   - installer (installer script)
   - documentation.pdf (notes on this part of the course, exercise)

 *CODES (in ./src folder)
   - Makefile          makefile
   - stochastic.f90    generates a sequence of exponentially correlated Gaussian
                       random numbers ; 
   - histogram.f90     calculates the probability distribution function ;
   - correlation.f90   calculates the correlation function ;
   - blocking.f90      performs the blocking analysis ;
 
 *FUNCTIONS
   - ran1.f            uniform pseudorandom numbers generator
   - gasdev.f          generates normally distributed deviate with zero mean and
                       unit variance           
 
 
--------------------------------------------------------------------------------
 Quick installation
--------------------------------------------------------------------------------

 To compile the programs from source codes type
 
   chmod +x installer
   ./installer
   
 A work directory is created in which you find a script file, launcher.sh.  Modify
 it according to your own needs and type "./launcher.sh" to perform statistical 
 analysis.
   

--------------------------------------------------------------------------------
 INPUT FILE
--------------------------------------------------------------------------------
 
 Use the same input file for stochastic, correlation and blocking programs 
 (see below an example).
 
  &stochastic
   nstep=10000                > number of data that will be generated (N)
   tau=10                     > correlation time (\tau)
  /
  &correlation
   inputfile="traj.dat"       > input file
   nstep=10000                > number of data in input (N) 
   dt=1                       > time step between two samples (use 1 -> internal unit)
   timemax=100                > max time for the autocorrelation function
   nsamp=1                    > steps between two samples of the autocorrelation function 
   it0=40                     > steps between two starting time
  /
  &blocking
   inputfile="traj.dat"       > input file
   nstep=10000                > number of data in input (N)
  /
  
################################################################################
################################################################################

