@@ -32,22 +32,22 @@ tissue1 = construct(Tissue, deepcopy([population, population2])) # Make a Tissue
32
32
tissue2 = construct (Tissue, deepcopy ([population2, population]))
33
33
embryo = construct (Embryo, deepcopy ([tissue1, tissue2])) # Make an embryo from Tissues
34
34
35
- cell_ode = function (t, cell, dcell )
35
+ cell_ode = function (dcell, cell,p,t )
36
36
m = mean (cell)
37
37
for i in eachindex (cell)
38
38
dcell[i] = - m* cell[i]
39
39
end
40
40
end
41
41
42
- f = function (t, embryo, dembryo )
42
+ f = function (dembryo, embryo,p,t )
43
43
for (cell, y, z) in LevelIterIdx (embryo, 2 )
44
- cell_ode (t, cell, @view dembryo[y: z])
44
+ cell_ode (@view ( dembryo[y: z]),cell,p,t )
45
45
end
46
46
end
47
47
48
48
tstop = [0.5 ]
49
49
50
- condition = function (t, u , integrator)
50
+ condition = function (u, t , integrator)
51
51
t ∈ tstop
52
52
end
53
53
@@ -82,7 +82,7 @@ sol = solve(prob, Rosenbrock23(), callback=shrinking_cb, tstops=tstop)
82
82
83
83
println (" Do the SDE Part" )
84
84
85
- g = function (t, u, du )
85
+ g = function (du,u,p,t )
86
86
for i in eachindex (u)
87
87
du[i] = 0.1 u[i]
88
88
end
0 commit comments