Bug in SEIRConn and incubation

UofUEpiBio · Jun 13, 2023 · 3da5e3f · 3da5e3f
1 parent 3b890bf
commit 3da5e3f
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/epiworld.hpp b/epiworld.hpp
@@ -10366,7 +10366,7 @@ inline epiworld_double Virus<TSeq>::get_incubation(
     if (incubation_fun)
         return incubation_fun(agent, *this, model);
 
-    return EPI_DEFAULT_VIRUS_INCUBATION;
+    return EPI_DEFAULT_INCUBATION_DAYS;
 
 }
 
@@ -16386,7 +16386,7 @@ inline ModelSEIRCONN<TSeq>::ModelSEIRCONN(
                 auto & v = p->get_virus(0u);
 
                 // Does the agent become infected?
-                if (m->runif() < 1.0/(v.get_incubation(m)))
+                if (m->runif() < 1.0/(v->get_incubation(m)))
                 {
 
                     p->change_state(m, ModelSEIRCONN<TSeq>::INFECTED);

diff --git a/include/epiworld/models/seirconnected.hpp b/include/epiworld/models/seirconnected.hpp
@@ -204,7 +204,7 @@ inline ModelSEIRCONN<TSeq>::ModelSEIRCONN(
                 auto & v = p->get_virus(0u);
 
                 // Does the agent become infected?
-                if (m->runif() < 1.0/(v.get_incubation(m)))
+                if (m->runif() < 1.0/(v->get_incubation(m)))
                 {
 
                     p->change_state(m, ModelSEIRCONN<TSeq>::INFECTED);