qsolve: Solve a linear system over the rationals

In dkahle/latter: Interface to 'LattE' and '4ti2'

Solve a linear system over the rationals

Description

qsolve runs 4ti2's qsolve program to compute the configuration matrix A corresponding to graphical statistical models given by a simplicial complex and levels on the nodes.

Usage

qsolve(mat, rel, sign, dir = tempdir(), quiet = TRUE, shell = FALSE, ...)

Arguments

mat	The A matrix (see the 4ti2 documentation or examples)
rel	A vector of "<" or ">" relations
sign	The signs of the individual
dir	Directory to place the files in, without an ending /
quiet	If FALSE, messages the 4ti2 output
shell	Messages the shell code used to do the computation
...	Additional arguments to pass to the function

Value

The configuration matrix of the model provided

Examples

if (has_4ti2()) {

# x + y > 0
# x + y < 0

mat <- rbind(
  c( 1,  1),
  c( 1,  1)
)
rel <- c(">", "<")
sign <- c(0, 0)

qsolve(mat, rel, sign, p = "arb")
qsolve(mat, rel, sign, p = "arb", quiet = FALSE)
qsolve(mat, rel, sign, p = "arb", shell = TRUE)

}

dkahle/latter documentation built on March 22, 2023, 9:05 p.m.