mrgsolve-tad.h

Go to the documentation of this file.

 
 
namespace mrgsolve {
class tadose {
public: 
  int cmt; 
  double told; 
  bool had_dose; 
  tadose(int cmt_);
  tadose();
  double tad(databox& self);
  void reset();
};
tadose::tadose(int cmt_) {
  cmt = cmt_;  
  told = -1e9;
}
tadose::tadose() {
  told = -1e9;  
}
 
double tadose::tad(databox& self) {
  if(self.newind <= 1) this->reset();
  if(self.evid==0) {
    return had_dose ? self.time - told : -1.0;
  }
  if((self.cmt == cmt) && (self.evid == 1 || self.evid == 4)) {
    told = self.time;
    had_dose = true;
    return 0.0;
  }
  return had_dose ? self.time - told : -1.0;
}
void tadose::reset() {
  had_dose = false;
  told = -1e9;
}
}