readbinarydisc: read disconnected loops from binary files

In hadron: Analysis Framework for Monte Carlo Simulation Data in Physics

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.

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.

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.

Author(s)

Carsten Urbach, [e-mail]

See Also

readcmidatafiles, readbinarycf, readcmidisc, readcmicor

Examples

## 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"))

hadron documentation built on Sept. 9, 2022, 5:06 p.m.