readcmidisc: reads disconnected loops in cmi format

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

reads disconnected loops in cmi format

Description

reads disconnected loops in cmi (Chris Michael) format from a list of files.

Usage

readcmidisc(files, obs = 9, ind.vec = c(2, 3, 4, 5, 6, 7, 8),
  excludelist = c(""), skip = 0, L, colClasses = c("integer", "integer",
  "integer", "integer", "numeric", "numeric", "numeric", "numeric"),
  debug = FALSE)

Arguments

files	list of filenames to be read. Can be created using getorderedfilelist.
obs	index of operator to parse from files
ind.vec	vector containing the index (column in file) of obs, t, samples, Re(local), Im(local, Re(smeared), Im(smeared).
excludelist	files to exclude from reading.
skip	lines to skip at beginning of each file.
L	the spatial lattice extent, set to Time/2 if missing.
colClasses	The column data type classes, the read.table.
debug	setting debug to TRUE makes the routine more verbose by spilling out separate filenames.

Value

returns a list with four arrays cf, icf scf and sicf containing real and imaginary parts of the local and smeared loops, respectively, and integers Time, nrStypes=2, nrSamples and nrObs=1. The four arrays have dimension c(Time, S, N), where S is the nubmer of samples, Time is the time extent and N is the number of measurements (gauges). Time is the time extent, nrStypes the number of smearing levels and nrObs the number of operators, which are currently fixed to 1 and 2, respectively. nrSamples is the number of samples.

Note that the arrays are normalised by 1/sqrt(L^2).

The routine expects that all files have identical content. Otherwise the routine will stop.

Author(s)

Carsten Urbach, [e-mail]

See Also

readcmidatafiles, readbinarycf, readbinarydisc, readcmicor

Examples

# a running toy example
hpath <- system.file(package="hadron")
files <- paste0(hpath, "/extdata/newdisc.0.1373.0.006.k0v4.10")
X <- readcmidisc(files=files)
X

## a more realistic example
## Not run: v4files <- character()
## Not run: for(i in seq(600,1744,8))
## Not run:   v4files <- 
## Not run:    c(v4files, paste("disc.0.163265.0.006.k0v4.", sprintf("%.04d", i), sep=""))
## Not run: v4data <- readcmidisc(v4files)

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