dMod (version 1.0.0)

symmetryDetection: Search for symmetries in the loaded model

Description

This function follows the method published in [1].

The function calls a python script via rPython. Usage problems might occur when different python versions are used. The script was written and tested for python 2.7.12, sympy 0.7.6.

Recently, users went into problems with RJSONIO when rPython was used. Unless a sound solution is available, please try to reinstall RJSONIO in these cases.

Usage

symmetryDetection(f, obsvect = NULL, prediction = NULL,
  initial = NULL, ansatz = "uni", pMax = 2, inputs = NULL,
  fixed = NULL, cores = 1, allTrafos = FALSE)

Arguments

f

object containing the ODE for which as.eqnvec() is defined

obsvect

vector of observation functions

prediction

vector containing prediction to be tested

initial

vector containing initial values

ansatz

type of infinitesimal ansatz used for the analysis (uni, par, multi)

pMax

maximal degree of infinitesimal ansatz

inputs

specify the input variables

fixed

variables to concider fixed

cores

maximal number of cores used for the analysis

allTrafos

do not remove transformations with a common parameter factor

References

[1] https://journals.aps.org/pre/abstract/10.1103/PhysRevE.92.012920

Examples

# NOT RUN {
eq <- NULL
eq <- addReaction(eq, "A", "B", "k1*A")
eq <- addReaction(eq, "B", "A", "k2*B")

observables <- eqnvec(Aobs = "alpha * A")

symmetryDetection(eq, observables)

# }