mrgsolve
|
#include <odeproblem.h>
Public Member Functions | |
odeproblem (Rcpp::NumericVector param, Rcpp::NumericVector init, Rcpp::List funs, int n_capture_) | |
~odeproblem () | |
void | do_init_calc (bool answer) |
void | advance (double tfrom, double tto, LSODA &solver) |
void | call_derivs (double *t, double *y, double *ydot) |
void | init (int pos, double value) |
double | init (int pos) |
void | init_call (const double &time) |
void | init_call_record (const double &time) |
void | init_derivs (double time) |
void | y_init (int pos, double value) |
void | y_init (Rcpp::NumericVector init) |
void | y_add (const unsigned int pos, const double &value) |
add value to compartment pos More... | |
void | table_call () |
Call $TABLE function. More... | |
void | table_init_call () |
void | config_call () |
Call $PREAMBLE function. More... | |
void | set_d (rec_ptr this_rec) |
void | omega (Rcpp::NumericMatrix &x) |
void | sigma (Rcpp::NumericMatrix &x) |
arma::mat | mv_omega (int n) |
arma::mat | mv_sigma (int n) |
void | pass_envir (Rcpp::Environment *x) |
bool | CFONSTOP () |
const std::vector< double > & | param () |
void | param (int pos, double value) |
void | rate (unsigned int pos, double value) |
double | rate (unsigned int pos) |
void | rate0 (unsigned int pos, double value) |
double | rate0 (unsigned int pos) |
int | rate_count (unsigned int pos) |
void | rate_add (unsigned int pos, const double &value) |
void | rate_rm (unsigned int pos, const double &value) |
void | rate_bump (const unsigned int pos, const double &value) |
void | rate_main (rec_ptr rec) |
void | rate_reset () |
Reset all infusion rates. More... | |
void | dur (unsigned int pos, double value) |
double | dur (unsigned int pos) |
void | reset_newid (const double id_) |
Reset odeproblem object for new individual. More... | |
void | eta (int pos, double value) |
void | eps (int pos, double value) |
unsigned short int | systemoff () |
void | on (unsigned short int cmt) |
void | off (unsigned short int cmt) |
int | is_on (unsigned int eq_n) |
void | fbio (unsigned int pos, double value) |
double | fbio (unsigned int pos) |
double | alag (int cmt) |
void | time (double time_) |
void | newind (unsigned int newind_) |
unsigned int | newind () |
void | advan (int x) |
int | advan () |
void | advan2 (const double &tfrom, const double &tto) |
void | advan4 (const double &tfrom, const double &tto) |
void | neta (const int n) |
set number of ETAs in the model More... | |
void | neps (const int n) |
set number of EPSs in the model More... | |
void | nid (int n) |
sets the number of IDs More... | |
void | nrow (int n) |
sets the number of data set rows More... | |
void | idn (int n) |
sets the current ID number More... | |
void | rown (int n) |
sets the currenw data set row number More... | |
dvec & | get_capture () |
double | capture (int i) |
void | copy_parin (const Rcpp::List &parin) |
copies items passed in through parin into the odeproblem object More... | |
void | copy_funs (const Rcpp::List &funs) |
bool | any_mtime () |
std::vector< mrgsolve::evdata > | mtimes () |
void | clear_mtime () |
void | y (const int pos, const double value) |
double | y (const int pos) |
int | istate () |
void | istate (int value) |
void | lsoda_init () |
int | npar () |
returns the number of parameters More... | |
int | neq () |
returns the number of state variables More... | |
void | tol (double atol, double rtol) |
sets the absolute and relative tolerances More... | |
Public Attributes | |
std::vector< double > | Y |
compartment amounts More... | |
std::vector< double > | Ydot |
dxdt values More... | |
std::vector< double > | Yout |
used to hold Y values during solving More... | |
std::vector< double > | Param |
parameter vector More... | |
std::vector< double > | Capture |
captured data items More... | |
double | Atol |
absolute tolerance used by ODE solver More... | |
double | Rtol |
relative tolerance used by ODE solver More... | |
double | ssAtol |
absolute tolerance when finding steady state More... | |
double | ssRtol |
int | Npar |
number of parameters More... | |
int | Neq |
number of equations More... | |
int | Istate |
istate value More... | |
bool | ss_fixed |
If true, then no warning is issued if SS not reached in ss_n doses. More... | |
int | ss_n |
Max number of doses during SS advance before warning is issued. More... | |
bool | ss_flag |
flag indicating when the system is advancing to SS More... | |
std::vector< int > | Ss_cmt |
vector of compartments to consider for SS More... | |
std::vector< double > | R0 |
acutal current infusion rate More... | |
std::vector< unsigned int > | infusion_count |
number of active infusions More... | |
std::vector< double > | R |
receive user input for infusion rate More... | |
std::vector< double > | D |
receive user input for infusion duration More... | |
std::vector< double > | Init_value |
initial conditions More... | |
std::vector< double > | Init_dummy |
initial conditions for user input More... | |
std::vector< double > | F |
bioavability More... | |
std::vector< double > | Alag |
dosing lag time More... | |
std::vector< int > | On |
compartment on/off indicator More... | |
databox | d |
various data passed to model functions More... | |
int | Advan |
simulation mode: 1/2/3/4 (PK models) or 13 (odes) More... | |
std::vector< double > | a |
used for advan 1/2/3/4 calculations More... | |
std::vector< double > | alpha |
used for advan 1/2/3/4 calculation More... | |
mrgsolve::resim | simeta |
functor for resimulating etas More... | |
mrgsolve::resim | simeps |
functor for resimulating epsilons More... | |
arma::mat | Omega |
variance/covariance matrix for between-subject variability More... | |
arma::mat | Sigma |
variance/covariance matrix for within-subject variability More... | |
std::vector< double > | pred |
brings clearances, volumes, and & for advan 1/2/3/4 More... | |
deriv_func | Derivs |
$ODE function More... | |
init_func | Inits |
$MAIN function More... | |
table_func | Table |
$TABLE function More... | |
config_func | Config |
$PREAMBLE function More... | |
bool | Do_Init_Calc |
Flag regulating whether or not initials are taken from $MAIN More... | |
odeproblem::odeproblem | ( | Rcpp::NumericVector | param, |
Rcpp::NumericVector | init, | ||
Rcpp::List | funs, | ||
int | n_capture_ | ||
) |
|
inline |
|
inline |
void odeproblem::advan | ( | int | x | ) |
void odeproblem::advan2 | ( | const double & | tfrom, |
const double & | tto | ||
) |
void odeproblem::advan4 | ( | const double & | tfrom, |
const double & | tto | ||
) |
void odeproblem::advance | ( | double | tfrom, |
double | tto, | ||
LSODA & | solver | ||
) |
double odeproblem::alag | ( | int | cmt | ) |
|
inline |
void odeproblem::call_derivs | ( | double * | t, |
double * | y, | ||
double * | ydot | ||
) |
|
inline |
|
inline |
|
inline |
void odeproblem::config_call | ( | ) |
Call $PREAMBLE
function.
void odeproblem::copy_funs | ( | const Rcpp::List & | funs | ) |
void odeproblem::copy_parin | ( | const Rcpp::List & | parin | ) |
copies items passed in through parin into the odeproblem object
|
inline |
|
inline |
|
inline |
|
inline |
|
inline |
double odeproblem::fbio | ( | unsigned int | pos | ) |
|
inline |
|
inline |
|
inline |
sets the current ID number
|
inline |
|
inline |
void odeproblem::init_call | ( | const double & | time | ) |
Call $MAIN
to get the initial conditions.
time | the time to assume for the calculation |
void odeproblem::init_call_record | ( | const double & | time | ) |
Call $MAIN
with the dummy initial condition vector.
time | the time to assume when making the call. |
void odeproblem::init_derivs | ( | double | time | ) |
|
inline |
|
inline |
|
inline |
|
inline |
|
inline |
arma::mat odeproblem::mv_omega | ( | int | n | ) |
arma::mat odeproblem::mv_sigma | ( | int | n | ) |
void odeproblem::neps | ( | const int | n | ) |
set number of EPSs
in the model
|
inline |
returns the number of state variables
void odeproblem::neta | ( | const int | n | ) |
set number of ETAs
in the model
|
inline |
|
inline |
|
inline |
sets the number of IDs
|
inline |
returns the number of parameters
|
inline |
sets the number of data set rows
void odeproblem::off | ( | unsigned short int | cmt | ) |
void odeproblem::omega | ( | Rcpp::NumericMatrix & | x | ) |
void odeproblem::on | ( | unsigned short int | cmt | ) |
|
inline |
|
inline |
|
inline |
|
inline |
|
inline |
|
inline |
|
inline |
void odeproblem::rate_add | ( | unsigned int | pos, |
const double & | value | ||
) |
void odeproblem::rate_bump | ( | const unsigned int | pos, |
const double & | value | ||
) |
|
inline |
void odeproblem::rate_main | ( | rec_ptr | rec | ) |
void odeproblem::rate_reset | ( | ) |
Reset all infusion rates.
void odeproblem::rate_rm | ( | unsigned int | pos, |
const double & | value | ||
) |
void odeproblem::reset_newid | ( | const double | id_ | ) |
Reset odeproblem
object for new individual.
|
inline |
sets the currenw data set row number
void odeproblem::set_d | ( | rec_ptr | this_rec | ) |
void odeproblem::sigma | ( | Rcpp::NumericMatrix & | x | ) |
|
inline |
void odeproblem::table_call | ( | ) |
Call $TABLE
function.
void odeproblem::table_init_call | ( | ) |
|
inline |
void odeproblem::tol | ( | double | atol, |
double | rtol | ||
) |
sets the absolute and relative tolerances
|
inline |
|
inline |
void odeproblem::y_add | ( | const unsigned int | pos, |
const double & | value | ||
) |
add value
to compartment pos
void odeproblem::y_init | ( | int | pos, |
double | value | ||
) |
Assigns values to both the compartment and the vector of initial conditions.
pos | the compartment number (C++ indexing) |
value | the value for the compartment |
void odeproblem::y_init | ( | Rcpp::NumericVector | init | ) |
std::vector<double> odeproblem::a |
used for advan 1/2/3/4 calculations
int odeproblem::Advan |
simulation mode: 1/2/3/4 (PK models) or 13 (odes)
std::vector<double> odeproblem::Alag |
dosing lag time
std::vector<double> odeproblem::alpha |
used for advan 1/2/3/4 calculation
double odeproblem::Atol |
absolute tolerance used by ODE solver
std::vector<double> odeproblem::Capture |
captured data items
config_func odeproblem::Config |
$PREAMBLE
function
std::vector<double> odeproblem::D |
receive user input for infusion duration
databox odeproblem::d |
various data passed to model functions
deriv_func odeproblem::Derivs |
$ODE
function
bool odeproblem::Do_Init_Calc |
Flag regulating whether or not initials are taken from $MAIN
std::vector<double> odeproblem::F |
bioavability
std::vector<unsigned int> odeproblem::infusion_count |
number of active infusions
std::vector<double> odeproblem::Init_dummy |
initial conditions for user input
std::vector<double> odeproblem::Init_value |
initial conditions
init_func odeproblem::Inits |
$MAIN
function
int odeproblem::Istate |
istate value
int odeproblem::Neq |
number of equations
int odeproblem::Npar |
number of parameters
arma::mat odeproblem::Omega |
variance/covariance matrix for between-subject variability
std::vector<int> odeproblem::On |
compartment on/off indicator
std::vector<double> odeproblem::Param |
parameter vector
std::vector<double> odeproblem::pred |
brings clearances, volumes, and & for advan 1/2/3/4
std::vector<double> odeproblem::R |
receive user input for infusion rate
std::vector<double> odeproblem::R0 |
acutal current infusion rate
double odeproblem::Rtol |
relative tolerance used by ODE solver
arma::mat odeproblem::Sigma |
variance/covariance matrix for within-subject variability
mrgsolve::resim odeproblem::simeps |
functor for resimulating epsilons
mrgsolve::resim odeproblem::simeta |
functor for resimulating etas
std::vector<int> odeproblem::Ss_cmt |
vector of compartments to consider for SS
bool odeproblem::ss_fixed |
If true, then no warning is issued if SS not reached in ss_n doses.
bool odeproblem::ss_flag |
flag indicating when the system is advancing to SS
int odeproblem::ss_n |
Max number of doses during SS advance before warning is issued.
double odeproblem::ssAtol |
absolute tolerance when finding steady state
double odeproblem::ssRtol |
table_func odeproblem::Table |
$TABLE
function
std::vector<double> odeproblem::Y |
compartment amounts
std::vector<double> odeproblem::Ydot |
dxdt values
std::vector<double> odeproblem::Yout |
used to hold Y values during solving