|
4 | 4 | #include "generics.h" |
5 | 5 | #include "internal.h" |
6 | 6 |
|
7 | | -static int host1 = 0; |
8 | | -static int host2 = 1; |
9 | | -static int outside = 2; |
| 7 | +static const int host1 = 0; |
| 8 | +static const int host2 = 1; |
| 9 | +static const int outside = 2; |
10 | 10 |
|
11 | 11 | //! S2I2R2 process state. |
12 | 12 | typedef struct { |
@@ -148,77 +148,77 @@ double s2i2r2_proc_t::event_rates (double *rate, int n) const { |
148 | 148 | template<> |
149 | 149 | void s2i2r2_genealogy_t::rinit (void) { |
150 | 150 | state.S1 = params.S1_0; |
151 | | -state.I1 = params.I1_0; |
152 | | -state.R1 = params.R1_0; |
153 | | -state.S2 = params.S2_0; |
154 | | -state.I2 = params.I2_0; |
155 | | -state.R2 = params.R2_0; |
156 | | -state.N1 = double(params.S1_0+params.I1_0+params.R1_0); |
157 | | -state.N2 = double(params.S2_0+params.I2_0+params.R2_0); |
158 | | -graft(host1,params.I1_0); |
159 | | -graft(host2,params.I2_0); |
| 151 | + state.I1 = params.I1_0; |
| 152 | + state.R1 = params.R1_0; |
| 153 | + state.S2 = params.S2_0; |
| 154 | + state.I2 = params.I2_0; |
| 155 | + state.R2 = params.R2_0; |
| 156 | + state.N1 = double(params.S1_0+params.I1_0+params.R1_0); |
| 157 | + state.N2 = double(params.S2_0+params.I2_0+params.R2_0); |
| 158 | + graft(host1,params.I1_0); |
| 159 | + graft(host2,params.I2_0); |
160 | 160 | } |
161 | 161 |
|
162 | 162 | template<> |
163 | 163 | void s2i2r2_genealogy_t::jump (int event) { |
164 | 164 | switch (event) { |
165 | 165 | case 0: |
166 | | - state.S1 -= 1; state.I1 += 1; birth(host1,host1); |
167 | | - break; |
168 | | - case 1: |
169 | | - state.S2 -= 1; state.I2 += 1; birth(host2,host2); |
170 | | - break; |
171 | | - case 2: |
172 | | - state.S1 -= 1; state.I1 += 1; birth(host2,host1); |
173 | | - break; |
174 | | - case 3: |
175 | | - state.I1 -= 1; state.R1 += 1; death(host1); |
176 | | - break; |
177 | | - case 4: |
178 | | - state.I2 -= 1; state.R2 += 1; death(host2); |
179 | | - break; |
180 | | - case 5: |
181 | | - state.R1 -= 1; state.S1 += 1; |
182 | | - break; |
183 | | - case 6: |
184 | | - state.R2 -= 1; state.S2 += 1; |
185 | | - break; |
186 | | - case 7: |
187 | | - sample(host1); |
188 | | - break; |
189 | | - case 8: |
190 | | - sample(host2); |
191 | | - break; |
192 | | - case 9: |
193 | | - state.S1 -= 1; state.I1 += 1; graft(outside); migrate(outside,host1); |
194 | | - break; |
195 | | - case 10: |
196 | | - state.S2 -= 1; state.I2 += 1; graft(outside); migrate(outside,host2); |
197 | | - break; |
198 | | - case 11: |
199 | | - state.S1 -= 1; state.N1 -= 1; |
200 | | - break; |
201 | | - case 12: |
202 | | - state.S2 -= 1; state.N2 -= 1; |
203 | | - break; |
204 | | - case 13: |
205 | | - state.I1 -= 1; state.N1 -= 1; death(host1); |
206 | | - break; |
207 | | - case 14: |
208 | | - state.I2 -= 1; state.N2 -= 1; death(host2); |
209 | | - break; |
210 | | - case 15: |
211 | | - state.R1 -= 1; state.N1 -= 1; |
212 | | - break; |
213 | | - case 16: |
214 | | - state.R2 -= 1; state.N2 -= 1; |
215 | | - break; |
216 | | - case 17: |
217 | | - state.S1 += 1; state.N1 += 1; |
218 | | - break; |
219 | | - case 18: |
220 | | - state.S2 += 1; state.N2 += 1; |
221 | | - break; |
| 166 | + state.S1 -= 1; state.I1 += 1; birth(host1,host1); |
| 167 | + break; |
| 168 | + case 1: |
| 169 | + state.S2 -= 1; state.I2 += 1; birth(host2,host2); |
| 170 | + break; |
| 171 | + case 2: |
| 172 | + state.S1 -= 1; state.I1 += 1; birth(host2,host1); |
| 173 | + break; |
| 174 | + case 3: |
| 175 | + state.I1 -= 1; state.R1 += 1; death(host1); |
| 176 | + break; |
| 177 | + case 4: |
| 178 | + state.I2 -= 1; state.R2 += 1; death(host2); |
| 179 | + break; |
| 180 | + case 5: |
| 181 | + state.R1 -= 1; state.S1 += 1; |
| 182 | + break; |
| 183 | + case 6: |
| 184 | + state.R2 -= 1; state.S2 += 1; |
| 185 | + break; |
| 186 | + case 7: |
| 187 | + sample(host1); |
| 188 | + break; |
| 189 | + case 8: |
| 190 | + sample(host2); |
| 191 | + break; |
| 192 | + case 9: |
| 193 | + state.S1 -= 1; state.I1 += 1; graft(outside); migrate(outside,host1); |
| 194 | + break; |
| 195 | + case 10: |
| 196 | + state.S2 -= 1; state.I2 += 1; graft(outside); migrate(outside,host2); |
| 197 | + break; |
| 198 | + case 11: |
| 199 | + state.S1 -= 1; state.N1 -= 1; |
| 200 | + break; |
| 201 | + case 12: |
| 202 | + state.S2 -= 1; state.N2 -= 1; |
| 203 | + break; |
| 204 | + case 13: |
| 205 | + state.I1 -= 1; state.N1 -= 1; death(host1); |
| 206 | + break; |
| 207 | + case 14: |
| 208 | + state.I2 -= 1; state.N2 -= 1; death(host2); |
| 209 | + break; |
| 210 | + case 15: |
| 211 | + state.R1 -= 1; state.N1 -= 1; |
| 212 | + break; |
| 213 | + case 16: |
| 214 | + state.R2 -= 1; state.N2 -= 1; |
| 215 | + break; |
| 216 | + case 17: |
| 217 | + state.S1 += 1; state.N1 += 1; |
| 218 | + break; |
| 219 | + case 18: |
| 220 | + state.S2 += 1; state.N2 += 1; |
| 221 | + break; |
222 | 222 | default: // #nocov |
223 | 223 | assert(0); // #nocov |
224 | 224 | break; // #nocov |
|
0 commit comments