hadron (version 3.1.2)

readcmidisc: 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.

See Also

readcmidatafiles, readbinarycf, readbinarydisc, readcmicor

Examples

# NOT RUN {
# 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)
# }