hadron (version 3.1.2)

readbinarydisc: read disconnected loops from binary files

Description

Reads disconnected loops from binary files.

Usage

readbinarydisc(files, Time = 48, obs = 5, endian = "little",
  excludelist = c(""), nrSamples = 1, path = "")

Arguments

files

list of filenames to be read. Can be created for instance using getorderedfilelist. The filelist is assumed to be ordered with number of samples running fastest, and the next to fastest nubmer of gauges.

Time

time extent of correlation functions.

obs

each file may contain Time*obs correlation functions. With obs one choses which observable to read in.

endian

the endianess of the binary file.

excludelist

files to exclude from reading.

nrSamples

the number of samples

path

path to be prepended to every filename.

Value

returns a list with two arrays cf and icf with real and imaginary parts of the loops, and integers Time, nrStypes=1, nrSamples and nrObs=1. Both of the arrays have dimension c(Time, N), where N is the number of measurements (gauges) and Time the time extent, nrStypes the number of smearing levels and nrObs the number of operators, both of which are currently fixed to 1.

Details

It is assumend that each file contains O*Time complex doubles, where Time is the time extent and O the number of observables in the file. It is assumed that complex is the fastest running index, next time and then observables. The different samples are assumend to be in different files. The file list is assumed to be ordered with number of samples running fastest, and then number of gauges.

See Also

readcmidatafiles, readbinarycf, readcmidisc, readcmicor

Examples

# NOT RUN {
## running toy example
file <- paste0(system.file("extdata", package = "hadron"), "/C2_pi0.dat")
X <- readbinarydisc(files=file, Time=64, obs=0)
X$cf

## more realistic example
# }
# NOT RUN {
files <- character()
# }
# NOT RUN {
for(i in seq(600,1744,8)) 
# }
# NOT RUN {
  files <- c(files, "C2_dis_u_conf", sprintf("%.04d", i), ".dat", sep="")
# }
# NOT RUN {
cf <- readbinarydisc(files, obs=4, excludelist=c("C2_pi0_conf0632.dat"))
# }
# NOT RUN {
# }