mrgsolve-tad.h

Go to the documentation of this file.

 
 
namespace mrgsolve {
class tadose {
public: 
  int cmt = -1e9; 
  double told = -1e9; 
  bool had_dose = false; 
  tadose(int cmt_);
  tadose(){};
  double tad(databox& self);
  void reset();
};

tadose::tadose(int cmt_) {
  cmt = cmt_;
}

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;
}
}