Create an odin model from a file, text string(s) or expression.
The odin_
version is a "standard evaluation" escape hatch.
odin(x, verbose = NULL, target = NULL, workdir = NULL, validate = NULL,
pretty = NULL, skip_cache = NULL, compiler_warnings = NULL,
debug_enable = NULL, no_check_unused_equations = NULL, options = NULL)
odin_(x, verbose = NULL, target = NULL, workdir = NULL,
validate = NULL, pretty = NULL, skip_cache = NULL,
compiler_warnings = NULL, debug_enable = NULL,
no_check_unused_equations = NULL, options = NULL)
Either the name of a file to read, a text string (if length is greater than 1 elements will be joined with newlines) or an expression.
Logical scalar indicating if the compilation should
be verbose. Defaults to the value of the option
odin.verbose
or FALSE
otherwise.
Compilation target. Options are "c", "r" or "js",
defaulting to the option odin.target
or "c" otherwise.
Directory to use for any generated files. This is
only relevant for the "c" target. Defaults to the value of the
option odin.workdir
or tempdir()
otherwise.
Validate the model's intermediate representation
against the included schema. Normally this is not needed and is
intended primarily for development use. Defaults to the value
of the option odin.validate
or FALSE
otherwise.
Pretty-print the model's intermediate
representation. Normally this is not needed and is intended
primarily for development use. Defaults to the value of the
option odin.pretty
or FALSE
otherwise.
Skip odin's cache. This might be useful if the
model appears not to compile when you would expect it to.
Hopefully this will not be needed often. Defaults to the option
odin.skip_cache
or FALSE
otherwise.
Previously this attempted detection of compiler warnings (with some degree of success), but is currently ignored. This may become supported again in a future version depending on underlying support in pkgbuild.
Enable debugging commands in generated code
(currently print()
). If TRUE
then these are generated by
odin targets that support them, and will generally make your
program slower.
If TRUE
, then don't print
messages about unused variables. Defaults to the option
odin.no_check_unused_equations
or FALSE
otherwise.
Named list of options. If provided, then all other options are ignored.
An odin_generator
object (an R6 class) which can be used
to create model instances.
Do not use odin::odin
in a package; you almost certainly want
to use odin_package instead.
A generated model can return information about itself; odin_ir
If the model accepts user parameters, then the parameter to the
constructor or the $set_user()
method can be used to control
the behaviour when unknown user actions are passed into the
model. Possible values are the strings stop
(throw an error),
warning
(issue a warning but keep going), message
(print a
message and keep going) or ignore
(do nothing). Defaults to
the option odin.unused_user_action
, or warning
otherwise.
When generating a model one must chose between using the
dde
package to solve the system or the default
deSolve
. Future versions may allow this to switch when
using run
, but for now this requires tweaking the generated
code to a point where one must decide at generation. dde
implements only the Dormand-Prince 5th order dense output solver,
with a delay equation solver that may perform better than the
solvers in deSolve. For non-delay equations, deSolve
is
very likely to outperform the simple solver implemented.
## Compile the model; exp_decay here is an R6ClassGenerator and will
## generate instances of a model of exponential decay:
exp_decay <- odin::odin({
deriv(y) <- -0.5 * y
initial(y) <- 1
}, target = "r")
#> Generating model in r
## Generate an instance; there are no parameters here so all instances
## are the same and this looks a bit pointless. But this step is
## required because in general you don't want to have to compile the
## model every time it is used (so the generator will go in a
## package).
mod <- exp_decay$new()
## Run the model for a series of times from 0 to 10:
t <- seq(0, 10, length.out = 101)
y <- mod$run(t)
plot(y, xlab = "Time", ylab = "y", main = "", las = 1)